The local superconvergence of the linear finite element method for the poisson problem

Assume that . In this study, the Richardson extrapolation for the tensor‐product block element and the linear finite element theory of the Green's function will be combined to study the local superconvergence of finite element methods for the Poisson equation in a bounded polytopic domain (polygonal or polyhedral domain for ), where a family of tensor‐product block partitions is not required or the solution need not have high global smoothness. We present a special family of partitions satisfying, for any , e is a tensor‐product block whenever where denotes the distance between e and . By the linear finite element theory of the Green's function and the Richardson extrapolation for the tensor‐product block element, we obtain the local superconvergence of the displacement for the linear finite element method over the special family of partitions . © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 930–946, 2014     

Analysis of expanded mixed finite element methods for the generalized forchheimer flows of slightly compressible fluids

The nonlinear Forchheimer equations are used to describe the dynamics of fluid flows in porous media when Darcy's law is not applicable. In this article, we consider the generalized Forchheimer flows for slightly compressible fluids, and then study the expanded mixed finite element method applied to the initial boundary value problem for the resulting degenerate parabolic equation for pressure. The bounds for the solutions, time derivative, and gradient of solutions are established. Utilizing the monotonicity properties of Forchheimer equation and boundedness of solutions, a priori error estimates for solution are obtained in ‐norm, ‐norm as well as for its gradient in ‐norm for all . Optimal ‐error estimates are shown for solutions under some additional regularity assumptions. Numerical results using the lowest order Raviart–Thomas mixed element confirm the theoretical analysis regarding convergence rates. Published 2015. Numer Methods Partial Differential Eq 32: 60–85, 2016     

A new nonconforming finite element with a conforming finite element approximation for a coupled continuum pipe‐flow/Darcy model in Karst aquifers

Numerical method is considered for a coupled continuum pipe‐flow/Darcy model describing flow in porous media with an embedded conduit pipe. A new nonconforming element is constructed to solve the Darcy equation on porous matrix. The existence and uniqueness of the approximation solution are deduced. Optimal error estimates are obtained in and norms. Some numerical examples show the accuracy and efficiency of the presented method. With the same number of nodal‐points and the same amount of computation, the results using the new nonconforming element are much better than those by both conforming element and Wilson nonconforming element on the same mesh. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 778–798, 2016     

Finite element approximation of a phase field model arising in nanostructure patterning

We consider a fully practical finite element approximation of the nonlinear parabolic Cahn–Hilliard system subject to an initial condition on the conserved order parameter , and mixed boundary conditions. Here, is the interfacial parameter, is the field strength parameter, is the obstacle potential, is the diffusion coefficient, and denotes differentiation with respect to the second argument. Furthermore, w is the chemical potential and is the electrostatic potential. The system, in the context of nanostructure patterning, has been proposed to model the manipulation of morphologies in organic solar cells with the help of an applied electric field. In the limit , it reduces to a sharp interface problem that models the evolution of an unstable interface between two dielectric media in the presence of a quasistatic electric field. On introducing a finite element approximation for the above Cahn–Hilliard system, we prove existence and stability of a discrete solution. Moreover, in the case of two space dimensions, we are able to prove convergence and hence existence of a solution to the considered system of partial differential equations. We demonstrate the practicality of our finite element approximation with several numerical simulations in two and three space dimensions. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1890–1924, 2015     

Rapid Solution of Minimal Riesz Energy Problems

In , , we compute the solution to both the unconstrained and constrained Gauss variational problem, considered for the Riesz kernel of order and a pair of compact, disjoint, boundaryless ‐dimensional ‐manifolds , , where , each being charged with Borel measures with the sign prescribed. Such variational problems over a cone of Borel measures can be formulated as minimization problems over the corresponding cone of surface distributions belonging to the Sobolev–Slobodetski space , where and (see Harbrecht et al., Math. Nachr. 287 (2014), 48–69). We thus approximate the sought density by piecewise constant boundary elements and apply the primal‐dual active set strategy to impose the desired inequality constraints. The boundary integral operator which is defined by the Riesz kernel under consideration is efficiently approximated by means of an ‐matrix approximation. This particularly enables the application of a preconditioner for the iterative solution of the first‐order optimality system. Numerical results in are given to demonstrate our approach. © 2016Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1535–1552, 2016     

Two unconditionally stable and convergent difference schemes with the extrapolation method for the one‐dimensional distributed‐order differential equations

The Grünwald formula is used to solve the one‐dimensional distributed‐order differential equations. Two difference schemes are derived. It is proved that the schemes are unconditionally stable and convergent with the convergence orders and in maximum norm, respectively, where and are step sizes in time, space and distributed order. The extrapolation method is applied to improve the approximate accuracy to the orders and respectively. An illustrative numerical example is given to confirm the theoretical results. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 591–615, 2016     

Rainbow Numbers for Cycles in Plane Triangulations

In the article, the existence of rainbow cycles in edge colored plane triangulations is studied. It is shown that the minimum number of colors that force the existence of a rainbow C₃ in any n‐vertex plane triangulation is equal to . For a lower bound and for an upper bound of the number is determined.     

Lipschitz extensions of definable p‐adic functions

In this paper, we prove a definable version of Kirszbraun's theorem in a non‐Archimedean setting for definable families of functions in one variable. More precisely, we prove that every definable function , where and , that is λ‐Lipschitz in the first variable, extends to a definable function that is λ‐Lipschitz in the first variable.     

Several iterative schemes for the stationary natural convection equations at different Rayleigh numbers

Several iterative schemes based on finite element discretization with triangulation for solving two‐dimensional natural convection equations are studied in this article. We establish some reference points for evaluation of the possible impact from three kinds of schemes with respect to Rayleigh numbers. In case of , all schemes are stable and convergent. Moreover, in case of , Schemes I and II can run well. Finally, in case of , only Scheme I is still stable and convergent. Numerical experiment is presented and discussed for testing of the performances of the proposed schemes, which confirms the theoretic analysis. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 761–776, 2015     

Unconditional convergence and error estimates for bounded numerical solutions of the barotropic Navier–Stokes system

We consider a mixed finite‐volume finite‐element method applied to the Navier–Stokes system of equations describing the motion of a compressible, barotropic, viscous fluid. We show convergence as well as error estimates for the family of numerical solutions on condition that: (a) the underlying physical domain as well as the data are smooth; (b) the time step and the parameter of the spatial discretization are proportional, ; and (c) the family of numerical densities remains bounded for . No a priori smoothness is required for the limit (exact) solution. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1208–1223, 2017     

Computational aspects of satisfiability in probability logic

We consider the complexity of satisfiability in ε‐logic, a probability logic. We show that for the relational fragment this problem is ‐complete for rational , answering a question by Terwijn. In contrast, we show that satisfiability in 0‐logic is decidable. The methods we employ to prove this fact also allow us to show that 0‐logic is compact, while it was previously shown that ε‐logic is not compact for .     

Electromagnetic eigenfields in checkerboard patterns

Here are described four solvers for time‐harmonic electromagnetic fields in checkerboard patterns. A pattern is built by four squares with constant permittivity, or . It is enclosed by conducting walls or is a unit cell of a periodic structure. The field is represented in two ways: by , the transverse component of the magnetic induction, and by , the magnetic vector potential in Lorenz gauge. and satisfy Helmholtz equations in each square as well as transmission and boundary conditions (BCs). These governing equations yield eigensolutions and , which are found to be at worst. Variational versions of the governing equations are introduced. The weak formulations for are standard, while those for are new. They imply that the derivative transmission and BCs are satisfied weakly on interfaces between regions with different permittivity. Eigenpairs are computed approximately by spectral element methods. They yield mutually consistent eigenpairs. However, only about half of the eigenpairs () correspond to eigenpairs (). For each set of BCs, the first few eigenfrequencies are given by tables, and some of the eigenfunctions are presented by contour plots. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 418–444, 2016     

Further Results on Large Sets of Resolvable Idempotent Latin Squares

An idempotent Latin square of order v is called resolvable and denoted by RILS(v) if the off‐diagonal cells can be resolved into disjoint transversals. A large set of resolvable idempotent Latin squares of order v, briefly LRILS(v), is a collection of RILS(v)s pairwise agreeing on only the main diagonal. In this paper, it is established that there exists an LRILS(v) for any positive integer , except for , and except possibly for .     

The Large‐Time Solution of Burgers' Equation with Time‐Dependent Coefficients. I. The Coefficients Are Exponential Functions

In this paper, we consider an initial‐value problem for Burgers' equation with variable coefficients where x and t represent dimensionless distance and time, respectively, and , are given functions of t. In particular, we consider the case when the initial data have algebraic decay as , with as and as . The constant states and are problem parameters. Two specific initial‐value problems are considered. In initial‐value problem 1 we consider the case when and , while in initial‐value problem 2 we consider the case when and . The method of matched asymptotic coordinate expansions is used to obtain the large‐t asymptotic structure of the solution to both initial‐value problems over all parameter values.     

Convergence analysis of a new mixed finite element method for Biot's consolidation model

In this article, we propose a mixed finite element method for the two‐dimensional Biot's consolidation model of poroelasticity. The new mixed formulation presented herein uses the total stress tensor and fluid flux as primary unknown variables as well as the displacement and pore pressure. This method is based on coupling two mixed finite element methods for each subproblem: the standard mixed finite element method for the flow subproblem and the Hellinger–Reissner formulation for the mechanical subproblem. Optimal a‐priori error estimates are proved for both semidiscrete and fully discrete problems when the Raviart–Thomas space for the flow problem and the Arnold–Winther space for the elasticity problem are used. In particular, optimality in the stress, displacement, and pressure has been proved in when the constrained‐specific storage coefficient is strictly positive and in the weaker norm when is nonnegative. We also present some of our numerical results.Copyright © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1189–1210, 2014     

Asymptotically exact a posteriori local discontinuous Galerkin error estimates for the one‐dimensional second‐order wave equation

In this article, we analyze a residual‐based a posteriori error estimates of the spatial errors for the semidiscrete local discontinuous Galerkin (LDG) method applied to the one‐dimensional second‐order wave equation. These error estimates are computationally simple and are obtained by solving a local steady problem with no boundary condition on each element. We apply the optimal L² error estimates and the superconvergence results of Part I of this work [Baccouch, Numer Methods Partial Differential Equations 30 (2014), 862–901] to prove that, for smooth solutions, these a posteriori LDG error estimates for the solution and its spatial derivative, at a fixed time, converge to the true spatial errors in the L²‐norm under mesh refinement. The order of convergence is proved to be , when p‐degree piecewise polynomials with are used. As a consequence, we prove that the LDG method combined with the a posteriori error estimation procedure yields both accurate error estimates and superconvergent solutions. Our computational results show higher convergence rate. We further prove that the global effectivity indices, for both the solution and its derivative, in the L²‐norm converge to unity at rate while numerically they exhibit and rates, respectively. Numerical experiments are shown to validate the theoretical results. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1461–1491, 2015     

A fourth‐order Hermitian box‐scheme with fast solver for the Poisson problem in a cube

We present a fourth‐order Hermitian box‐scheme (HB‐scheme) for the Poisson problem in a cube. A single‐nonstaggered regular grid is used supporting the discrete unknowns u and . The scheme is fourth‐order accurate for u and in norm. The fast numerical resolution uses a matrix capacitance method, resulting in a computational complexity of . Numerical results are reported on several examples including nonseparable problems. The present scheme is the extension to the three‐dimensional case of the HB‐scheme presented in Abbas and Croisille [J Sci Comp 49 (2011), 239–267]. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 609–629, 2015     

Least‐squares spectral method for velocity‐vorticity‐pressure form of the Stokes equations

The aim of this article is to present and analyze first‐order system least‐squares spectral method for the Stokes equations in two‐dimensional spaces. The Stokes equations are transformed into a first‐order system of equations by introducing vorticity as a new variable. The least‐squares functional is then defined by summing up the ‐ and ‐norms of the residual equations. The ‐norm in the least‐squares functional is replaced by suitable operator. Continuous and discrete homogeneous least‐squares functionals are shown to be equivalent to ‐norm of velocity and ‐norm of vorticity and pressure for spectral Galerkin and pseudospectral method. The spectral convergence of the proposed methods are given and the theory is validated by numerical experiment. Mass conservation is also briefly investigated. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 661–680, 2016     

Carleson measures and balayage for Bergman spaces of strongly pseudoconvex domains

Given a bounded strongly pseudoconvex domain D in with smooth boundary, we characterize ‐Bergman Carleson measures for , , and . As an application, we show that the Bergman space version of the balayage of a Bergman Carleson measure on D belongs to BMO in the Kobayashi metric.     

Large time asymptotics of Feynman–Kac functionals for symmetric stable processes

Let μ be a positive Kato measure on associated with the Green kernel of the transient symmetric α‐stable process, the Markov process with generator (). Let be the positive continuous additive functional in the Revuz correspondence with μ. If, in addition, μ is of compact support, we give exact large time asymptotics of the expectation of the Feynman–Kac functional, .