University of GuilanJournal of Mathematical Modeling2345-394X12120240301Convergence analysis of compact finite difference method for the solution of anti-periodic boundary value problems115719710.22124/jmm.2023.25342.2251ENAbdol BaseerSaqibDepartment of Mathematical Sciences, Yazd University, Yazd, Iran0000-0002-8311-3267GhasemBarid LoghmaniDepartment of Mathematical Sciences, Yazd University, Yazd, IranMohammadHeydariDepartment of Mathematical Sciences, Yazd University, Yazd, IranJournal Article20230822The main objective of this paper is to introduce the fourth and sixth-order compact finite difference methods for solving anti-periodic boundary value problems. Compact finite difference formulas can approximate the derivatives of a function more accurately than the standard finite difference formulas for the same number of grid points. The convergence analysis of the proposed method is also investigated. This analysis shows how the error between the approximate and exact solutions decreases as the grid space is reduced. To validate the proposed method's accuracy and efficiency, some computational experiments are provided. Moreover, a comparison is performed between the standard and compact finite difference methods. The experiments indicate that the compact finite difference method is more accurate and efficient than the standard one.https://jmm.guilan.ac.ir/article_7197_e4973decc257e78ed4c45b66d552f215.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X12120240301Finding a time-dependent reaction coefficient of a nonlinear heat source in an inverse heat conduction problem1732720910.22124/jmm.2023.24413.2186ENKamalRashediDepartment of Mathematics, University of Science and Technology of Mazandaran, Behshahr, IranJournal Article20230429In this article, we propose an approximate technique for reconstructing a time-dependent reaction coefficient together with the surface heat flux histories and temperature distribution in a nonlinear inverse heat conduction problem (IHCP). We assume that the initial condition and the transient heat flux on the accessible boundary along with the temperature measured at specified interior locations in the domain of the problem are given as the input data. By applying the given measurements in a transformation, the main problem is reformulated as a certain parabolic problem and later a procedure based upon deploying the Ritz approximation along with the collocation method is applied which converts the problem to a nonlinear system of algebraic equations. Accurate numerical results in dealing with the exact initial and boundary data are obtained and regarding the perturbed boundary data, the regularization method based on cubic spline approximation is used, which results in obtaining stable numerical derivatives.https://jmm.guilan.ac.ir/article_7209_c9b0c293836ad4b8a2d46d10b4ffa046.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X12120240301A fitted mesh method for a class of two-parameter nonlinear singular perturbation problems3349725910.22124/jmm.2023.25201.2237ENManikandanMariappanDepartment of Mathematics, School of Engineering, Presidency University, Bengaluru-560 064, Karnataka, India0000-0001-8268-7800Journal Article20230808A class of two-parameter singularly perturbed nonlinear second order ordinary differential equations is considered in this article. A fitted mesh method which is a combination of finite difference scheme and a Shishkin mesh is developed to solve the problems. The method is proved to be essentially first order parameter independent convergent. Numerical experiments support the established theoretical results.https://jmm.guilan.ac.ir/article_7259_62f9f41e724e64d8e599552af1f42a96.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X12120240301Analysis of a coupled pair of Cahn-Hilliard equations with nondegenerate mobility5170726010.22124/jmm.2023.25558.2272ENGhufran A.Al-MusawiDepartment of Mathematics, College of Sciences, University of Basrah, Basrah, IraqAkil JassimHarfashDepartment of Mathematics, College of Sciences, University of Basrah, Basrah, Iraq0000-0002-3738-4242Journal Article20230915A mathematical analysis is performed for a system consisting of two coupled Cahn-Hilliard equations. These equations incorporate a diffusional mobility that depends on concentration. This modeling approach is often used to describe the process of phase separation in a thin layer of a binary liquid mixture covering a substrate, particularly when one of the components wets the substrate. The analysis establishes the existence of a weak formulation for this problem, which is supported by the use of a Lyapunov functional. Additionally, the analysis provides insights into the regularity properties of the weak formulation.https://jmm.guilan.ac.ir/article_7260_4fe3fd7d7504363e428bea9bbc71f036.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X12120240301A new numerical method for discretization of the nonlinear Klein-Gordon model arising in light waves7184728610.22124/jmm.2023.25115.2230ENHamidMesgaraniDepartment of Mathematics, Faculty of Science, Shahid Rajaee Teacher Training University, Tehran, 16785 -136, I. R. IranYonesEsmaeelzade AghdamDepartment of Mathematics, Faculty of Science, Shahid Rajaee Teacher Training University, Tehran, 16785 -136, I. R. Iran0000-0001-5109-1561EzzatollahDarabiDepartment of Mathematics, Faculty of Science, Shahid Rajaee Teacher Training University, Tehran, 16785 -136, I. R. IranJournal Article20230729Due to the importance of the generalized nonlinear Klein-Gordon equation (NL-KGE) in describing the behavior of light waves and nonlinear optical materials, including phenomena such as optical switching by manipulating the dispersion and nonlinearity of optical fibers and stable solitons, we explain the approximation of this model by evaluating different classical and fractional terms in this paper. To estimate the fundamental function, we use a first-order finite difference approach in the temporal direction and a collocation method based on Gegenbauer polynomials (GP) in the spatial direction to solve the NL-KGE model. Moreover, the stability and convergence analysis is proved by examining the order of the new method in the time direction as $\mathcal{O}( \delta t )$. To demonstrate the efficiency of this design, we presented numerical examples and made comparisons with other methods in the literature.https://jmm.guilan.ac.ir/article_7286_19b9d1f4d7853b95a673faf5b32e24aa.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X12120240301An improved extended block Arnoldi method for solving low-rank Lyapunov equation8598731810.22124/jmm.2023.25670.2281ENIliasAbdaouiENSA Oujda, Equipe MSN, Lab. LM2N, Université Mohammed Premier, Oujda, MoroccoJournal Article20230928We are interested in the numerical solution of the continuous-time Lyapunov equation. Generally, classical Krylov subspace methods for solving matrix equations use the Petrov-Galerkin condition to obtain projected equations from the original ones. The projected problems involves the restrictions of the coefficient matrices to a Krylov subspace. Alternatively, we propose a scheme based on the extended block Krylov subspace that leads to a smaller-scale equation, which also incorporates the restriction of the inverse of the Lyapunov equation's square coefficient. The effectiveness of this approach is experimentally confirmed, particularly in terms of the required CPU time.https://jmm.guilan.ac.ir/article_7318_dff6b430202ad958d6708c49eeebba7f.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X12120240301A fifth-order symmetrical weighted hybrid ENO-flux limiter scheme for traffic flow model on networks99115733910.22124/jmm.2023.24976.2220ENRooholahAbedianSchool of Engineering Science, College of Engineering, University of Tehran, Iran0000-0002-1739-5964Journal Article20230714In this research work, a fifth-order weighted essentially non-oscillatory (WENO) scheme is created for traffic flow problems on networks. Street systems can be numerically demonstrated as a graph, whose edges are a limited number of streets that connect at intersections. A scalar hyperbolic conservation law can portray the advancement on each street, and traffic distribution matrices are considered to define coupling conditions at the network intersections. In this paper, numerical results for road networks with rich solution structures will be presented. These numerical results show that the new proposed scheme in this paper can generate essentially non-oscillatory and high resolution solutions.https://jmm.guilan.ac.ir/article_7339_876197e1acbec9c8d3b252caa44fcbcf.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X12120240301Inverse eigenvalue problem of nonnegative matrices via unit lower triangular matrices (Part I)117130735210.22124/jmm.2023.21759.2092ENAlimohammadNazariDepartment of Mathematics, Arak University, P.O. Box 38156-8-8943, Arak, Iran0000-0002-3231-0340AtiyehNezamiDepartment of Mathematics, Arak University, P.O. Box 38156-8-8943, Arak, IranJournal Article20221215This paper uses unit lower triangular matrices to solve the nonnegative inverse eigenvalue problem for various sets of real numbers. This problem has remained unsolved for many years for $n \geq 5.$ The inverse of the unit lower triangular matrices can be easily calculated and the matrix similarities are also helpful to be able to solve this important problem to a considerable extent. It is assumed that in the given set of eigenvalues, the number of positive eigenvalues is less than or equal to the number of nonpositive eigenvalues to find a nonnegative matrix such that the given set is its spectrum.https://jmm.guilan.ac.ir/article_7352_3f3b6276cb29703ff071720a7048095a.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X12120240301A new $(I+P)$-like preconditioner for the SOR method for solving multi-linear systems with $ \mathcal{M} $-tensors131144736810.22124/jmm.2023.25368.2253ENAfsanehHasanpourDepartment of Mathematics, University of Sistan and Baluchestan, Zahedan, IranMaryamMojarrabDepartment of Mathematics, University of Sistan and Baluchestan, Zahedan, Iran0000-0002-9619-211XJournal Article20230823The use of preconditioning techniques has been shown to offer significant advantages in solving multi-linear systems involving nonsingular $\mathcal{M}$-tensors. In this paper, we introduce a new preconditioner that employs $(I+P)$-like preconditioning techniques, and give the proof of its convergence. We also present numerical examples and comparison results that demonstrate the superior efficiency of our preconditioner compared to both the original SOR method and the previously proposed preconditioned SOR method.https://jmm.guilan.ac.ir/article_7368_c3cef72c8af15df5f57d71122d50d240.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X12120240301Moore-Penrose inverse of an interval matrix and its application145155738210.22124/jmm.2023.24972.2219ENMarziehDehghani-MadisehDepartment of Mathematics, Faculty of Mathematical Sciences and Computer, Shahid Chamran University of Ahvaz, Ahvaz, IranJournal Article20230714 In this paper, we introduce the concept of Moore-Penrose inverse of a rectangular interval matrix based on a modified interval arithmetic. We determine the Moore-Penrose inverse in such a way that it satisfies all the four criteria similar to the real case. Also, we use the Moore-Penrose inverse for solving rectangular interval linear systems, algebraically.https://jmm.guilan.ac.ir/article_7382_d8dfa86be5aa3de5459e398d7491fa2a.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X12120240301A uniformly convergent computational method for singularly perturbed parabolic partial differential equation with integral boundary condition157175739410.22124/jmm.2023.24590.2198ENGetu MekonnenWondimuMathematics Department, Adama Raya University, Maichew, Ethiopia0000-0003-3730-4780Gemechis FileDuressaDepartment of Mathematics, Jimma University, Jimma, Ethiopia0000-0003-1889-4690Mesfin MekuriaWoldaregayDepartment of Applied Mathematics, Adama Science and Technology University, Adama, Ethiopia0000-0002-6555-7534Tekle GemechuDinkaApplied Mathematics, Adama Science and Technology University, Adama, Ethiopia0000-0003-0592-9508Journal Article20230524This paper presents a numerical method for a class of singularly perturbed parabolic partial differential equations with integral boundary conditions (IBC). The solution to the considered problem exhibits pronounced boundary layers on both the left and right sides of the spatial domain. To address this challenging problem, we propose the use of the implicit Euler method for time discretization and a finite difference method on a well-designed piecewise uniform Shishkin mesh for spatial discretization. The integral boundary condition is approximated using Simpson's $\frac{1}{3}$ rule. The presented method demonstrates almost second-order uniform convergence in the discretization of the spatial derivative and first-order convergence in the discretization of the time derivative. To validate the applicability and accuracy of the proposed method, two illustrative examples are employed. The computational results not only accurately reflect the theoretical estimations but also highlight the method's effectiveness in capturing the intricate features of singularly perturbed parabolic partial differential equations with integral boundary conditions.https://jmm.guilan.ac.ir/article_7394_1c073111c73b1a08462b46083684e3e6.pdfUniversity of GuilanJournal of Mathematical Modeling2345-394X12120240301An algorithm for multi-objective fuzzy linear programming problem with interval type-2 fuzzy numbers and ambiguity in parameters177197743210.22124/jmm.2023.25400.2258ENHassanMishmast NehiFaculty of Mathematics, University of Sistan and Baluchestan, Zahedan, IranShokouhSargolzaeiFaculty of Mathematics, University of Sistan and Baluchestan, Zahedan, IranJournal Article20230827Considering multiple criteria and objectives simultaneously in a single real-world problem and the fuzzy nature of this type of problem is of particular importance and application. As a result, multi-objective interval type-2 fuzzy linear programming problems have received much attention. However, there are few and a limited number of methods available for solving multi-objective interval type-2 triangular fuzzy linear programming problems with ambiguity-type imprecision (interval type-2 triangular fuzzy numbers) in almost all of the problem parameters. This research first considers a multi-objective interval type-2 fuzzy linear programming problem with ambiguity in all coefficients, in which, all problem coefficients are interval triangular fuzzy numbers. In addition, using the weighted sum method and the concept of nearest interval approximation, the problem is solved, and an example is provided.https://jmm.guilan.ac.ir/article_7432_499a93695db7ad5fe996b02050aed9ac.pdf