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.展开更多
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.展开更多
For any given positive integer n ≥ 1, the Euler function φ（n） is defined to be the number of positive integers not exceeding n which are relatively prime to n. w（n） is defined to be the number of different prime...For any given positive integer n ≥ 1, the Euler function φ（n） is defined to be the number of positive integers not exceeding n which are relatively prime to n. w（n） is defined to be the number of different prime divisors of n. Some kind of equations involving Euler＇s function is studied in the paper.展开更多
The depth from extreme points(DEXP)method can be used for estimating source depths and providing a rough image as a starting model for inversion.However,the application of the DEXP method is limited by the lack of pri...The depth from extreme points(DEXP)method can be used for estimating source depths and providing a rough image as a starting model for inversion.However,the application of the DEXP method is limited by the lack of prior information regarding the structural index.Herein,we describe an automatic DEXP method derived from Euler’s Homogeneity equation,and we call it the Euler–DEXP method.We prove that its scaling field is independent of structural indices,and the scaling exponent is a constant for any potential field or its derivative.Therefore,we can simultaneously estimate source depths with diff erent geometries in one DEXP image.The implementation of the Euler–DEXP method is fully automatic.The structural index can be subsequently determined by utilizing the estimated depth.This method has been tested using synthetic cases with single and multiple sources.All estimated solutions are in accordance with theoretical source parameters.We demonstrate the practicability of the Euler–DEXP method with the gravity field data of the Hastings Salt Dome.The results ultimately represent a better understanding of the geometry and depth of the salt dome.展开更多
Solutions to the differential equation in Smith’s Prize Examination taken by Maxwell are discussed. It was a competitive examination using which skill full students were identified and James Clerk Maxwell was one of ...Solutions to the differential equation in Smith’s Prize Examination taken by Maxwell are discussed. It was a competitive examination using which skill full students were identified and James Clerk Maxwell was one of them. He later formulated the theory of Electromagnetism and predicted the light speed & its value was subsequently confirmed by experiments. Light travel in a direction perpendicular to oscillating electric and magnetic field through a vacuum from sun. In the same exam paper, Maxwell answered the question related to Stokes Theorem of vector calculus which was used in the formalism of Electromagnetic theory.展开更多
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.展开更多
We consider optimal control problems for the flow of gas in a pipe network. The equations of motions are taken to be represented by a semi-linear model derived from the fully nonlinear isothermal Euler gas equations. ...We consider optimal control problems for the flow of gas in a pipe network. The equations of motions are taken to be represented by a semi-linear model derived from the fully nonlinear isothermal Euler gas equations. We formulate an optimal control problem on a given network and introduce a time discretization thereof. We then study the well-posedness of the corresponding time-discrete optimal control problem. In order to further reduce the complexity, we consider an instantaneous control strategy. The main part of the paper is concerned with a non-overlapping domain decomposition of the semi-linear elliptic optimal control problem on the graph into local problems on a small part of the network, ultimately on a single edge.展开更多
The viscous dissipation limit of weak solutions is considered for the Navier-Stokes equations of compressible isentropic flows confined in a bounded domain.We establish a Kato-type criterion for the validity of the in...The viscous dissipation limit of weak solutions is considered for the Navier-Stokes equations of compressible isentropic flows confined in a bounded domain.We establish a Kato-type criterion for the validity of the inviscid limit for the weak solutions of the Navier-Stokes equations in a function space with the regularity index close to Onsager’s critical threshold.In particular,we prove that under such a regularity assumption,if the viscous energy dissipation rate vanishes in a boundary layer of thickness in the order of the viscosity,then the weak solutions of the Navier-Stokes equations converge to a weak admissible solution of the Euler equations.Our approach is based on the commutator estimates and a subtle foliation technique near the boundary of the domain.展开更多
As is well known,the definitions of fractional sum and fractional difference of f(z)on non-uniform lattices x(z)=c1z^(2)+c2z+c3 or x(z)=c1q^(z)+c2q^(-z)+c3 are more difficult and complicated.In this article,for the fi...As is well known,the definitions of fractional sum and fractional difference of f(z)on non-uniform lattices x(z)=c1z^(2)+c2z+c3 or x(z)=c1q^(z)+c2q^(-z)+c3 are more difficult and complicated.In this article,for the first time we propose the definitions of the fractional sum and fractional difference on non-uniform lattices by two different ways.The analogue of Euler’s Beta formula,Cauchy’Beta formula on non-uniform lattices are established,and some fundamental theorems of fractional calculas,the solution of the generalized Abel equation on non-uniform lattices are obtained etc.展开更多
This paper introduces computational fluid used in aerospace engineering, to deal with surface physical and mathematical foundations of CFD, this traffic problems such as queue/platoon distribution, dynamics （CFD）, a...This paper introduces computational fluid used in aerospace engineering, to deal with surface physical and mathematical foundations of CFD, this traffic problems such as queue/platoon distribution, dynamics （CFD）, a numerical traffic flow related problems. approach widely and successfully After a brief introduction of the paper develops CFD implementation methodology for modeling shockwave propagation, and prediction of system performance. Some theoretical and practical applications are discussed in this paper to illustrate the implementation methodology. It is found that CFD approach can facilitate a superior insight into the formation and propagation of congestion, thereby supporting more effective methods to alleviate congestion. In addition, CFD approach is found capable of assessing freeway system performance using less ITS detectors, and enhancing the coverage and reliability of a traffic detection system.展开更多
We present an adaptive lattice Boltzmann model to simulate super- sonic flows.The particle velocities are determined by the mean velocity and internal energy.The adaptive nature of particle velocities permits the mean...We present an adaptive lattice Boltzmann model to simulate super- sonic flows.The particle velocities are determined by the mean velocity and internal energy.The adaptive nature of particle velocities permits the mean flow to have high Mach number.A particle potential energy is introduced so that the model is suitable for the perfect gas with arbitrary specific heat ratio.The Navier-Stokes equations are derived by the Chapman-Enskog method from the BGK Boltzmann equation. As preliminary tests,two kinds of simulations have been performed on hexagonal lattices.One is the one-dimensional simulation for sinusoidal velocity distributions. The velocity distributions are compared with the analytical solution and the mea- sured viscosity is compared with the theoretical values.The agreements are basically good.However,the discretion error may cause some non-isotropic effects.The other simulation is the 29 degree shock reflection.展开更多
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.展开更多
文摘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.
文摘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.
基金Foundation item： Supported by the National Natural Science Foundation of China（10671056）
文摘For any given positive integer n ≥ 1, the Euler function φ（n） is defined to be the number of positive integers not exceeding n which are relatively prime to n. w（n） is defined to be the number of different prime divisors of n. Some kind of equations involving Euler＇s function is studied in the paper.
基金supported by the National Natural Science Foundation of China (Grant No.42176186).
文摘The depth from extreme points(DEXP)method can be used for estimating source depths and providing a rough image as a starting model for inversion.However,the application of the DEXP method is limited by the lack of prior information regarding the structural index.Herein,we describe an automatic DEXP method derived from Euler’s Homogeneity equation,and we call it the Euler–DEXP method.We prove that its scaling field is independent of structural indices,and the scaling exponent is a constant for any potential field or its derivative.Therefore,we can simultaneously estimate source depths with diff erent geometries in one DEXP image.The implementation of the Euler–DEXP method is fully automatic.The structural index can be subsequently determined by utilizing the estimated depth.This method has been tested using synthetic cases with single and multiple sources.All estimated solutions are in accordance with theoretical source parameters.We demonstrate the practicability of the Euler–DEXP method with the gravity field data of the Hastings Salt Dome.The results ultimately represent a better understanding of the geometry and depth of the salt dome.
文摘Solutions to the differential equation in Smith’s Prize Examination taken by Maxwell are discussed. It was a competitive examination using which skill full students were identified and James Clerk Maxwell was one of them. He later formulated the theory of Electromagnetism and predicted the light speed & its value was subsequently confirmed by experiments. Light travel in a direction perpendicular to oscillating electric and magnetic field through a vacuum from sun. In the same exam paper, Maxwell answered the question related to Stokes Theorem of vector calculus which was used in the formalism of Electromagnetic theory.
文摘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.
文摘We consider optimal control problems for the flow of gas in a pipe network. The equations of motions are taken to be represented by a semi-linear model derived from the fully nonlinear isothermal Euler gas equations. We formulate an optimal control problem on a given network and introduce a time discretization thereof. We then study the well-posedness of the corresponding time-discrete optimal control problem. In order to further reduce the complexity, we consider an instantaneous control strategy. The main part of the paper is concerned with a non-overlapping domain decomposition of the semi-linear elliptic optimal control problem on the graph into local problems on a small part of the network, ultimately on a single edge.
基金supported by National Science Foundation of USA(Grant No.DMS-1907584)supported by the Fundamental Research Funds for the Central Universities(Grant No.JBK 2202045)+1 种基金supported by National Science Foundation of USA(Grant Nos.DMS-1907519 and DMS-2219384)supported by National Natural Science Foundation of China(Grant No.12271122)。
文摘The viscous dissipation limit of weak solutions is considered for the Navier-Stokes equations of compressible isentropic flows confined in a bounded domain.We establish a Kato-type criterion for the validity of the inviscid limit for the weak solutions of the Navier-Stokes equations in a function space with the regularity index close to Onsager’s critical threshold.In particular,we prove that under such a regularity assumption,if the viscous energy dissipation rate vanishes in a boundary layer of thickness in the order of the viscosity,then the weak solutions of the Navier-Stokes equations converge to a weak admissible solution of the Euler equations.Our approach is based on the commutator estimates and a subtle foliation technique near the boundary of the domain.
基金Supported by the National Natural Science Foundation Fujian province of China(2016J01032).
文摘As is well known,the definitions of fractional sum and fractional difference of f(z)on non-uniform lattices x(z)=c1z^(2)+c2z+c3 or x(z)=c1q^(z)+c2q^(-z)+c3 are more difficult and complicated.In this article,for the first time we propose the definitions of the fractional sum and fractional difference on non-uniform lattices by two different ways.The analogue of Euler’s Beta formula,Cauchy’Beta formula on non-uniform lattices are established,and some fundamental theorems of fractional calculas,the solution of the generalized Abel equation on non-uniform lattices are obtained etc.
文摘This paper introduces computational fluid used in aerospace engineering, to deal with surface physical and mathematical foundations of CFD, this traffic problems such as queue/platoon distribution, dynamics （CFD）, a numerical traffic flow related problems. approach widely and successfully After a brief introduction of the paper develops CFD implementation methodology for modeling shockwave propagation, and prediction of system performance. Some theoretical and practical applications are discussed in this paper to illustrate the implementation methodology. It is found that CFD approach can facilitate a superior insight into the formation and propagation of congestion, thereby supporting more effective methods to alleviate congestion. In addition, CFD approach is found capable of assessing freeway system performance using less ITS detectors, and enhancing the coverage and reliability of a traffic detection system.
基金The project supported by the National Natural Science Foundation of China(Grant Nos.19672030 and 19972037)and by the Scientific Research Foundation for Returned Overseas Chinese ScholarsState Education Ministry
文摘We present an adaptive lattice Boltzmann model to simulate super- sonic flows.The particle velocities are determined by the mean velocity and internal energy.The adaptive nature of particle velocities permits the mean flow to have high Mach number.A particle potential energy is introduced so that the model is suitable for the perfect gas with arbitrary specific heat ratio.The Navier-Stokes equations are derived by the Chapman-Enskog method from the BGK Boltzmann equation. As preliminary tests,two kinds of simulations have been performed on hexagonal lattices.One is the one-dimensional simulation for sinusoidal velocity distributions. The velocity distributions are compared with the analytical solution and the mea- sured viscosity is compared with the theoretical values.The agreements are basically good.However,the discretion error may cause some non-isotropic effects.The other simulation is the 29 degree shock reflection.
文摘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.
