Source author record

Waixiang Cao

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

Catalog footprint

Research overview

8works
2topics
4close collaborators

Published work

8 publications

preprint2022arXiv

An LpL^p- Primal-Dual Weak Galerkin method for div-curl Systems

This paper presents a new LpL^p-primal-dual weak Galerkin (PDWG) finite element method for the div-curl system with the normal boundary condition for p>1p>1. Two crucial features for the proposed LpL^p-PDWG finite element scheme are as follows: (1) it offers an accurate and reliable numerical solution to the div-curl system under the low Wα,pW^{α, p}-regularity (α>0α>0) assumption for the exact solution; (2) it offers an effective approximation of the normal harmonic vector fields on domains with complex topology. An optimal order error estimate is established in the LqL^q-norm for the primal variable where 1p+1q=1\frac{1}{p}+\frac{1}{q}=1. A series of numerical experiments are presented to demonstrate the performance of the proposed LpL^p-PDWG algorithm.

preprint2021arXiv

Unconditionally optimal convergence of an energy-conserving and linearly implicit scheme for nonlinear wave equations

In this paper, we present and analyze an energy-conserving and linearly implicit scheme for solving the nonlinear wave equations. Optimal error estimates in time and superconvergent error estimates in space are established without time-step dependent on the spatial mesh size. The key is to estimate directly the solution bounds in the H2H^2-norm for both the nonlinear wave equation and the corresponding fully discrete scheme, while the previous investigations rely on the temporal-spatial error splitting approach. Numerical examples are presented to confirm energy-conserving properties, unconditional convergence, and optimal error estimates, respectively, of the proposed fully discrete schemes.

preprint2020arXiv

A C1C^1 Petrov-Galerkin method and Gauss collocation method for 1D general elliptic problems and superconvergence

In this paper, we present and study C1C^1 Petrov-Galerkin and Gauss collocation methods with arbitrary polynomial degree kk (≥3\ge 3) for one-dimensional elliptic equations. We prove that, the solution and its derivative approximations converge with rate 2k−22k-2 at all grid points; and the solution approximation is superconvergent at all interior roots of a special Jacobi polynomial of degree k+1k+1 in each element, the first-order derivative approximation is superconvergent at all interior k−2k-2 Lobatto points, and the second-order derivative approximation is superconvergent at k−1k-1 Gauss points, with an order of k+2k+2, k+1k+1, and kk, respectively. As a by-product, we prove that both the Petrov-Galerkin solution and the Gauss collocation solution are superconvergent towards a particular Jacobi projection of the exact solution in H2H^2, H1H^1, and L2L^2 norms. All theoretical findings are confirmed by numerical experiments.

preprint2018arXiv

An Accurate and Efficient Algorithm for The Time-fractional Molecular Beam Epitaxy Model with Slope Selection

In this paper, we propose a time-fractional molecular beam epitaxy (MBE) model with slope selection and its efficient, accurate, full discrete, linear numerical approximation. The numerical scheme utilizes the fast algorithm for the Caputo fractional derivative operator in time discretization and Fourier spectral method in spatial discretization. Refinement tests are conducted to verify the 2−α2-α order of time convergence, with α∈(0,1]α\in (0, 1] the fractional order of derivative. Several numerical simulations are presented to demonstrate the accuracy and efficiency of our newly proposed scheme. By exploring the fast algorithm calculating the Caputo fractional derivative, our numerical scheme makes it practice for long time simulation of MBE coarsening, which is essential for MBE model in practice. With the proposed fractional MBE model, we observe that the scaling law for the energy decays as O(t−α3) O(t^{-\fracα{3}}) and the roughness increases as O(tα3)O(t^{\fracα{3}}), during the coarsening dynamics with random initial condition. That is to say, the coarsening rate of MBE model could be manipulated by the fractional order αα, and it is linearly proportional to αα. This is the first time in literature to report/discover such scaling correlation. It provides a potential application field for fractional differential equations. Besides, the numerical approximation strategy proposed in this paper can be readily applied to study many classes of time-fractional and high dimensional phase field models.

preprint2014arXiv

Is 2k2k-Conjecture valid for finite volume methods?

This paper is concerned with superconvergence properties of a class of finite volume methods of arbitrary order over rectangular meshes. Our main result is to prove {\it 2k-conjecture}: at each vertex of the underlying rectangular mesh, the bi-kk degree finite volume solution approximates the exact solution with an order O(h2k) O(h^{2k}), where hh is the mesh size. As byproducts, superconvergence properties for finite volume discretization errors at Lobatto and Gauss points are also obtained. All theoretical findings are confirmed by numerical experiments.

preprint2014arXiv

Superconvergence of Local Discontinuous Galerkin method for one-dimensional linear parabolic equations

In this paper, we study superconvergence properties of the local discontinuous Galerkin method for one-dimensional linear parabolic equations when alternating fluxes are used. We prove, for any polynomial degree kk, that the numerical fluxes converge at a rate of 2k+12k+1 (or 2k+1/22k+1/2) for all mesh nodes and the domain average under some suitable initial discretization. We further prove a k+1k+1th superconvergence rate for the derivative approximation and a k+2k+2th superconvergence rate for the function value approximation at the Radau points. Numerical experiments demonstrate that in most cases, our error estimates are optimal, i.e., the error bounds are sharp.

preprint2012arXiv

Any order superconvergence finite volume schemes for 1D general elliptic equations

We present and analyze a finite volume scheme of arbitrary order for elliptic equations in the one-dimensional setting. In this scheme, the control volumes are constructed by using the Gauss points in subintervals of the underlying mesh. We provide a unified proof for the inf-sup condition, and show that our finite volume scheme has optimal convergence rate under the energy and L2L^2 norms of the approximate error. Furthermore, we prove that the derivative error is superconvergent at all Gauss points and in some special case, the convergence rate can reach h2rh^{2r}, where rr is the polynomial degree of the trial space. All theoretical results are justified by numerical tests.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.