Search

Purpose: The purpose of this study is to develop a new numerical algorithm to simulate the phase-field model. Design/methodology/approach: First, the derivative of the temporal direction is discretized by a second-order linearized finite difference scheme where it conserves the energy stability of the mathematical model. Then, the isogeometric collocation (IGC) method is used to approximate the derivative of spacial direction. The IGC procedure can be applied on irregular physical domains. The IGC method is constructed based upon the nonuniform rational B-splines (NURBS). Each curve and surface can be approximated by the NURBS. Also, a map will be defined to project the physical domain to a simple computational domain. In this procedure, the partial derivatives will be transformed to the new domain by the Jacobian and Hessian matrices. According to the mentioned procedure, the first- and second-order differential matrices are built. Furthermore, the pseudo-spectral algorithm is used to derive the first- and second-order nodal differential matrices. In the end, the Greville Abscissae points are used to the collocation method. Findings: In the numerical experiments, the efficiency and accuracy of the proposed method are assessed through two examples, demonstrating its performance on both rectangular and nonrectangular domains. Originality/value: This research work introduces the IGC method as a simulation technique for the phase-field crystal model. [ABSTRACT FROM AUTHOR]

In this manuscript, a new localized meshfree (meshless) approximation, i.e., a combination of the reproducing kernel particle (RKP) shape functions with the cardinal functions of the scalable polyharmonic spline radial kernel with polynomial augmentation (PHS+poly) is introduced. It is called the RKP+PHS+poly approximation, and its convergence rate is of order O (h m + 1) , where m is the total degree of polynomials. We apply this method to construct the spaces of trial and test functions in a Galerkin scheme of an extended version of the biological mathematical model in two dimensions describing the interactions between endothelial cells, fibronectin, angiogenic growth factors, and fasentin concentrations. By considering the row-sum method, the eigenvalue stability is also numerically carried out for the discrete equations corresponding to the obtained weak formulation. Accordingly, a semi-implicit form of the backward difference method of order 1 (SBDF1) has been utilized to approximate the weak form in time. We complete our numerical algorithm by solving the obtained full-discretized problem overtime via the biconjugate gradient stabilized (BiCGSTAB) solver with a proper preconditioner. Some simulation results are investigated by estimating the maximum velocity parameter of an enzymatic reaction from the experimental dose–response curve of fasentin to demonstrate the effect of using fasentin drug. [ABSTRACT FROM AUTHOR]

In this paper, we develop a new meshless numerical procedure for simulating the combustion model. To that end, we employ a local meshless collocation method according to the moving Taylor polynomial (MTP) approximation. The space derivative is approximated by using the local approach and then the Crank–Nicolson algorithm is utilized to approximate the time derivative. The stability and convergence of the time-discrete formulation are discussed, analytically and numerically. The Broyden method is applied to solve this nonlinear system. Since the size of the physical domain is large, we employ the non-overlapping domain decomposition method (DDM) to obtain a faster numerical algorithm. The local meshless approaches are efficient numerical techniques to simulate models in the fluid flow. The obtained results show that the proposed numerical formulation has efficient results for solving this mathematical model. [ABSTRACT FROM AUTHOR]