Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridyna...Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.展开更多
The paper presents the comparative study on numerical methods of Euler method,Improved Euler method and fourth-order Runge-Kutta method for solving the engineering problems and applications.The three proposed methods ...The paper presents the comparative study on numerical methods of Euler method,Improved Euler method and fourth-order Runge-Kutta method for solving the engineering problems and applications.The three proposed methods are quite efficient and practically well suited for solving the unknown engineering problems.This paper aims to enhance the teaching and learning quality of teachers and students for various levels.At each point of the interval,the value of y is calculated and compared with its exact value at that point.The next interesting point is the observation of error from those methods.Error in the value of y is the difference between calculated and exact value.A mathematical equation which relates various functions with its derivatives is known as a differential equation.It is a popular field of mathematics because of its application to real-world problems.To calculate the exact values,the approximate values and the errors,the numerical tool such as MATLAB is appropriate for observing the results.This paper mainly concentrates on identifying the method which provides more accurate results.Then the analytical results and calculates their corresponding error were compared in details.The minimum error directly reflected to realize the best method from different numerical methods.According to the analyses from those three approaches,we observed that only the error is nominal for the fourth-order Runge-Kutta method.展开更多
Our purpose of this paper is to apply the improved Kudryashov method for solving various types of nonlinear fractional partial differential equations. As an application, the time-space fractional Korteweg-de Vries-Bur...Our purpose of this paper is to apply the improved Kudryashov method for solving various types of nonlinear fractional partial differential equations. As an application, the time-space fractional Korteweg-de Vries-Burger (KdV-Burger) equation is solved using this method and we get some new travelling wave solutions. To acquire our purpose a complex transformation has been also used to reduce nonlinear fractional partial differential equations to nonlinear ordinary differential equations of integer order, in the sense of the Jumarie’s modified Riemann-Liouville derivative. Afterwards, the improved Kudryashov method is implemented and we get our required reliable solutions where the results are justified by mathematical software Maple-13.展开更多
Laplace transform is one of the powerful tools for solving differential equations in engineering and other science subjects.Using the Laplace transform for solving differential equations,however,sometimes leads to sol...Laplace transform is one of the powerful tools for solving differential equations in engineering and other science subjects.Using the Laplace transform for solving differential equations,however,sometimes leads to solutions in the Laplace domain that are not readily invertible to the real domain by analyticalmeans.Thus,we need numerical inversionmethods to convert the obtained solution fromLaplace domain to a real domain.In this paper,we propose a numerical scheme based on Laplace transform and numerical inverse Laplace transform for the approximate solution of fractal-fractional differential equations with orderα,β.Our proposed numerical scheme is based on three main steps.First,we convert the given fractal-fractional differential equation to fractional-differential equation in Riemann-Liouville sense,and then into Caputo sense.Secondly,we transformthe fractional differential equation in Caputo sense to an equivalent equation in Laplace space.Then the solution of the transformed equation is obtained in Laplace domain.Finally,the solution is converted into the real domain using numerical inversion of Laplace transform.Three inversion methods are evaluated in this paper,and their convergence is also discussed.Three test problems are used to validate the inversion methods.We demonstrate our results with the help of tables and figures.The obtained results show that Euler’s and Talbot’s methods performed better than Stehfest’s method.展开更多
To maximize energy profit with the participation of electricity,natural gas,and district heating networks in the day-ahead market,stochastic scheduling of energy hubs taking into account the uncertainty of photovoltai...To maximize energy profit with the participation of electricity,natural gas,and district heating networks in the day-ahead market,stochastic scheduling of energy hubs taking into account the uncertainty of photovoltaic and wind resources,has been carried out.This has been done using a new meta-heuristic algorithm,improved artificial rabbits optimization(IARO).In this study,the uncertainty of solar and wind energy sources is modeled using Hang’s two-point estimating method(TPEM).The IARO algorithm is applied to calculate the best capacity of hub energy equipment,such as solar and wind renewable energy sources,combined heat and power(CHP)systems,steamboilers,energy storage,and electric cars in the day-aheadmarket.The standard ARO algorithmis developed to mimic the foraging behavior of rabbits,and in this work,the algorithm’s effectiveness in avoiding premature convergence is improved by using the dystudynamic inertia weight technique.The proposed IARO-based scheduling framework’s performance is evaluated against that of traditional ARO,particle swarm optimization(PSO),and salp swarm algorithm(SSA).The findings show that,in comparison to previous approaches,the suggested meta-heuristic scheduling framework based on the IARO has increased energy profit in day-ahead electricity,gas,and heating markets by satisfying the operational and energy hub limitations.Additionally,the results show that TPEM approach dependability consideration decreased hub energy’s profit by 8.995%as compared to deterministic planning.展开更多
In this work, we apply the Zhou’s method [1] or differential transformation method (DTM) for solving the Euler equidimensional equation. The Zhou’s method may be considered as alternative and efficient for finding t...In this work, we apply the Zhou’s method [1] or differential transformation method (DTM) for solving the Euler equidimensional equation. The Zhou’s method may be considered as alternative and efficient for finding the approximate solutions of initial values problems. We prove superiority of this method by applying them on the some Euler type equation, in this case of order 2 and 3 [2]. The power series solution of the reduced equation transforms into an approximate implicit solution of the original equations. The results agreed with the exact solution obtained via transformation to a constant coefficient equation.展开更多
Energy methods and the principle of virtual work are commonly used for obtaining solutions of boundary value problems (BVPs) and initial value problems (IVPs) associated with homogeneous, isotropic and non-homogeneous...Energy methods and the principle of virtual work are commonly used for obtaining solutions of boundary value problems (BVPs) and initial value problems (IVPs) associated with homogeneous, isotropic and non-homogeneous, non-isotropic matter without using (or in the absence of) the mathematical models of the BVPs and the IVPs. These methods are also used for deriving mathematical models for BVPs and IVPs associated with isotropic, homogeneous as well as non-homogeneous, non-isotropic continuous matter. In energy methods when applied to IVPs, one constructs energy functional (<i>I</i>) consisting of kinetic energy, strain energy and the potential energy of loads. The first variation of this energy functional (<em>δI</em>) set to zero is a necessary condition for an extremum of <i>I</i>. In this approach one could use <i>δI</i> = 0 directly in constructing computational processes such as the finite element method or could derive Euler’s equations (differential or partial differential equations) from <i>δI</i> = 0, which is also satisfied by a solution obtained from <i>δI</i> = 0. The Euler’s equations obtained from <i>δI</i> = 0 indeed are the mathematical model associated with the energy functional <i>I</i>. In case of BVPs we follow the same approach except in this case, the energy functional <i>I</i> consists of strain energy and the potential energy of loads. In using the principle of virtual work for BVPs and the IVPs, we can also accomplish the same as described above using energy methods. In this paper we investigate consistency and validity of the mathematical models for isotropic, homogeneous and non-isotropic, non-homogeneous continuous matter for BVPs that are derived using energy functional consisting of strain energy and the potential energy of loads. Similar investigation is also presented for IVPs using energy functional consisting of kinetic energy, strain energy and the potential energy of loads. The computational approaches for BVPs and the IVPs designed using energy functional and principle of virtual work, their consistency and validity are also investigated. Classical continuum mechanics (CCM) principles <i>i.e.</i> conservation and balance laws of CCM with consistent constitutive theories and the elements of calculus of variations are employed in the investigations presented in this paper.展开更多
An inexact Halley's method-Halley-PCG(preconditioned conjugate gradient) method is proposed for solving the systems of linear equations for improved Halley method either by Cholesky factorization exactly or by prec...An inexact Halley's method-Halley-PCG(preconditioned conjugate gradient) method is proposed for solving the systems of linear equations for improved Halley method either by Cholesky factorization exactly or by preconditioned conjugate gradient method approximately. The convergence result is given and the efficiency of the method compared to the improved Halley's method is shown.展开更多
This work considers a special case of Lotka-Volterra equations, which means that in the system of two ordinary differential equations, we take the four parameters equal to one. The reason is that we want just to illus...This work considers a special case of Lotka-Volterra equations, which means that in the system of two ordinary differential equations, we take the four parameters equal to one. The reason is that we want just to illustrate the procedure to reduce that system to only one ordinary differential equation, such that we know its analytical solution. This idea will be applied to study the relations between a system of three ordinary differential equations, and a couple of partial differential equations. Lotka-Volterra equations are solved numerically by a fourth-order predictor-corrector method, which is initialized by an improved Euler method with a rather small time step because it is only a second-order algorithm. Then, it is proposed a model with three species, defined by ordinary differential equations.展开更多
文摘Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.
文摘The paper presents the comparative study on numerical methods of Euler method,Improved Euler method and fourth-order Runge-Kutta method for solving the engineering problems and applications.The three proposed methods are quite efficient and practically well suited for solving the unknown engineering problems.This paper aims to enhance the teaching and learning quality of teachers and students for various levels.At each point of the interval,the value of y is calculated and compared with its exact value at that point.The next interesting point is the observation of error from those methods.Error in the value of y is the difference between calculated and exact value.A mathematical equation which relates various functions with its derivatives is known as a differential equation.It is a popular field of mathematics because of its application to real-world problems.To calculate the exact values,the approximate values and the errors,the numerical tool such as MATLAB is appropriate for observing the results.This paper mainly concentrates on identifying the method which provides more accurate results.Then the analytical results and calculates their corresponding error were compared in details.The minimum error directly reflected to realize the best method from different numerical methods.According to the analyses from those three approaches,we observed that only the error is nominal for the fourth-order Runge-Kutta method.
文摘Our purpose of this paper is to apply the improved Kudryashov method for solving various types of nonlinear fractional partial differential equations. As an application, the time-space fractional Korteweg-de Vries-Burger (KdV-Burger) equation is solved using this method and we get some new travelling wave solutions. To acquire our purpose a complex transformation has been also used to reduce nonlinear fractional partial differential equations to nonlinear ordinary differential equations of integer order, in the sense of the Jumarie’s modified Riemann-Liouville derivative. Afterwards, the improved Kudryashov method is implemented and we get our required reliable solutions where the results are justified by mathematical software Maple-13.
文摘Laplace transform is one of the powerful tools for solving differential equations in engineering and other science subjects.Using the Laplace transform for solving differential equations,however,sometimes leads to solutions in the Laplace domain that are not readily invertible to the real domain by analyticalmeans.Thus,we need numerical inversionmethods to convert the obtained solution fromLaplace domain to a real domain.In this paper,we propose a numerical scheme based on Laplace transform and numerical inverse Laplace transform for the approximate solution of fractal-fractional differential equations with orderα,β.Our proposed numerical scheme is based on three main steps.First,we convert the given fractal-fractional differential equation to fractional-differential equation in Riemann-Liouville sense,and then into Caputo sense.Secondly,we transformthe fractional differential equation in Caputo sense to an equivalent equation in Laplace space.Then the solution of the transformed equation is obtained in Laplace domain.Finally,the solution is converted into the real domain using numerical inversion of Laplace transform.Three inversion methods are evaluated in this paper,and their convergence is also discussed.Three test problems are used to validate the inversion methods.We demonstrate our results with the help of tables and figures.The obtained results show that Euler’s and Talbot’s methods performed better than Stehfest’s method.
基金This research is supported by the Deputyship forResearch&Innovation,Ministry of Education in Saudi Arabia under Project Number(IFP-2022-35).
文摘To maximize energy profit with the participation of electricity,natural gas,and district heating networks in the day-ahead market,stochastic scheduling of energy hubs taking into account the uncertainty of photovoltaic and wind resources,has been carried out.This has been done using a new meta-heuristic algorithm,improved artificial rabbits optimization(IARO).In this study,the uncertainty of solar and wind energy sources is modeled using Hang’s two-point estimating method(TPEM).The IARO algorithm is applied to calculate the best capacity of hub energy equipment,such as solar and wind renewable energy sources,combined heat and power(CHP)systems,steamboilers,energy storage,and electric cars in the day-aheadmarket.The standard ARO algorithmis developed to mimic the foraging behavior of rabbits,and in this work,the algorithm’s effectiveness in avoiding premature convergence is improved by using the dystudynamic inertia weight technique.The proposed IARO-based scheduling framework’s performance is evaluated against that of traditional ARO,particle swarm optimization(PSO),and salp swarm algorithm(SSA).The findings show that,in comparison to previous approaches,the suggested meta-heuristic scheduling framework based on the IARO has increased energy profit in day-ahead electricity,gas,and heating markets by satisfying the operational and energy hub limitations.Additionally,the results show that TPEM approach dependability consideration decreased hub energy’s profit by 8.995%as compared to deterministic planning.
文摘In this work, we apply the Zhou’s method [1] or differential transformation method (DTM) for solving the Euler equidimensional equation. The Zhou’s method may be considered as alternative and efficient for finding the approximate solutions of initial values problems. We prove superiority of this method by applying them on the some Euler type equation, in this case of order 2 and 3 [2]. The power series solution of the reduced equation transforms into an approximate implicit solution of the original equations. The results agreed with the exact solution obtained via transformation to a constant coefficient equation.
文摘Energy methods and the principle of virtual work are commonly used for obtaining solutions of boundary value problems (BVPs) and initial value problems (IVPs) associated with homogeneous, isotropic and non-homogeneous, non-isotropic matter without using (or in the absence of) the mathematical models of the BVPs and the IVPs. These methods are also used for deriving mathematical models for BVPs and IVPs associated with isotropic, homogeneous as well as non-homogeneous, non-isotropic continuous matter. In energy methods when applied to IVPs, one constructs energy functional (<i>I</i>) consisting of kinetic energy, strain energy and the potential energy of loads. The first variation of this energy functional (<em>δI</em>) set to zero is a necessary condition for an extremum of <i>I</i>. In this approach one could use <i>δI</i> = 0 directly in constructing computational processes such as the finite element method or could derive Euler’s equations (differential or partial differential equations) from <i>δI</i> = 0, which is also satisfied by a solution obtained from <i>δI</i> = 0. The Euler’s equations obtained from <i>δI</i> = 0 indeed are the mathematical model associated with the energy functional <i>I</i>. In case of BVPs we follow the same approach except in this case, the energy functional <i>I</i> consists of strain energy and the potential energy of loads. In using the principle of virtual work for BVPs and the IVPs, we can also accomplish the same as described above using energy methods. In this paper we investigate consistency and validity of the mathematical models for isotropic, homogeneous and non-isotropic, non-homogeneous continuous matter for BVPs that are derived using energy functional consisting of strain energy and the potential energy of loads. Similar investigation is also presented for IVPs using energy functional consisting of kinetic energy, strain energy and the potential energy of loads. The computational approaches for BVPs and the IVPs designed using energy functional and principle of virtual work, their consistency and validity are also investigated. Classical continuum mechanics (CCM) principles <i>i.e.</i> conservation and balance laws of CCM with consistent constitutive theories and the elements of calculus of variations are employed in the investigations presented in this paper.
文摘An inexact Halley's method-Halley-PCG(preconditioned conjugate gradient) method is proposed for solving the systems of linear equations for improved Halley method either by Cholesky factorization exactly or by preconditioned conjugate gradient method approximately. The convergence result is given and the efficiency of the method compared to the improved Halley's method is shown.
文摘This work considers a special case of Lotka-Volterra equations, which means that in the system of two ordinary differential equations, we take the four parameters equal to one. The reason is that we want just to illustrate the procedure to reduce that system to only one ordinary differential equation, such that we know its analytical solution. This idea will be applied to study the relations between a system of three ordinary differential equations, and a couple of partial differential equations. Lotka-Volterra equations are solved numerically by a fourth-order predictor-corrector method, which is initialized by an improved Euler method with a rather small time step because it is only a second-order algorithm. Then, it is proposed a model with three species, defined by ordinary differential equations.