Source author record

Kejia Pan

Kejia Pan appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

6works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

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.

Published work

6 published item(s)

preprint2022arXiv

Convergence study of IB methods for Stokes equations with no-slip boundary conditions

Peskin's Immersed Boundary (IB) model and method are among one of the most important modeling tools and numerical methods. The IB method has been known to be first order accurate in the velocity. However, almost no rigorous theoretical proof can be found in the literature for Stokes equations with a prescribed velocity boundary condition. In this paper, it has been shown that the pressure of the Stokes equation has a convergence order $O(\sqrt{h} |\log h| )$ in the $L^2$ norm while the velocity has an $O(h |\log h| )$ convergence order in the infinity norm in two-space dimensions. The proofs are based on splitting the singular source terms, discrete Green functions on finite lattices with homogeneous and Neumann boundary conditions, a new discovered simplest $L^2$ discrete delta function, and the convergence proof of the IB method for elliptic interface problems \cite{li:mathcom}. The conclusion in this paper can apply to problems with different boundary conditions as long as the problems are wellposed. The proof process also provides an efficient way to decouple the system into three Helmholtz/Poisson equations without affecting the order of convergence. A non-trivial numerical example is also provided to confirm the theoretical analysis and the simple new discrete delta function.

preprint2022arXiv

High order compact schemes for flux type BCs

In this paper new innovative fourth order compact schemes for Robin and Neumann boundary conditions have been developed for boundary value problems of elliptic PDEs in two and three dimensions. Different from traditional finite difference operator approach, which may not work for flux type of boundary conditions, carefully designed undetermined coefficient methods are utilized in developing high order compact (HOC) schemes. The new methods not only can be utilized to design HOC schemes for flux type of boundary conditions but can also be applied to general elliptic PDEs including Poisson, Helmholtz, diffusion-advection, and anisotropic equations with linear boundary conditions. In the new developed HOC methods, the coefficient matrices are generally M-matrices, which guarantee the discrete maximum principle for well-posed problems, so the convergence of the HOC methods. The developed HOC methods are versatile and can cover most of high order compact schemes in the literature. The HOC methods for Robin boundary conditions and for anisotropic diffusion and advection equations with Robin or even Dirichlet boundary conditions are likely the first ones that have ever been developed. With the help of pseudo-inverse, or SVD solutions, we have also observed that the developed HOC methods usually have smaller error constants compared with traditional HOC methods when applicable. Non-trivial examples with large wave numbers and oscillatory solutions are presented to confirm the performance of the new HOC methods.

preprint2021arXiv

A multigrid-reduction-in-time solver with a new two-level convergence for unsteady fractional Laplacian problems

The multigrid-reduction-in-time (MGRIT) technique has proven to be successful in achieving higher run-time speedup by exploiting parallelism in time. The goal of this article is to develop and analyze a MGRIT algorithm, using FCF-relaxation with time-dependent time-grid propagators, to seek the finite element approximations of unsteady fractional Laplacian problems. The multigrid with line smoother proposed in [L. Chen, R. H. Nochetto, E. Ot{á}rola, A. J. Salgado, Math. Comp. 85 (2016) 2583--2607] is chosen to be the spatial solver. Motivated by [B. S. Southworth, SIAM J. Matrix Anal. Appl. 40 (2019) 564--608], we provide a new temporal eigenvalue approximation property and then deduce a generalized two-level convergence theory which removes the previous unitary diagonalization assumption on the fine and coarse time-grid propagators required in [X. Q. Yue, S. Shu, X. W. Xu, W. P. Bu, K. J. Pan, Comput. Math. Appl. 78 (2019) 3471--3484]. Numerical computations are included to confirm the theoretical predictions and demonstrate the sharpness of the derived convergence upper bound.

preprint2016arXiv

An extrapolation cascadic multigrid method combined with a fourth order compact scheme for 3D poisson equation

In this paper, we develop an EXCMG method to solve the three-dimensional Poisson equation on rectangular domains by using the compact finite difference (FD) method with unequal meshsizes in different coordinate directions. The resulting linear system from compact FD discretization is solved by the conjugate gradient (CG) method with a relative residual stopping criterion. By combining the Richardson extrapolation and tri-quartic Lagrange interpolation for the numerical solutions from two-level of grids (current and previous grids), we are able to produce an extremely accurate approximation of the actual numerical solution on the next finer grid, which can greatly reduce the number of relaxation sweeps needed. Additionally, a simple method based on the midpoint extrapolation formula is used for the fourth-order FD solutions on two-level of grids to achieve sixth-order accuracy on the entire fine grid cheaply and directly. The gradient of the numerical solution can also be easily obtained through solving a series of tridiagonal linear systems resulting from the fourth-order compact FD discretizations. Numerical results show that our EXCMG method is much more efficient than the classical V-cycle and W-cycle multigrid methods. Moreover, only few CG iterations are required on the finest grid to achieve full fourth-order accuracy in both the $L^2$-norm and $L^{\infty}$-norm for the solution and its gradient when the exact solution belongs to $C^6$. Finally, numerical result shows that our EXCMG method is still effective when the exact solution has a lower regularity, which widens the scope of applicability of our EXCMG method.

preprint2015arXiv

A linearly implicit conservative difference scheme for the generalized Rosenau-Kawahara-RLW equation

This paper concerns the numerical study for the generalized Rosenau-Kawahara-RLW equation obtained by coupling the generalized Rosenau-RLW equation and the generalized Rosenau-Kawahara equation. We first derive the energy conservation law of the equation, and then develop a three-level linearly implicit difference scheme for solving the equation. We prove that the proposed scheme is energy-conserved, unconditionally stable and second-order accurate both in time and space variables. Finally, numerical experiments are carried out to confirm the energy conservation, the convergence rates of the scheme and effectiveness for long-time simulation.

preprint2015arXiv

Nanostructures imaging via numerical solution of a 3-d inverse scattering problem without the phase information

Inverse scattering problems without the phase information arise in imaging of nanostructures whose sizes are hundreds of nanometers as well as in imaging of biological cells. The governing equation is the 3-d generalized Helmholtz equation with the unknown coefficient, which represents the spatially distributed dielectric constant. It is assumed that only the modulus of the complex valued wave field is measured on a frequency range. The phase is not measured. Two rigorous numerical methods are tested and their performances are compared for realistic ranges of parameters.