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

05 June 2026, Volume 44 Issue 4
    

  • Select all
    |
  • A. J. A. Ramos, L. G. M. Rosário, B. Feng, M. M. Freitas, L. S. Veras
    Journal of Computational Mathematics. 2026, 44(4): 891-916. https://doi.org/10.4208/jcm.2505-m2024-0233
    Abstract ( ) Download PDF   Knowledge map   Save
    In this work, we investigate a thermoviscoelastic system governed by the Guyer-Krumhansl model. The system is free from the paradox of infinite heat propagation speed and, furthermore, is more suitable for modeling complex problems involving heterogeneous materials on a macroscale and at room temperature. Firstly, we establish the well-posedness of the system using the theory of semigroups of linear operators, and then we prove uniform exponential decay with respect to a given physical parameter using the multiplier method. Subsequently, we discretize the system and propose a monotone and consistent numerical scheme using finite differences. The convergence of the numerical solution is proven by the Lax equivalence theorem. Finally, we present numerical experiments using MATLAB to demonstrate the accuracy and efficiency of the scheme that reproduce the theoretical results.
  • Wenjia Liu, Huilan Zeng, Shuo Zhang
    Journal of Computational Mathematics. 2026, 44(4): 917-947. https://doi.org/10.4208/jcm.2509-m2024-0165
    Abstract ( ) Download PDF   Knowledge map   Save
    In this paper, two O(h2)-accurate conservative finite element schemes with low-degree polynomials for the incompressible Stokes equations are presented. The schemes use respective H(div) finite element spaces, namely the third-order Brezzi-Douglas-Marini space and Brezzi-Douglas-Fortin-Marini space, with enhanced smoothness for the velocity and piecewise quadratic polynomials for the pressure, and are denoted as sBDM3-P2 and sBDFM3-P2 schemes, respectively. The discrete Korn inequality holds for both sBDM3 and sBDFM3 finite element spaces. For the sBDM3-P2 scheme, the inf-sup condition holds on general triangulations, and for the sBDFM3-P2 scheme, the inf-sup condition holds on triangulations with mild restriction. Both schemes achieve an energy norm of velocity errors of O(h2) order and an L2-norm of pressure errors of O(h2) order. Numerical experiments support the theoretical constructions.
  • Chaobao Huang, Yujie Yu, Na An, Hu Chen
    Journal of Computational Mathematics. 2026, 44(4): 948-966. https://doi.org/10.4208/jcm.2505-m2024-0260
    Abstract ( ) Download PDF   Knowledge map   Save
    In this work, a nonlinear subdiffusion initial-boundary value problem with a variable exponent is considered, whose solution behaves a weak singularity at initial time. By adding a corrected term to the nonuniform L1 scheme, a novel scheme is investigated to approximate the time-fractional Caputo derivative with a variable exponent. This scheme allows us to use a smaller grading parameter r to obtain a similar level of accuracy as that of the L1 method. Combining the proposed scheme with the finite element method in space and the Newton linearization for the nonlinear term, a fully discrete scheme is constructed. To obtain the unconditional optimal error estimate, the temporal-spatial splitting technique is adopted to derive the boundedness of the computed solution Uhn in L-norm. With the help of this bound and the discrete fractional Gronwall inequality, the optimal error analysis without certain temporal restrictions dependent on the spatial mesh size is derived. Furthermore, by using a simple postprocessing technique of the computed solution, the convergence order in the spatial direction is improved. Finally, numerical experiments are presented to verify the theoretical findings.
  • Qiang Han, Shihao Lan, Quanxin Zhu
    Journal of Computational Mathematics. 2026, 44(4): 967-990. https://doi.org/10.4208/jcm.2505-m2024-0206
    Abstract ( ) Download PDF   Knowledge map   Save
    In this paper, we design a novel explicit second order scheme with one step for forward backward stochastic differential equations, and the Crank-Nicolson scheme is a specific case of our proposed framework. We first establish a rigorous stability result, and then we derive precise error estimates. Moreover, we confirm that the proposed novel scheme is second order convergent. The theoretical results for the proposed methods are supported by numerical experiments.
  • Fenglong Qu, Mengyue Wang, Yanli Cui
    Journal of Computational Mathematics. 2026, 44(4): 991-1012. https://doi.org/10.4208/jcm.2505-m2024-0234
    Abstract ( ) Download PDF   Knowledge map   Save
    We consider an inverse problem of scattering by a mixed-type scatterer consisting of an inhomogeneous penetrable conductive medium and an impenetrable obstacle with generalized oblique derivative boundary condition induced by incident plane waves scattering. Relying on the well-posedness of the direct problem which can be proved directly by a variational method, we are interested in studying the inverse problem of developing a modified factorization method to simultaneously reconstruct the shape and location of the mixed-type scatterer. The complex refractive index and the generalized oblique derivative boundary condition may bring new challenges since the factorization method is closely related to the refractive index and the mixed boundary conditions. Finally, some numerical examples are given to show the effectiveness and feasibility of the inversion algorithm.
  • Wenhui Liu, Anhua Wan
    Journal of Computational Mathematics. 2026, 44(4): 1013-1048. https://doi.org/10.4208/jcm.2505-m2024-0057
    Abstract ( ) Download PDF   Knowledge map   Save
    Recovery of block sparse signals with partially-known block support information is of particular importance in compressed sensing. A uniform sufficient condition guaranteeing stable recovery of non-strictly block k-sparse signals is established via the weighted l2,pl2,q nonconvex minimization method, and the reconstruction error is precisely bounded in terms of the residual of block-sparsity and the measurement error. Furthermore, a series of contrastive numerical experiments reveal that exploiting the approximate block-sparsity characteristic and the nonuniform prior block support estimate substantially promotes the performance of reconstruction for block-structural signals.
  • Hongjin He, Kai Wang, Jintao Yu
    Journal of Computational Mathematics. 2026, 44(4): 1049-1082. https://doi.org/10.4208/jcm.2505-m2024-0095
    Abstract ( ) Download PDF   Knowledge map   Save
    In this paper, we propose a new primal-dual algorithmic framework for a class of convexconcave saddle point problems frequently arising from image processing and machine learning. Our algorithmic framework updates the primal variable between the twice calculations of the dual variable, thereby appearing a symmetric iterative scheme, which is accordingly called the symmetric primal-dual algorithm (SPIDA). It is noteworthy that the subproblems of our SPIDA are equipped with Bregman proximal regularization terms, which make SPIDA versatile in the sense that it enjoys an algorithmic framework to understand the iterative schemes of some existing algorithms, such as the classical augmented Lagrangian method (ALM), linearized ALM, and Jacobian splitting algorithms for linearly constrained optimization problems. Besides, our algorithmic framework allows us to derive some customized versions so that SPIDA works as efficiently as possible for structured optimization problems. Theoretically, under some mild conditions, we prove the global convergence of SPIDA and estimate the linear convergence rate under a generalized error bound condition defined by Bregman distance. Finally, a series of numerical experiments on the basis pursuit, robust principal component analysis, and image restoration demonstrate that our SPIDA works well on synthetic and real-world datasets.
  • Qiang-Gui Jin, Yao-Hui Ma
    Journal of Computational Mathematics. 2026, 44(4): 1083-1104. https://doi.org/10.4208/jcm.2505-m2024-0124
    Abstract ( ) Download PDF   Knowledge map   Save
    The process of direct method to solve large sparse linear equation mainly includes reordering, symbolic factorization, numerical factorization and triangular solving. Traditional symbolic factorization predicts the pattern of L based on single column and single row index. We propose to directly partition supernodes based on characteristics of matrix reordered by METIS, and then perform parallel symbolic factorization based on supernodes and row index fragments. A parallel block supernode numerical factorization strategy is proposed based on the concept of task pool here. In triangular solving stage, unlike traditional algorithms based on DAXPY and DDOT operations, we propose a new parallel triangular solving algorithm based on DGEMM and DTRSM operations. We name the parallel solver as finite element analysis direct solver (FEADS) and compare it with the advanced MKL PARDISO and MUMPS. The stiffness equations of 394770 and 719871 dimensions are solved using the solvers on two different computers. On the first computer, the solving efficiency of FEADS and MKL PARDISO is comparable, while MUMPS is relatively backward. On the second computer, FEADS performs especially well. For solving the case with 394770 dimensions, FEADS leads MKL PARDISO and MUMPS by 21.92% and 42.35%, respectively. For solving the case with 719871 dimensions, FEADS leads 34.75% and 38.38% respectively.
  • Liying Zhang, Qi Zhang
    Journal of Computational Mathematics. 2026, 44(4): 1105-1125. https://doi.org/10.4208/jcm.2505-m2024-0262
    Abstract ( ) Download PDF   Knowledge map   Save
    In this paper, we propose the parareal algorithms for stochastic Maxwell equations with the damping term driven by additive noise. The proposed algorithms proceed as twolevel temporal parallelizable integrators with the stochastic exponential integrator as the coarse G-propagator and both the exact solution integrator and the stochastic exponential integrator as the fine F-propagator. The mean-square convergence order of the proposed algorithms consistently increases to k, regardless of whether the exact solution integrator or the stochastic exponential integrator is chosen as the fine F-propagator. Several numerical experiments are illustrated in order to verify our theoretical findings for different choices of the iteration number k and the damping coefficient σ.
  • Ruimin Gao, Dongfang Li, Hongyu Qin
    Journal of Computational Mathematics. 2026, 44(4): 1126-1150. https://doi.org/10.4208/jcm.2505-m2024-0212
    Abstract ( ) Download PDF   Knowledge map   Save
    This paper proposes an energy-dissipative scheme for solving two- and three-dimensional time-fractional Navier-Stokes equations. The numerical scheme is constructed, using nonuniform L2-1σ approximation in the temporal direction and the Fourier spectral method in the spatial direction. It is shown that the numerical scheme can keep discrete energy stable and the numerical solutions are uniformly bounded without any restriction on step sizes. Error estimates of the fully-discrete scheme are presented. Moreover, a fast algorithm is applied to accelerate the computation. Numerical results in long time intervals are presented to confirm the effectiveness and high efficiency of the scheme.
  • Kaifang Liu
    Journal of Computational Mathematics. 2026, 44(4): 1151-1163. https://doi.org/10.4208/jcm.2506-m2025-0066
    Abstract ( ) Download PDF   Knowledge map   Save
    In this work, we develop a low-regularity error analysis for the interior-penalty discontinuous Galerkin (IPDG) method, incorporating numerical fluxes originally proposed by Brezzi et al. [Numer. Methods Partial Differential Equations, 16 (2000)]. Our analysis specifically addresses elliptic problems with solutions residing in the low-regularity space Hs, where 0 ≤ s < 1/2. Notably, our error estimates hold under two critical settings: discontinuous coefficients and general Lipschitz domains, precisely capturing the essential features of practical applications. We establish error estimates in the energy norm and the L2-norm, providing a complete theoretical framework for the IPDG method in low-regularity limitation. To systematically verify the theoretical results, we conduct some numerical experiments incorporating precision-controlled parameters that directly correspond to the analytical model’s constraints.
  • Fei Zhao, Zhiqiang Sheng, Guangwei Yuan
    Journal of Computational Mathematics. 2026, 44(4): 1164-1190. https://doi.org/10.4208/jcm.2502-m2024-0081
    Abstract ( ) Download PDF   Knowledge map   Save
    In this paper, we introduce a nonlinear finite volume scheme preserving discrete strong extremum principle (DSEP) for diffusion equations on tetrahedral meshes. In the construction of our nonlinear scheme, the key is to reformulate a discrete normal flux with local extremum principle structure, which is based on a modification of a second order linear scheme. In the construction of existing cell-centered finite volume schemes that maintain the discrete maximum principle, it is required to represent auxiliary unknowns as convex combinations of primary unknowns, which results in strong constraints on the smoothness of the mesh and diffusion coefficient. By contrast, our new scheme avoids this kind of constraints. Moreover, we will prove that there holds the DSEP for any solution of our scheme and there exists at least one solution preserving DSEP for our scheme. Furthermore, a modified Picard iteration with the Anderson acceleration (mP-AA) for solving the nonlinear scheme is proposed, and the nonlinear convergence of the modified Picard iteration is also proved. Finally, numerical examples are presented to show that the new scheme preserves DSEP and obtains second order accuracy, as well as the mP-AA method is effective.
  • Zhihao Ge, Wenlong He
    Journal of Computational Mathematics. 2026, 44(4): 1191-1218. https://doi.org/10.4208/jcm.2412-m2024-0188
    Abstract ( ) Download PDF   Knowledge map   Save
    In this paper, we propose and analyze a multiphysics finite element method for a nonlinear poroelasticity model. To more effectively capture the deformation and diffusion processes, we reformulate the original nonlinear fluid-solid coupling problem into a fluid-fluid coupling problem using a multiphysics approach. We then establish the growth, coercivity, and monotonicity properties of the nonlinear stress-strain relation, derive energy estimates, and use Schaefer’s fixed point theorem to prove the existence and uniqueness of the weak solution. Furthermore, we design a fully discrete time-stepping scheme – multiphysics finite element method with P2-P1-P1 elements for spatial variables and the backward Euler method for time variable. To handle nonlinearity, we employ the Newton iterative method, establish discrete energy laws and give the optimal convergence order error estimates. Finally, we show some numerical examples to verify the rationality of theoretical analysis and the proposed method has no “locking phenomenon”.