中国科学院数学与系统科学研究院期刊网
JCM

15 May 2025, Volume 43 Issue 3
    

  • Select all
    |
  • NUMERICAL ENERGY DISSIPATION FOR TIME-FRACTIONAL PHASE-FIELD EQUATIONS
    Chaoyu Quan, Tao Tang, Jiang Yang
    Journal of Computational Mathematics. 2025, 43(3): 515-539. https://doi.org/10.4208/jcm.2311-m2021-0199
    Abstract ( ) Download PDF   Knowledge map   Save
    The numerical integration of phase-field equations is a delicate task which needs to recover at the discrete level intrinsic properties of the solution such as energy dissipation and maximum principle. Although the theory of energy dissipation for classical phase field models is well established, the corresponding theory for time-fractional phase-field models is still incomplete. In this article, we study certain nonlocal-in-time energies using the first-order stabilized semi-implicit L1 scheme. In particular, we will establish a discrete fractional energy law and a discrete weighted energy law. The extension for a (2-α)-order L1 scalar auxiliary variable scheme will be investigated. Moreover, we demonstrate that the energy bound is preserved for the L1 schemes with nonuniform time steps. Several numerical experiments are carried to verify our theoretical analysis.
  • PROXIMAL ADMM APPROACH FOR IMAGE RESTORATION WITH MIXED POISSON-GAUSSIAN NOISE
    Miao Chen, Yuchao Tang, Jie Zhang, Tieyong Zeng
    Journal of Computational Mathematics. 2025, 43(3): 540-568. https://doi.org/10.4208/jcm.2212-m2022-0122
    Abstract ( ) Download PDF   Knowledge map   Save
    Image restoration based on total variation has been widely studied owing to its edgepreservation properties. In this study, we consider the total variation infimal convolution (TV-IC) image restoration model for eliminating mixed Poisson-Gaussian noise. Based on the alternating direction method of multipliers (ADMM), we propose a complete splitting proximal bilinear constraint ADMM algorithm to solve the TV-IC model. We prove the convergence of the proposed algorithm under mild conditions. In contrast with other algorithms used for solving the TV-IC model, the proposed algorithm does not involve any inner iterations, and each subproblem has a closed-form solution. Finally, numerical experimental results demonstrate the efficiency and effectiveness of the proposed algorithm.
  • NUMERICAL METHODS FOR APPROXIMATING STOCHASTIC SEMILINEAR TIME-FRACTIONAL RAYLEIGH-STOKES EQUATIONS
    Mariam Al-Maskari
    Journal of Computational Mathematics. 2025, 43(3): 569-587. https://doi.org/10.4208/jcm.2311-m2023-0047
    Abstract ( ) Download PDF   Knowledge map   Save
    This paper investigates a semilinear stochastic fractional Rayleigh-Stokes equation featuring a Riemann-Liouville fractional derivative of order α ∈ (0, 1) in time and a fractional time-integral noise. The study begins with an examination of the solution’s existence, uniqueness, and regularity. The spatial discretization is then carried out using a finite element method, and the error estimate is analyzed. A convolution quadrature method generated by the backward Euler method is employed for the time discretization resulting in a fully discrete scheme. The error estimate for the fully discrete solution is considered based on the regularity of the solution, and a strong convergence rate is established. The paper concludes with numerical tests to validate the theoretical findings.
  • CENTRAL LIMIT THEOREM FOR TEMPORAL AVERAGE OF BACKWARD EULER-MARUYAMA METHOD
    Diancong Jin
    Journal of Computational Mathematics. 2025, 43(3): 588-614. https://doi.org/10.4208/jcm.2311-m2023-0147
    Abstract ( ) Download PDF   Knowledge map   Save
    This work focuses on the temporal average of the backward Euler-Maruyama (BEM) method, which is used to approximate the ergodic limit of stochastic ordinary differential equations (SODEs). We give the central limit theorem (CLT) of the temporal average of the BEM method, which characterizes its asymptotics in distribution. When the deviation order is smaller than the optimal strong order, we directly derive the CLT of the temporal average through that of original equations and the uniform strong order of the BEM method. For the case that the deviation order equals to the optimal strong order, the CLT is established via the Poisson equation associated with the generator of original equations. Numerical experiments are performed to illustrate the theoretical results. The main contribution of this work is to generalize the existing CLT of the temporal average of numerical methods to that for SODEs with super-linearly growing drift coefficients.
  • A HIGH ORDER SCHEME FOR FRACTIONAL DIFFERENTIAL EQUATIONS WITH THE CAPUTO-HADAMARD DERIVATIVE
    Xingyang Ye, Junying Cao, Chuanju Xu
    Journal of Computational Mathematics. 2025, 43(3): 615-640. https://doi.org/10.4208/jcm.2312-m2023-0098
    Abstract ( ) Download PDF   Knowledge map   Save
    In this paper, we consider numerical solutions of the fractional diffusion equation with the α order time fractional derivative defined in the Caputo-Hadamard sense. A high order time-stepping scheme is constructed, analyzed, and numerically validated. The contribution of the paper is twofold: 1) regularity of the solution to the underlying equation is investigated, 2) a rigorous stability and convergence analysis for the proposed scheme is performed, which shows that the proposed scheme is 3 + α order accurate. Several numerical examples are provided to verify the theoretical statement.
  • CRANK-NICOLSON GALERKIN APPROXIMATIONS FOR LOGARITHMIC KLEIN-GORDON EQUATION
    Fang Chen, Meng Li, Yanmin Zhao
    Journal of Computational Mathematics. 2025, 43(3): 641-672. https://doi.org/10.4208/jcm.2312-m2023-0185
    Abstract ( ) Download PDF   Knowledge map   Save
    This paper presents three regularized models for the logarithmic Klein-Gordon equation. By using a modified Crank-Nicolson method in time and the Galerkin finite element method (FEM) in space, a fully implicit energy-conservative numerical scheme is constructed for the local energy regularized model that is regarded as the best one among the three regularized models. Then, the cut-off function technique and the time-space error splitting technique are innovatively combined to rigorously analyze the unconditionally optimal and high-accuracy convergence results of the numerical scheme without any coupling condition between the temporal step size and the spatial mesh width. The theoretical framework is uniform for the other two regularized models. Finally, numerical experiments are provided to verify our theoretical results. The analytical techniques in this work are not limited in the FEM, and can be directly extended into other numerical methods. More importantly, this work closes the gap for the unconditional error/stability analysis of the numerical methods for the logarithmic systems in higher dimensional spaces.
  • AN ITERATIVE TWO-GRID METHOD FOR STRONGLY NONLINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
    Jiajun Zhan, Lei Yang, Xiaoqing Xing, Liuqiang Zhong
    Journal of Computational Mathematics. 2025, 43(3): 673-689. https://doi.org/10.4208/jcm.2305-m2023-0088
    Abstract ( ) Download PDF   Knowledge map   Save
    We design and analyze an iterative two-grid algorithm for the finite element discretizations of strongly nonlinear elliptic boundary value problems in this paper. We propose an iterative two-grid algorithm, in which a nonlinear problem is first solved on the coarse space, and then a symmetric positive definite problem is solved on the fine space. The main contribution in this paper is to establish a first convergence analysis, which requires dealing with four coupled error estimates, for the iterative two-grid methods. We also present some numerical experiments to confirm the efficiency of the proposed algorithm.
  • ERROR ANALYSIS OF FRACTIONAL COLLOCATION METHODS FOR VOLTERRA INTEGRO-DIFFERENTIAL EQUATIONS WITH NONCOMPACT OPERATORS
    Zheng Ma, Chengming Huang, Anatoly A. Alikhanov
    Journal of Computational Mathematics. 2025, 43(3): 690-707. https://doi.org/10.4208/jcm.2401-m2023-0196
    Abstract ( ) Download PDF   Knowledge map   Save
    This paper is concerned with the numerical solution of Volterra integro-differential equations with noncompact operators. The focus is on the problems with weakly singular solutions. To handle the initial weak singularity of the solution, a fractional collocation method is applied. A rigorous hp-version error analysis of the numerical method under a weighted H1-norm is carried out. The result shows that the method can achieve high order convergence for such equations. Numerical experiments are also presented to confirm the effectiveness of the proposed method.
  • A DECOUPLED, LINEARLY IMPLICIT AND UNCONDITIONALLY ENERGY STABLE SCHEME FOR THE COUPLED CAHN-HILLIARD SYSTEMS
    Dan Zhao, Dongfang Li, Yanbin Tang, Jinming Wen
    Journal of Computational Mathematics. 2025, 43(3): 708-730. https://doi.org/10.4208/jcm.2402-m2023-0079
    Abstract ( ) Download PDF   Knowledge map   Save
    We present a decoupled, linearly implicit numerical scheme with energy stability and mass conservation for solving the coupled Cahn-Hilliard system. The time-discretization is done by leap-frog method with the scalar auxiliary variable (SAV) approach. It only needs to solve three linear equations at each time step, where each unknown variable can be solved independently. It is shown that the semi-discrete scheme has second-order accuracy in the temporal direction. Such convergence results are proved by a rigorous analysis of the boundedness of the numerical solution and the error estimates at different time-level. Numerical examples are presented to further confirm the validity of the methods.
  • ERROR ANALYSIS OF VIRTUAL ELEMENT METHODS FOR THE TIME-DEPENDENT POISSON-NERNST-PLANCK EQUATIONS
    Ying Yang, Ya Liu, Yang Liu, Shi Shu
    Journal of Computational Mathematics. 2025, 43(3): 731-770. https://doi.org/10.4208/jcm.2401-m2023-0130
    Abstract ( ) Download PDF   Knowledge map   Save
    We discuss and analyze the virtual element method on general polygonal meshes for the time-dependent Poisson-Nernst-Planck (PNP) equations, which are a nonlinear coupled system widely used in semiconductors and ion channels. After presenting the semi-discrete scheme, the optimal H1 norm error estimates are presented for the time-dependent PNP equations, which are based on some error estimates of a virtual element energy projection. The Gummel iteration is used to decouple and linearize the PNP equations and the error analysis is also given for the iteration of fully discrete virtual element approximation. The numerical experiment on different polygonal meshes verifies the theoretical convergence results and shows the efficiency of the virtual element method.