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.展开更多
Several efficient analytical methods have been developed to solve the solid-state diffusion problem, for constant diffusion coefficient problems. However, these methods cannot be applied for concentration-dependent di...Several efficient analytical methods have been developed to solve the solid-state diffusion problem, for constant diffusion coefficient problems. However, these methods cannot be applied for concentration-dependent diffusion coefficient problems and numerical methods are used instead. Herein, grid-based numerical methods derived from the control volume discretization are presented to resolve the characteristic nonlinear system of partial differential equations. A novel hybrid backward Euler control volume (HBECV) method is presented which requires only one iteration to reach an implicit solution. The HBECV results are shown to be stable and accurate for a moderate number of grid points. The computational speed and accuracy of the HBECV, justify its use in battery simulations, in which the solid-state diffusion coefficient is a strong function of the concentration.展开更多
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.展开更多
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.展开更多
The article explores the issue of designing a new design of a loading cylinder with a casing filled with vulcanized rubber for pneumomechanical spinning machines. The theoretical calculation of the deformed state of a...The article explores the issue of designing a new design of a loading cylinder with a casing filled with vulcanized rubber for pneumomechanical spinning machines. The theoretical calculation of the deformed state of a cylindrical shell filled with vulcanized rubber is given. Deflections and stresses in the rubber layer are determined, which we use approximately for the Ritz methods. The theory of the radial and axial moving rubber layer was analyzed. The specific energy of deformation of a cylindrical layer of a compound cylinder is determined. The statics of the case and the loading cylinder of spinning machines are thoroughly studied.展开更多
文摘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.
文摘Several efficient analytical methods have been developed to solve the solid-state diffusion problem, for constant diffusion coefficient problems. However, these methods cannot be applied for concentration-dependent diffusion coefficient problems and numerical methods are used instead. Herein, grid-based numerical methods derived from the control volume discretization are presented to resolve the characteristic nonlinear system of partial differential equations. A novel hybrid backward Euler control volume (HBECV) method is presented which requires only one iteration to reach an implicit solution. The HBECV results are shown to be stable and accurate for a moderate number of grid points. The computational speed and accuracy of the HBECV, justify its use in battery simulations, in which the solid-state diffusion coefficient is a strong function of the concentration.
文摘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.
文摘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.
文摘The article explores the issue of designing a new design of a loading cylinder with a casing filled with vulcanized rubber for pneumomechanical spinning machines. The theoretical calculation of the deformed state of a cylindrical shell filled with vulcanized rubber is given. Deflections and stresses in the rubber layer are determined, which we use approximately for the Ritz methods. The theory of the radial and axial moving rubber layer was analyzed. The specific energy of deformation of a cylindrical layer of a compound cylinder is determined. The statics of the case and the loading cylinder of spinning machines are thoroughly studied.
