On the numerical approximation of p‐biharmonic and ∞‐biharmonic functions

The ∞ ‐Bilaplacian is a third‐order fully nonlinear PDE given by (1) In this work, we build a numerical method aimed at quantifying the nature of solutions to this problem, which we call ∞ ‐biharmonic functions. For fixed p we design a mixed finite element scheme for the prelimiting equation, the p‐Bilaplacian (2) We prove convergence of the numerical solution to the weak solution of and show that we are able to pass to the limit p → ∞ . We perform various tests aimed at understanding the nature of solutions of and we prove convergence of our discretization to an appropriate weak solution concept of this problem that of ‐solutions.     

Fast high order difference schemes for the time fractional telegraph equation

In this paper, a fast high order difference scheme is first proposed to solve the time fractional telegraph equation based on the ℱℒ 2-1_σ formula for the Caputo fractional derivative, which reduces the storage and computational cost for calculation. A compact scheme is then presented to improve the convergence order in space. The unconditional stability and convergence in maximum norm are proved for both schemes, with the accuracy order and , respectively. Difficulty arising from the two Caputo fractional derivatives is overcome by some detailed analysis. Finally, we carry out numerical experiments to show the efficiency and accuracy, by comparing with the ℒ 2-1_σ method.     

Analysis and efficient implementation of alternating direction implicit finite volume method for Riesz space‐fractional diffusion equations in two space dimensions

In this article, we develop a Crank–Nicolson alternating direction implicit finite volume method for time‐dependent Riesz space‐fractional diffusion equation in two space dimensions. Norm‐based stability and convergence analysis are given to show that the developed method is unconditionally stable and of second‐order accuracy both in space and time. Furthermore, we develop a lossless matrix‐free fast conjugate gradient method for the implementation of the numerical scheme, which only has memory requirement and computational complexity per iteration with N being the total number of spatial unknowns. Several numerical experiments are presented to demonstrate the effectiveness and efficiency of the proposed scheme for large‐scale modeling and simulations.     

Uniformly convergent scheme for two‐parameter singularly perturbed problems with non‐smooth data

A numerical scheme is constructed for the problems in which the diffusion and convection parameters (?₁ and ?₂ , respectively) both are small, and the convection and source terms have a jump discontinuity in the domain of consideration. Depending on the magnitude of the ratios , and two different cases have been considered separately. Through rigorous analysis, the theoretical error bounds on the singular and regular components of the solution are obtained separately, which shows that in both cases the method is convergent uniformly irrespective of the size of the parameters ?₁, ?₂ . Two test problems are included to validate the theoretical results.     

Ultraconvergence of finite element method by Richardson extrapolation for elliptic problems with inhomogeneous boundary conditions

In this article, Richardson extrapolation technique is employed to investigate the local ultraconvergence properties of Lagrange finite element method using piecewise polynomials of degrees () for the second order elliptic problem with inhomogeneous boundary. A sequence of special graded partition are proposed and a new interpolation operator is introduced to achieve order local ultraconvergence for the displacement and derivative.     

H1-superconvergence of finite difference method based on Q1-element on quasi-uniform mesh for the 3D Poisson equation

In this paper, the finite difference (FD) method is considered for the 3D Poisson equation by using the Q₁-element on a quasi-uniform mesh. First, under the regularity assumption of , the H¹-superconvergence of the FD solution u_h based on the Q₁-element to the first-order interpolation function is obtained. Next, the H¹-superconvergence of the second-order interpolation postprocessing function based on the FD solution u_h to u is provided. Finally, numerical tests are presented to show the H¹-superconvergence result of the FD postprocessing function to u if .     

Half‐inverse problems for the quadratic pencil of the Sturm–Liouville equations with impulse

In this paper, we consider the inverse spectral problem for the impulsive Sturm–Liouville differential pencils on [0, π] with the Robin boundary conditions and the jump conditions at the point . We prove that two potentials functions on the whole interval and the parameters in the boundary and jump conditions can be determined from a set of eigenvalues for two cases: (i) the potentials given on and (ii) the potentials given on , where 0 < α < 1 , respectively. Inverse spectral problems, Sturm–Liouville operator, spectrum, uniqueness.     

A robust scheme based on novel-operational matrices for some classes of time-fractional nonlinear problems arising in mechanics and mathematical physics

In this paper, we present a novel approach based on shifted Gegenbauer wavelets to attain approximate solutions of some classed of time-fractional nonlinear problems. First, we present the approximation of a function of two variables u(x,t) with help of shifted Gegenbauer wavelets and then some novel operational matrices are proposed with the help of piecewise functions to investigate the positive integer derivative ( D _x and D _t), fractional-order derivative ( and ), fractional-order integration ( and ) and delay terms ( and ) of approximated function u(x,t). In order to transform the discussed nonlinear problem into linear problem Picard iterative scheme has been adopt. The current scheme converts the discussed highly nonlinear time-fractional problem into system of linear algebraic equation the help of developed operational matrices and Picard idea. Analysis on the error bound and convergence to authenticate the mathematical formulation of the computational algorithm. We solve various test problems, such as the van der Pol oscillator model, generalized Burger–Huxley, neutral delay parabolic differential equations, sine-Gordon, parabolic integro-differential equation and nonlinear Schrödinger equations to show the efficiency and accuracy of linearized shifted Gegenbauer wavelets method. A comprehensive comparative examination shows the credibility, accuracy, and reliability of the presently proposed computational approach. Also, this scheme can be extended conveniently to other multi-dimensional physical problems of highly nonlinear fractional or variable order of complex nature.     

Superconvergence analysis of an energy stable scheme with three step backward differential formula-finite element method for nonlinear reaction–diffusion equation

A three step backward differential formula scheme is proposed for nonlinear reaction–diffusion equation and superconvergence results are studied with Galerkin finite element method unconditionally. Energy stability is testified for the constructed scheme with an artificial term. Splitting technique is utilized to get rid of the ratio between the time step size and the subdivision parameter . Temporal error estimate in H²-norm is derived, which leads to the boundedness of the solutions of the time-discrete equations. Unconditional spatial error estimate in L²-norm is deduced which help bound the numerical solutions in L^∞-norm. Superconvergent property of in H¹-norm with order is obtained by taking difference between two time levels of the error equations unconditionally. The global superconvergent property is deduced through the above results. Two numerical examples show the validity of the theoretical analysis.     

Numerical analysis for a system coupling curve evolution attached orthogonally to a fixed boundary,to a reaction–diffusion equation on the curve

We consider a semidiscrete finite element approximation for a system consisting of the evolution of a planar curve evolving by forced curve shortening flow inside a given bounded domain , such that the curve meets the boundary orthogonally, and the forcing is a function of the solution of a reaction–diffusion equation that holds on the evolving curve. We prove optimal order error bounds for the resulting approximation and present numerical experiments.     

Gray codes extending quadratic matchings

Is it true that every matching in the n-dimensional hypercube can be extended to a Gray code? More than two decades have passed since Ruskey and Savage asked this question and the problem still remains open. A solution is known only in some special cases, including perfect matchings or matchings of linear size. This article shows that the answer to the Ruskey–Savage problem is affirmative for every matching of size at most . The proof is based on an inductive construction that extends balanced matchings in the completion of the hypercube by edges of into a Hamilton cycle of . On the other hand, we show that for every there is a balanced matching in of size that cannot be extended in this way.     

A mixed hp FEM for the approximation of fourth‐order singularly perturbed problems on smooth domains

We consider fourth‐order singularly perturbed problems posed on smooth domains and the approximation of their solution by a mixed Finite Element Method on the so‐called Spectral Boundary Layer Mesh. We show that the method converges uniformly, with respect to the singular perturbation parameter, at an exponential rate when the error is measured in the energy norm. Numerical examples illustrate our theoretical findings.     

Well‐posedness of fractional integro‐differential equations in vector‐valued functional spaces

We study the well‐posedness of the fractional differential equations with infinite delay on Lebesgue–Bochner spaces and Besov spaces , where A and B are closed linear operators on a Banach space X satisfying ,  and . Under suitable assumptions on the kernels a and b, we completely characterize the well‐posedness of in the above vector‐valued function spaces on by using known operator‐valued Fourier multiplier theorems. We also give concrete examples where our abstract results may be applied.     

The local ultraconvergence of high‐order finite element method for second‐order elliptic problems with constant coefficients over a rectangular partition

In this article, we will discuss the local ultraconvergence of high‐degree finite element method based on a rectangular partition for the second‐degree elliptic problem with constant coefficients in Ω ? ?² , u( y ) = 0 on ?Ω . Based on suitable regularity, ultraconvergence of the displacement of the extrapolated kth (k ≥ 3) degree finite element solution has been obtained by an extrapolation technique. Finally, numerical experiments are applied to demonstrate our theoretical findings.     

Swendsen‐Wang algorithm on the mean‐field Potts model

We study the q‐state ferromagnetic Potts model on the n‐vertex complete graph known as the mean‐field (Curie‐Weiss) model. We analyze the Swendsen‐Wang algorithm which is a Markov chain that utilizes the random cluster representation for the ferromagnetic Potts model to recolor large sets of vertices in one step and potentially overcomes obstacles that inhibit single‐site Glauber dynamics. Long et al. studied the case q = 2, the Swendsen‐Wang algorithm for the mean‐field ferromagnetic Ising model, and showed that the mixing time satisfies: (i) for , (ii) for , (iii) for , where β_c is the critical temperature for the ordered/disordered phase transition. In contrast, for there are two critical temperatures that are relevant. We prove that the mixing time of the Swendsen‐Wang algorithm for the ferromagnetic Potts model on the n‐vertex complete graph satisfies: (i) for , (ii) for , (iii) for , and (iv) for . These results complement refined results of Cuff et al. on the mixing time of the Glauber dynamics for the ferromagnetic Potts model.     

Numerical study of a conservative weighted compact difference scheme for the symmetric regularized long wave equations

In this article, a new weighted and compact conservative difference scheme for the symmetric regularized long wave (SRLW) equations is considered. The new scheme is decoupled and linearized in practical computation, that is, at each time step only two tridiagonal systems of linear algebraic equations need to be solved. It is proved by the discrete energy method that the compact scheme is uniquely solvable, the convergence and stability of the difference scheme are obtained, and its numerical convergence order is in the ‐norm. Numerical experiment results show that the scheme is efficient and reliable.     

An efficient difference scheme for the coupled nonlinear fractional Ginzburg–Landau equations with the fractional Laplacian

In this article, an efficient difference scheme for the coupled fractional Ginzburg–Landau equations with the fractional Laplacian is studied. We construct the discrete scheme based on the implicit midpoint method in time and a weighted and shifted Grünwald difference method in space. Then, we prove that the scheme is uniquely solvable, and the numerical solutions are bounded and unconditionally convergent in the norm. Finally, numerical tests are given to confirm the theoretical results and show the effectiveness of the scheme.     

Numerical method for generalized time fractional KdV-type equation

In this article, an efficient numerical method for linearized and nonlinear generalized time-fractional KdV-type equations is proposed by combining the finite difference scheme and Petrov–Galerkin spectral method. The scale and weight functions involved in generalized fractional derivative cause too much difficulty in discretization and numerical analysis. Fortunately, motivated by finite difference method for fractional differential equation on graded mesh, the stability and convergence of the constructed method are established rigorously. It is proved that the full discretization schemes of generalized time-fractional KdV-type equation is unconditionally stable in linear case. While for nonlinear case, it is stable under a CFL condition and for not small ϵ, coefficient of the high-order spatial differential term. In addition, the full discretization schemes with respect to linear and nonlinear cases respectively converge to the associated exact solutions with orders and , where τ, α, N and m accordingly indicate the time step size, the order of the fractional derivative, polynomial degree, and regularity of the exact solution. Numerical experiments are carried out to support the theoretical results.     

Decomposition of integral self‐affine multi‐tiles

In contrast to the situation with self‐affine tiles, the representation of self‐affine multi‐tiles may not be unique (for a fixed dilation matrix). Let be an integral self‐affine multi‐tile associated with an integral, expansive matrix B and let K tile by translates of . In this work, we propose a stepwise method to decompose K into measure disjoint pieces  satisfying in such a way that the collection of sets forms an integral self‐affine collection associated with the matrix B and this with a minimum number of pieces . When used on a given measurable subset K which tiles by translates of , this decomposition terminates after finitely many steps if and only if the set K is an integral self‐affine multi‐tile. Furthermore, we show that the minimal decomposition we provide is unique.     

Fractional powers of the noncommutative Fourier's law by the S‐spectrum approach

Let e_?, for ? = 1,2,3, be orthogonal unit vectors in and let be a bounded open set with smooth boundary ?Ω. Denoting by a point in Ω, the heat equation, for nonhomogeneous materials, is obtained replacing the Fourier law, given by the following: into the conservation of energy law, here a, b, are given functions. With the S‐spectrum approach to fractional diffusion processes we determine, in a suitable way, the fractional powers of T. Then, roughly speaking, we replace the fractional powers of T into the conservation of energy law to obtain the fractional evolution equation. This method is important for nonhomogeneous materials where the Fourier law is not simply the negative gradient. In this paper, we determine under which conditions on the coefficients a, b, the fractional powers of T exist in the sense of the S‐spectrum approach. More in general, this theory allows to compute the fractional powers of vector operators that arise in different fields of science and technology. This paper is devoted to researchers working in fractional diffusion and fractional evolution problems, partial differential equations, and noncommutative operator theory.