Some abstract error estimates of a finite volume scheme for a nonstationary heat equation on general nonconforming multidimensional spatial meshes

A general class of nonconforming meshes has been recently studied for stationary anisotropic heterogeneous diffusion problems, see Eymard et al. (IMA J. Numer. Anal. 30 (2010), 1009–1043). Thanks to the basic ideas developed in the stated reference for stationary problems, we derive a new discretization scheme in order to approximate the nonstationary heat problem. The unknowns of this scheme are the values at the centre of the control volumes, at some internal interfaces, and at the mesh points of the time discretization. We derive error estimates in discrete norms L ^∞(0, T;H ₀ ¹ (Ω)) and W ^1,∞(0, T;L ²(Ω)), and an error estimate for an approximation of the gradient, in a general framework in which the discrete bilinear form involved in the finite volume scheme satisfies some ellipticity condition.     

On a wave equation with supercritical interior and boundary sources and damping terms

This article addresses nonlinear wave equations with supercritical interior and boundary sources, and subject to interior and boundary damping. The presence of a nonlinear boundary source alone is known to pose a significant difficulty since the linear Neumann problem for the wave equation is not, in general, well‐posed in the finite‐energy space H¹(Ω) × L₂(?Ω) with boundary data in L₂ due to the failure of the uniform Lopatinskii condition. Further challenges stem from the fact that both sources are non‐dissipative and are not locally Lipschitz operators from H¹(Ω) into L₂(Ω), or L₂(?Ω). With some restrictions on the parameters in the model and with careful analysis involving the Nehari Manifold, we obtain global existence of a unique weak solution, and establish exponential and algebraic uniform decay rates of the finite energy (depending on the behavior of the dissipation terms). Moreover, we prove a blow up result for weak solutions with nonnegative initial energy.     

Riesz potential and Sobolev embeddings on generalized Lebesgue and Sobolev spaces Lp(·) and Wk,p(·)

We study the Riesz potentials I_αf on the generalized Lebesgue spaces L^p(·)(?^d), where 0 < α < d and I_αf(x) ? ∫equation/tex2gif-inf-3.gif |f(y)| |x – y|^{α – d} dy. Under the assumptions that p locally satisfies |p(x) – p(x)| ≤ C/(– ln |x – y|) and is constant outside some large ball, we prove that I_α : L^p(·)(?^d) → L^p?(·)(?^d), where . If p is given only on a bounded domain Ω with Lipschitz boundary we show how to extend p to on ?^d such that there exists a bounded linear extension operator ? : W^1,p(·)(Ω) ? (?^d), while the bounds and the continuity condition of p are preserved. As an application of Riesz potentials we prove the optimal Sobolev embeddings W^k,p(·)(?^d) ?L^p*(·)(R^d) with and W^1,p(·)(Ω) ? L^p*(·)(Ω) for k = 1. We show compactness of the embeddings W^1,p(·)(Ω) ? L^q(·)(Ω), whenever q(x) ≤ p*(x) – ε for some ε > 0. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The discrete duality finite volume method for a class of quasi‐Newtonian Stokes flows

In this paper, we propose a discrete duality finite volume (DDFV) scheme for the incompressible quasi‐Newtonian Stokes equation. The DDFV method is based on the use of discrete differential operators which satisfy some duality properties analogous to their continuous counterparts in a discrete sense. The DDFV method has a great ability to handle general geometries and meshes. In addition, every component of the velocity gradient can be reconstructed directly, which makes it suitable to deal with the nonlinear terms in the quasi‐Newtonian Stokes equation. We prove that the proposed DDFV scheme is uniquely solvable and of first‐order convergence in the discrete L²‐norms for the velocity, the strain rate tensor, and the pressure, respectively. Ample numerical tests are provided to highlight the performance of the proposed DDFV scheme and to validate the theoretical error analysis, in particular on locally refined nonconforming and polygonal meshes.     

Stabilized DDFV schemes for stokes problem with variable viscosity on general 2D meshes

“Discrete Duality Finite Volume” schemes (DDFV for short) on general meshes are studied here for Stokes problems with variable viscosity with Dirichlet boundary conditions. The aim of this work is to analyze the well‐posedness of the scheme and its convergence properties. The DDFV method requires a staggered scheme, the discrete unknowns, the components of the velocity and the pressure, are located on different nodes. The scheme is stabilized using a finite volume analogue to Brezzi‐Pitkäranta techniques. This scheme is proved to be well‐posed on general meshes and to be first order convergent in a discrete H¹ ‐norm and a discrete L² ‐norm for respectively the velocity and the pressure. Finally, numerical experiments confirm the theoretical prediction, in particular on locally refined non conformal meshes. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1666–1706, 2011     

The steady states and convergence to equilibria for a 1‐D chemotaxis model with volume‐filling effect

We consider a chemotaxis model with volume‐filling effect introduced by Hillen and Painter. They also proved the existence of global solutions for a compact Riemannian manifold without boundary. Moreover, the existence of a global attractor in W^{1, p}(Ω??ⁿ), p>n, p?2, was proved by Wrzosek. He also proved that the ω‐limit set consists of regular stationary solutions. In this paper, we prove that the 1‐D stationary problem has at most an infinitely countable number of regular solutions. Furthermore, we show that as t→∞ the solution of the 1‐D evolution problem converges to an equilibrium in W^{1, p}, p?2. Copyright © 2009 John Wiley & Sons, Ltd.     

Superconvergence of the direct discontinuous Galerkin method for a time‐fractional initial‐boundary value problem

A time‐fractional reaction–diffusion initial‐boundary value problem with periodic boundary condition is considered on Q ? Ω × [0, T] , where Ω is the interval [0, l] . Typical solutions of such problem have a weak singularity at the initial time t = 0. The numerical method of the paper uses a direct discontinuous Galerkin (DDG) finite element method in space on a uniform mesh, with piecewise polynomials of degree k ≥ 2 . In the temporal direction we use the L1 approximation of the Caputo derivative on a suitably graded mesh. We prove that at each time level of the mesh, our L1‐DDG solution is superconvergent of order k + 2 in L²(Ω) to a particular projection of the exact solution. Moreover, the L1‐DDG solution achieves superconvergence of order (k + 2) in a discrete L²(Q) norm computed at the Lobatto points, and order (k + 1) superconvergence in a discrete H¹(Q) seminorm at the Gauss points; numerical results show that these estimates are sharp.     

A space–time discontinuous Galerkin method for linear convection-dominated Sobolev equations

This article presents a space–time discontinuous Galerkin (DG) finite element method for linear convection-dominated Sobolev equations. The finite element method has basis functions that are continuous in space and discontinuous in time, and variable spatial meshes and time steps are allowed. In the discrete intervals of time, using properties of the Radau quadrature rule, eliminates the restriction to space–time meshes of convectional space–time Galerkin methods. The existence and uniqueness of the approximate solution are proved. An optimal priori error estimate in L^∞(H¹) is derived. Numerical experiments are presented to confirm theoretical results.     

Supraconvergence and Supercloseness of a Scheme for Elliptic Equations on Nonuniform Grids

In this paper, we study the convergence of a finite difference scheme on nonuniform grids for the solution of second-order elliptic equations with mixed derivatives and variable coefficients in polygonal domains subjected to Dirichlet boundary conditions. We show that the scheme is equivalent to a fully discrete linear finite element approximation with quadrature. It exhibits the phenomenon of supraconvergence, more precisely, for s ∈ [1,2] order O(h ^s )-convergence of the finite difference solution, and its gradient is shown if the exact solution is in the Sobolev space H ^1+s (Ω). In the case of an equation with mixed derivatives in a domain containing oblique boundary sections, the convergence order is reduced to O(h ^3/2?ε) with ε > 0 if u ∈ H ³(Ω). The second-order accuracy of the finite difference gradient is in the finite element context nothing else than the supercloseness of the gradient. For s ∈ {1,2}, the given error estimates are strictly local.     

Optimal error estimates of the local discontinuous Galerkin method for nonlinear second‐order elliptic problems on Cartesian grids

In this paper, we study the local discontinuous Galerkin (LDG) methods for two‐dimensional nonlinear second‐order elliptic problems of the type u_xx + u_yy = f(x, y, u, u_x, u_y) , in a rectangular region Ω with classical boundary conditions on the boundary of Ω . Convergence properties for the solution and for the auxiliary variable that approximates its gradient are established. More specifically, we use the duality argument to prove that the errors between the LDG solutions and the exact solutions in the L² norm achieve optimal (p + 1)th order convergence, when tensor product polynomials of degree at most p are used. Moreover, we prove that the gradient of the LDG solution is superclose with order p + 1 toward the gradient of Gauss–Radau projection of the exact solution. The results are valid in two space dimensions on Cartesian meshes using tensor product polynomials of degree p ≥ 1 , and for both mixed Dirichlet–Neumann and periodic boundary conditions. Preliminary numerical experiments indicate that our theoretical findings are optimal.     

Global existence and decay of energy to systems of wave equations with damping and supercritical sources

This paper is concerned with a system of nonlinear wave equations with supercritical interior and boundary sources and subject to interior and boundary damping terms. It is well-known that the presence of a nonlinear boundary source causes significant difficulties since the linear Neumann problem for the single wave equation is not, in general, well-posed in the finite-energy space H ¹(Ω) × L ²(?Ω) with boundary data from L ²(?Ω) (due to the failure of the uniform Lopatinskii condition). Additional challenges stem from the fact that the sources considered in this article are non-dissipative and are not locally Lipschitz from H ¹(Ω) into L ²(Ω) or L ²(?Ω). With some restrictions on the parameters in the system and with careful analysis involving the Nehari Manifold, we obtain global existence of a unique weak solution and establish (depending on the behavior of the dissipation in the system) exponential and algebraic uniform decay rates of energy. Moreover, we prove a blow-up result for weak solutions with nonnegative initial energy.     

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.     

Rapid error reduction for block Gauss–Seidel based on p‐hierarchical basis

We consider a two‐level block Gauss–Seidel iteration for solving systems arising from finite element discretizations employing higher‐order elements. A p‐hierarchical basis is used to induce this block structure. Using superconvergence results normally employed in the analysis of gradient recovery schemes, we argue that a massive reduction of the H¹‐error occurs in the first iterate, so that the discrete solution is adequately resolved in very few iterates—sometimes a single iteration is sufficient. Numerical experiments on uniform and adapted meshes support these claims. Copyright © 2012 John Wiley & Sons, Ltd.     

Ridged Domains,Embedding Theorems and Poincaré Inequalities

We study the embeddings E : W(X(Ω), Y(Ω)) ↪ Z(Ω), where X(Ω), Y(Ω) and Z(Ω) are rearrangement–invariant Banach function spaces (BFS) defined on a generalized ridged domain Ω, and W denotes a first–order Sobolev–type space. We obtain two–sided estimates for the measure of non–compactness of E when Z(Ω) = X(Ω) and, in turn, necessary and sufficient conditions for a Poincaré–type inequality to be valid and also for E to be compact. The results are used to analyse the example of a trumpet–shaped domain Ω in Lorentz spaces. We consider the problem of determining the range of possible target spaces Z(Ω), in which case we prove that the problem is equivalent to an analogue on the generalized ridge Γ of Ω. The range of target spaces Z(Ω) is determined amongst a scale of (weighted) Lebesgue spaces for “rooms and passages” and trumpet–shaped domains.     

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     

Elliptic boundary value problems with λ-dependent boundary conditions

In this paper second order elliptic boundary value problems on bounded domains Ω⊂Rn with boundary conditions on ∂Ω depending nonlinearly on the spectral parameter are investigated in an operator theoretic framework. For a general class of locally meromorphic functions in the boundary condition a solution operator of the boundary value problem is constructed with the help of a linearization procedure. In the special case of rational Nevanlinna or Riesz-Herglotz functions on the boundary the solution operator is obtained in an explicit form in the product Hilbert space L²(Ω)⊕(L²m(∂Ω)), which is a natural generalization of known results on λ-linear elliptic boundary value problems and λ-rational boundary value problems for ordinary second order differential equations.     

K4−‐factor in a graph

Let G be a graph of order 4k and let δ(G) denote the minimum degree of G. Let F be a given connected graph. Suppose that |V(G)| is a multiple of |V(F)|. A spanning subgraph of G is called an F‐factor if its components are all isomorphic to F. In this paper, we prove that if δ(G)≥5/2k, then G contains a K₄^?‐factor (K₄^? is the graph obtained from K₄ by deleting just one edge). The condition on the minimum degree is best possible in a sense. In addition, the proof can be made algorithmic. © 2002 John Wiley & Sons, Inc. J Graph Theory 39: 111–128, 2002     

Boundary blow-up solutions for a class of quasilinear elliptic equations

For a given bounded domain Ω in R ^N with smooth boundary ?Ω, we give sufficient conditions on f so that the m-Laplacian equation △ _m u = f(x, u, ?u) admits a boundary blow-up solution u ∈ W ^1,p (Ω). Our main results are new and extend the results in J.V. Concalves and Angelo Roncalli [Boundary blow-up solutions for a class of elliptic equations on a bounded domain, Appl. Math. Comput. 182 (2006), pp. 13–23]. Our approach employs the method of lower–upper solution theorem, fixed point theory and weak comparison principle.     

On the Uniform Poincaré Inequality

We give a proof of the Poincaré inequality in W ^{1, p} (Ω) with a constant that is independent of Ω ? , where  is a set of uniformly bounded and uniformly Lipschitz domains in ? ⁿ . As a byproduct, we obtain the following: The first non vanishing eigenvalues λ₂(Ω) of the standard Neumann (variational) boundary value problem on Ω for the Laplace operator are bounded below by a positive constant if the domains Ω vary and remain uniformly bounded and uniformly Lipschitz regular.     

Local well‐posedness for a system of quadratic nonlinear Schrödinger equations in one or two dimensions

In this article, the local well‐posedness of Cauchy's problem is explored for a system of quadratic nonlinear Schrödinger equations in the space L^p( R ⁿ). In a special case of mass resonant 2 × 2 system, it is well known that this problem is well posed in H^s(s≥0) and ill posed in H^s(s < 0) in two‐space dimensions. By translation on a linear semigroup, we show that the general system becomes locally well posed in L^p( R ²) for 1 < p < 2, for which p can arbitrarily be close to the scaling limit p_c=1. In one‐dimensional case, we show that the problem is locally well posed in L¹( R ); moreover, it has a measure valued solution if the initial data are a Dirac function. Copyright © 2016 John Wiley & Sons, Ltd.