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

14 September 2025, Volume 46 Issue 3
    

  • Select all
    |
  • Lv Minrui, Xu Xianmin, Lu Benzhuo
    Journal on Numerica Methods and Computer Applications. 2025, 46(3): 165-188. https://doi.org/10.12288/szjs.s2024-0990
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The effect of an oscillating external electric field on ion transport in nanoscale channels is investigated within the framework of the classical Poisson-Nernst-Planck (PNP) model. Three cases are analyzed by multiscale method: (1) an externally applied time-oscillating electric field along the channel direction; (2) a time-oscillating electric field coupled with a spatially periodic electric field along the channel; (3) an oscillating electric field applied along the channel direction in the presence of periodic surface charge distributions on the channel walls. An effective model for high oscillation frequencies is derived by considering the leading order approximation, which shows that the ion distribution and average transport properties within the channel depend only on the time and space averaged properties of the external oscillating electric field. Numerical simulations on a two-dimensional single-channel nanoscale model confirm the validity of these analytical results.
  • Nian Chenyu, Bao Wendi, Deng Shuaihao, Liu Shuaidong, Wang Dongrui
    Journal on Numerica Methods and Computer Applications. 2025, 46(3): 189-202. https://doi.org/10.12288/szjs.s2024-0976
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In this paper, firstly, inspired by the ideas of two algorithms: the greedy randomized Kaczmarz algorithm and the geometric probability randomized Kaczmarz algorithm for solving linear equations, two novel greedy randomized algorithms for solving matrix equations are proposed. Then, the convergence of these two methods in matrix equations is proved based on important inequalities. Finally, the numerical experiments are implemented to verify the effectiveness of the proposed methods. The numerical results show that the geometric probability randomized Kaczmarz algorithm outperforms the greedy randomized Kaczmarz algorithm for large-scale systems.
  • Chen Xiwen, Xiao Lifen, Ke Yifen, Wen Shuhong
    Journal on Numerica Methods and Computer Applications. 2025, 46(3): 203-213. https://doi.org/10.12288/szjs.s2024-0985
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    This paper presents a modulus-based matrix splitting iteration method for solving a class of vertical linear complementarity problems. By reformulating the vertical linear complementarity problem as an equivalent nonlinear system of equations, a new class of modulus-based matrix splitting iteration methods is created, and the convergence of the algorithm is proved under certain conditions. Finally, two numerical examples are provided to demonstrate the effectiveness of the proposed algorithm.
  • Su Ziyao, Zhang Yan, Zhu Jun
    Journal on Numerica Methods and Computer Applications. 2025, 46(3): 214-234. https://doi.org/10.12288/szjs.s2024-0991
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    In this paper, a new hybridization strategy for the sixth order unequal-sized weighted essentially non-oscillatory (HUS-WENO) scheme is designed for solving the nonlinear degenerate parabolic equations on structured meshes. This HUS-WENO scheme only uses the information defined on three unequal-sized stencils to obtain the optimal sixth order accuracy in smooth regions while maintaining stable, non-oscillatory and sharp discontinuity transitions. In addition, the linear weights in the proposed scheme are artificially set to be any random positive numbers with the only requirement that their sum equals one. To reduce CPU time, a hybridization strategy is designed based on the reconstruction polynomial of the six-point stencil in US-WENO scheme, which can accurately, efficiently, and automatically identify the troubled cells, and dose not contain any artificial parameters related to the problems. Finally, some numerical examples are used to verify the performance of the presented WENO scheme in various aspects, such as computational efficiency, low dissipation characteristics, shock capture ability, etc.
  • Wang Tao, Zeng Ling, Ren Wuyue
    Journal on Numerica Methods and Computer Applications. 2025, 46(3): 235-249. https://doi.org/10.12288/szjs.j2025-0994
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    The accurate simulation of Fluid-Structure Interaction (FSI) is important for the safety design and operation of nuclear reactors. This paper investigates the simulation of FSI by using a two-way coupled FSI solver based on the high-order discontinuous Galerkin method, which is implemented by the open-source software ExaDG. The solver is verified through the examination of a cantilever beam vibrating in the quiescent water, where the simulated vibration frequencies and damping ratios are compared with theoretical predictions. Furthermore, a numerical example involving turbulent flow is presented. This example employs the Large Eddy Simulation (LES) method to resolve the large scale eddies and utilizes the Synthetic Eddy Method (SEM) to generate inlet boundary conditions, showcasing the potential for applications in engineering and scientific research.
  • Liu Kaiyang, Yan Fuyou, Li Yonghui
    Journal on Numerica Methods and Computer Applications. 2025, 46(3): 250-262. https://doi.org/10.12288/szjs.s2025-0995
    Abstract ( ) Download PDF ( )   Knowledge map   Save
    Based on the classical spectral collocation method for solving boundary value problems of partial differential equations, the integral and differential operators of fractional order Chebyshev polynomials with integer and fractional order are deduced in vector form, and the solution function that satisfy the consolidation equation and its boundary conditions is established. After substituting the function into the consolidation equation, a system of algebraic equations is deduced with the help of the spectral collocation method. Finally, the numerical algorithm for solving the one-dimensional fractional viscoelastic consolidation problem is presented. It is a direct method for solving the consolidation equation in the time domain, and the excess pore water pressure and the effective stress, as well as settlement are the explicit functions of time, those are composed of finite term series, and it is a useful mathematical tool for solving problems of fractional viscoelastic consolidation with variable loading or multi stage loading. Some numerical examples are shown to illustrate the accuracy and effective of the proposed numerical algorithm.