Source author record

Buyang Li

Buyang Li 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

29works
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

29 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.

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]$.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

preprint2016arXiv

Convergence of a decoupled mixed FEM for the dynamic Ginzburg--Landau equations in nonsmooth domains with incompatible initial data

In this paper, we propose a fully discrete mixed finite element method for solving the time-dependent Ginzburg--Landau equations, and prove the convergence of the finite element solutions in general curved polyhedra, possibly nonconvex and multi-connected, without assumptions on the regularity of the solution. Global existence and uniqueness of weak solutions for the PDE problem are also obtained in the meantime. A decoupled time-stepping scheme is introduced, which guarantees that the discrete solution has bounded discrete energy, and the finite element spaces are chosen to be compatible with the nonlinear structure of the equations. Based on the boundedness of the discrete energy, we prove the convergence of the finite element solutions by utilizing a uniform $L^{3+δ}$ regularity of the discrete harmonic vector fields, establishing a discrete Sobolev embedding inequality for the Nédélec finite element space, and introducing a $\ell^2(W^{1,3+δ})$ estimate for fully discrete solutions of parabolic equations. The numerical example shows that the constructed mixed finite element solution converges to the true solution of the PDE problem in a nonsmooth and multi-connected domain, while the standard Galerkin finite element solution does not converge.

preprint2016arXiv

Maximal $\bf L^p$ analysis of finite element solutions for parabolic equations with nonsmooth coefficients in convex polyhedra

The paper is concerned with Galerkin finite element solutions for parabolic equations in a convex polygon or polyhehron with a diffusion coefficient in $W^{1,N+ε}$ for some $ε>0$, where $N$ denotes the dimension of the domain. We prove the analyticity of the semigroup generated by the discrete elliptic operator, the discrete maximal $L^p$ regularity and the optimal $L^p$ error estimate of the finite element solution for the parabolic equation.

preprint2016arXiv

Maximum norm analysis of implicit-explicit backward difference formulas for nonlinear parabolic equations

We establish optimal order a priori error estimates for implicit-explicit BDF methods for abstract semilinear parabolic equations with time-dependent operators in a complex Banach space settings, under a sharp condition on the non-self-adjointness of the linear operator. Our approach relies on the discrete maximal parabolic regularity of implicit BDF schemes for autonomous linear parabolic equations, recently established in [20], and on ideas from [7]. We illustrate the applicability of our results to four initial and boundary value problems, namely two for second order, one for fractional order, and one for fourth order, namely the Cahn-Hilliard, parabolic equations.

preprint2015arXiv

A-stable time discretizations preserve maximal parabolic regularity

It is shown that for a parabolic problem with maximal $L^p$-regularity (for $1<p<\infty$), the time discretization by a linear multistep method or Runge--Kutta method has maximal $\ell^p$-regularity uniformly in the stepsize if the method is A-stable (and satisfies minor additional conditions). In particular, the implicit Euler method, the Crank-Nicolson method, the second-order backward difference formula (BDF), and the Radau IIA and Gauss Runge--Kutta methods of all orders preserve maximal regularity. The proof uses Weis' characterization of maximal $L^p$-regularity in terms of $R$-boundedness of the resolvent, a discrete operator-valued Fourier multiplier theorem by Blunck, and generating function techniques that have been familiar in the stability analysis of time discretization methods since the work of Dahlquist. The A($α$)-stable higher-order BDF methods have maximal $\ell^p$-regularity under an $R$-boundedness condition in a larger sector. As an illustration of the use of maximal regularity in the error analysis of discretized nonlinear parabolic equations, it is shown how error bounds are obtained without using any growth condition on the nonlinearity or for nonlinearities having singularities.

preprint2014arXiv

A new approach for numerical simulation of the time-dependent Ginzburg-Landau equations

We introduce a new approach for finite element simulations of the time-dependent Ginzburg-Landau equations (TDGL) in a general curved polygon, possibly with reentrant corners. Specifically, we reformulate the TDGL into an equivalent system of equations by decomposing the magnetic potential to the sum of its divergence-free and curl-free parts, respectively. Numerical simulations of vortex dynamics show that, in a domain with reentrant corners, the new approach is much more stable and accurate than the old approaches of solving the TDGL directly (under either the temporal gauge or the Lorentz gauge); in a convex domain, the new approach gives comparably accurate solutions as the old approaches.

preprint2014arXiv

Convergence of a decoupled mixed FEM for miscible displacement in interfacial porous media

In this paper, we study the stability and convergence of a decoupled and linearized mixed finite element method (FEM) for incompressible miscible displacement in a porous media whose permeability and porosity are discontinuous across some interfaces. We show that the proposed scheme has optimal-order convergence rate unconditionally, without restriction on the grid ratio (between the time-step size and spatial mesh size). Previous works all required certain restrictions on the grid ratio except for the problem with globally smooth permeability and porosity. Our idea is to introduce an intermediate system of elliptic interface problems, whose solution is uniformly regular in each subdomain separated by the interfaces and its finite element solution coincides with the fully discrete solution of the original problem. In order to prove the boundedness of the fully discrete solution, we study the finite element discretization of the intermediate system of elliptic interface problems.

preprint2014arXiv

Global well-posedness of the time-dependent Ginzburg-Landau superconductivity model in curved polyhedra

We study the time-dependent Ginzburg--Landau equations in a three-dimensional curved polyhedron (possibly nonconvex). Compared with the previous works, we prove existence and uniqueness of a global weak solution based on weaker regularity of the solution in the presence of edges or corners, where the magnetic potential may not be in $L^2(0,T;H^1(Ω)^3)$.

preprint2014arXiv

Linearized FE approximations to a nonlinear gradient flow

We study fully discrete linearized Galerkin finite element approximations to a nonlinear gradient flow, applications of which can be found in many areas. Due to the strong nonlinearity of the equation, existing analyses for implicit schemes require certain restrictions on the time step and no analysis has been explored for linearized schemes. This paper focuses on the unconditionally optimal $L^2$ error estimate of a linearized scheme. The key to our analysis is an iterated sequence of time-discrete elliptic equations and a rigorous analysis of its solution. We prove the $W^{1,\infty}$ boundedness of the solution of the time-discrete system and the corresponding FE solution, based on a more precise estimate of elliptic PDEs in $W^{2,2+ε}$ and a physical feature of the gradient-dependent diffusion coefficient. Numerical examples are provided to support our theoretical analysis.

preprint2014arXiv

Mathematical and numerical analysis of time-dependent Ginzburg--Landau equations in nonconvex polygons based on Hodge decomposition

We prove well-posedness of time-dependent Ginzburg--Landau system in a nonconvex polygonal domain, and decompose the solution as a regular part plus a singular part. We see that the magnetic potential is not in $H^1$ in general, and the finite element method (FEM) may give incorrect solutions. To remedy this situation, we reformulate the equations into an equivalent system of elliptic and parabolic equations based on the Hodge decomposition, which avoids direct calculation of the magnetic potential. The essential unknowns of the reformulated system admit $H^1$ solutions and can be solved correctly by the FEMs. We then propose a decoupled and linearized FEM to solve the reformulated equations and present error estimates based on proved regularity of the solution. Numerical examples are provided to support our theoretical analysis and show the efficiency of the method.

preprint2014arXiv

Maximum-norm stability and maximal L^p regularity of FEMs for parabolic equations with Lipschitz continuous coefficients

In this paper, we study the semi-discrete Galerkin finite element method for parabolic equations with Lipschitz continuous coefficients. We prove the maximum-norm stability of the semigroup generated by the corresponding elliptic finite element operator, and prove the space-time stability of the parabolic projection onto the finite element space in $L^\infty(Ω_T)$ and $L^p((0,T);L^q(Ω))$, $1<p,q<\infty$. The maximal $L^p$ regularity of the parabolic finite element equation is also established.

preprint2014arXiv

Regularity of the diffusion-dispersion tensor and error analysis of Galerkin FEMs for a porous media flow

We study Galerkin finite element methods for an incompressible miscible flow in porous media with the commonly-used Bear-Scheidegger diffusion-dispersion tensor $D({\bf u}) = Φd_m I + |{\bf u}| \big ( α_T I + (α_L - α_T) \frac{{\bf u} \otimes {\bf u}}{|{\bf u}|^2}\big)$. The traditional approach to optimal $L^\infty((0,T);L^2)$ error estimates is based on an elliptic Ritz projection, which usually requires the regularity of $\nabla_x\partial_tD({\bf u}(x,t)) \in L^p(Ω_T)$. However, the Bear-Scheidegger diffusion-dispersion tensor may not satisfy the regularity condition even for a smooth velocity field ${\bf u}$. A new approach is presented in this paper, in terms of a parabolic projection, which only requires the Lipschitz continuity of $D({\bf u})$. With the new approach, we establish optimal $L^p((0,T);L^q)$ error estimates and an almost optimal $L^\infty((0,T);L^\infty)$ error estimate.

preprint2013arXiv

Error analysis of linearized semi-implicit Galerkin finite element methods for nonlinear parabolic equations

This paper is concerned with the time-step condition of commonly-used linearized semi-implicit schemes for nonlinear parabolic PDEs with Galerkin finite element approximations. In particular, we study the time-dependent nonlinear Joule heating equations. We present optimal error estimates of the semi-implicit Euler scheme in both the $L^2$ norm and the $H^1$ norm without any time-step restriction. Theoretical analysis is based on a new splitting of the error and precise analysis of a corresponding time-discrete system. The method used in this paper can be applied to more general nonlinear parabolic systems and many other linearized (semi)-implicit time discretizations for which previous works often require certain restriction on the time-step size $τ$.

preprint2013arXiv

Unconditional convergence and optimal error estimates of a Galerkin-mixed FEM for incompressible miscible flow in porous media

In this paper, we study the unconditional convergence and error estimates of a Galerkin-mixed FEM with the linearized semi-implicit Euler time-discrete scheme for the equations of incompressible miscible flow in porous media. We prove that the optimal $L^2$ error estimates hold without any time-step (convergence) condition, while all previous works require certain time-step condition. Our theoretical results provide a new understanding on commonly-used linearized schemes for nonlinear parabolic equations. The proof is based on a splitting of the error function into two parts: the error from the time discretization of the PDEs and the error from the finite element discretization of corresponding time-discrete PDEs. The approach used in this paper is applicable for more general nonlinear parabolic systems and many other linearized (semi)-implicit time discretizations.

preprint2013arXiv

Unconditionally optimal error analysis of fully discrete Galerkin methods for general nonlinear parabolic equations

The paper focuses on unconditionally optimal error analysis of the fully discrete Galerkin finite element methods for a general nonlinear parabolic system in $\R^d$ with $d=2,3$. In terms of a corresponding time-discrete system of PDEs as proposed in \cite{LS1}, we split the error function into two parts, one from the temporal discretization and one the spatial discretization. We prove that the latter is $τ$-independent and the numerical solution is bounded in the $L^{\infty}$ and $W^{1,\infty}$ norms by the inverse inequalities. With the boundedness of the numerical solution, optimal error estimates can be obtained unconditionally in a routine way. Several numerical examples in two and three dimensional spaces are given to support our theoretical analysis.

preprint2013arXiv

Uniform BMO estimate of parabolic equations and global well-posedness of the thermistor problem

Global well-posedness of the time-dependent (degenerate) thermistor problem remains open for many years. In this paper, we solve the problem by establishing a uniform-in-time BMO estimate of inhomogeneous parabolic equations. Applying this estimate to the temperature equation, we derive a BMO bound of the temperature uniform with respect to time, which implies that the electric conductivity is a $A_2$ weight. The Hölder continuity of the electric potential is then proved by applying the De Giorgi--Nash--Moser estimate for degenerate elliptic equations with $A_2$ coefficient. Uniqueness of solution is proved based on the established regularity of the weak solution. Our results also imply the existence of a global classical solution when the initial and boundary data are smooth.

preprint2012arXiv

Unconditionally optimal error estimates of a Crank--Nicolson Galerkin method for the nonlinear thermistor equations

This paper focuses on unconditionally optimal error analysis of an uncoupled and linearized Crank--Nicolson Galerkin finite element method for the time-dependent nonlinear thermistor equations in $d$-dimensional space, $d=2,3$. We split the error function into two parts, one from the spatial discretization and one from the temporal discretization, by introducing a corresponding time-discrete (elliptic) system. We present a rigorous analysis for the regularity of the solution of the time-discrete system and error estimates of the time discretization. With these estimates and the proved regularity, optimal error estimates of the fully discrete Crank--Nicolson Galerkin method are obtained unconditionally. Numerical results confirm our analysis and show the efficiency of the method.