A unifying theory of a posteriori error control for nonconforming finite element methods

Residual-based a posteriori error estimates were derived within one unifying framework for lowest-order conforming, nonconforming, and mixed finite element schemes in Carstensen [Numer Math 100:617–637, 2005]. Therein, the key assumption is that the conforming first-order finite element space annulates the linear and bounded residual ℓ written . That excludes particular nonconforming finite element methods (NCFEMs) on parallelograms in that . The present paper generalises the aforementioned theory to more general situations to deduce new a posteriori error estimates, also for mortar and discontinuous Galerkin methods. The key assumption is the existence of some bounded linear operator with some elementary properties. It is conjectured that the more general hypothesis (H1)–(H3) can be established for all known NCFEMs. Applications on various nonstandard finite element schemes for the Laplace, Stokes, and Navier–Lamé equations illustrate the presented unifying theory of a posteriori error control for NCFEM. Supported by DFG Research Center MATHEON “Mathematics for key technologies” in Berlin and the German Indian Project DST-DAAD (PPP-05). J. Hu was partially supported by National Science Foundation of China under Grant No.10601003.     

<Emphasis Type="Italic">hp</Emphasis>-Discontinuous Galerkin methods for strongly nonlinear elliptic boundary value problems

In this paper, we have analyzed a one parameter family of hp-discontinuous Galerkin methods for strongly nonlinear elliptic boundary value problems with Dirichlet boundary conditions. These methods depend on the values of the parameter , where θ = + 1 corresponds to the nonsymmetric and θ = −1 corresponds to the symmetric interior penalty methods when and f(u,∇u) = −f, that is, for the Poisson problem. The error estimate in the broken H ¹ norm, which is optimal in h (mesh size) and suboptimal in p (degree of approximation) is derived using piecewise polynomials of degree p ≥ 2, when the solution . In the case of linear elliptic problems also, this estimate is optimal in h and suboptimal in p. Further, optimal error estimate in the L ² norm when θ = −1 is derived. Numerical experiments are presented to illustrate the theoretical results. Supported by DST-DAAD (PPP-05) project.     

Multilevel frames for sparse tensor product spaces

For Au = f with an elliptic differential operator and stochastic data f, the m-point correlation function of the random solution u satisfies a deterministic equation with the m-fold tensor product operator A ^(m) of A. Sparse tensor products of hierarchic FE-spaces in are known to allow for approximations to which converge at essentially the rate as in the case m = 1, i.e. for the deterministic problem. They can be realized by wavelet-type FE bases (von Petersdorff and Schwab in Appl Math 51(2):145–180, 2006; Schwab and Todor in Computing 71:43–63, 2003). If wavelet bases are not available, we show here how to achieve the fast computation of sparse approximations of for Galerkin discretizations of A by multilevel frames such as BPX or other multilevel preconditioners of any standard FEM approximation for A. Numerical examples illustrate feasibility and scope of the method.     

Analysis of FETI methods for multiscale PDEs

In this paper, we study a variant of the finite element tearing and interconnecting (FETI) method which is suitable for elliptic PDEs with highly heterogeneous (multiscale) coefficients α(x); in particular, coefficients with strong variation within subdomains and/or jumps that are not aligned with the subdomain interfaces. Using energy minimisation and cut-off arguments we can show rigorously that for an arbitrary (positive) coefficient function the condition number of the preconditioned FETI system can be bounded by C(α) (1 + log(H/h))² where H is the subdomain diameter and h is the mesh size, and where the function C(α) depends only on the coefficient variation in the vicinity of subdomain interfaces. In particular, if varies only mildly in a layer Ω _i,η of width η near the boundary of each of the subdomains Ω _i , then , independent of the variation of α in the remainder Ω _i \Ω _i,η of each subdomain and independent of any jumps of α across subdomain interfaces. The quadratic dependence of C(α) on H/η can be relaxed to a linear dependence under stronger assumptions on the behaviour of α in the interior of the subdomains. Our theoretical findings are confirmed in numerical tests. C. Pechstein was supported by the Austrian Science Funds (FWF) under grant F1306.     

Vectorspacelike representation of absolute planes

The pointset E of an absolute plane can be provided with a binary operation "+" such that (E, +) becomes a loop and for each a E \ {o} the line [a] through o and a is a commutative subgroup of (E, +). Two elements a, b E \ {o} are called independent if [a] ∩ [b] = {o} and the absolute plane is called vectorspacelike if for any two independent elements we have E = [a] + [b] := {x + y | x [a], y [b]}. If is singular then (E, +) is a commutative group and is vectorspacelike iff is Euclidean. If is a hyperbolic plane then is vectorspacelike and in the continous case if a, b are independent, each point p has a unique representation as a quasilinear combination p = α · a + μ · b where α · a [a]and β · b [b] are points, α, β real numbers such that λ (o, λ · a) = |λ|· λ (o, a) and λ (o, μ · b) = |μ|. λ(o, b) and λ is the distance function. This work was partially supported by the Research Project of MIUR (Italian Ministery of Education and University) “Geometria combinatoria e sue applicazioni” and by the research group GNSAGA of INDAM. Dedicated to Walter Benz on the occasion of his 75 ^th birthday, in friendship     

Multigrid methods for the computation of singular solutions and stress intensity factors II: Crack singularities

We consider the Poisson equation −Δu=f with homogeneous Dirichlet boundary condition on a two-dimensional polygonal domain Ω with cracks. Multigrid methods for the computation of singular solutions and stress intensity factors using piecewise linear functions are analyzed. The convergence rate for the stress intensity factors is whenfεL ²(Ω) and whenfεH ¹(Ω). The convergence rate in the energy norm is in the first case and in the second case. The costs of these multigrid methods are proportional to the number of elements in the triangulation. The general case wherefεH ^m (Ω) is also discussed. The work of the first author was partially supported by NSF under grant DMS-96-00133     

The method of mothers for non-overlapping non-matching DDM

In this paper we introduce a variant of the three-field formulation where we use only two sets of variables. Considering, to fix the ideas, the homogeneous Dirichlet problem for in , our variables are i) the approximations of u in each sub-domain (each on its own grid), and ii) an approximation of u on the skeleton (the union of the interfaces of the sub-domains) on an independent grid (that could often be uniform). The novelty is in the way to derive, from , the values of each trace of on the boundary of each . We do it by solving an auxiliary problem on each that resembles the mortar method but is more flexible. Under suitable assumptions, quasi-optimal error estimates are proved, uniformly with respect to the number and size of the subdomains. A preliminary version of the method and of its theoretical analysis has been presented in Bertoluzza et al. (15th international conference on domain decomposition methods, 2002).     

Discrete spectrum in the gaps of a perturbed pseudorelativistic hamiltonian

The pseudorelativistic Hamiltonian is considered under wide conditions on potentials A(x), W(x). It is assumed that a real point λ is regular for G_1/2. Let G_1/2(α)=G_1/2−αV, where α>0, V(x)≥0, and V ∈L _d(ℝ^d). Denote by N(λ, α) the number of eigenvalues of G_1/2(t) that cross the point λ as t increases from 0 to α. A Weyl-type asymptotics is obtained for N(λ, α) as α→∞. Bibliography: 5 titles. To O. A. Ladyzhenskaya Translated fromZapiski Nauchnykh Seminarov POMI, Vol. 249, 1997. pp. 102–117. Translated by A. B. Pushnitskii.     

On the estimation of the regression coefficients of a continuous parameter process with stationary residual

The asymptotic expressions of the covariance matrices for both the least square estimates _L α _T and Markov (best linear) estimates are obtained, based on a sample in a finite interval (0, T) of the regression co-efficients α = (α ₁, …, α _{m 0})′ of a parameter-continuous process with a stationary residual. We assume that the regression variables φ _ν(t), t ⩾ 0, ν = 1, …, m ₀, are continuous in t, and satisfy conditions (3.1)–(3.3). For the residual, we assume that it is a stationary process that possesses a bounded continuous spectral density f(λ). Under these assumptions, it is proven that

On a regularity criterion in terms of the gradient of pressure for the Navier-Stokes equations in \mathbb{R}^{N}

In this paper we establish a Serrin’s type regularity criterion on the gradient of pressure for weak solutions to the Navier–Stokes equations in It is proved that if the gradient of pressure belongs to L^{α, γ} with then the weak solution actually is regular and unique. Received: May 4, 2004     

Bounds and monotonicity for the generalized Robin problem

We consider the principal eigenvalue λ ₁^Ω(α) corresponding to Δu = λ (α) u in on ∂Ω, with α a fixed real, and a C ^0,1 bounded domain. If α > 0 and small, we derive bounds for λ ₁^Ω(α) in terms of a Stekloff-type eigenvalue; while for α > 0 large we study the behavior of its growth in terms of maximum curvature. We analyze how domain monotonicity of the principal eigenvalue depends on the geometry of the domain, and prove that domains which exhibit domain monotonicity for every α are calibrable. We conjecture that a domain has the domain monotonicity property for some α if and only if it is calibrable. Robert Smits: This author was partially supported by a grant of the National Security Agency, grant #H98230-05-1-0060.     

Finite elements approximation of second order linear elliptic equations in divergence form with right-hand side in L 1

In this paper we consider, in dimension d≥ 2, the standard finite elements approximation of the second order linear elliptic equation in divergence form with coefficients in L ^∞(Ω) which generalizes Laplace’s equation. We assume that the family of triangulations is regular and that it satisfies an hypothesis close to the classical hypothesis which implies the discrete maximum principle. When the right-hand side belongs to L ¹(Ω), we prove that the unique solution of the discrete problem converges in (for every q with ) to the unique renormalized solution of the problem. We obtain a weaker result when the right-hand side is a bounded Radon measure. In the case where the dimension is d = 2 or d = 3 and where the coefficients are smooth, we give an error estimate in when the right-hand side belongs to L ^r (Ω) for some r > 1.     

Asymptotics and weighted estimates of meixner polynomials orthogonal on the gird {0, δ, 2δ, ...}

Suppose that 0<δ≤1,N=1/δ, and α, ga≥0, is an integer. For the classical Meixner polynomials orthonormal on the gird {0, δ, 2δ, ...} with weight ρ(x)=(1-e ^−δ)^αг(Nx+α+ 1)/г(Nx+1), the following asymptotic formula is obtained: . The remainderv _n,N ^α (z) forn≤λN satisfies the estimate

Multiplicity of Positive Solutions for a Class of Inhomogeneous Neumann Problems Involving the p(x)-Laplacian

We study the existence and multiplicity of positive solutions for the inhomogeneous Neumann boundary value problems involving the p(x)-Laplacian of the form

Existence and Multiplicity Results for a Degenerate Elliptic Equation

The goal of this paper is to study the multiplicity result of positive solutions of a class of degenerate elliptic equations. On the basis of the mountain pass theorems and the sub- and supersolutions argument for p-Laplacian operators, under suitable conditions on the nonlinearity f（x, s）, we show the following problem：-△pu=λu^α-a（x）u^q in Ω,u│δΩ=0 possesses at least two positive solutions for large λ, where Ω is a bounded open subset of R^N, N ≥ 2, with C^2 boundary, λ is a positive parameter, Ap is the p-Laplacian operator with p 〉 1, α, q are given constants such that p - 1 〈α 〈 q, and a（x） is a continuous positive function in Ω^-.     

Approximation on the limit continuum of a degenerate Kleinian group

Let Γ be a degenerate Kleinian group with limit conlinuum K. Then the linear combinations of the fractions , are dense in C(K) and λ^α(K). Bibliography: 6 titles. Translated fromZapiski Nauchnykh Seminarov POMI, Vol. 206, 1993, pp. 137–150. Translated by S. V. Kislyakov.     

Free boundary problem for the equations of spherically symmetric motion of compressible gas with density-dependent viscosity

We consider a free boundary problem for the equations of spherically symmetric motion of a isentropic gas with a density-dependent viscosity , where and λ are positive constants. We prove that the problem admits a weak solution provided that 0 < λ < 1/4.      

Sharp estimates for intrinsic ultracontractivity on C 1,α-domains

Let be the first Dirichlet eigenfunction on a connected bounded C ^1,α-domain in and the corresponding Dirichlet heat kernel. It is proved that where λ₂ > λ₁ are the first two Dirichlet eigenvalues. This estimate is sharp for both short and long times. Bounded Lipschitz domains, elliptic operators on manifolds, and a general framework are also discussed. Supported in part by Creative Research Group Fund of the National Foundation of China (no. 10121101), the 973-Project in China and RFDP(20040027009).     

$${\mathcal{R}}$$ -boundedness,pseudodifferential operators,and maximal regularity for some classes of partial differential operators

It is shown that an elliptic scattering operator A on a compact manifold with boundary with operator valued coefficients in the morphisms of a bundle of Banach spaces of class () and Pisier’s property (α) has maximal regularity (up to a spectral shift), provided that the spectrum of the principal symbol of A on the scattering cotangent bundle avoids the right half-plane. This is accomplished by representing the resolvent in terms of pseudodifferential operators with -bounded symbols, yielding by an iteration argument the -boundedness of λ(A−λ)⁻¹ in for some . To this end, elements of a symbolic and operator calculus of pseudodifferential operators with -bounded symbols are introduced. The significance of this method for proving maximal regularity results for partial differential operators is underscored by considering also a more elementary situation of anisotropic elliptic operators on with operator valued coefficients.     

Dispersion of Lagrangian trajectories in a random large-scale velocity field

We study the distribution of the distance R(t) between two Lagrangian trajectories in a spatially smooth turbulent velocity field with an arbitrary correlation time and a non-Gaussian distribution. There are two dimensionless parameters, the degree of deviation from the Gaussian distribution α and β=τD, where τ is the velocity correlation time and D is a characteristic velocity gradient. Asymptotically, R(t) has a lognormal distribution characterized by the mean runaway velocity and the dispersion Δ. We use the method of higher space dimensions d to estimate and Δ for different values of α and β. It was shown previously that for β≪ 1 and for β≫ 1. The estimate of Δ is then nonuniversal and depends on details of the two-point velocity correlator. Higher-order velocity correlators give an additional contribution to Δ estimated as αD²τ for β≪1 and αβ/τ for β≫1. For α above some critical value σ_cr, the values of and Δ are determined by higher irreducible correlators of the velocity gradient, and our approach loses its applicability. This critical value can be estimated as for β≪1 and for β≫1. Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 122, No. 3, pp. 456–467, March, 2000.