The Time Fractional Burger equation was solved in this study using the Mabel software and the Variational Iteration approach. where a number of instances of the Time Fractional Burger Equation were handled using this ...The Time Fractional Burger equation was solved in this study using the Mabel software and the Variational Iteration approach. where a number of instances of the Time Fractional Burger Equation were handled using this technique. Tables and images were used to present the collected numerical results. The difference between the exact and numerical solutions demonstrates the effectiveness of the Mabel program’s solution, as well as the accuracy and closeness of the results this method produced. It also demonstrates the Mabel program’s ability to quickly and effectively produce the numerical solution.展开更多
In this paper, the Adomian decomposition method was used to solve the Time Fractional Burger equation using Mabel program. This method was applied to a number of examples of the Time Fractional Burger Equation. The ob...In this paper, the Adomian decomposition method was used to solve the Time Fractional Burger equation using Mabel program. This method was applied to a number of examples of the Time Fractional Burger Equation. The obtained numerical results were presented in the form of tables and graphics. The difference between the exact solutions and the numerical solutions shows us the effectiveness of the solution using the Mabel program and that this method gave accurate results and was close to the exact solution, in addition to its ability to obtain the numerical solution quickly and efficiently using the Mabel program.展开更多
In this paper, our objective is to explore novel solitary wave solutions of the Burgers-Fisher equation, which characterizes the interplay between diffusion and reaction phenomena. Understanding this equation is cruci...In this paper, our objective is to explore novel solitary wave solutions of the Burgers-Fisher equation, which characterizes the interplay between diffusion and reaction phenomena. Understanding this equation is crucial for addressing challenges in fluid, chemical kinetics and population dynamics. We tackle this task by employing the Riccati equation and employing various function transformations to solve the Burgers-Fisher equation. By adopting different coefficients in the Riccati equation, we obtain a wide range of exact solutions, many of which have not been previously documented. These abundant solitary wave solutions serve as valuable tools for comprehending the Burgers-Fisher equation and contribute to expanding our knowledge in this field.展开更多
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.展开更多
In this work,we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion,the instantaneous nonlinear cubic ...In this work,we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion,the instantaneous nonlinear cubic Kerr response,and the nonlinear delayed Raman molecular vibrational response.Unlike the first-order PDE-ODE governing equations considered previously in Bokil et al.(J Comput Phys 350:420–452,2017)and Lyu et al.(J Sci Comput 89:1–42,2021),a model of mixed-order form is adopted here that consists of the first-order PDE part for Maxwell’s equations coupled with the second-order ODE part(i.e.,the auxiliary differential equations)modeling the linear and nonlinear dispersion in the material.The main contribution is a new numerical strategy to treat the Kerr and Raman nonlinearities to achieve provable energy stability property within a second-order temporal discretization.A nodal discontinuous Galerkin(DG)method is further applied in space for efficiently handling nonlinear terms at the algebraic level,while preserving the energy stability and achieving high-order accuracy.Indeed with d_(E)as the number of the components of the electric field,only a d_(E)×d_(E)nonlinear algebraic system needs to be solved at each interpolation node,and more importantly,all these small nonlinear systems are completely decoupled over one time step,rendering very high parallel efficiency.We evaluate the proposed schemes by comparing them with the methods in Bokil et al.(2017)and Lyu et al.(2021)(implemented in nodal form)regarding the accuracy,computational efficiency,and energy stability,by a parallel scalability study,and also through the simulations of the soliton-like wave propagation in one dimension,as well as the spatial-soliton propagation and two-beam interactions modeled by the two-dimensional transverse electric(TE)mode of the equations.展开更多
The solution of Poisson’s Equation plays an important role in many areas, including modeling high-intensity and high-brightness beams in particle accelerators. For the computational domain with a large aspect ratio, ...The solution of Poisson’s Equation plays an important role in many areas, including modeling high-intensity and high-brightness beams in particle accelerators. For the computational domain with a large aspect ratio, the integrated Green’s function method has been adopted to solve the 3D Poisson equation subject to open boundary conditions. In this paper, we report on the efficient implementation of this method, which can save more than a factor of 50 computing time compared with the direct brute force implementation and its improvement under certain extreme conditions.展开更多
Our study identifies a subtle deviation from Newton’s third law in the derivation of the ideal rocket equation, also known as the Tsiolkovsky Rocket Equation (TRE). TRE can be derived using a 1D elastic collision mod...Our study identifies a subtle deviation from Newton’s third law in the derivation of the ideal rocket equation, also known as the Tsiolkovsky Rocket Equation (TRE). TRE can be derived using a 1D elastic collision model of the momentum exchange between the differential propellant mass element (dm) and the rocket final mass (m1), in which dm initially travels forward to collide with m1 and rebounds to exit through the exhaust nozzle with a velocity that is known as the effective exhaust velocity ve. We observe that such a model does not explain how dm was able to acquire its initial forward velocity without the support of a reactive mass traveling in the opposite direction. We show instead that the initial kinetic energy of dm is generated from dm itself by a process of self-combustion and expansion. In our ideal rocket with a single particle dm confined inside a hollow tube with one closed end, we show that the process of self-combustion and expansion of dm will result in a pair of differential particles each with a mass dm/2, and each traveling away from one another along the tube axis, from the center of combustion. These two identical particles represent the active and reactive sub-components of dm, co-generated in compliance with Newton’s third law of equal action and reaction. Building on this model, we derive a linear momentum ODE of the system, the solution of which yields what we call the Revised Tsiolkovsky Rocket Equation (RTRE). We show that RTRE has a mathematical form that is similar to TRE, with the exception of the effective exhaust velocity (ve) term. The ve term in TRE is replaced in RTRE by the average of two distinct exhaust velocities that we refer to as fast-jet, vx<sub>1</sub>, and slow-jet, vx<sub>2</sub>. These two velocities correspond, respectively, to the velocities of the detonation pressure wave that is vectored directly towards the exhaust nozzle, and the retonation wave that is initially vectored in the direction of rocket propagation, but subsequently becomes reflected from the thrust surface of the combustion chamber to exit through the exhaust nozzle with a time lag behind the detonation wave. The detonation-retonation phenomenon is supported by experimental evidence in the published literature. Finally, we use a convolution model to simulate the composite exhaust pressure wave, highlighting the frequency spectrum of the pressure perturbations that are generated by the mutual interference between the fast-jet and slow-jet components. Our analysis offers insights into the origin of combustion oscillations in rocket engines, with possible extensions beyond rocket engineering into other fields of combustion engineering.展开更多
We applied a spatial high-order finite-difference-time-domain (HO-FDTD) scheme to solve 2D Maxwell’s equations in order to develop a fluid model employed to study the production of terahertz radiation by the filament...We applied a spatial high-order finite-difference-time-domain (HO-FDTD) scheme to solve 2D Maxwell’s equations in order to develop a fluid model employed to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma. We examined the performance of the applied scheme, in this context, we implemented the developed model to study selected phenomena in terahertz radiation production, such as the excitation energy and conversion efficiency of the produced THz radiation, in addition to the influence of the pulse chirping on properties of the produced radiation. The obtained numerical results have clarified that the applied HO-FDTD scheme is precisely accurate to solve Maxwell’s equations and sufficiently valid to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma.展开更多
传统的数值求解方法面临维数灾难和效率与精度平衡问题,而基于数据驱动的神经网络求解方法又存在训练量冗余和不可解释性问题。针对此问题,物理信息神经网络(Physical Information Neural Networks,PINNs)关注了训练数据中隐含的物理先...传统的数值求解方法面临维数灾难和效率与精度平衡问题,而基于数据驱动的神经网络求解方法又存在训练量冗余和不可解释性问题。针对此问题,物理信息神经网络(Physical Information Neural Networks,PINNs)关注了训练数据中隐含的物理先验知识,融合了神经网络拟合复杂变量的能力,赋予了传统神经网络所缺乏的物理可解释性。应用该算法模型,提出了一种基于PINN的Burgers方程求解模型,该算法模型在训练中施加物理信息约束,因此能用少量的训练样本学习预测到分布在时空域上的偏微分方程模型。实验结果表明,在1+1维Burgers方程算例下,所提方法相比于经典的机器学习算法能有效捕抓到方程的变化并进行精确模拟,相比于有限差分法,可以大幅度缩短模拟时间。通过对不同的网络参数进行比较实验,所提方法在10%的噪声破坏下能产生合理的识别准确度,网络逼近方程的待定系数误差在0.001以内。展开更多
Burgers equation in random environment is studied. In order to give the exact solutions of random Burgers equation, we only consider the Wick-type stochastic Burgers equation which is the perturbation of the Burgers e...Burgers equation in random environment is studied. In order to give the exact solutions of random Burgers equation, we only consider the Wick-type stochastic Burgers equation which is the perturbation of the Burgers equation with variable coefficients by white noise W(t)=Bt, where Bt is a Brown motion. The auto-Baecklund transformation and stochastic soliton solutions of the Wick-type stochastic Burgers equation are shown by the homogeneous balance and Hermite transform. The generalization of the Wick-type stochastic Burgers equation is also studied.展开更多
A class of generalized complex Burgers equation is considered. First, a set of equations of the complex value functions are solved by using the homotopic mapping method. The approximate solution for the original gener...A class of generalized complex Burgers equation is considered. First, a set of equations of the complex value functions are solved by using the homotopic mapping method. The approximate solution for the original generalized complex Burgers equation is obtained. This method can find the approximation of arbitrary order of precision simply and reliably.展开更多
This paper is concerned with the stability of the rarefaction wave for the Burgers equationwhere 0 ≤ a < 1/4p (q is determined by (2.2)). Roughly speaking, under the assumption that u_ < u+, the authors prove t...This paper is concerned with the stability of the rarefaction wave for the Burgers equationwhere 0 ≤ a < 1/4p (q is determined by (2.2)). Roughly speaking, under the assumption that u_ < u+, the authors prove the existence of the global smooth solution to the Cauchy problem (I), also find the solution u(x, t) to the Cauchy problem (I) satisfying sup |u(x, t) -uR(x/t)| → 0 as t → ∞, where uR(x/t) is the rarefaction wave of the non-viscous Burgersequation ut + f(u)x = 0 with Riemann initial data u(x, 0) =展开更多
It is obtained the existence of the weak solution for a degenerate generalized Burgers equation under the restriction u0 ∈ L∞. The main method is to add viscosity perturbation and obtain some estimates in L1 norm. M...It is obtained the existence of the weak solution for a degenerate generalized Burgers equation under the restriction u0 ∈ L∞. The main method is to add viscosity perturbation and obtain some estimates in L1 norm. Meanwhile it is obtained the solution is exponential decay when the initial data has compact support.展开更多
In this work we generate the numerical solutions of the Burgers’ equation by applying the Crank-Nicolson method directly to the Burgers’ equation, i.e., we do not use Hopf-Cole transformation to reduce Burgers’ equ...In this work we generate the numerical solutions of the Burgers’ equation by applying the Crank-Nicolson method directly to the Burgers’ equation, i.e., we do not use Hopf-Cole transformation to reduce Burgers’ equation into the linear heat equation. Absolute error of the present method is compared to the absolute error of the two existing methods for two test problems. The method is also analyzed for a third test problem, nu-merical solutions as well as exact solutions for different values of viscosity are calculated and we find that the numerical solutions are very close to exact solution.展开更多
In the current work, we extend the local discontinuous Galerkin method to a more general application system. The Burgers and coupled Burgers equations are solved by the local discontinuous Galerkin method. Numerical e...In the current work, we extend the local discontinuous Galerkin method to a more general application system. The Burgers and coupled Burgers equations are solved by the local discontinuous Galerkin method. Numerical experiments are given to verify the efficiency and accuracy of our method. Moreover the numerical results show that the method can approximate sharp fronts accurately with minimal oscillation.展开更多
This paper studies the coupled Burgers equation and the high-order Boussinesq-Burgers equation. The Hirota bilinear method is applied to show that the two equations are completely integrable. Multiple-kink (soliton)...This paper studies the coupled Burgers equation and the high-order Boussinesq-Burgers equation. The Hirota bilinear method is applied to show that the two equations are completely integrable. Multiple-kink (soliton) solutions and multiple-singular-kink (soliton) solutions are derived for the two equations.展开更多
文摘The Time Fractional Burger equation was solved in this study using the Mabel software and the Variational Iteration approach. where a number of instances of the Time Fractional Burger Equation were handled using this technique. Tables and images were used to present the collected numerical results. The difference between the exact and numerical solutions demonstrates the effectiveness of the Mabel program’s solution, as well as the accuracy and closeness of the results this method produced. It also demonstrates the Mabel program’s ability to quickly and effectively produce the numerical solution.
文摘In this paper, the Adomian decomposition method was used to solve the Time Fractional Burger equation using Mabel program. This method was applied to a number of examples of the Time Fractional Burger Equation. The obtained numerical results were presented in the form of tables and graphics. The difference between the exact solutions and the numerical solutions shows us the effectiveness of the solution using the Mabel program and that this method gave accurate results and was close to the exact solution, in addition to its ability to obtain the numerical solution quickly and efficiently using the Mabel program.
文摘In this paper, our objective is to explore novel solitary wave solutions of the Burgers-Fisher equation, which characterizes the interplay between diffusion and reaction phenomena. Understanding this equation is crucial for addressing challenges in fluid, chemical kinetics and population dynamics. We tackle this task by employing the Riccati equation and employing various function transformations to solve the Burgers-Fisher equation. By adopting different coefficients in the Riccati equation, we obtain a wide range of exact solutions, many of which have not been previously documented. These abundant solitary wave solutions serve as valuable tools for comprehending the Burgers-Fisher equation and contribute to expanding our knowledge in this field.
文摘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.
基金supported by China Postdoctoral Science Foundation grant 2020TQ0344the NSFC grants 11871139 and 12101597the NSF grants DMS-1720116,DMS-2012882,DMS-2011838,DMS-1719942,DMS-1913072.
文摘In this work,we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion,the instantaneous nonlinear cubic Kerr response,and the nonlinear delayed Raman molecular vibrational response.Unlike the first-order PDE-ODE governing equations considered previously in Bokil et al.(J Comput Phys 350:420–452,2017)and Lyu et al.(J Sci Comput 89:1–42,2021),a model of mixed-order form is adopted here that consists of the first-order PDE part for Maxwell’s equations coupled with the second-order ODE part(i.e.,the auxiliary differential equations)modeling the linear and nonlinear dispersion in the material.The main contribution is a new numerical strategy to treat the Kerr and Raman nonlinearities to achieve provable energy stability property within a second-order temporal discretization.A nodal discontinuous Galerkin(DG)method is further applied in space for efficiently handling nonlinear terms at the algebraic level,while preserving the energy stability and achieving high-order accuracy.Indeed with d_(E)as the number of the components of the electric field,only a d_(E)×d_(E)nonlinear algebraic system needs to be solved at each interpolation node,and more importantly,all these small nonlinear systems are completely decoupled over one time step,rendering very high parallel efficiency.We evaluate the proposed schemes by comparing them with the methods in Bokil et al.(2017)and Lyu et al.(2021)(implemented in nodal form)regarding the accuracy,computational efficiency,and energy stability,by a parallel scalability study,and also through the simulations of the soliton-like wave propagation in one dimension,as well as the spatial-soliton propagation and two-beam interactions modeled by the two-dimensional transverse electric(TE)mode of the equations.
文摘The solution of Poisson’s Equation plays an important role in many areas, including modeling high-intensity and high-brightness beams in particle accelerators. For the computational domain with a large aspect ratio, the integrated Green’s function method has been adopted to solve the 3D Poisson equation subject to open boundary conditions. In this paper, we report on the efficient implementation of this method, which can save more than a factor of 50 computing time compared with the direct brute force implementation and its improvement under certain extreme conditions.
文摘Our study identifies a subtle deviation from Newton’s third law in the derivation of the ideal rocket equation, also known as the Tsiolkovsky Rocket Equation (TRE). TRE can be derived using a 1D elastic collision model of the momentum exchange between the differential propellant mass element (dm) and the rocket final mass (m1), in which dm initially travels forward to collide with m1 and rebounds to exit through the exhaust nozzle with a velocity that is known as the effective exhaust velocity ve. We observe that such a model does not explain how dm was able to acquire its initial forward velocity without the support of a reactive mass traveling in the opposite direction. We show instead that the initial kinetic energy of dm is generated from dm itself by a process of self-combustion and expansion. In our ideal rocket with a single particle dm confined inside a hollow tube with one closed end, we show that the process of self-combustion and expansion of dm will result in a pair of differential particles each with a mass dm/2, and each traveling away from one another along the tube axis, from the center of combustion. These two identical particles represent the active and reactive sub-components of dm, co-generated in compliance with Newton’s third law of equal action and reaction. Building on this model, we derive a linear momentum ODE of the system, the solution of which yields what we call the Revised Tsiolkovsky Rocket Equation (RTRE). We show that RTRE has a mathematical form that is similar to TRE, with the exception of the effective exhaust velocity (ve) term. The ve term in TRE is replaced in RTRE by the average of two distinct exhaust velocities that we refer to as fast-jet, vx<sub>1</sub>, and slow-jet, vx<sub>2</sub>. These two velocities correspond, respectively, to the velocities of the detonation pressure wave that is vectored directly towards the exhaust nozzle, and the retonation wave that is initially vectored in the direction of rocket propagation, but subsequently becomes reflected from the thrust surface of the combustion chamber to exit through the exhaust nozzle with a time lag behind the detonation wave. The detonation-retonation phenomenon is supported by experimental evidence in the published literature. Finally, we use a convolution model to simulate the composite exhaust pressure wave, highlighting the frequency spectrum of the pressure perturbations that are generated by the mutual interference between the fast-jet and slow-jet components. Our analysis offers insights into the origin of combustion oscillations in rocket engines, with possible extensions beyond rocket engineering into other fields of combustion engineering.
文摘We applied a spatial high-order finite-difference-time-domain (HO-FDTD) scheme to solve 2D Maxwell’s equations in order to develop a fluid model employed to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma. We examined the performance of the applied scheme, in this context, we implemented the developed model to study selected phenomena in terahertz radiation production, such as the excitation energy and conversion efficiency of the produced THz radiation, in addition to the influence of the pulse chirping on properties of the produced radiation. The obtained numerical results have clarified that the applied HO-FDTD scheme is precisely accurate to solve Maxwell’s equations and sufficiently valid to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma.
文摘传统的数值求解方法面临维数灾难和效率与精度平衡问题,而基于数据驱动的神经网络求解方法又存在训练量冗余和不可解释性问题。针对此问题,物理信息神经网络(Physical Information Neural Networks,PINNs)关注了训练数据中隐含的物理先验知识,融合了神经网络拟合复杂变量的能力,赋予了传统神经网络所缺乏的物理可解释性。应用该算法模型,提出了一种基于PINN的Burgers方程求解模型,该算法模型在训练中施加物理信息约束,因此能用少量的训练样本学习预测到分布在时空域上的偏微分方程模型。实验结果表明,在1+1维Burgers方程算例下,所提方法相比于经典的机器学习算法能有效捕抓到方程的变化并进行精确模拟,相比于有限差分法,可以大幅度缩短模拟时间。通过对不同的网络参数进行比较实验,所提方法在10%的噪声破坏下能产生合理的识别准确度,网络逼近方程的待定系数误差在0.001以内。
文摘Burgers equation in random environment is studied. In order to give the exact solutions of random Burgers equation, we only consider the Wick-type stochastic Burgers equation which is the perturbation of the Burgers equation with variable coefficients by white noise W(t)=Bt, where Bt is a Brown motion. The auto-Baecklund transformation and stochastic soliton solutions of the Wick-type stochastic Burgers equation are shown by the homogeneous balance and Hermite transform. The generalization of the Wick-type stochastic Burgers equation is also studied.
基金Project supported by the National Natural Science Foundation of China (Grant No. 40876010), the Main Direction Program of the Knowledge Innovation Project of the Chinese Academy of Sciences (Grant No. KZCX2-YW-Q03-08), the R & D Special Fund for Public Welfare Industry (Meteorology) (Grant No. GYHY200806010), the LASG State Key Laboratory Special Fund, the Foundation of E-Institutes of Shanghai Municipal Education Commission (Crant No. E03004) and the Natural Science Foundation of Zhejiang Province of China (Grant No. Y6090164).
文摘A class of generalized complex Burgers equation is considered. First, a set of equations of the complex value functions are solved by using the homotopic mapping method. The approximate solution for the original generalized complex Burgers equation is obtained. This method can find the approximation of arbitrary order of precision simply and reliably.
文摘This paper is concerned with the stability of the rarefaction wave for the Burgers equationwhere 0 ≤ a < 1/4p (q is determined by (2.2)). Roughly speaking, under the assumption that u_ < u+, the authors prove the existence of the global smooth solution to the Cauchy problem (I), also find the solution u(x, t) to the Cauchy problem (I) satisfying sup |u(x, t) -uR(x/t)| → 0 as t → ∞, where uR(x/t) is the rarefaction wave of the non-viscous Burgersequation ut + f(u)x = 0 with Riemann initial data u(x, 0) =
文摘It is obtained the existence of the weak solution for a degenerate generalized Burgers equation under the restriction u0 ∈ L∞. The main method is to add viscosity perturbation and obtain some estimates in L1 norm. Meanwhile it is obtained the solution is exponential decay when the initial data has compact support.
文摘In this work we generate the numerical solutions of the Burgers’ equation by applying the Crank-Nicolson method directly to the Burgers’ equation, i.e., we do not use Hopf-Cole transformation to reduce Burgers’ equation into the linear heat equation. Absolute error of the present method is compared to the absolute error of the two existing methods for two test problems. The method is also analyzed for a third test problem, nu-merical solutions as well as exact solutions for different values of viscosity are calculated and we find that the numerical solutions are very close to exact solution.
基金supported by the National Natural Science Foundation of China(Grant No.11171038)
文摘In the current work, we extend the local discontinuous Galerkin method to a more general application system. The Burgers and coupled Burgers equations are solved by the local discontinuous Galerkin method. Numerical experiments are given to verify the efficiency and accuracy of our method. Moreover the numerical results show that the method can approximate sharp fronts accurately with minimal oscillation.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10871117 and 10571110)
文摘This paper studies the coupled Burgers equation and the high-order Boussinesq-Burgers equation. The Hirota bilinear method is applied to show that the two equations are completely integrable. Multiple-kink (soliton) solutions and multiple-singular-kink (soliton) solutions are derived for the two equations.