Local mesh refinement for the discretization of Neumann boundary control problems on polyhedra

This paper deals with optimal control problems constrained by linear elliptic partial differential equations. The case where the right‐hand side of the Neumann boundary is controlled, is studied. The variational discretization concept for these problems is applied, and discretization error estimates are derived. On polyhedral domains, one has to deal with edge and corner singularities, which reduce the convergence rate of the discrete solutions, that is, one cannot expect convergence order two for linear finite elements on quasi‐uniform meshes in general. As a remedy, a local mesh refinement strategy is presented, and a priori bounds for the refinement parameters are derived such that convergence with optimal rate is guaranteed. As a by‐product, finite element error estimates in the H¹(Ω)‐norm, L²(Ω)‐norm and L²(Γ)‐norm for the boundary value problem are obtained, where the latter one turned out to be the main challenge. Copyright © 2015 John Wiley & Sons, Ltd.     

A posteriori error estimates of stabilized finite volume method for the Stokes equations

In this work, the residual‐type posteriori error estimates of stabilized finite volume method are studied for the steady Stokes problem based on two local Gauss integrations. By using the residuals between the source term and numerical solutions, the computable global upper and local lower bounds for the errors of velocity in H¹ norm and pressure in L² norm are derived. Furthermore, a global upper bound of u ? u_h in L²‐norm is also derived. Finally, some numerical experiments are provided to verify the performances of the established error estimators. Copyright © 2015 John Wiley & Sons, Ltd.     

Polynomial approximation on the sphere using scattered data

We consider the problem of approximately reconstructing a function f defined on the surface of the unit sphere in the Euclidean space ℝ^{q +1} by using samples of f at scattered sites. A central role is played by the construction of a new operator for polynomial approximation, which is a uniformly bounded quasi‐projection in the de la Vallée Poussin style, i.e. it reproduces spherical polynomials up to a certain degree and has uniformly bounded L^p operator norm for 1 ≤ p ≤ ∞. Using certain positive quadrature rules for scattered sites due to Mhaskar, Narcowich and Ward, we discretize this operator obtaining a polynomial approximation of the target function which can be computed from scattered data and provides the same approximation degree of the best polynomial approximation. To establish the error estimates we use Marcinkiewicz–Zygmund inequalities, which we derive from our continuous approximating operator. We give concrete bounds for all constants in the Marcinkiewicz–Zygmund inequalities as well as in the error estimates. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Stabilized mixed finite element methods for the Navier‐Stokes equations with damping

In this paper, the stabilized mixed finite element methods are presented for the Navier‐Stokes equations with damping. The existence and uniqueness of the weak solutions are proven by use of the Brouwer fixed‐point theorem. Then, optimal error estimates for the H¹‐norm and L²‐norm of the velocity and the L²‐norm of the pressure are derived. Moreover, on the basis of the optimal L²‐norm error estimate of the velocity, a stabilized two‐step method is proposed, which is more efficient than the usual stabilized methods. Finally, two numerical examples are implemented to confirm the theoretical analysis.     

On the finite element approximation of elliptic QVIs with noncoercive operators

In this paper, we extend the approach developed by the author for the standard finite element method in the L^∞‐norm of the noncoercive variational inequalities (VI) (Numer Funct Anal Optim.2015;36:1107‐1121.) to impulse control quasi‐variational inequality (QVI). We derive the optimal error estimate, combining the so‐called Bensoussan‐Lions Algorithm and the concept of subsolutions for VIs.     

A Posteriori Error Estimates for the Approximations of the Stresses in the Hencky Plasticity Problem

In this article, we derive a posteriori error estimates for the Hencky plasticity problem. These estimates are formulated in terms of the stresses and present guaranteed and computable bounds of the difference between the exact stress field and any approximation of it from the energy space of the dual variational problem. They consist of quantities that can be considered as penalties for the violations of the equilibrium equations, the yield condition and the constitutive relations that must hold for the exact stresses and strains. It is proved that the upper bound tends to zero for any sequence of stresses that tends to the exact solution of the Haar–Karman variational problem. An important ingredient of our analysis is a collection of Poincaré type inequalities involving the L ¹ norms of the tensors of small deformation. Estimates of this form are not new, however we will present computable upper bounds for the constants being involved even for rather complicated domains.     

A posteriori error estimates for finite element approximations to the wave equation with discontinuous coefficients

We derive residual‐based a posteriori error estimates of finite element method for linear wave equation with discontinuous coefficients in a two‐dimensional convex polygonal domain. A posteriori error estimates for both the space‐discrete case and for implicit fully discrete scheme are discussed in L^∞(L²) norm. The main ingredients used in deriving a posteriori estimates are new Clément type interpolation estimates in conjunction with appropriate adaption of the elliptic reconstruction technique of continuous and discrete solutions. We use only an energy argument to establish a posteriori error estimates with optimal order convergence in the L^∞(L²) norm.     

On the best constants in inequalities of the Markov and Wirtinger types for polynomials on the half-line

The main topic of the paper is best constants in Markov-type inequalities between the norms of higher derivatives of polynomials and the norms of the polynomials themselves. The norm is the L² norm with Laguerre weight. The leading term of the asymptotics of the constants is determined and tight bounds for the principal coefficient in this term, which is the operator norm of a Volterra operator, are given. For best constants in inequalities of the Wirtinger type, the limit is computed and an asymptotic formula for the error term is presented.     

Discontinuous Galerkin finite element methods for variational inequalities of first and second kinds

We develop the error analysis for the h‐version of the discontinuous Galerkin finite element discretization for variational inequalities of first and second kinds. We establish an a priori error estimate for the method which is of optimal order in a mesh dependant as well as L²‐norm.© 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2007     

A posteriori error estimates of edge residual type of finite element method for nonmonotone quasi‐linear elliptic problems

In this article, we study the edge residual‐based a posteriori error estimates of conforming linear finite element method for nonmonotone quasi‐linear elliptic problems. It is proven that edge residuals dominate a posteriori error estimates. Up to higher order perturbations, edge residuals can act as a posteriori error estimators. The global reliability and local efficiency bounds are established both in H ¹‐norm and L ²‐norm. Numerical experiments are provided to illustrate the performance of the proposed error estimators. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 813–837, 2014     

Error bounds and variational methods for nonlinear boundary value problems

Complementary variational principles are used to deriveL ² andL error bounds for approximate solutions of nonlinear boundary value problems for ordinary differential equations.     

Sharp Integral Inequalities of the Hermite–Hadamard Type

We consider a family of two-point quadrature formulae and establish sharp estimates for the remainders under various regularity conditions. Improved forms of certain integral inequalities due to Hermite and Hadamard, Iyengar, Milovanovi and Pecari , and others are obtained as special cases. Our results can also be interpreted as analogues to a theorem of Ostrowski on the deviation of a function from its averages. Furthermore, we establish a generalization of a result of Fink concerning L^p estimates for the remainder of the trapezoidal rule and present the best constants in the error bounds.     

Estimates of deviations from exact solutions of variational inequalities based upon Payne–Weinberger inequality

A new method for obtaining computable estimates for the difference between exact solutions of elliptic variational inequalities and arbitrary functions in the respective energy space is suggested. The estimates are obtained by transforming the corresponding variational inequality without the use of variational duality arguments. These estimates are valid for any function in the energy class and contain no constants depending on the mesh used to find an approximate solution. This method for linear elliptic and parabolic problems was earlier suggested by the author. The guaranteed error bounds we derive can be of two types. Estimates of the first type contain only one global constant, which is a constant in the Friedrichs type inequality. Estimates of the second type are based on the decomposition of Ω into convex subdomains and the Payne–Weinberger inequalities for these subdomains. Bibliography: 20 titles. Translated from Problems in Mathematical Analysis 39 February, 2009, pp. 81–90.     

Parallel Galerkin domain decomposition procedures for parabolic equation on general domain

Parallel Galerkin domain decomposition procedures for parabolic equation on general domain are given. These procedures use implicit Galerkin method in the subdomains and simple explicit flux calculation on the interdomain boundaries by integral mean method or extrapolation method to predict the inner‐boundary conditions. Thus, the parallelism can be achieved by these procedures. These procedures are conservative both in the subdomains and across interboundaries. The explicit nature of the flux prediction induces a time‐step limitation that is necessary to preserve stability, but this constraint is less severe than that for a fully explicit method. L²‐norm error estimates are derived for these procedures. Compared with the work of Dawson and Dupont [Math Comp 58 (1992), 21–35], these L²‐norm error estimates avoid the loss of H^?1/2 factor. Experimental results are presented to confirm the theoretical results. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

Local discontinuous galerkin approximation of convection‐dominated diffusion optimal control problems with control constraints

In this article, we investigate local discontinuous Galerkin approximation of stationary convection‐dominated diffusion optimal control problems with distributed control constraints. The state variable and adjoint state variable are approximated by piecewise linear polynomials without continuity requirement, whereas the control variable is discretized by variational discretization concept. The discrete first‐order optimality condition is derived. We show that optimization and discretization are commutative for the local discontinuous Galerkin approximation. Because the solutions to convection‐dominated diffusion equations often admit interior or boundary layers, residual type a posteriori error estimate in L² norm is proved, which can be used to guide mesh refinement. Finally, numerical examples are presented to illustrate the theoretical findings. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 339–360, 2014     

Analysis of finite element methods for second order boundary value problems using mesh dependent norms

This paper presents a new approach to the analysis of finite element methods based onC ⁰-finite elements for the approximate solution of 2nd order boundary value problems in which error estimates are derived directly in terms of two mesh dependent norms that are closely ralated to theL ₂ norm and to the 2nd order Sobolev norm, respectively, and in which there is no assumption of quasi-uniformity on the mesh family. This is in contrast to the usual analysis in which error estimates are first derived in the 1st order Sobolev norm and subsequently are derived in theL ₂ norm and in the 2nd order Sobolev norm — the 2nd order Sobolev norm estimates being obtained under the assumption that the functions in the underlying approximating subspaces lie in the 2nd order Sobolev space and that the mesh family is quasi-uniform.     

Semidiscrete Galerkin method for equations of motion arising in Kelvin‐Voigt model of viscoelastic fluid flow

Finite element Galerkin method is applied to equations of motion arising in the Kelvin–Voigt model of viscoelastic fluids for spatial discretization. Some new a priori bounds which reflect the exponential decay property are obtained for the exact solution. For optimal L^∞( L ²) estimate in the velocity, a new auxiliary operator which is based on a modification of the Stokes operator is introduced and analyzed. Finally, optimal error bounds for the velocity in L^∞( L ²) as well as in L^∞( H )‐norms and the pressure in L^∞(L²)‐norm are derived which again preserves the exponential decay property. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Propagation of and Maxwellian weighted bounds for derivatives of solutions to the homogeneous elastic Boltzmann equation

We consider the n-dimensional space homogeneous Boltzmann equation for elastic collisions for variable hard potentials with Grad (angular) cutoff. We prove sharp moment inequalities, the propagation of L¹-Maxwellian weighted estimates, and consequently, the propagation L^∞-Maxwellian weighted estimates to all derivatives of the initial value problem associated to the afore mentioned equation. More specifically, we extend to all derivatives of the initial value problem associated to this class of Boltzmann equations corresponding sharp moment (Povzner) inequalities and time propagation of L¹-Maxwellian weighted estimates as originally developed Bobylev [A.V. Bobylev, Moment inequalities for the Boltzmann equation and applications to spatially homogeneous problems, J. Statist. Phys. 88 (1997) 1183–1214] in the case of hard spheres in 3 dimensions. To achieve this goal we implement the program presented in Bobylev–Gamba–Panferov [A.V. Bobylev, I.M. Gamba, V. Panferov, Moment inequalities and high-energy tails for Boltzmann equation with inelastic interactions, J. Statist. Phys. 116 (5–6) (2004) 1651–1682], which includes a full analysis of the moments by means of sharp moment inequalities and the control of L¹-exponential bounds, in the case of stationary states for different inelastic Boltzmann related problems with ‘heating’ sources where high energy tail decay rates depend on the inelasticity coefficient and the type of ‘heating’ source. More recently, this work was extended to variable hard potentials with angular cutoff by Gamba–Panferov–Villani [I.M. Gamba, V. Panferov, C. Villani, Upper Maxwellian bounds for the spatially homogeneous Boltzmann equation, ARMA (2008), in press] in the elastic case collision case where the L¹-Maxwellian weighted norm was shown to propagate if initial states have such property. In addition, we also extend to all derivatives the propagation of L^∞-Maxwellian weighted estimates, proven in [I.M. Gamba, V. Panferov, C. Villani, Upper Maxwellian bounds for the spatially homogeneous Boltzmann equation, ARMA (2008), in press], to solutions of the initial value problem to the Boltzmann equations for elastic collisions for variable hard potentials with Grad (angular) cutoff.     

Global Well‐Posedness of the Three‐Dimensional Primitive Equations with Only Horizontal Viscosity and Diffusion

In this paper, we consider the initial boundary value problem of the three‐dimensional primitive equations for planetary oceanic and atmospheric dynamics with only horizontal eddy viscosity in the horizontal momentum equations and only horizontal diffusion in the temperature equation. Global well‐posedness of the strong solution is established for any H² initial data. An N‐dimensional logarithmic Sobolev embedding inequality, which bounds the L^∞‐norm in terms of the L^q‐norms up to a logarithm of the L^p‐norm for p > N of the first‐order derivatives, and a system version of the classic Grönwall inequality are exploited to establish the required a~priori H² estimates for global regularity.© 2016 Wiley Periodicals, Inc.     

Fully discrete interior penalty discontinuous Galerkin methods for nonlinear parabolic equations

In this article, we investigate interior penalty discontinuous Galerkin (IPDG) methods for solving a class of two‐dimensional nonlinear parabolic equations. For semi‐discrete IPDG schemes on a quasi‐uniform family of meshes, we obtain a priori bounds on solutions measured in the L² norm and in the broken Sobolev norm. The fully discrete IPDG schemes considered are based on the approximation by forward Euler difference in time and broken Sobolev space. Under a restriction related to the mesh size and time step, an hp ‐version of an a priori l^∞(L²) and l²(H¹) error estimate is derived and numerical experiments are presented.© 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 288–311, 2012