Methods for the numerical solution of the Benjamin‐Bona‐Mahony‐Burgers equation

L^∞‐error estimates for finite element for Galerkin solutions for the Benjamin‐Bona‐Mahony‐Burgers (BBMB) equation are considered. A priori bound and the semidiscrete Galerkin scheme are studied using appropriate projections. For fully discrete Galerkin schemes, we consider the backward Euler method and analyze the corresponding error estimates. For a second order accuracy in time, we propose a three‐level backward method. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Fully discrete second-order backward difference method for Kelvin-Voigt fluid flow model

In this article, based on a second-order backward difference method, a completely discrete scheme is discussed for a Kelvin-Voigt viscoelastic fluid flow model with nonzero forcing function, which is either independent of time or in L ^∞ (L ²). After deriving some a priori bounds for the solution of a semidiscrete Galerkin finite element scheme, a second-order backward difference method is applied for temporal discretization. Then, a priori estimates in Dirichlet norm are derived, which are valid uniformly in time using a combination of discrete Gronwall’s lemma and Stolz-Cesaro’s classical result on sequences. Moreover, an existence of a discrete global attractor for the discrete problem is established. Further, optimal a priori error estimates are obtained, whose bounds may depend exponentially in time. Under uniqueness condition, these estimates are shown to be uniform in time. Finally, several numerical experiments are conducted to confirm our theoretical findings.     

A posteriori error estimates for H 1-Galerkin mixed finite-element method for parabolic problems

The purpose of this article is to derive a posteriori error estimates for the H ¹-Galerkin mixed finite element method for parabolic problems. We study both semidiscrete and fully discrete a posteriori error analyses using standard energy argument. A fully discrete a posteriori error analysis based on the backward Euler method is analysed and upper bounds for the errors are derived. The estimators yield upper bounds for the errors which are global in space and time. Our analysis is based on residual approach and the estimators are free from edge residuals.     

and Norms Error Estimates in Finite Element Method for Linear Parabolic Interface Problems

We derive new a priori error estimates for linear parabolic equations with discontinuous coefficients. Due to low global regularity of the solutions the error analysis of the standard finite element method for parabolic problems is difficult to adopt for parabolic interface problems. A finite element procedure is, therefore, proposed and analyzed in this paper. We are able to show that the standard energy technique of finite element method for non-interface parabolic problems can be extended to parabolic interface problems if we allow interface triangles to be curved triangles. Optimal pointwise-in-time error estimates in the L ²(Ω) and H ¹(Ω) norms are shown to hold for the semidiscrete scheme. A fully discrete scheme based on backward Euler method is analyzed and pointwise-in-time error estimates are derived. The interfaces are assumed to be arbitrary shape but smooth for our purpose.     

On Kirchhoff's Model of Parabolic Type

In this article, the existence of a global strong solution for all finite time is derived for the Kirchhoff's model of parabolic type. Based on exponential weight function, some new regularity results which reflect the exponential decay property are obtained for the exact solution. For the related dynamics, the existence of a global attractor is shown to hold for the problem when the non-homogeneous forcing function is either independent of time or in L^∞(L²). With the finite element Galerkin method applied in spatial direction keeping time variable continuous, a semidiscrete scheme is analyzed, and it is also established that the semidiscrete system has a global discrete attractor. Optimal error estimates in L^∞(H¹) norm are derived which are valid uniformly in time. Further, based on a backward Euler method, a completely discrete scheme is analyzed and error estimates are derived. It is also further, observed that in cases where f = 0 or f = O(e^?γ₀t) with γ₀ > 0, the discrete solutions and error estimates decay exponentially in time. Finally, some numerical experiments are discussed which confirm our theoretical findings.     

Higher order fully discrete scheme combined withH 1-Galerkin mixed finite element method for semilinear reaction-diffusion equations

We first apply a first order splitting to a semilinear reaction-diffusion equation and then discretize the resulting system by anH ¹-Galerkin mixed finite element method in space. This semidiscrete method yields a system of differential algebraic equations (DAEs) ofindex one. Apriori error estimates for semidiscrete scheme are derived for both differential as well as algebraic components. For fully discretization, an implicit Runge-Kutta (IRK) methods is applied to the temporal direction and the error estimates are discussed for both components. Finally, we conclude the paper with a numerical example.     

A finite element method for the Sivashinsky equation

The Sivashinsky equation is a nonlinear evolutionary equation of fourth order in space. In this paper we have analyzed a semidiscrete finite element method and completely discrete scheme based on the backward Euler method and Crank–Nicolson–Galerkin scheme. A linearized backward Euler method have been developed and error bounds are derived for an L² projection.     

Discontinuous finite volume element method for parabolic problems

In this article, we consider the semidiscrete and the backward Euler fully discrete discontinuous finite volume element methods for the second‐order parabolic problems and obtain the optimal order error estimates in a mesh dependent norm and in the L²‐norm. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Finite element analysis of parabolic integro-differential equations of Kirchhoff type

The aim of this paper is to study parabolic integro-differential equations of Kirchhoff type. We prove the existence and uniqueness of the solution for this problem via Galerkin method. Semidiscrete formulation for this problem is presented using conforming finite element method. As a consequence of the Ritz–Volterra projection, we derive error estimates for both semidiscrete solution and its time derivative. To find the numerical solution of this class of equations, we develop two different types of numerical schemes, which are based on backward Euler–Galerkin method and Crank–Nicolson–Galerkin method. A priori bounds and convergence estimates in spatial as well as temporal direction of the proposed schemes are established. Finally, we conclude this work by implementing some numerical experiments to confirm our theoretical results.     

Finite element methods for semilinear parabolic interface problems

In this exposition we study the finite element methods for second-order semilinear parabolic interface problems in two dimensional convex polygonal domains with smooth interface. Both semidiscrete and fully discrete schemes are analyzed. Optimal order error estimates in the L²(0, T; H¹(Ω))-norm are established. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Low‐curvature image simplifiers: Global regularity of smooth solutions and Laplacian limiting schemes

We consider a class of fourth‐order nonlinear diffusion equations motivated by Tumblin and Turk's “low‐curvature image simplifiers” for image denoising and segmentation. The PDE for the image intensity u is of the form where g(s) = k²/(k² + s²) is a “curvature” threshold and λ denotes a fidelity‐matching parameter. We derive a priori bounds for Δu that allow us to prove global regularity of smooth solutions in one space dimension, and a geometric constraint for finite‐time singularities from smooth initial data in two space dimensions. This is in sharp contrast to the second‐order Perona‐Malik equation (an ill‐posed problem), on which the original LCIS method is modeled. The estimates also allow us to design a finite difference scheme that satisfies discrete versions of the estimates, in particular, a priori bounds on the smoothness estimator in both one and two space dimensions. We present computational results that show the effectiveness of such algorithms. Our results are connected to recent results for fourth‐order lubrication‐type equations and the design of positivity‐preserving schemes for such equations. This connection also has relevance for other related fourth‐order imaging equations. © 2004 Wiley Periodicals, Inc.     

Finite element approximation of nonlocal parabolic problem

In this article, we present a finite element scheme combined with backward Euler method to solve a nonlocal parabolic problem. An important issue in the numerical solution of nonlocal problems while using Newton's method is related to its structure. In fact differently from the local case where the Jacobian matrix is sparse and banded, in the nonlocal case the Jacobian matrix is dense and computations are much more onerous compared to that for differential equations. In order to avoid this difficulty, we use the technique given by Gudi (SIAM J Numer Anal 50 (2012), 657–668) for elliptic nonlocal problem of Kirchhoff type. We discuss the well‐posedness of the weak formulation at continuous as well as at discrete levels. We also derive a priori error estimates for semidiscrete and fully discrete formulations in L² and H¹ norms. Results based on the usual finite element method are provided to confirm the theoretical estimates. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 786–813, 2017     

An Alternate Approach to Optimal L 2-Error Analysis of Semidiscrete Galerkin Methods for Linear Parabolic Problems with Nonsmooth Initial Data

In this article, we propose and analyze an alternate proof of a priori error estimates for semidiscrete Galerkin approximations to a general second order linear parabolic initial and boundary value problem with rough initial data. Our analysis is based on energy arguments without using parabolic duality. Further, it follows the spirit of the proof technique used for deriving optimal error estimates for finite element approximations to parabolic problems with smooth initial data and hence, it unifies both theories, that is, one for smooth initial data and other for nonsmooth data. Moreover, the proposed technique is also extended to a semidiscrete mixed method for linear parabolic problems. In both cases, optimal L ²-error estimates are derived, when the initial data is in L ². A superconvergence phenomenon is also observed, which is then used to prove L ^∞-estimates for linear parabolic problems defined on two-dimensional spatial domain again with rough initial data.     

W m,2 -Fehlerabschätzungen für Galerkin-Näherungen schwach elliptischer quasilinearer Randwertprobleme 2m-ter Ordnung

Summary For weak elliptic quasilinear boundary value problems of order 2m inn dimensionsW ^m,2 -error estimates for the Galerkin method are established, in whichL -norms of certain derivatives of the Galerkin approximations still occur. The order of these derivatives depends on several conditions on the coefficients of the differential operator. With the help of appropriate a priori bounds for the discrete solutions asymptotic error estimates for the finite element method may be obtained from this. This procedure yields quasioptimal results in several cases. Finally some examples are discussed.

     

Sharp estimates for finite element approximations to parabolic problems with Neumann boundary data of low regularity

Consider a homogeneous parabolic problem on a smooth bounded domain in ℝ ^N but with initial data and Neumann boundary data of low regularity. Sharp interior maximum norm error estimates are given for a semidiscrete C ⁰ finite element approximation to this problem. These estimates are obtained by first establishing a new localized L ^∞ estimate for semidiscrete finite element approximations on interior subdomains. Numerical examples illustrate the findings. AMS subject classification (2000) 65N30     

Error estimates of expanded mixed methods for optimal control problems governed by hyperbolic integro‐differential equations

In this article, we investigate the L^∞(L²) ‐error estimates of the semidiscrete expanded mixed finite element methods for quadratic optimal control problems governed by hyperbolic integrodifferential equations. The state and the costate are discretized by the order k Raviart‐Thomas mixed finite element spaces, and the control is approximated by piecewise polynomials of order k(k ≥ 0). We derive error estimates for both the state and the control approximation. Numerical experiments are presented to test the theoretical results. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

Semidiscrete and asymptotic approximations for the nonstationary radiative–conductive heat transfer problem in a periodic system of grey heat shields

We consider semidiscrete and asymptotic approximations to a solution to the nonstationary nonlinear initial-boundary-value problem governing the radiative–conductive heat transfer in a periodic system consisting of n grey parallel plate heat shields of width ε = 1/n, separated by vacuum interlayers. We study properties of special semidiscrete and homogenized problems whose solutions approximate the solution to the problem under consideration. We establish the unique solvability of the problem and deduce a priori estimates for the solutions. We obtain error estimates of order O( ?{e} ) O\left( {\sqrt {\varepsilon } } \right) and O(ε) for semidiscrete approximations and error estimates of order O( ?{e} ) O\left( {\sqrt {\varepsilon } } \right) and O(ε ^3/4) for asymptotic approximations. Bibliography: 9 titles.     

Mortar element methods for parabolic problems

In this article a standard mortar finite element method and a mortar element method with Lagrange multiplier are used for spatial discretization of a class of parabolic initial‐boundary value problems. Optimal error estimates in L^∞(L²) and L^∞(H¹)‐norms for semidiscrete methods for both the cases are established. The key feature that we have adopted here is to introduce a modified elliptic projection. In the standard mortar element method, a completely discrete scheme using backward Euler scheme is discussed and optimal error estimates are derived. The results of numerical experiments support the theoretical results obtained in this article. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

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 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