BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Buyang Li

Buyang Li contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math.NANumerical Analysismath.AP

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

12 published item(s)

preprint2022arXiv

A second-order low-regularity correction of Lie splitting for the semilinear Klein--Gordon equation

The numerical approximation of the semilinear Klein--Gordon equation in the $d$-dimensional space, with $d=1,2,3$, is studied by analyzing the consistency errors in approximating the solution. By discovering and utilizing a new cancellation structure in the semilinear Klein--Gordon equation, a low-regularity correction of the Lie splitting method is constructed, which can have second-order convergence in the energy space under the regularity condition $(u,\partial_tu)\in L^\infty(0,T;H^{1+\frac{d}{4}}\times H^{\frac{d}{4}})$, where $d=1,2,3$ denotes the dimension of space. In one dimension, the proposed method is shown to have a convergence order arbitrarily close to $\frac53$ in the energy space for solutions in the same space, i.e. no additional regularity in the solution is required. Rigorous error estimates are presented for a fully discrete spectral method with the proposed low-regularity time-stepping scheme. Numerical examples are provided to support the theoretical analysis and to illustrate the performance of the proposed method in approximating both nonsmooth and smooth solutions of the semilinear Klein--Gordon equation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An unfiltered low-regularity integrator for the KdV equation with solutions below ${\bf H^1}$

This article is concerned with the construction and analysis of new time discretizations for the KdV equation on a torus for low-regularity solutions below $H^1$. New harmonic analysis tools, including new averaging approximations to the exponential phase functions, new frequency decomposition techniques, and new trilinear estimates of the KdV operator, are established for the construction and analysis of time discretizations with higher convergence orders under low-regularity conditions. In addition, new techniques are introduced to establish stability estimates of time discretizations under low-regularity conditions without using filters when the energy techniques fail. The proposed method is proved to be convergent with order $γ$ (up to a logarithmic factor) in $L^2$ under the regularity condition $u\in C([0,T];H^γ)$ for $γ\in(0,1]$.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Maximal regularity of backward difference time discretization for evolving surface PDEs and its application to nonlinear problems

Maximal parabolic $L^p$-regularity of linear parabolic equations on an evolving surface is shown by pulling back the problem to the initial surface and studying the maximal $L^p$-regularity on a fixed surface. By freezing the coefficients in the parabolic equations at a fixed time and utilizing a perturbation argument around the freezed time, it is shown that backward difference time discretizations of linear parabolic equations on an evolving surface along characteristic trajectories can preserve maximal $L^p$-regularity in the discrete setting. The result is applied to prove the stability and convergence of time discretizations of nonlinear parabolic equations on an evolving surface, with linearly implicit backward differentiation formulae characteristic trajectories of the surface, for general locally Lipschitz nonlinearities. The discrete maximal $L^p$-regularity is used to prove the boundedness and stability of numerical solutions in the $L^\infty(0,T;W^{1,\infty})$ norm, which is used to bound the nonlinear terms in the stability analysis. Optimal-order error estimates of time discretizations in the $L^\infty(0,T;W^{1,\infty})$ norm is obtained by combining the stability analysis with the consistency estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A fully discrete low-regularity integrator for the 1D periodic cubic nonlinear Schrödinger equation

A fully discrete and fully explicit low-regularity integrator is constructed for the one-dimensional periodic cubic nonlinear Schrödinger equation. The method can be implemented by using fast Fourier transform with $O(N\ln N)$ operations at every time level, and is proved to have an $L^2$-norm error bound of $O(τ\sqrt{\ln(1/τ)}+N^{-1})$ for $H^1$ initial data, without requiring any CFL condition, where $τ$ and $N$ denote the temporal stepsize and the degree of freedoms in the spatial discretisation, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Analysis of fully discrete finite element methods for 2D Navier--Stokes equations with critical initial data

First-order convergence in time and space is proved for a fully discrete semi-implicit finite element method for the two-dimensional Navier--Stokes equations with $L^2$ initial data in convex polygonal domains, without extra regularity assumptions or grid-ratio conditions. The proof utilises the smoothing properties of the Navier--Stokes equations, an appropriate duality argument, and the smallness of the numerical solution in the discrete $L^2(0,t_m;H^1)$ norm when $t_m$ is smaller than some constant. Numerical examples are provided to support the theoretical analysis.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Convergence of Dziuk's semidiscrete finite element method for mean curvature flow of closed surfaces with high-order finite elements

Dziuk's surface finite element method for mean curvature flow has had significant impact on the development of parametric and evolving surface finite element methods for surface evolution equations and curvature flows. However, Dziuk's surface finite element method for mean curvature flow of closed surfaces remains still open since it was proposed in 1990. In this article, we prove convergence of Dziuk's semidiscrete surface finite element method with high-order finite elements for mean curvature flow of closed surfaces by utilizing the matrix-vector formulation of evolving surface finite element methods and a monotone structure of the nonlinear discrete surface Laplacian proved in this paper.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

High-order mass- and energy-conserving SAV-Gauss collocation finite element methods for the nonlinear Schrödinger equation

A family of arbitrarily high-order fully discrete space-time finite element methods are proposed for the nonlinear Schrödinger equation based on the scalar auxiliary variable formulation, which consists of a Gauss collocation temporal discretization and the finite element spatial discretization. The proposed methods are proved to be well-posed and conserving both mass and energy at the discrete level. An error bound of the form $O(h^p+τ^{k+1})$ in the $L^\infty(0,T;H^1)$-norm is established, where $h$ and $τ$ denote the spatial and temporal mesh sizes, respectively, and $(p,k)$ is the degree of the space-time finite elements. Numerical experiments are provided to validate the theoretical results on the convergence rates and conservation properties. The effectiveness of the proposed methods in preserving the shape of a soliton wave is also demonstrated by numerical results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A convergent algorithm for forced mean curvature flow driven by diffusion on the surface

The evolution of a closed two-dimensional surface driven by both mean curvature flow and a reaction--diffusion process on the surface is formulated into a system, which couples the velocity law not only to the surface partial differential equation but also to the evolution equations for the geometric quantities, namely the normal vector and the mean curvature on the surface. Two algorithms are considered for the obtained system. Both methods combine surface finite elements as a space discretisation and linearly implicit backward difference formulae for time integration. Based on our recent results for mean curvature flow, one of the algorithms directly admits a convergence proof for its full discretisation in the case of finite elements of polynomial degree at least two and backward difference formulae of orders two to five. Numerical examples are provided to support and complement the theoretical convergence results (demonstrating the convergence properties of the method without error estimate), and demonstrate the effectiveness of the methods in simulating a three-dimensional tumour growth model.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A convergent evolving finite element algorithm for Willmore flow of closed surfaces

A proof of convergence is given for a novel evolving surface finite element semi-discretization of Willmore flow of closed two-dimensional surfaces, and also of surface diffusion flow. The numerical method proposed and studied here discretizes fourth-order evolution equations for the normal vector and mean curvature, reformulated as a system of second-order equations, and uses these evolving geometric quantities in the velocity law interpolated to the finite element space. This numerical method admits a convergence analysis in the case of continuous finite elements of polynomial degree at least two. The error analysis combines stability estimates and consistency estimates to yield optimal-order $H^1$-norm error bounds for the computed surface position, velocity, normal vector and mean curvature. The stability analysis is based on the matrix--vector formulation of the finite element method and does not use geometric arguments. The geometry enters only into the consistency estimates. Numerical experiments illustrate and complement the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Maximal regularity of multistep fully discrete finite element methods for parabolic equations

This article extends the semidiscrete maximal $L^p$-regularity results in [27] to multistep fully discrete finite element methods for parabolic equations with more general diffusion coefficients in $W^{1,d+β}$, where $d$ is the dimension of space and $β>0$. The maximal angles of $R$-boundedness are characterized for the analytic semigroup $e^{zA_h}$ and the resolvent operator $z(z-A_h)^{-1}$, respectively, associated to an elliptic finite element operator $A_h$. Maximal $L^p$-regularity, optimal $\ell^p(L^q)$ error estimate, and $\ell^p(W^{1,q})$ estimate are established for fully discrete finite element methods with multistep backward differentiation formula.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Subdiffusion with Time-Dependent Coefficients: Improved Regularity and Second-Order Time Stepping

This article concerns second-order time discretization of subdiffusion equations with time-dependent diffusion coefficients. High-order differentiability and regularity estimates are established for subdiffusion equations with time-dependent coefficients. Using these regularity results and a perturbation argument of freezing the diffusion coefficient, we prove that the convolution quadrature generated by the second-order backward differentiation formula, with proper correction at the first time step, can achieve second-order convergence for both nonsmooth initial data and incompatible source term. Numerical experiments are consistent with the theoretical results.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Weak discrete maximum principle of finite element methods in convex polyhedra

We prove that the Galerkin finite element solution $u_h$ of the Laplace equation in a convex polyhedron $\varOmega$, with a quasi-uniform tetrahedral partition of the domain and with finite elements of polynomial degree $r\ge 1$, satisfies the following weak maximum principle: \begin{align*} \left\|u_{h}\right\|_{L^{\infty}(\varOmega)} \le C\left\|u_{h}\right\|_{L^{\infty}(\partial \varOmega)} , \end{align*} with a constant $C$ independent of the mesh size $h$. By using this result, we show that the Ritz projection operator $R_h$ is stable in $L^\infty$ norm uniformly in $h$ for $r\geq 2$, i.e. \begin{align*} \|R_hu\|_{L^{\infty}(\varOmega)} \le C\|u\|_{L^{\infty}(\varOmega)} . \end{align*} Thus we remove a logarithmic factor appearing in the previous results for convex polyhedral domains.

0 citations0 reviewsNo code linkedOpen access