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

15 July 2023, Volume 41 Issue 4
    

  • Select all
    |
  • Yongxia Hao, Ting Li
    Journal of Computational Mathematics. 2023, 41(4): 551-568. https://doi.org/10.4208/jcm.2106-m2021-0050
    Abstract ( ) Download PDF   Knowledge map   Save
    In this paper, we present a method for generating Bézier surfaces from the boundary information based on a general second order functional and a third order functional associated with the triharmonic equation. By solving simple linear equations, the internal control points of the resulting Bézier surface can be obtained as linear combinations of the given boundary control points. This is a generalization of previous works on Plateau-Bézier problem, harmonic, biharmonic and quasi-harmonic Bézier surfaces. Some representative examples show the effectiveness of the presented method.
  • Zhiyun Yu, Dongyang Shi, Huiqing Zhu
    Journal of Computational Mathematics. 2023, 41(4): 569-587. https://doi.org/10.4208/jcm.2107-m2021-0114
    Abstract ( ) Download PDF   Knowledge map   Save
    A low order nonconforming mixed finite element method (FEM) is established for the fully coupled non-stationary incompressible magnetohydrodynamics (MHD) problem in a bounded domain in 3D. The lowest order finite elements on tetrahedra or hexahedra are chosen to approximate the pressure, the velocity field and the magnetic field, in which the hydrodynamic unknowns are approximated by inf-sup stable finite element pairs and the magnetic field by H1(?)-conforming finite elements, respectively. The existence and uniqueness of the approximate solutions are shown. Optimal order error estimates of L2(H1)-norm for the velocity field, L2(L2)-norm for the pressure and the broken L2(H1)-norm for the magnetic field are derived.
  • Weina Wang, Nannan Tian, Chunlin Wu
    Journal of Computational Mathematics. 2023, 41(4): 588-622. https://doi.org/10.4208/jcm.2108-m2021-0057
    Abstract ( ) Download PDF   Knowledge map   Save
    Two-phase image segmentation is a fundamental task to partition an image into foreground and background. In this paper, two types of nonconvex and nonsmooth regularization models are proposed for basic two-phase segmentation. They extend the convex regularization on the characteristic function on the image domain to the nonconvex case, which are able to better obtain piecewise constant regions with neat boundaries. By analyzing the proposed non-Lipschitz model, we combine the proximal alternating minimization framework with support shrinkage and linearization strategies to design our algorithm. This leads to two alternating strongly convex subproblems which can be easily solved. Similarly, we present an algorithm without support shrinkage operation for the nonconvex Lipschitz case. Using the Kurdyka-Lojasiewicz property of the objective function, we prove that the limit point of the generated sequence is a critical point of the original nonconvex nonsmooth problem. Numerical experiments and comparisons illustrate the effectiveness of our method in two-phase image segmentation.
  • Xiuhui Guo, Lulu Tian, Yang Yang, Hui Guo
    Journal of Computational Mathematics. 2023, 41(4): 623-642. https://doi.org/10.4208/jcm.2108-m2021-0143
    Abstract ( ) Download PDF   Knowledge map   Save
    In this paper, we apply local discontinuous Galerkin (LDG) methods for pattern formation dynamical model in polymerizing actin flocks. There are two main difficulties in designing effective numerical solvers. First of all, the density function is non-negative, and zero is an unstable equilibrium solution. Therefore, negative density values may yield blow-up solutions. To obtain positive numerical approximations, we apply the positivitypreserving (PP) techniques. Secondly, the model may contain stiff source. The most commonly used time integration for the PP technique is the strong-stability-preserving Runge-Kutta method. However, for problems with stiff source, such time discretizations may require strictly limited time step sizes, leading to large computational cost. Moreover, the stiff source any trigger spurious filament polarization, leading to wrong numerical approximations on coarse meshes. In this paper, we combine the PP LDG methods with the semi-implicit Runge-Kutta methods. Numerical experiments demonstrate that the proposed method can yield accurate numerical approximations with relatively large time steps.
  • Xiaoqiang Yan, Xu Qian, Hong Zhang, Songhe Song, Xiujun Cheng
    Journal of Computational Mathematics. 2023, 41(4): 643-662. https://doi.org/10.4208/jcm.2109-m2021-0020
    Abstract ( ) Download PDF   Knowledge map   Save
    Block boundary value methods (BBVMs) are extended in this paper to obtain the numerical solutions of nonlinear delay-differential-algebraic equations with singular perturbation (DDAESP). It is proved that the extended BBVMs in some suitable conditions are globally stable and can obtain a unique exact solution of the DDAESP. Besides, whenever the classic Lipschitz conditions are satisfied, the extended BBVMs are preconsistent and pth order consistent. Moreover, through some numerical examples, the correctness of the theoretical results and computational validity of the extended BBVMs is further confirmed.
  • Yidan Geng, Minghui Song, Mingzhu Liu
    Journal of Computational Mathematics. 2023, 41(4): 663-682. https://doi.org/10.4208/jcm.2109-m2021-0116
    Abstract ( ) Download PDF   Knowledge map   Save
    In this paper, we consider the stochastic differential equations with piecewise continuous arguments (SDEPCAs) in which the drift coefficient satisfies the generalized one-sided Lipschitz condition and the diffusion coefficient satisfies the linear growth condition. Since the delay term t-[t] of SDEPCAs is not continuous and differentiable, the variable substitution method is not suitable. To overcome this difficulty, we adopt new techniques to prove the boundedness of the exact solution and the numerical solution. It is proved that the truncated Euler-Maruyama method is strongly convergent to SDEPCAs in the sense of Lq(q ≥ 2). We obtain the convergence order with some additional conditions. An example is presented to illustrate the analytical theory.
  • Ziang Chen, Andre Milzarek, Zaiwen Wen
    Journal of Computational Mathematics. 2023, 41(4): 683-716. https://doi.org/10.4208/jcm.2110-m2020-0317
    Abstract ( ) Download PDF   Knowledge map   Save
    We propose a trust-region type method for a class of nonsmooth nonconvex optimization problems where the objective function is a summation of a (probably nonconvex) smooth function and a (probably nonsmooth) convex function. The model function of our trust-region subproblem is always quadratic and the linear term of the model is generated using abstract descent directions. Therefore, the trust-region subproblems can be easily constructed as well as efficiently solved by cheap and standard methods. When the accuracy of the model function at the solution of the subproblem is not sufficient, we add a safeguard on the stepsizes for improving the accuracy. For a class of functions that can be "truncated", an additional truncation step is defined and a stepsize modification strategy is designed. The overall scheme converges globally and we establish fast local convergence under suitable assumptions. In particular, using a connection with a smooth Riemannian trust-region method, we prove local quadratic convergence for partly smooth functions under a strict complementary condition. Preliminary numerical results on a family of ${\ell _1}$-optimization problems are reported and demonstrate the efficiency of our approach.
  • Yazid Dendani, Radouen Ghanem
    Journal of Computational Mathematics. 2023, 41(4): 717-740. https://doi.org/10.4208/jcm.2110-m2021-0131
    Abstract ( ) Download PDF   Knowledge map   Save
    In this paper we deal with the convergence analysis of the finite element method for an elliptic penalized unilateral obstacle optimal control problem where the control and the obstacle coincide. Error estimates are established for both state and control variables. We apply a fixed point type iteration method to solve the discretized problem.
    To corroborate our error estimations and the efficiency of our algorithms, the convergence results and numerical experiments are illustrated by concrete examples.
  • Ming-Jun Lai, Jiaxin Xie, Zhiqiang Xu
    Journal of Computational Mathematics. 2023, 41(4): 741-770. https://doi.org/10.4208/jcm.2201-m2021-0130
    Abstract ( ) Download PDF   Knowledge map   Save
    Graph sparsification is to approximate an arbitrary graph by a sparse graph and is useful in many applications, such as simplification of social networks, least squares problems, and numerical solution of symmetric positive definite linear systems. In this paper, inspired by the well-known sparse signal recovery algorithm called orthogonal matching pursuit (OMP), we introduce a deterministic, greedy edge selection algorithm, which is called the universal greedy approach (UGA) for the graph sparsification problem. For a general spectral sparsification problem, e.g., the positive subset selection problem from a set of m vectors in $\mathbb{R}{^n}$, we propose a nonnegative UGA algorithm which needs O(mn2 + n3/∈2) time to find a $\frac{{1 + \in/\beta }}{{1-\in/\beta }}$-spectral sparsifier with positive coefficients with sparsity at most $\left[{\frac{n}{{{ \in ^2}}}} \right]$, where β is the ratio between the smallest length and largest length of the vectors. The convergence of the nonnegative UGA algorithm is established. For the graph sparsification problem, another UGA algorithm is proposed which can output a $\frac{{1 + O/(\in)}}{{1-O/(\in)}}$-spectral sparsifier with $\left[{\frac{n}{{{ \in ^2}}}} \right]$ edges in O(m+n2/∈2) time from a graph with m edges and n vertices under some mild assumptions. This is a linear time algorithm in terms of the number of edges that the community of graph sparsification is looking for. The best result in the literature to the knowledge of the authors is the existence of a deterministic algorithm which is almost linear, i.e. O(m1+o(1)) for some o(1)=O($\frac{{{{(\log \log (m))}^{2/3}}}}{{{{\log }^{1/3}}(m)}}$). Finally, extensive experimental results, including applications to graph clustering and least squares regression, show the effectiveness of proposed approaches.
  • Weizhu Bao, Quan Zhao
    Journal of Computational Mathematics. 2023, 41(4): 771-796. https://doi.org/10.4208/jcm.2205-m2021-0237
    Abstract ( ) Download PDF   Knowledge map   Save
    We propose an accurate and energy-stable parametric finite element method for solving the sharp-interface continuum model of solid-state dewetting in three-dimensional space. The model describes the motion of the film/vapor interface with contact line migration and is governed by the surface diffusion equation with proper boundary conditions at the contact line. We present a weak formulation for the problem, in which the contact angle condition is weakly enforced. By using piecewise linear elements in space and backward Euler method in time, we then discretize the formulation to obtain a parametric finite element approximation, where the interface and its contact line are evolved simultaneously. The resulting numerical method is shown to be well-posed and unconditionally energystable. Furthermore, the numerical method is generalized to the case of anisotropic surface energies in the Riemannian metric form. Numerical results are reported to show the convergence and efficiency of the proposed numerical method as well as the anisotropic effects on the morphological evolution of thin films in solid-state dewetting.