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

14 September 2026, Volume 47 Issue 3
    

  • Select all
    |
  • Fan Jie, Huang Zhongyi, Wu Hao, Yu Bin, Zhang Weiming
    Journal on Numerica Methods and Computer Applications. 2026, 47(3): 299-321. https://doi.org/10.12288/szjs.s2025-1025
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In this paper, we propose the optimal production demand transport model and the optimal production demand transport with fixed mass model. These models extend the classical optimal transport model and the optimal production transport model to address the total cost minimization problem when both production and demand quantities can be adjusted within certain ranges. For example, these models can be applied to power generation, transport and demand scheduling in electricity markets. To handle the inequality constraints, we introduce slack variables and employ a multiple regularization method. By leveraging Lagrangian duality, we develop a generalized alternating Sinkhorn algorithm to solve the proposed models. Numerical experiments demonstrate that the generalized alternating Sinkhorn algorithm outperforms traditional linear programming methods and the iterative Bregman projections in terms of computational efficiency, accuracy, and memory usage.
  • Zhang Chaobing, Luo Renkai, Tang Qinglin
    Journal on Numerica Methods and Computer Applications. 2026, 47(3): 322-334. https://doi.org/10.12288/szjs.s2025-1068
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    This paper investigates high-order numerical methods for solving the nonlinear Schrödinger equation (NLSE) in the case of bright solitons. Using Dirichlet boundary conditions to truncate the bright-soliton domain, we apply a compact finite-difference scheme for spatial discretization and combine the Crank-Nicolson and IMEX Runge-Kutta methods for time discretization, constructing schemes of order $O(\tau^2 +h^4)$ and $O(\tau^4 +h^4)$. Numerical experiments confirm the expected convergence rates and successfully reproduce interaction dynamics of bright solitons.
  • Wang Xinyue, Chen Yumei
    Journal on Numerica Methods and Computer Applications. 2026, 47(3): 335-350. https://doi.org/10.12288/szjs.s2025-1066
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    To numerically solve singularly perturbed Volterra-Fredholm integro-differential equations, a hybrid numerical method is proposed by combining the central difference scheme, central upwind scheme, and Legendre wavelet method. First, the computational domain is discretized via a Shishkin-type non-uniform mesh, which partitions the domain into a boundary layer region and a non-boundary layer region. Subsequently, a region-adaptive discretization strategy is applied to the differential terms: the central difference scheme is adopted for the differential component in the boundary layer region, whereas the central upwind scheme is utilized for that in the non-boundary layer region. Meanwhile, the Fredholm and Volterra integral terms in the governing equation are uniformly discretized using the Legendre wavelet method, leading to the fully discrete formulation of the singularly perturbed Volterra-Fredholm integro-differential equation. The numerical solution of the original problem is then obtained by solving this discrete system. In addition, a rigorous uniform error estimate is derived for the proposed hybrid method, and this estimate is proven to be independent of the perturbation parameter. Finally, numerical experiments are conducted to validate the effectiveness of the method. Comparative studies against existing hybrid methods reported in the literature demonstrate that the proposed method achieves higher computational accuracy and a faster convergence rate as the number of mesh nodes increases.
  • Liu Jichuan
    Journal on Numerica Methods and Computer Applications. 2026, 47(3): 351-366. https://doi.org/10.12288/szjs.s2025-1007
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In this paper, an inverse source problem of the Poisson equation is investigated. We want to seek reconstruction algorithms to detect the hidden source within a body from Cauchy data on the boundary. Our goal is to detect the location, the size and the shape of the hidden source. This problem is ill-posed, regularization techniques should be employed to obtain the regularized solution. Numerical examples show that our proposed algorithms are valid and effective.
  • He Jingjin, Lu Changna, Miao Yuxing
    Journal on Numerica Methods and Computer Applications. 2026, 47(3): 367-382. https://doi.org/10.12288/szjs.s2025-1009
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    This paper proposes a Convolutional Neural Network (CNN)-based WENO scheme, named WENO-CNN, which aims to optimize the nonlinear weights in the WENO framework through machine learning techniques to enhance the accuracy and robustness of numerical simulations. Traditional WENO-JS schemes address discontinuities by weighting different reconstruction polynomials. To enable weight adaptation to diverse flow conditions, WENO-CNN employs a CNN to predict weights at identified troubled cells, allowing dynamic weight adjustments based on flow field characteristics. This approach improves numerical performance in scenarios involving strong discontinuities and complex flow structures. A specialized CNN architecture is developed for flow field problems, incorporating residual layers with skip connections. The network design includes an input layer that captures essential flow field features, ReLU activation functions, a customized loss function, and the Adam optimizer for training. Extensive validation through five classical test cases governed by the Euler equations demonstrates that the WENO-CNN scheme achieves higher shock wave resolution with refined detail capture while maintaining stability in smooth regions. The scheme effectively suppresses numerical oscillations, showing significant improvements over conventional methods in both shock and discontinuity flow simulations.
  • Piao Yuhao, Li Yushan, Wang Huimin
    Journal on Numerica Methods and Computer Applications. 2026, 47(3): 383-401. https://doi.org/10.12288/szjs.s2025-1010
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Based on terminal data, this paper identifies the spatial source term in the distributed-order time-space fractional diffusion equation. First, the existence and uniqueness of the solution to the forward problem are established, and the uniqueness of the solution to the inverse problem is proven by using the regularity of the solution to the forward problem. Subsequently, the inverse problem is transformed into a variational problem using the Tikhonov regularization method, and the optimal perturbation method is employed to solve the minimizer of the variational problem, thereby obtaining an approximate solution depending solely on the spatial source term. In the iterative process of solving the forward problem, the matrix transformation technique combined with the finite difference algorithm is adopted. Finally, numerical examples in one and two dimensions are provided to verify the effectiveness and stability of the proposed algorithm in determining the spatial source term.
  • Wang Zhiyuan, Tang Lingyan
    Journal on Numerica Methods and Computer Applications. 2026, 47(3): 402-419. https://doi.org/10.12288/szjs.s2025-1011
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The rotating shallow water equations with Coriolis force serve as a fundamental mathematical model in geophysical fluid dynamics for simulating large-scale oceanic and atmospheric circulation. The development of high-order well-balanced numerical methods for these equations is of critical importance. In this study, we propose a Coriolis potential function to equivalently reformulate the Coriolis force term into a bottom topography term. Through a source term splitting technique, we establish a strict correspondence between the source term and flux derivatives under equilibrium conditions. By modifying numerical fluxes and implementing identical nonlinear reconstruction operators for conserved variables, bottom topography, and the Coriolis potential function, we develop a well-balanced weighted compact nonlinear scheme (WCNS) with fifth-order spatial accuracy. Rigorous theoretical analysis demonstrates the scheme's exact preservation of geostrophic balance for both 1D and 2D shallow water equations. Comprehensive numerical experiments validate its high-order accuracy, well-balanced properties, and ability of capturing small perturbations.
  • Yin Sile, Yan Yun, Zhang Fangmei
    Journal on Numerica Methods and Computer Applications. 2026, 47(3): 420-430. https://doi.org/10.12288/szjs.s2025-1014
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    For the one-dimensional Stefan problem, a linear finite difference scheme with second-order accuracy in both time and space is constructed using the Landau transformation and second-order explicit extrapolation techniques. The scheme's unique solvability and stability are proven by the discrete energy method and the freezing coefficient method. Numerical experimental results are consistent with theoretical analysis.
  • Tao Yanling, Xiong Jiang, Chen Wanshun, Liu Jinkui
    Journal on Numerica Methods and Computer Applications. 2026, 47(3): 431-445. https://doi.org/10.12288/szjs.s2025-1020
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Based on the hybrid Liu-Storey-Dai-Yuan (HLSDY) method proposed by Liu and Li[9] in Euclidean space, this paper extends it to Riemannian manifolds by using retraction and vector transport mapping, and proposes a novel Riemannian Liu-Storey-Dai-Yuan (RLSDY) method. By incorporating the Powell restart condition, the proposed method has the sufficient descent property and the global convergence under Riemannian strong Wolfe conditions. Numerical experiments demonstrate its high efficiency in solving classical Riemannian optimization problems.
  • Yang Yuehan, Long Xiaoxiao, Yu Jiajia
    Journal on Numerica Methods and Computer Applications. 2026, 47(3): 446-462. https://doi.org/10.12288/szjs.s2025-1027
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Based on the improved BFGS methods proposed by Yuan and Wei[15] and Andrei, Neculaii[7], this paper proposes an improved single-parameter BFGS method. By introducing an adjustable parameter, the update formula proposed by Yuan is adaptively scaled to reduce the influence of large eigenvalues of the Hessian approximation matrix of the objective function, and a new update formula is constructed based on the quasi-Newton equation. The theoretical properties of this method for convex functions are analyzed. Numerical experiments are carried out using 73 unconstrained optimization test functions to verify the effectiveness of the algorithm, and its potential value in the field of hydrological engineering is also explored.