Finite Element Approximation of a Rigid Punch Indenting a Membrane

Optimal order H¹ and L error bounds are obtained for a continuouspiecewise linear finite element approximation of an obstacleproblem, where the obstacle's height as well as the contactzone, _c, are a priori unknown. The problem models the indentationof a membrane by a rigid punch. For R², given ,g R⁺ and an obstacle defined over E we consider the minimization of |v|²_1,+gµover (v, µ) H¹₀() x R subject to v+µ on E. In additionwe show under certain nondegeneracy conditions that dist (_c,^h_c)Ch ln 1/h, where ^h_c is the finite element approximation to_c. Finally we show that the resulting algebraic problem canbe solved using a projected SOR algorithm.     

Total Flux Estimates for a Finite-Element Approximation of Elliptic Equations

An elliptic boundary-value problem on a domain with prescribedDirichlet data on _I is approximated using a finite-elementspace of approximation power h^K in the L² norm. It is shownthat the total flux across _I can be approximated with an errorof O(h^K) when is a curved domain in Rⁿ (n = 2 or 3) and isoparametricelements are used. When is a polyhedron, an O(h^2K–2)approximation is given. We use these results to study the finite-elementapproximation of elliptic equations when the prescribed boundarydata on _I is the total flux. Present address: School of Mathematical and Physical Sciences,University of Sussex, Brighton, Sussex BN1 9QH.     

Finite-Element Approximation of a Plasma Equilibrium Problem

The plasma problem studied is: given R⁺ find (, d, u) R ?R ? H¹() such that Let ₁ < ₂ be the first two eigenvalues of the associatedlinear eigenvalue problem: find $$\left(\lambda ,\phi \right)\in\mathrm{R;}\times {\hbox{ H }}_{0}^{1}\left(\Omega \right)$$such that For ₀(0,₂) it is well known that there exists a unique solution(₀, d₀, u₀) to the above problem. We show that the standard continuous piecewise linear Galerkinfinite-element approximatinon $$\left({\lambda }_{0},{\hbox{d }}_{0}^{k},{u}_{0}^{h}\right)$$, for ₀(0,₂), converges atthe optimal rate in the H¹, L², and L norms as h, the mesh length,tends to 0. In addition, we show that dist (, ^h)Ch² ln 1/h,where $${\Gamma }^{\left(h\right)}=\left\{x\in \Omega :{u}_{0}^{\left(h\right)}\left(x\right)=0\right\}$$.Finally we consider a more practical approximation involvingnumerical integration.     

The P-version Penalty Finite Element Method

A p-version penalty finite element method is used to solve themodel problem –u=f in , u=g on . Error estimates are derivedin H¹-norm. The p-version penalty method with extrapolationyields an approximate solution which converges at the optimalrate. Numerical results show the effectiveness of the p-versionpenalty method with extrapolation.     

A note on least-squares mixed finite elements in relation to standard and mixed finite elements

^** Email: brandts{at}science.uva.nl The least-squares mixed finite-element method for second-orderelliptic problems yields an approximation u_h V_h H₀¹() of thepotential u together with an approximation p_{h h} H(div ; )of the vector field p = – Au. Comparing u_h with the standardfinite-element approximation of u in V_h, and p_h with the mixedfinite-element approximation of p, it turns out that they arehigher-order perturbations of each other. In other words, theyare ‘superclose’. Refined a priori bounds and superconvergenceresults can now be proved. Also, the local mass conservationerror is of higher order than could be concluded from the standarda priori analysis.     

The Model Plasma Problem: Existence of a Solution Branch

We are interested in the model plasma problem –u = u⁺in ,u = –d on , _au⁺ dx=j where is a bounded domain in with boundary ; here, j isa given positive number, the function u and the positive number are the unknowns of the problem, and d is a real parameter.Using a variant of the implicit function theorem, we can provethe existence of a global solution branch parametrized by d.The method has the advantage that it can be used for analysingthe approximation of the above problem by a finite-element method.     

Discontinuous Galerkin finite element approximation of quasilinear elliptic boundary value problems I: the scalar case

^** Email: Paul.Houston{at}mcs.le.ac.uk^*** Email: Janice.Robson{at}comlab.ox.ac.uk^**** Email: Endre.Suli{at}comlab.ox.ac.uk We develop a one-parameter family of hp-version discontinuousGalerkin finite element methods, parameterised by [–1,1], for the numerical solution of quasilinear elliptic equationsin divergence form on a bounded open set ^d, d 2. In particular,we consider the analysis of the family for the equation –·{µ(x, |u|)u} = f(x) subject to mixed Dirichlet–Neumannboundary conditions on . It is assumed that µ is a real-valuedfunction, µ C( x [0, )), and thereexist positive constants m_µ and M_µ such that m_µ(t– s) µ(x, t)t – µ(x, s)s M_µ(t– s) for t s 0 and all x . Using a result from the theory of monotone operators for any valueof [–1, 1], the corresponding method is shown to havea unique solution u_DG in the finite element space. If u C¹() H^k(), k 2, then with discontinuous piecewise polynomials ofdegree p 1, the error between u and u_DG, measured in the brokenH¹()-norm, is (h^s–1/p^k–3/2), where 1 s min {p+ 1, k}.     

Finite-Element Approximation of Elliptic Equations with a Neumann or Robin Condition on a Curved Boundary

This paper considers a finite-element approximation of a second-orderself adjoint elliptic equation in a region Rⁿ (with n=2 or 3)having a curved boundary on which a Neumann or Robin conditionis prescribed. If the finite-element space defined over , a union of elements, has approximation power h^kin the L² norm, and if the region of integration is approximatedby ^h with dist (, ^h)Ch^k, then it is shown that one retains optimalrates of convergence for the error in the H¹ and L² norms, whetherQ^h is fitted or unfitted , provided that the numerical integration scheme has sufficientaccuracy.     

Rate of convergence of numerical approximations to homoclinic bifurcation points

For x=f (x, ), x Rⁿ, R, having a hyperbolic or semihyperbolicequilibrium p(), we study the numerical approximation of parametervalues ^* at which there is an orbit homoclinic to p(). We approximate^* by solving a finite-interval boundary value problem on J=[T_–,T₊], T_–<0<T₊, with boundary conditions that sayx(T_–) and x(T₊) are in approximations to appropriate invariantmanifolds of p(). A phase condition is also necessary to makethe solution unique. Using a lemma of Xiao-Biao Lin, we improve,for certain phase conditions, existing estimates on the rateof convergence of the computed homoclinic bifurcation parametervalue , to the true value ^*. The estimates we obtain agree withthe rates of convergence observed in numerical experiments.Unfortunately, the phase condition most commonly used in numericalwork is not covered by our results.     

On the Behaviour of the First Eigenfunction of the p-Laplacian Near its Critical Points

In this paper, the behaviour of the positive eigenfunction of in u_| = 0, p > 1, isstudied near its critical points. Under some convexity and symmetryassumptions on , is seen to have a unique critical point atx = 0; also, the behaviour of both and is determined nearby.Positive solutions u to some general problems –_pu = f(u)in , u_| = 0, are also considered, with some convexity restrictionson u. 2000 Mathematics Subject Classification 35B05 (primary),35J65, 35J70 (secondary).     

Finite Element Methods for a Model for Full Waveform Acoustic Logging

A model is defined to simulate the propagation of waves in aradially symmetric, isotropic, composite system consisting ofa fluid-filled well bore ^f through a fluid-saturated poroussolid ^p. Biot's equations of motion are chosen to describe thepropagation of waves in ^p, while the standard equation of motionfor compressible inviscid fluids is used for ^f, with appropriateboundary conditions at the contact surface between ^f and ^p.Also, absorbing boundary conditions for the artificial boundariesof ^p are derived for the model, their effect being to make themtransparent for waves arriving normally First, results on the existence and uniqueness of the solutionof the differential problem are given and then a discrete-time,explicit finite element procedure is defined and analysed, withfinite element spaces suited for radially symmetric problemsbeing used for the spatial discretisation.     

Approximation of Ground State Eigenvalues and Eigenfunctions of Dirichlet Laplacians

Let be a bounded connected open set in R^N, N 2, and let –0be the Dirichlet Laplacian defined in L²(). Let > 0 be thesmallest eigenvalue of –, and let > 0 be its correspondingeigenfunction, normalized by ||||₂ = 1. For sufficiently small>0 we let R() be a connected open subset of satisfying Let – 0 be the Dirichlet Laplacian on R(), and let >0and >0 be its ground state eigenvalue and ground state eigenfunction,respectively, normalized by ||||₂=1. For functions f definedon , we let Sf denote the restriction of f to R(). For functionsg defined on R(), we let Tg be the extension of g to satisfying 1991 Mathematics SubjectClassification 47F05.     

Central L-Values and Toric Periods for GL(2)

Let be a cusp form on GL(2) over a number field F and let Ebe a quadratic extension of F. Denote by _E the base change of to E and by a unitary character of A^x_E/ E^x. We use the relativetrace formula to give an explicit formula for L(1/2, _E ) interms of period integrals of Gross–Prasad test vectors.We give an application of this formula to equidistribution ofgeodesics on a hyperbolic 3-fold.     

Discrete methods in the study of an inverse problem for Laplace's equation

Let u be harmonic in the interior of a rectangle and satisfythe third-kind boundary condition u_n + yu = where 0, y 0with supports included in the bottom and in the top side of, respectively. Recovering y from a knowledge of and of thetrace of u on the bottom is a nonlinear inverse problem ofinterest in the field of nondestructive evaluation. A convergentGalerkin method for approximating y is proposed and tested innumerical experiments.     

Logarithmic Convexity for Supremum Norms of Harmonic Functions

We prove the following convexity property for supremum normsof harmonic functions. Let be a domain in Rⁿ, ₀ and E a subdomainand a compact sebset of ,respectively. Then there exists a constant = (E, ₀, ) (0, 1) such that for all harmonic functions u on, the inequality is valid.The case of concentric balls E plays a key role in the proof.For positive harmonic funcitons ono osuch balls, we determinethe sharp constant in the inequlity.     

On hearing the shape of a bounded domain with Robin boundary conditions

The asymptotic expansions of the trace of the heat kernel (t)= [sum ]_j=1 exp(-t_j) for small positive t, where {_j} _j=1 arethe eigenvalues of the negative Laplacian -_n = -[sum ]ⁿ_k=1 (/x^k)²in Rⁿ (n = 2 or 3), are studied for a general multiply connectedbounded domain which is surrounded by simply connected boundeddomains _i with smooth boundaries _i (i = 1,...,m), where smoothfunctions Y_i (i = 1,...,m) are assuming the Robin boundary conditions(n_i + Y_i) = 0 on _i. Here /n_i denote differentiations along theinward-pointing normals to _i (i = 1,...,m). Some applicationsof an ideal gas enclosed in the multiply connected bounded containerwith Neumann or Robin boundary conditions are given.     

On Restrictions and Extensions of the Besov and Triebel-Lizorkin Spaces with Respect to Lipschitz Domains

The restrictions B^s_pq() and F^s_pq() of the Besov and Triebel–Lizorkinspaces of tempered distributions B^s_pq(Rⁿ) and F^s_pq(Rⁿ) to Lipschitzdomains Rⁿ are studied. For general values of parameters (sR,p>0, q>0) a ‘universal’ linear bounded extensionoperator from B^s_pq() and F^s_pq() into the corresponding spaceson Rⁿ is constructed. The construction is based on a new variantof the Calderón reproducing formula with kernels supportedin a fixed cone. Explicit characterizations of the elementsof B^s_pq() and F^s_pq() in terms of their values in are also obtained.     

Convergence rates for sparse chaos approximations of elliptic problems with stochastic coefficients

^** Email: todor{at}math.ethz.ch^*** Corresponding author. Email: schwab{at}math.ethz.ch A scalar, elliptic boundary-value problem in divergence formwith stochastic diffusion coefficient a(x, ) in a bounded domainD ^d is reformulated as a deterministic, infinite-dimensional,parametric problem by separation of deterministic (x D) andstochastic ( ) variables in a(x, ) via Karhúnen–Loèveor Legendre expansions of the diffusion coefficient. Deterministic,approximate solvers are obtained by projection of this probleminto a product probability space of finite dimension M and sparsediscretizations of the resulting M-dimensional parametric problem.Both Galerkin and collocation approximations are considered.Under regularity assumptions on the fluctuation of a(x, ) inthe deterministic variable x, the convergence rate of the deterministicsolution algorithm is analysed in terms of the number N of deterministicproblems to be solved as both the chaos dimension M and themultiresolution level of the sparse discretization resp. thepolynomial degree of the chaos expansion increase simultaneously.     

A bifurcation problem governed by the boundary condition II

In this work we consider the problem u = a(x)u^p in on , where is a smooth bounded domain, isthe outward unit normal to , is regarded as a parameter and0 < p < 1. We consider both cases where a(x) > 0 in or a(x) is allowed to vanish in a whole subdomain ₀ of . Ourmain results include existence of non-negative non-trivial solutionsin the range 0 < < ₁, where ₁ is characterized by meansof an eigenvalue problem, uniqueness and bifurcation from infinityof such solutions for small , and the appearance of dead coresfor large enough .     

Base Sizes and Regular Orbits for Coprime Affine Permutation Groups

Let G be a permutation group on a finite set . A sequence B=(₁,..., _b) of points in is called a base if its pointwise stabilizerin G is the identity. Bases are of fundamental importance incomputational algorithms for permutation groups. For both practicaland theoretical reasons, one is interested in the minimal basesize for (G, ), For a nonredundant base B, the elementary inequality2^|B||G|||^|B| holds; in particular, |B|log|G|/log||. In the casewhen G is primitive on , Pyber [8, p. 207] has conjectured thatthe minimal base size is less than Clog|G|/log|| for some (large)universal constant C. It appears that the hardest case of Pyber's conjecture is thatof primitive affine groups. Let H=GV be a primitive affine group;here the point stabilizer G acts faithfully and irreduciblyon the elementary abelian regular normal subgroup V of H, andwe may assume that =V. For positive integers m, let mV denotethe direct sum of m copies of V. If (v₁, ..., v_m)mV belongsto a regular G-orbit, then (0, v₁, ..., v_m) is a base for theprimitive affine group H. Conversely, a base (₁, ..., _b) forH which contains 0V= gives rise to a regular G-orbit on (b–1)V. Thus Pyber's conjecture for affine groups can be viewed asa regular orbit problem for G-modules, and it is therefore aspecial case of an important problem in group representationtheory. For a related result on regular orbits for quasisimplegroups, see [4, Theorem 6].