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

15 July 2024, Volume 42 Issue 4
    

  • Select all
    |
  • Yuping Zeng, Mingchao Cai, Liuqiang Zhong
    Journal of Computational Mathematics. 2024, 42(4): 911-931. https://doi.org/10.4208/jcm.2212-m2021-0231
    CSTR: 32031.14.jcm.2212-m2021-0231
    Abstract ( ) Download PDF   Knowledge map   Save
    A mixed finite element method is presented for the Biot consolidation problem in poroelasticity. More precisely, the displacement is approximated by using the Crouzeix-Raviart nonconforming finite elements, while the fluid pressure is approximated by using the node conforming finite elements. The well-posedness of the fully discrete scheme is established, and a corresponding priori error estimate with optimal order in the energy norm is also derived. Numerical experiments are provided to validate the theoretical results.
  • Debora Cores, Johanna Figueroa
    Journal of Computational Mathematics. 2024, 42(4): 932-954. https://doi.org/10.4208/jcm.2301-m2021-0313
    CSTR: 32031.14.jcm.2301-m2021-0313
    Abstract ( ) Download PDF   Knowledge map   Save
    Recently, the authors proposed a low-cost approach, named Optimization Approach for Linear Systems (OPALS) for solving any kind of a consistent linear system regarding the structure, characteristics, and dimension of the coefficient matrix A. The results obtained by this approach for matrices with no structure and with indefinite symmetric part were encouraging when compare with other recent and well-known techniques. In this work, we proposed to extend the OPALS approach for solving the Linear Least-Squares Problem (LLSP) and the Minimum Norm Linear System Problem (MNLSP) using any iterative lowcost gradient-type method, avoiding the construction of the matrices ATA or AAT, and taking full advantage of the structure and form of the gradient of the proposed nonlinear objective function in the gradient direction. The combination of those conditions together with the choice of the initial iterate allow us to produce a novel and efficient low-cost numerical scheme for solving both problems. Moreover, the scheme presented in this work can also be used and extended for the weighted minimum norm linear systems and minimum norm linear least-squares problems. We include encouraging numerical results to illustrate the practical behavior of the proposed schemes.
  • Chunxiong Zheng, Xianwei Wen, Jinyu Zhang, Zhenya Zhou
    Journal of Computational Mathematics. 2024, 42(4): 955-978. https://doi.org/10.4208/jcm.2301-m2022-0208
    CSTR: 32031.14.jcm.2301-m2022-0208
    Abstract ( ) Download PDF   Knowledge map   Save
    Asymptotic theory for the circuit envelope analysis is developed in this paper. A typical feature of circuit envelope analysis is the existence of two significantly distinct timescales: one is the fast timescale of carrier wave, and the other is the slow timescale of modulation signal. We first perform pro forma asymptotic analysis for both the driven and autonomous systems. Then resorting to the Floquet theory of periodic operators, we make a rigorous justification for first-order asymptotic approximations. It turns out that these asymptotic results are valid at least on the slow timescale. To speed up the computation of asymptotic approximations, we propose a periodization technique, which renders the possibility of utilizing the NUFFT algorithm. Numerical experiments are presented, and the results validate the theoretical findings.
  • Dongyang Shi, Houchao Zhang
    Journal of Computational Mathematics. 2024, 42(4): 979-998. https://doi.org/10.4208/jcm.2210-m2021-0337
    CSTR: 32031.14.jcm.2210-m2021-0337
    Abstract ( ) Download PDF   Knowledge map   Save
    The focus of this paper is on two novel linearized Crank-Nicolson schemes with nonconforming quadrilateral finite element methods (FEMs) for the nonlinear coupled Schrödinger-Helmholtz equations. Optimal L2 and H1 estimates of orders $\mathcal{O}$(h2 +τ2) and $\mathcal{O}$(h+τ2) are derived respectively without any grid-ratio condition through the following two keys. One is that a time-discrete system is introduced to split the error into the temporal error and the spatial error, which leads to optimal temporal error estimates of order $\mathcal{O}$(τ2) in L2 and the broken H1- norms, as well as the uniform boundness of numerical solutions in L- norm. The other is that a novel projection is utilized, which can iron out the difficulty of the existence of the consistency errors. This leads to derive optimal spatial error estimates of orders $\mathcal{O}$(h2) in L2-norm and $\mathcal{O}$(h) in the broken H1-norm under the H2 regularity of the solutions for the time-discrete system. At last, two numerical examples are provided to confirm the theoretical analysis. Here, h is the subdivision parameter, and τ is the time step.
  • Jiaofen Li, Lingchang Kong, Xuefeng Duan, Xuelin Zhou, Qilun Luo
    Journal of Computational Mathematics. 2024, 42(4): 999-1031. https://doi.org/10.4208/jcm.2211-m2021-0043
    CSTR: 32031.14.jcm.2211-m2021-0043
    Abstract ( ) Download PDF   Knowledge map   Save
    The truncated singular value decomposition has been widely used in many areas of science including engineering, and statistics, etc. In this paper, the original truncated complex singular value decomposition problem is formulated as a Riemannian optimization problem on a product of two complex Stiefel manifolds, a practical algorithm based on the generic Riemannian trust-region method of Absil et al. is presented to solve the underlying problem, which enjoys the global convergence and local superlinear convergence rate. Numerical experiments are provided to illustrate the efficiency of the proposed method. Comparisons with some classical Riemannian gradient-type methods, the existing Riemannian version of limited-memory BFGS algorithms in the MATLAB toolbox Manopt and the Riemannian manifold optimization library ROPTLIB, and some latest infeasible methods for solving manifold optimization problems, are also provided to show the merits of the proposed approach.
  • Suna Ma, Huiyuan Li, Zhimin Zhang, Hu Chen, Lizhen Chen
    Journal of Computational Mathematics. 2024, 42(4): 1032-1062. https://doi.org/10.4208/jcm.2304-m2022-0243
    CSTR: 32031.14.jcm.2304-m2022-0243
    Abstract ( ) Download PDF   Knowledge map   Save
    An efficient spectral-Galerkin method for eigenvalue problems of the integral fractional Laplacian on a unit ball of any dimension is proposed in this paper. The symmetric positive definite linear system is retained explicitly which plays an important role in the numerical analysis. And a sharp estimate on the algebraic system’s condition number is established which behaves as N4s with respect to the polynomial degree N, where 2s is the fractional derivative order. The regularity estimate of solutions to source problems of the fractional Laplacian in arbitrary dimensions is firstly investigated in weighted Sobolev spaces. Then the regularity of eigenfunctions of the fractional Laplacian eigenvalue problem is readily derived. Meanwhile, rigorous error estimates of the eigenvalues and eigenvectors are obtained. Numerical experiments are presented to demonstrate the accuracy and efficiency and to validate the theoretical results.
  • Yayun Fu, Dongdong Hu, Wenjun Cai, Yushun Wang
    Journal of Computational Mathematics. 2024, 42(4): 1063-1079. https://doi.org/10.4208/jcm.2302-m2020-0279
    CSTR: 32031.14.jcm.2302-m2020-0279
    Abstract ( ) Download PDF   Knowledge map   Save
    In the paper, we propose a novel linearly implicit structure-preserving algorithm, which is derived by combing the invariant energy quadratization approach with the exponential time differencing method, to construct efficient and accurate time discretization scheme for a large class of Hamiltonian partial differential equations (PDEs). The proposed scheme is a linear system, and can be solved more efficient than the original energy-preserving exponential integrator scheme which usually needs nonlinear iterations. Various experiments are performed to verify the conservation, efficiency and good performance at relatively large time step in long time computations.
  • Jian Lu, Huaxuan Hu, Yuru Zou, Zhaosong Lu, Xiaoxia Liu, Keke Zu, Lin Li
    Journal of Computational Mathematics. 2024, 42(4): 1080-1108. https://doi.org/10.4208/jcm.2301-m2022-0091
    CSTR: 32031.14.jcm.2301-m2022-0091
    Abstract ( ) Download PDF   Knowledge map   Save
    Low-dose computed tomography (LDCT) contains the mixed noise of Poisson and Gaussian, which makes the image reconstruction a challenging task. In order to describe the statistical characteristics of the mixed noise, we adopt the sinogram preprocessing as a standard maximum a posteriori (MAP). Based on the fact that the sinogram of LDCT has nonlocal self-similarity property, it exhibits low-rank characteristics. The conventional way of solving the low-rank problem is implemented in matrix forms, and ignores the correlations among similar patch groups. To avoid this issue, we make use of a nonlocal KroneckerBasis-Representation (KBR) method to depict the low-rank problem. A new denoising model, which consists of the sinogram preprocessing for data fidelity and the nonlocal KBR term, is developed in this work. The proposed denoising model can better illustrate the generative mechanism of the mixed noise and the prior knowledge of the LDCT. Numerical results show that the proposed denoising model outperforms the state-of-the-art algorithms in terms of peak-signal-to-noise ratio (PSNR), feature similarity (FSIM), and normalized mean square error (NMSE).
  • Guangqiang Lan, Yu Jiang
    Journal of Computational Mathematics. 2024, 42(4): 1109-1123. https://doi.org/10.4208/jcm.2302-m2022-0246
    CSTR: 32031.14.jcm.2302-m2022-0246
    Abstract ( ) Download PDF   Knowledge map   Save
    Convergence of modified truncated Euler-Maruyama (MTEM) method for stochastic differential equations (SDEs) with (1/2+α)-Hölder continuous diffusion coefficients are investigated in this paper. We prove that the MTEM method for SDE converges to the exact solution in Lq sense under given conditions. Two examples are provided to support our conclusions.
  • Hanzhang Hu, Yanping Chen, Jianwei Zhou
    Journal of Computational Mathematics. 2024, 42(4): 1124-1144. https://doi.org/10.4208/jcm.2302-m2022-0033
    CSTR: 32031.14.jcm.2302-m2022-0033
    Abstract ( ) Download PDF   Knowledge map   Save
    A two-grid finite element method with L1 scheme is presented for solving two-dimensional time-fractional nonlinear Schrödinger equation. The finite element solution in the L-norm are proved bounded without any time-step size conditions (dependent on spatialstep size). The classical L1 scheme is considered in the time direction, and the two-grid finite element method is applied in spatial direction. The optimal order error estimations of the two-grid solution in the Lp-norm is proved without any time-step size conditions. It is shown, both theoretically and numerically, that the coarse space can be extremely coarse, with no loss in the order of accuracy.
  • Pierluigi Amodio, Luigi Brugnano, Gianluca Frasca-Caccia, Felice Iavernaro
    Journal of Computational Mathematics. 2024, 42(4): 1145-1171. https://doi.org/10.4208/jcm.2301-m2022-0065
    CSTR: 32031.14.jcm.2301-m2022-0065
    Abstract ( ) Download PDF   Knowledge map   Save
    In this paper, we define arbitrarily high-order energy-conserving methods for Hamiltonian systems with quadratic holonomic constraints. The derivation of the methods is made within the so-called line integral framework. Numerical tests to illustrate the theoretical findings are presented.
  • Jiayin Zhai, Xiao Li, Zhijun Shen
    Journal of Computational Mathematics. 2024, 42(4): 1172-1196. https://doi.org/10.4208/jcm.2301-m2022-0177
    CSTR: 32031.14.jcm.2301-m2022-0177
    Abstract ( ) Download PDF   Knowledge map   Save
    This paper presents a cell-centered Godunov method based on staggered data distribution in Eulerian framework. The motivation is to reduce the intrinsic entropy dissipation of classical Godunov methods in the calculation of an isentropic or rarefaction flow. At the same time, the property of accurate shock capturing is also retained. By analyzing the factors that cause nonphysical entropy in the conventional Godunov methods, we introduce two velocities rather than a single velocity in a cell to reduce kinetic energy dissipation. A series of redistribution strategies are adopted to update subcell quantities in order to improve accuracy. Numerical examples validate that the present method can dramatically reduce nonphysical entropy increase.