Finite volume transport schemes

We analyze arbitrary order linear finite volume transport schemes and show asymptotic stability in L ¹ and L ^∞ for odd order schemes in dimension one. It gives sharp fractional order estimates of convergence for BV solutions. It shows odd order finite volume advection schemes are better than even order finite volume schemes. Therefore the Gibbs phenomena is controlled for odd order finite volume schemes. Numerical experiments sustain the theoretical analysis. In particular the oscillations of the Lax–Wendroff scheme for small Courant numbers are correlated with its non stability in L ¹. A scheme of order three is proved to be stable in L ¹ and L ^∞.     

On the quadrature error in operational quadrature methods for convolutions

Summary We analyze the quadrature error associated with operational quadrature methods for convolution equations. The assumptions are that the convolution kernel is inL ¹ and that its Laplace transform is analytic and bounded in an obtuse sector of the complex plane. Under these circumstances the Laplace transform has a slow variation property which admits a Fourier analysis of the quadrature error. We provide generalL ^p error estimates assuming suitable smoothness conditions on the function under convolution.     

The finite element method with nodes moving along the characteristics for convection-diffusion equations

Summary We examine the convergence properties of the finite element method with nodes moving along the characteristics for one-dimensional convection-diffusion equations. For linear elements, we demonstrate optimal rates of convergence in theL ²,H ¹ andL norms. Both linear and nonlinear problems are considered.This work forms part of the research programme of the Oxford/Reading Institute for Computational Dynamics.     

Analysis of a BDF–DGFE scheme for nonlinear convection–diffusion problems

We deal with the numerical solution of a scalar nonstationary nonlinear convection–diffusion equation. We employ a combination of the discontinuous Galerkin finite element method for the space semi-discretization and the k-step backward difference formula for the time discretization. The diffusive and stabilization terms are treated implicitly whereas the nonlinear convective term is treated by a higher order explicit extrapolation method, which leads to the necessity to solve only a linear algebraic problem at each time step. We analyse this scheme and derive a priori asymptotic error estimates in the discrete L ^∞(L ²)-norm and the L ²(H ¹)-seminorm with respect to the mesh size h and time step τ for k = 2,3. Numerical examples verifying the theoretical results are presented. This work is a part of the research project MSM 0021620839 financed by the Ministry of Education of the Czech Republic and was partly supported by the Grant No. 316/2006/B-MAT/MFF of the Grant Agency of the Charles University Prague. The research of M. Vlasák was supported by the project LC06052 of the Ministry of Education of the Czech Republic (Jindřich Nečas Center for Mathematical Modelling).     

On hypersurfaces inR n+1 with integral bounds on curvature

We show that the L ^p norm of the second fundamental form of hypersurfaces in R ⁿ⁺¹ is very coercive in the GMT setting of Gauss graphs currents, when p exceeds the dimension n. A compactness result for immersed hypersurfaces and its application to a variational problem are provided.     

ELLIPTIC H^2-VOLTERRA PROJECTION AND THE H^1-GALERKIN METHODS FOR THE INTEGRO-DIFFERENTIAL EQUATIONS OF EVOLUTION

This paper studies H^1-Galerkin methods for the integro-differential equations of evolution. The elliptic H^2-Volterra projection is induced and then used in the derivations of error estimates for semi-discrete and full-discrete H^1-Galerkin methods.The optimal L^2, H^1 and H^2 norm error estimates are obtained.     

Besov regularity and new error estimates for finite volume approximations of the p-laplacian

Summary. In [1], we have constructed a family of finite volume schemes on rectangular meshes for the p-laplacian and we proved error estimates in case the exact solution lies in W^2,p. Actually, W^2,p is not a natural space for solutions of the p-laplacian in the case p>2. Indeed, for general L^p data it can be shown that the solution only belongs to the Besov space In this paper, we prove Besov kind a priori estimates on the approximate solution for any data in L^p. We then obtain new error estimates for such solutions in the case of uniform meshes     

A θ-stable approximations of abstract Cauchy problems

Summary We study the approximation of linear parabolic Cauchy problems by means of Galerkin methods in space andA -stable multistep schemes of arbitrary order in time. The error is evaluated in the norm ofL _t ² (H _x ¹ ) L _t (L _x ² ).     

A posteriori error estimates for elliptic problems with Dirac delta source terms

The aim of this paper is to introduce residual type a posteriori error estimators for a Poisson problem with a Dirac delta source term, in L ^p norm and W^1,p seminorm. The estimators are proved to yield global upper and local lower bounds for the corresponding norms of the error. They are used to guide adaptive procedures, which are experimentally shown to lead to optimal orders of convergence.     

Cubic Spline Prewavelets on the Four-Directional Mesh

Dedicated to Professor M. J. D. Powell on the occasion of his sixty-fifth birthday and his retirement. In this paper, we design differentiable, two-dimensional, piecewise polynomial cubic prewavelets of particularly small compact support. They are given in closed form, and provide stable, orthogonal decompositions of L ² (R ² ) . In particular, the splines we use in our prewavelet constructions give rise to stable bases of spline spaces that contain all cubic polynomials, whereas the more familiar box spline constructions cannot reproduce all cubic polynomials, unless resorting to a box spline of higher polynomial degree.     

Operator-splitting galerkin methods for one kind of oin reaction model for the pollution in groundwater

Abstract. The operator-splitting methods for the mathematic model of one kind of oin reactionsfor the problem of groundwater are considered. Optimal error estimates in Lz and Hl norm areobtained for the approximation solution.     

L p estimates of boundary integral equations for some nonlinear boundary value problems

Summary Recently, Galerkin and collocation methods have been analyzed for boundary integral equation formulations of some potential problems in the plane with nonlinear boundary conditions, and stability results and error estimates in theH ^1/2-norm have been proved (Ruotsalainen and Wendland, and Ruotsalainen and Saranen). We show that these results extend toL ^p setting without any extra conditions. These extensions are proved by studying the uniform boundedness of the inverses of the linearized integral operators, and then considering the nonlinear equations. The fact that inH ^1/2 setting the nonlinear operator is a homeomorphism with Lipschitz continuous inverse plays a crucial role. Optimal error estimates for the Galerkin and collocation method inL ^p space then follow.This research was performed while the second author was visiting professor at the University of Delaware, spring 1989     

Error estimates for the numerical solution of a time-periodic linear parabolic problem

In this paper we consider the numerical solution of a time-periodic linear parabolic problem. We derive optimal order error estimates inL ₂() for approximate solutions obtained by discretizing in space by a Galerkin finite-element method and in time by single-step finite difference methods, using known estimates for the associated initial value problem. We generalize this approach and obtain error estimates for more general discretization methods in the norm of a Banach spaceB L ₂(), e.g.,B=H ₀ ¹ () orL (). Finally, we consider some computational aspects and give a numerical example.     

Walsh Multipliers for Dyadic Hardy Spaces

In this article we are interested in conditions on the coefficients of a Walsh multiplier operator that imply the operator is bounded on certain dyadic Hardy spaces H ^p , 0 < p < ∞. In particular, we consider two classical coefficient conditions, originally introduced for the trigonometric case, the Marcinkiewicz and the Hörmander–Mihlin conditions. They are known to be sufficient in the spaces L ^p , 1 < p < ∞. Here we study the corresponding problem on dyadic Hardy spaces, and find the values of p for which these conditions are sufficient. Then, we consider the cases of H ¹ and L ¹ which are of special interest. Finally, based on a recent integrability condition for Walsh series, a new condition is provided that implies that the multiplier operator is bounded from L ¹ to L ¹, and from H ¹ to H ¹. We note that existing multiplier theorems for Hardy spaces give growth conditions on the dyadic blocks of the Walsh series of the kernel, but these growth are not computable directly in terms of the coefficients.     

Mappings of finite distortion: the degree of regularity

This paper investigates the self-improving integrability properties of the so-called mappings of finite distortion. Let K(x)1 be a measurable function defined on a domain ΩRⁿ, n2, and such that exp(βK(x))L_loc¹(Ω), β>0. We show that there exist two universal constants c₁(n),c₂(n) with the following property: Let f be a mapping in W_loc^1,1(Ω,Rⁿ) with |Df(x)|ⁿK(x)J(x,f) for a.e. xΩ and such that the Jacobian determinant J(x,f) is locally in L¹ log^−c₁(n)βL. Then automatically J(x,f) is locally in L¹ log^c₂(n)βL(Ω). This result constitutes the appropriate analog for the self-improving regularity of quasiregular mappings and clarifies many other interesting properties of mappings of finite distortion. Namely, we obtain novel results on the size of removable singularities for bounded mappings of finite distortion, and on the area distortion under this class of mappings.     

Galerkin methods for nonlinear Sobolev equations

Summary We study Galerkin approximations to the solution of nonlinear Sobolev equations with homogeneous Dirichlet boundary condition in two spatial dimensions and derive optimalL ² error estimates for the continuous Crank — Nicolson and Extrapolated Crank — Nicolson approximations.     

Invariant measure for a class of parabolic degenerate equations

In this paper we consider a stochastic flow in Rⁿ which leaves a closed convex set K invariant. By using a recent characterization of the invariance, involving the distance function, we study the corresponding transition semigroup P_t and its infinitesimal generator N. Due to the invariance property, N is a degenerate elliptic operator. We study existence of an invariant measure ν of P_t and the realization of N in L² (H, ν).     

Two composition methods for solving certain systems of linear equations

Summary This note is concerned with the following problem: Given a systemA·x=b of linear equations and knowing that certains of its subsystemsA ¹·x ¹=b ¹, ...,A ^m ·x ^m =b ^m can be solved uniquely what can be said about the regularity ofA and how to find the solutionx fromx ¹, ...,x ^m ? This question is of particular interest for establishing methods computing certain linear or quasilinear sequence transformations recursively [7, 13, 15].Work performed under NATO Research Grant 027-81     

Optimal L∞-error estimates for nonconforming and mixed finite element methods of lowest order

Summary For second order linear elliptic problems, it is proved that theP ₁-nonconforming finite element method has the sameL -asymptotic accuracy as theP ₁-conforming one. This result is applied to derive optimalL -error estimates for both the displacement and the stress fields of the lowest order Raviart-Thomas mixed finite element method, and a superconvergence result at the barycenter of each element.Performed in the research program of Istituto di Analisi Numerica of C.N.R. of PaviaPartially supported by MPI, GNIM of CNR, ItalySupported by Consejo Nacional de Investigaciones Cientificas y Técnicas, Argentina     

Finite element analysis of exponentially fitted Lumped schemes for time-dependent convection-diffusion problems

Summary In the paper we consider a singularly perturbed linear parabolic initialboundary value problem in one space variable. Two exponential fitted schemes are derived for the problem using Petrov-Galerkin finite element methods with various choices of trial and test spaces. On rectangular meshes which are either arbitrary or slightly restricted, we derive global energy norm andL ² norm and localL error bounds which are uniform in the diffusion parameter. Numerical results are also persented.