University of TabrizComputational Methods for Differential Equations2345-398211420230701Mathematical modeling of a nonlinear two-phase flow in a porous medium and the inflow of volatile oil to a well taking into account inertial effects6646751606810.22034/cmde.2023.48570.2032ENFikret AhmadaliAlievInstitute of Applied Mathematics, Baku State University, Baku, Azerbaijan.Institute of Information Technologies of the National Academy of
Sciences of Azerbaijan, Baku, Azerbaijan.MahammadJamalbayovSocar "Oil gas scientific research project" Institute, AZ1122, Hasan
Bey Zardabi, Baku, Azerbaijan.IlyasHasanovSocar "Oil gas scientific research project" Institute, AZ1122, Hasan
Bey Zardabi, Baku, Azerbaijan.NazimValiyevSocar Head OfficeNazileHajiyevaInstitute of Applied Mathematics, Baku State University, Baku, Azerbaijan.Journal Article20211026This paper discusses a semi-analytic solution for the volatile oil influx into the well on the base of Forchheimer flow law. The solution is developed employing the concept of binary model for the two-phase petroleum hydrocarbon system in view of phase transformations and interphase mass transfer. Algorithms are developed for calculating the volatile oil reservoir key performance indicators by applying the material balance equations, which take into account the compaction behavior of rocks. A computer simulator for the volatile oil reservoir is modeled, proceeding from these algorithms. The inertial effects on the development process of a volatile oil reservoir, the rocks of which are exposed to elastic deformation, are studied by this simulator. In regard thereto, the reservoir development process is simulated in two variants in conformity with the constant depression: in the first case, it is assumed that the filtration occurs according to Darcy's law, while in the second one, the process is considered on the base of Forchheimer equation. A comparison of the results of these options made it possible to demonstrate the nature of the inertial effects on the volatile oil reservoir key performance indicators.https://cmde.tabrizu.ac.ir/article_16068_690fec16072c4d6c289b4948131598c4.pdfUniversity of TabrizComputational Methods for Differential Equations2345-398211420230701An explicit split-step truncated Milstein method for stochastic differential equations6766951563910.22034/cmde.2022.51249.2132ENAmirHaghighiDepartment of Mathematics, Faculty of Science, Razi University, Kermanshah 67149, Iran.0000-0001-8928-9557Journal Article20220422In this paper, we propose an explicit split-step truncated Milstein method for stochastic differential equations (SDEs) with commutative noise. We discuss the mean-square convergence properties of the new method for numerical solutions of a class of highly nonlinear SDEs in a finite time interval. As a result, we show that the strong convergence rate of the new method can be arbitrarily close to one under some additional conditions. Finally, we use an illustrative example to highlight the advantages of our new findings in terms of both stability and accuracy compared to the results in Guo et al. (2018).https://cmde.tabrizu.ac.ir/article_15639_0374edec83e9f9ace24d55d461e071b5.pdfUniversity of TabrizComputational Methods for Differential Equations2345-398211420230701A mathematical study on non-linear ordinary differential equation for Magnetohydrodynamic flow of the Darcy-Forchheimer nanofluid6967151616610.22034/cmde.2023.55330.2302ENSathyamoorthySivasankariResearch Scholar, Research Centre and PG Department of Mathematics, The Madura College (Affiliated to Madurai Kamaraj University), Madurai,
Tamil Nadu, India.VembuAnanthaswamyResearch Scholar, Research Centre and PG Department of Mathematics, The Madura College (Affiliated to Madurai Kamaraj University), Madurai,
Tamil Nadu, India.0000-0002-2938-8745Journal Article20230208An analytical study is carried out to obtain the approximate solution for the Magnetohydrodynamic (MHD) flow issue of Darcy-Forchheimer nanofluid containing motile microorganisms having viscous dissipation effect through a non-linear extended sheet employing a new approximate analytical method namely Ananthaswamy-Sivasankari Method (ASM) and also Modified Homotopy Analysis method (MHAM). The derived analytical solution is given in explicit form and is compared with the numerical solution. The graphical results are interlined to reflect the effects of various physical parameters involved in the problem. The numerical computation of the Nusselt number, the local skin friction parameter, and the Sherwood number are compared and shown in the table. Faster convergence is acquired using this strategy. The solution obtained by this method is closer to the exact solution. Also, the solution is in the simplest and most explicit form. It is applicable for all initial and boundary value problems with non-zero boundary conditions. This method can be easily extended to solve other non-linear higher order boundary value problems in physical, chemical, and biological sciences.https://cmde.tabrizu.ac.ir/article_16166_6aff29b24ce2e5a1a4a91612eafa4143.pdfUniversity of TabrizComputational Methods for Differential Equations2345-398211420230701A numerical approximation for the solution of a time-fractional telegraph equation by the moving least squares approach7167261610210.22034/cmde.2023.55070.2286ENHaniyeHajinezhadDepartment of Mathematics, Payame Noor University, Tehran, Iran.Ali R.SoheiliDepartment of Applied Mathematics, Faculty of Mathematical Sciences, Ferdowsi University of Mashhad, Mashhad, Iran.Journal Article20230123This paper focuses on the numerical solution of the time-fractional telegraph equation in Caputo sense with $1 < \beta < 2$. The time-fractional telegraph equation models neutron transport inside the core of a nuclear reactor. The proposed numerical solution consists of two stages. First, the time-discretized scheme of this equation is obtained by the Crank-Nicolson method. The stability and convergence of results from the semi-discretized scheme are presented. In the second stage, the numerical approximation of the unknown function at specific points is achieved through the collocation method using the moving least square method. The numerical experiments analyze the impact of some parameters of the proposed method.https://cmde.tabrizu.ac.ir/article_16102_2ee8e45582c982ffdc5c8403a898f19e.pdfUniversity of TabrizComputational Methods for Differential Equations2345-398211420230701Multi-Soliton Solutions to the Generalized Boussinesq Equation of Tenth Order7277371626810.22034/cmde.2023.55261.2297ENBharathaKalegowdaPG Department of Mathematics, PG Studies and Research Centre, St. Philomena's College, Mysuru-570 015, India.RangarajanRaghavacharDepartment of Studies in Mathematics,University of Mysore, Manasagangothri, Mysuru-570 006, India.Journal Article20230206In the recent literature, many researchers are interested to apply standard computational methods for exact or numerical solutions of many classical nonlinear partial differential equations. Some leading methods are based on Lie group analysis, Painleve Analysis, G0/G expansion techniques, homotopy perturbation methods, and so on. The equations include complicated Navier-Stokes equation, Schrodinger equation, KdV-like equations, and so on. As a result, the glory of nonlinear dynamics can be witnessed through its applications in many fields namely: ocean engineering, plasma physics, optical communications, fluid dynamics, and much more. One of the significant observations is that whatever may be the order of nonlinear PDE, as far as the soliton and multisolitons of KdV like equation or Boussinesq equation are concerned Hirota’s method and tanh − coth method play a crucial role. The main result of the paper demonstrates that the above novel theme works well with the generalized Boussinesq equation of 10th order. In this paper, the Boussinesq equation of order ten is derived and its multi-soliton solutions are deduced by the Hirota’s method. The one soliton solution is reconfirmed using the tanh method. https://cmde.tabrizu.ac.ir/article_16268_00efb3e429268179e0ce9c0692f17616.pdfUniversity of TabrizComputational Methods for Differential Equations2345-398211420230701Extrapolated triangular splitting method to interval system of linear algebraic equations7387521634410.22034/cmde.2023.53101.2240ENGhodratEbadiFaculty of Mathematical Science, University of Tabriz, 51666-14766, Tabriz, Iran.BehnamMoradiFaculty of Mathematical Science, University of Tabriz, 51666-14766, Tabriz, Iran.Sergey P.SharyFederal Research Center for Information and Computational Technologies, (previously Institute of Computational Technologies SB RAS), Novosi-birsk, Russia.Journal Article20220820This paper presents an extrapolated triangular splitting method (\texttt{ETrnSplit}) to find the formal solution for the interval system of linear algebraic equations in which this method uses Kaucher interval arithmetic. Some numerical experiments are given to show the efficiency of this method.https://cmde.tabrizu.ac.ir/article_16344_9442f863525765a20f105f1fe66cc417.pdfUniversity of TabrizComputational Methods for Differential Equations2345-398211420230701Infinitely smooth multiquadric RBFs combined high-resolution compact discretization for nonlinear 2D elliptic PDEs on a scattered grid network7537751592610.22034/cmde.2023.51339.2136ENNavnitJhaDepartment of Mathematics, South Asian University, Maidan Garhi, New Delhi-110068, India.0000-0001-7884-8640ShikhaVermaDepartment of Mathematics, South Asian University, Maidan Garhi, New Delhi-110068, India.0000-0003-4700-4052Journal Article20220428Multiquadric radial basis functions combined with compact discretization to estimate solutions of two dimensions nonlinear elliptic type partial differential equations are presented. The scattered grid network with continuously varying step sizes helps tune the solution accuracies depending upon the location of high oscillation. The radial basis functions employing a nine-point grid network are used to improve the functional evaluations by compact formulation, and it saves memory space and computing time. A detailed description of convergence theory is presented to estimate the error bounds. The analysis is based on a strongly connected graph of the Jacobian matrix, and their monotonicity occurred in the scheme. It is shown that the present strategy improves the approximate solution values for the elliptic equations exhibiting a sharp changing character in a thin zone. Numerical simulations for the convection-diffusion equation, Graetz-Nusselt equation, Schr¨odinger equation, Burgers equation, and Gelfand-Bratu equation are reported to illustrate the utility of the new algorithm.https://cmde.tabrizu.ac.ir/article_15926_057663f18e1d38eba69b306cc69f4fcc.pdfUniversity of TabrizComputational Methods for Differential Equations2345-398211420230701Numerical solution of one-phase Stefan problem for the non-classical heat equation with a convective condition7767841631810.22034/cmde.2023.54177.2266ENSeyed ShojaeddinPishnamaz MohammadiDepartment of Applied Mathematics, Faculty of Mathematical Sciences, University of Tabriz, Tabriz, Iran.KarimIvazDepartment of Applied Mathematics, Faculty of Mathematical Sciences, University of Tabriz, Tabriz, Iran.Journal Article20221119A numerical technique for the solution of the one-phase Stefan problem for the non-classical heat equation with a convective condition is discussed. This approach is based on a scheme introduced in [16]. The compatibility and convergence of the method is proven. Numerical examples round out the discussion.https://cmde.tabrizu.ac.ir/article_16318_aee31b57f9ec89f35b13252f8a53d95b.pdfUniversity of TabrizComputational Methods for Differential Equations2345-398211420230701Existence, uniqueness, and finite-time stability of solutions for Ψ-Caputo fractional differential equations with time delay7858021607010.22034/cmde.2023.52613.2209ENNaoufelHatimeLMACS Laboratory, Sultan Moulay Slimane University, Beni Mellal, Morocco.SaidMellianiLMACS Laboratory, Sultan Moulay Slimane University, Beni Mellal, Morocco.AliEl MfadelLMACS Laboratory, Sultan Moulay Slimane University, Beni Mellal, Morocco.Superior School of Technology, Sultan Moulay Slimane University, Khenifra, Morocco.MhamedElomariLMACS Laboratory, Sultan Moulay Slimane University, Beni Mellal, Morocco.Journal Article20220720In this paper, we study the existence, uniqueness, and finite-time stability results for fractional delayed Newton cooling law equation involving Ψ-Caputo fractional derivatives of order α ∈ (0, 1). By using Banach fixed point theorem, Henry Gronwall type retarded integral inequalities, and some techniques of Ψ-Caputo fractional calculus, we establish the existence and uniqueness of solutions for our proposed model. Based on the heat transfer model, a new criterion for finite time stability and some estimated results of solutions with time delay are derived. In addition, we give some specific examples with graphs and numerical experiments to illustrate the obtained results. More importantly, the comparison of model predictions versus experimental data, classical model, and non-delayed model shows the effectiveness of our proposed model with a reasonable precision.https://cmde.tabrizu.ac.ir/article_16070_711a8407b7decab41a212c75c1af24ee.pdfUniversity of TabrizComputational Methods for Differential Equations2345-398211420230701A finite generating set of differential invariants for Lie symmetry group of the fifth-order KdV types8038101606910.22034/cmde.2023.53967.2262ENGhorbanaliHaghighatdoostDepartment of Mathematics, Azarbaijan Shahid Madani University, Tabriz, Iran.0000-0002-7979-6122MostafaBazghandiDepartment of Mathematics, Azarbaijan Shahid Madani University, Tabriz, Iran.FiroozPashaeiDepartment of Mathematics, University of Maragheh, P.O.Box 55181-83111, Maragheh, Iran.Journal Article20221030In this paper, we study the algebraic structure of differential invariants of a fifth-order KdV equation, known as the Kawahara KdV equation. Using the moving frames method, we locate a finite generating set of differential invariants, recurrence relations, and syzygies among the differential invariants generators of the equation. We prove that the differential invariant algebra of the equation can be generated by two first-order differential invariants.https://cmde.tabrizu.ac.ir/article_16069_e55e1b2c66684193891fbb52b83d4fd6.pdfUniversity of TabrizComputational Methods for Differential Equations2345-398211420230701On approximating eigenvalues and eigenfunctions of fractional order Sturm-Liouville problems8118211524010.22034/cmde.2022.52790.2221ENArezuAghazadehDepartment of Mathematics, Tabriz Branch, Islamic Azad University, Tabriz, Iran.0000-0003-2409-6740YaghoubMahmoudiDepartment of Mathematics, Tabriz Branch, Islamic Azad University, Tabriz, Iran.Journal Article20220728In this paper, the eigenvalues and corresponding eigenfunctions of a fractional order Sturm-Liouville problem (FSLP) are approximated by using the fractional differential transform method (FDTM), which is a generalization of the differential transform method (DTM). FDTM reduces the proposed fourth-order FSLP to a system of algebraic equations. The resulting coefficient matrix defines a characteristic polynomial which its roots correspond to the eigenvalues of FSLP. The obtained numerical results which are compared with the results of other papers confirm the efficiency of the method.https://cmde.tabrizu.ac.ir/article_15240_6c72b4827140cc917d3d318f2d9aadd2.pdfUniversity of TabrizComputational Methods for Differential Equations2345-398211420230701A second order numerical scheme for solving mixed type boundary value problems involving singular perturbation8228331607110.22034/cmde.2023.52094.2178ENSubal RanjanSahuDepartment of Mathematics, Larambha college, Bargarh, Orissa - 768102, India.JugalMohapatraDepartment of Mathematics, National Institute of Technology Rourkela, Odisha, India.0000-0001-5118-3933LoluguGovindaraoDepartment of Mathematics Amrita School of Engineering,
Amrita Vishwa Vidyapeetham, Coimbatore- 641112, India.Journal Article20220614A class of singularly perturbed mixed type boundary value problems is considered here in this work. The domain is partitioned into two subdomains. Convection-diffusion and reaction-diffusion problems are posed on the first and second subdomain, respectively. To approximate the problem, a hybrid scheme which consists of a second-order central difference scheme and a midpoint upwind scheme is constructed on Shishkin-type meshes. We have shown that the proposed scheme is second-order convergent in the maximum norm which is independent of the perturbation parameter. Numerical results are illustrated to support the theoretical findings.https://cmde.tabrizu.ac.ir/article_16071_07e7d00f6439054ecf01ce7711824570.pdfUniversity of TabrizComputational Methods for Differential Equations2345-398211420230701A Legendre Tau method for numerical solution of multi-order fractional mathematical model for COVID-19 disease8348501610310.22034/cmde.2023.53231.2245ENMarjanBidarianDepartment of Mathematics, Kerman Branch, Islamic Azad University, Kerman, Iran.HabibollahSaeediDepartment of Applied Mathematics, Faculty of Mathematics and Computer,
Shahid Bahonar University of Kerman, Kerman, Iran.0000-0001-5627-9678Mohammad RezaBaloochshahryariDepartment of Mathematics, Kerman Branch, Islamic Azad University, Kerman, Iran.Journal Article20220830In this paper, we describe a spectral Tau approach for approximating the solutions of a system of multi-order fractional differential equations which resulted from coronavirus disease mathematical modeling (COVID-19). The non-singular fractional derivative with a Mittag-Leffler kernel serves as the foundation for the fractional derivatives. Also, the operational matrix of fractional differentiation on the domain [0, a] is presented. Then, the convergence analysis of the proposed approximate approach is established and the error bounds are determined in a weighted L2 norm. Finally, by applying the Tau method, some of the important parameters in the model’s impact on the dynamics of the disease are graphically displayed for various values of the non-integer order of the ABC-derivative. https://cmde.tabrizu.ac.ir/article_16103_fe39a69d4cfb24087cc37cbe5fb9b065.pdfUniversity of TabrizComputational Methods for Differential Equations2345-398211420230701Existence of nonoscillatory solutions of second-order differential equations with mixed neutral term8518641615810.22034/cmde.2023.55524.2310ENOrhanOzdemirDepartment of Mathematics, Faculty of Arts and Sciences, Tokat
Gaziosmanpa\c{s}a University, 60240, Tokat, Turkey.0000-0003-1294-5346DemetBinbasiogluDepartment of Mathematics, Faculty of Arts and Sciences, Tokat
Gaziosmanpa\c{s}a University, 60240, Tokat, Turkey.0000-0001-7041-5277Journal Article20230221In this study, we aim to contribute to the increasing interest in functional differential equations by obtaining new existence theorems for non-oscillatory solutions of second-order neutral differential equations involving positive and negative terms which have not been performed in previous studies. We consider different cases for the ranges of the neutral coefficients, by utilizing the Banach contraction mapping principle. The applicability of the results is illustrated by several examples in the last section.https://cmde.tabrizu.ac.ir/article_16158_5a39b123d6c6b2d8ba4fd81f367cb154.pdfUniversity of TabrizComputational Methods for Differential Equations2345-398211420230701Solving Abel’s equations with the shifted Legendre polynomials8658751597510.22034/cmde.2023.52786.2222ENMaryamShahsavariDepartment of Mathematics, Faculty of Mathematics, Statistics and Computer Sciences, Semnan University, P. O. Box 35195-363, Semnan, Iran.LeilaTorkzadehDepartment of Mathematics, Faculty of Mathematics, Statistics and Computer Sciences, Semnan University, P. O. Box 35195-363, Semnan, Iran.0000-0002-2504-4048KazemNouriDepartment of Mathematics, Faculty of Mathematics, Statistics and Computer Sciences, Semnan University, P. O. Box 35195-363, Semnan, Iran.0000-0002-7922-5848Journal Article20220729In this article, a numerical method is presented to solve Abel’s equations. In the given method, the solution of the equation is found as a finite expansion of the shifted Legendre polynomials. To this end, the integral and differential parts of the equation are converted to vector-matrix representations. Therefore, the equation is converted to an algebraic system of the equations and by solving it, the solution of the equation is obtained. Further, the numerical example is given to illustrate the method’s efficiency.https://cmde.tabrizu.ac.ir/article_15975_cd59407bcfd2b0c8e2fa282de1e91a2a.pdf