Missile interception problem can be regarded as a two-person zero-sum differential games problem,which depends on the solution of Hamilton-Jacobi-Isaacs(HJI)equa-tion.It has been proved impossible to obtain a closed-f...Missile interception problem can be regarded as a two-person zero-sum differential games problem,which depends on the solution of Hamilton-Jacobi-Isaacs(HJI)equa-tion.It has been proved impossible to obtain a closed-form solu-tion due to the nonlinearity of HJI equation,and many iterative algorithms are proposed to solve the HJI equation.Simultane-ous policy updating algorithm(SPUA)is an effective algorithm for solving HJI equation,but it is an on-policy integral reinforce-ment learning(IRL).For online implementation of SPUA,the dis-turbance signals need to be adjustable,which is unrealistic.In this paper,an off-policy IRL algorithm based on SPUA is pro-posed without making use of any knowledge of the systems dynamics.Then,a neural-network based online adaptive critic implementation scheme of the off-policy IRL algorithm is pre-sented.Based on the online off-policy IRL method,a computa-tional intelligence interception guidance(CIIG)law is developed for intercepting high-maneuvering target.As a model-free method,intercepting targets can be achieved through measur-ing system data online.The effectiveness of the CIIG is verified through two missile and target engagement scenarios.展开更多
This paper centers on the integrated learning of English and law in China.Firstly,it outlines the importance of English in the solution of the ever increasing legal disputes between China and the outside world,which i...This paper centers on the integrated learning of English and law in China.Firstly,it outlines the importance of English in the solution of the ever increasing legal disputes between China and the outside world,which inevitably involves an integrated learning of English and law.Secondly,it points out that the content of legal English reflects a combination of legal knowledge and English skills.Thirdly,it expounds on the difficulties that Chinese English majors are facing in the process of learning English and law simultaneously and furnishes some practical suggestions.展开更多
In this paper, we present a new method to construct the conservation laws for relativistic mechanical systems by finding corresponding integrating factors. First, the Lagrange equations of relativistic mechanical syst...In this paper, we present a new method to construct the conservation laws for relativistic mechanical systems by finding corresponding integrating factors. First, the Lagrange equations of relativistic mechanical systems are established, and the definition of integrating factors of the systems is given; second, the necessary conditions for the existence of conserved quantities of the relativistic mechanical systems are studied in detail, and the relation between the conservation laws and the integrating factors of the systems is obtained and the generalized Killing equations for the determination of the integrating factors are given; finally, the conservation theorem and its inverse for the systems are established, and an example is given to illustrate the application of the results.展开更多
In this paper, the conservation laws of generalized Birkhoff system in event space are studied by using the method of integrating factors. Firstly, the generalized Pfaff-Birkhoff principle and the generalized Birkhoff...In this paper, the conservation laws of generalized Birkhoff system in event space are studied by using the method of integrating factors. Firstly, the generalized Pfaff-Birkhoff principle and the generalized Birkhoff equations are established, and the definition of the integrating factors for the system is given. Secondly, based on the concept of integrating factors, the conservation theorems and their inverse for the generalized Birkhoff system in the event space are presented in detail, and the relation between the conservation laws and the integrating factors of the system is obtained and the generalized Killing equations for the determination of the integrating factors are given. Finally, an example is given to illustrate the application of the results.展开更多
A numerical method is proposed to simulate the transverse vibrations of a viscoelastic moving string constituted by an integral law. In the numerical computation, the Galerkin method based on the Hermite functions is ...A numerical method is proposed to simulate the transverse vibrations of a viscoelastic moving string constituted by an integral law. In the numerical computation, the Galerkin method based on the Hermite functions is applied to discretize the state variables, and the Runge- Kutta method is applied to solve the resulting differential-integral equation system. A linear iterative process is designed to compute the integral terms at each time step, which makes the numerical method more efficient and accurate. As examples, nonlinear parametric vibrations of an axially moving viscoelastic string are analyzed.展开更多
This paper is concerned with the generalized nonlinear second-order equation.By the direct construction method,all of the first-order multipliers of the equation are obtained,and the corresponding complete conservatio...This paper is concerned with the generalized nonlinear second-order equation.By the direct construction method,all of the first-order multipliers of the equation are obtained,and the corresponding complete conservation laws(CLs) of such equations are provided.Furthermore,the integrability of the equation is considered in terms of the conservation laws.In addition,the relationship of multipliers and symmetries of the equations is investigated.展开更多
We formulate a “Moore’s law” for photonic integrated circuits (PICs) and their spatial integration density using two methods. One is decomposing the integrated photonics devices of diverse types into equivalent bas...We formulate a “Moore’s law” for photonic integrated circuits (PICs) and their spatial integration density using two methods. One is decomposing the integrated photonics devices of diverse types into equivalent basic elements, which makes a comparison with the generic elements of electronic integrated circuits more meaningful. The other is making a complex compo- nent equivalent to a series of basic elements of the same functionality, which is used to calculate the integration density for func- tional components realized with different structures. The results serve as a benchmark of the evolution of PICs and we can con- clude that the density of integration measured in this way roughly increases by a factor of 2 per year. The prospects for a continued increase of spatial integration density are discussed.展开更多
We construct a nonlinear integrable coupling of discrete soliton hierarchy, and establish the infinite conservation laws (CLs) for the nonlinear integrable coupling of the lattice hierarchy. As an explicit applicati...We construct a nonlinear integrable coupling of discrete soliton hierarchy, and establish the infinite conservation laws (CLs) for the nonlinear integrable coupling of the lattice hierarchy. As an explicit application of the method proposed in the paper, the infinite conservation laws of the nonlinear integrable coupling of the Volterra lattice hierarchy are presented.展开更多
Under investigation in this paper are two coupled integrable dispersionless (CID) equations modelingthe dynamics of the current-fed string within an external magnetic field.Through a set of the dependent variabletrans...Under investigation in this paper are two coupled integrable dispersionless (CID) equations modelingthe dynamics of the current-fed string within an external magnetic field.Through a set of the dependent variabletransformations, the bilinear forms for the CID equations are derived.Based on the Hirota method and symboliccomputation, the analytic N-soliton solutions are presented.Infinitely many conservation laws for the CID equationsare given through the known spectral problem.Propagation characteristics and interaction behaviors of the solitons areanalyzed graphically.展开更多
Two hierarchies of nonlinear integrable positive and negative lattice equations are derived from a discrete spectrak problem. The two lattice hierarchies are proved to have discrete zero curvature representations asso...Two hierarchies of nonlinear integrable positive and negative lattice equations are derived from a discrete spectrak problem. The two lattice hierarchies are proved to have discrete zero curvature representations associated with a discrete spectral problem, which also shows that the positive and negative hierarchies correspond to positive and negative power expansions of Lax operators with respect to the spectral parameter, respectively. Moreover, the integrable lattice models in the positive hierarchy are of polynomial type, and the integrable lattice models in the negative hierarchy are of rational type. Further, we construct infinite conservation laws about the positive hierarchy.展开更多
With the aid of Lenard recursion equations, an integrable hierarchy of nonlinear evolution equations associated with a 2 × 2 matrix spectral problem is proposed, in which the first nontrivial member in the positi...With the aid of Lenard recursion equations, an integrable hierarchy of nonlinear evolution equations associated with a 2 × 2 matrix spectral problem is proposed, in which the first nontrivial member in the positive flows can be reduced to a new generalization of the Wadati–Konno–Ichikawa(WKI) equation. Further, a new generalization of the Fokas–Lenells(FL) equation is derived from the negative flows. Resorting to these two Lax pairs and Riccati-type equations, the infinite conservation laws of these two corresponding equations are obtained.展开更多
The conservation law of J-integral in two-media with a crack paralleling to the interface of the two media was firstly proved by analytical and numerical finite element method. Then a schedule model was established th...The conservation law of J-integral in two-media with a crack paralleling to the interface of the two media was firstly proved by analytical and numerical finite element method. Then a schedule model was established that an interface crack is inserted in four media. According to the J-integral conservation law on multi-media, the energy release ratio of Ⅰ-type crack was considered to be conservation when the middle medium layers are very thin. And the conservation law was also convinced by numerical method. By means of the dimension analysis on the model, the asymptotic results and formula calculating the energy release ratio and complex stress intensity factor are presented.展开更多
As the international community place high value on the ocean,protecting the marine environment and promoting the sustainable development of the ocean have become the consensus of all countries around the world.Integra...As the international community place high value on the ocean,protecting the marine environment and promoting the sustainable development of the ocean have become the consensus of all countries around the world.Integrated ocean management,as an important part of the ocean governance mechanism within the framework of the United Nations,have approached to the scene view of human gradually.Taking the United Nations Convention on the Law of the Sea as the framework,the current international integrated ocean management is summarized and the future development trend of ocean management is analyzed,providing reference for China to build a maritime power and a maritime community with shared future.展开更多
In modem society, the rule or law is an important pattern of national governance. Realizing human rights is closely related to the rule of law. What the author intends to discuss is how to use integral thinking to und...In modem society, the rule or law is an important pattern of national governance. Realizing human rights is closely related to the rule of law. What the author intends to discuss is how to use integral thinking to understand human rights as well as the relationship between human rights and development and how to comprehend the rule of law and the promotion of human rights through the rule of law in the process of facilitating the cause of human rights.展开更多
A discrete matrix spectral problem and the associated hierarchy of Lax integrable lattice equations are presented, and it is shown that the resulting Lax integrable lattice equations are all Liouville integrable discr...A discrete matrix spectral problem and the associated hierarchy of Lax integrable lattice equations are presented, and it is shown that the resulting Lax integrable lattice equations are all Liouville integrable discrete Hamiltonian systems. A new integrable symplectic map is given by binary Bargmann constraint of the resulting hierarchy. Finally, an infinite set of conservation laws is given for the resulting hierarchy.展开更多
A novel hierarchy of integrable nonlinear evolution equations related to the combined Ablowitz–Kaup–Newell–Segur(AKNS) and Wadati–Konno–Ichikawa(WKI) spectral problems is proposed,from which the Lax pair for ...A novel hierarchy of integrable nonlinear evolution equations related to the combined Ablowitz–Kaup–Newell–Segur(AKNS) and Wadati–Konno–Ichikawa(WKI) spectral problems is proposed,from which the Lax pair for a corresponding negative flow and its infinite many conservation laws are obtained.Furthermore,a reduction of this hierarchy is discussed,by which a generalized sinh-Gordon equation is derived on the basis of its negative flow.展开更多
The classical power law relaxation, i.e. relaxation of current with inverse of power of time for a step-voltage excitation to dielectric—as popularly known as Curie-von Schweidler law is empirically derived and is ob...The classical power law relaxation, i.e. relaxation of current with inverse of power of time for a step-voltage excitation to dielectric—as popularly known as Curie-von Schweidler law is empirically derived and is observed in several relaxation experiments on various dielectrics studies since late 19th Century. This relaxation law is also regarded as “universal-law” for dielectric relaxations;and is also termed as power law. This empirical Curie-von Schewidler relaxation law is then used to derive fractional differential equations describing constituent expression for capacitor. In this paper, we give simple mathematical treatment to derive the distribution of relaxation rates of this Curie-von Schweidler law, and show that the relaxation rate follows Zipf’s power law distribution. We also show the method developed here give Zipfian power law distribution for relaxing time constants. Then we will show however mathematically correct this may be, but physical interpretation from the obtained time constants distribution are contradictory to the Zipfian rate relaxation distribution. In this paper, we develop possible explanation that as to why Zipfian distribution of relaxation rates appears for Curie-von Schweidler Law, and relate this law to time variant rate of relaxation. In this paper, we derive appearance of fractional derivative while using Zipfian power law distribution that gives notion of scale dependent relaxation rate function for Curie-von Schweidler relaxation phenomena. This paper gives analytical approach to get insight of a non-Debye relaxation and gives a new treatment to especially much used empirical Curie-von Schweidler (universal) relaxation law.展开更多
A new active control method was proposed, in which the analytical control law was deduced by using a step by step integral method to differential equation of motion under the condition of static error being zero. This...A new active control method was proposed, in which the analytical control law was deduced by using a step by step integral method to differential equation of motion under the condition of static error being zero. This control law is terse in mathematical expression and convenient for practical use. The simulation results demonstrate that the proposed method can provide much more remarkable peak response reduction of seismically excited structures than the classical LQR method.展开更多
In this paper, we propose a variational integrator for nonlinear Schrodinger equations with variable coefficients. It is shown that our variational integrator is naturally multi-symplectic. The discrete multi-symplect...In this paper, we propose a variational integrator for nonlinear Schrodinger equations with variable coefficients. It is shown that our variational integrator is naturally multi-symplectic. The discrete multi-symplectic structure of the integrator is presented by a multi-symplectic form formula that can be derived from the discrete Lagrangian boundary function. As two examples of nonlinear Schrodinger equations with variable coefficients, cubic nonlinear Schrodinger equations and Gross-Pitaevskii equations are extensively studied by the proposed integrator. Our numerical simulations demonstrate that the integrator is capable of preserving the mass, momentum, and energy conservation during time evolutions. Convergence tests are presented to verify that our integrator has second-order accuracy both in time and space.展开更多
文摘Missile interception problem can be regarded as a two-person zero-sum differential games problem,which depends on the solution of Hamilton-Jacobi-Isaacs(HJI)equa-tion.It has been proved impossible to obtain a closed-form solu-tion due to the nonlinearity of HJI equation,and many iterative algorithms are proposed to solve the HJI equation.Simultane-ous policy updating algorithm(SPUA)is an effective algorithm for solving HJI equation,but it is an on-policy integral reinforce-ment learning(IRL).For online implementation of SPUA,the dis-turbance signals need to be adjustable,which is unrealistic.In this paper,an off-policy IRL algorithm based on SPUA is pro-posed without making use of any knowledge of the systems dynamics.Then,a neural-network based online adaptive critic implementation scheme of the off-policy IRL algorithm is pre-sented.Based on the online off-policy IRL method,a computa-tional intelligence interception guidance(CIIG)law is developed for intercepting high-maneuvering target.As a model-free method,intercepting targets can be achieved through measur-ing system data online.The effectiveness of the CIIG is verified through two missile and target engagement scenarios.
文摘This paper centers on the integrated learning of English and law in China.Firstly,it outlines the importance of English in the solution of the ever increasing legal disputes between China and the outside world,which inevitably involves an integrated learning of English and law.Secondly,it points out that the content of legal English reflects a combination of legal knowledge and English skills.Thirdly,it expounds on the difficulties that Chinese English majors are facing in the process of learning English and law simultaneously and furnishes some practical suggestions.
基金Natural Science Foundation of High Education of Jiangsu Province of China,"Qing Lan" Project Foundation of Jiangsu Province
文摘In this paper, we present a new method to construct the conservation laws for relativistic mechanical systems by finding corresponding integrating factors. First, the Lagrange equations of relativistic mechanical systems are established, and the definition of integrating factors of the systems is given; second, the necessary conditions for the existence of conserved quantities of the relativistic mechanical systems are studied in detail, and the relation between the conservation laws and the integrating factors of the systems is obtained and the generalized Killing equations for the determination of the integrating factors are given; finally, the conservation theorem and its inverse for the systems are established, and an example is given to illustrate the application of the results.
基金supported by National Natural Science Foundation of China under Grant No. 10572021
文摘In this paper, the conservation laws of generalized Birkhoff system in event space are studied by using the method of integrating factors. Firstly, the generalized Pfaff-Birkhoff principle and the generalized Birkhoff equations are established, and the definition of the integrating factors for the system is given. Secondly, based on the concept of integrating factors, the conservation theorems and their inverse for the generalized Birkhoff system in the event space are presented in detail, and the relation between the conservation laws and the integrating factors of the system is obtained and the generalized Killing equations for the determination of the integrating factors are given. Finally, an example is given to illustrate the application of the results.
基金supported by the National Outstanding Young Scientists Fund of China (No. 10725209)the National ScienceFoundation of China (No. 10672092)+1 种基金Shanghai Municipal Education Commission Scientific Research Project (No. 07ZZ07)Shanghai Leading Academic Discipline Project (No. Y0103).
文摘A numerical method is proposed to simulate the transverse vibrations of a viscoelastic moving string constituted by an integral law. In the numerical computation, the Galerkin method based on the Hermite functions is applied to discretize the state variables, and the Runge- Kutta method is applied to solve the resulting differential-integral equation system. A linear iterative process is designed to compute the integral terms at each time step, which makes the numerical method more efficient and accurate. As examples, nonlinear parametric vibrations of an axially moving viscoelastic string are analyzed.
基金Supported by the National Natural Science Foundation of China under Grant Nos.11171041 and 10971018the Natural Science Foundation of Shandong Province under Grant No.ZR2010AM029+1 种基金the Promotive Research Fund for Young and Middle-Aged Scientists of Shandong Province under Grant No.BS2010SF001 the Doctoral Foundation of Binzhou University under Grant No.2009Y01
文摘This paper is concerned with the generalized nonlinear second-order equation.By the direct construction method,all of the first-order multipliers of the equation are obtained,and the corresponding complete conservation laws(CLs) of such equations are provided.Furthermore,the integrability of the equation is considered in terms of the conservation laws.In addition,the relationship of multipliers and symmetries of the equations is investigated.
文摘We formulate a “Moore’s law” for photonic integrated circuits (PICs) and their spatial integration density using two methods. One is decomposing the integrated photonics devices of diverse types into equivalent basic elements, which makes a comparison with the generic elements of electronic integrated circuits more meaningful. The other is making a complex compo- nent equivalent to a series of basic elements of the same functionality, which is used to calculate the integration density for func- tional components realized with different structures. The results serve as a benchmark of the evolution of PICs and we can con- clude that the density of integration measured in this way roughly increases by a factor of 2 per year. The prospects for a continued increase of spatial integration density are discussed.
基金Project supported by the Postdoctoral Science Foundation of China (Grant No. 2011M500404 )the Program for Liaoning Excellent Talents in University,China (Grant No. LJQ2011119)
文摘We construct a nonlinear integrable coupling of discrete soliton hierarchy, and establish the infinite conservation laws (CLs) for the nonlinear integrable coupling of the lattice hierarchy. As an explicit application of the method proposed in the paper, the infinite conservation laws of the nonlinear integrable coupling of the Volterra lattice hierarchy are presented.
基金Supported by the National Natural Science Foundation of China under Grant No.60772023the Open Fund No.BUAA-SKLSDE-09KF-04+2 种基金Supported Project No.SKLSDE-2010ZX-07 of the State Key Laboratory of Software Development Environment,Beijing University of Aeronautics and Astronauticsthe National Basic Research Program of China (973 Program) under Grant No.2005CB321901 the Specialized Research Fund for the Doctoral Program of Higher Education under Grant No.200800130006,Chinese Ministry of Education
文摘Under investigation in this paper are two coupled integrable dispersionless (CID) equations modelingthe dynamics of the current-fed string within an external magnetic field.Through a set of the dependent variabletransformations, the bilinear forms for the CID equations are derived.Based on the Hirota method and symboliccomputation, the analytic N-soliton solutions are presented.Infinitely many conservation laws for the CID equationsare given through the known spectral problem.Propagation characteristics and interaction behaviors of the solitons areanalyzed graphically.
基金supported by the "Chunlei" Project of Shandong University of Science and Technology of China under Grant No. 2008BWZ070
文摘Two hierarchies of nonlinear integrable positive and negative lattice equations are derived from a discrete spectrak problem. The two lattice hierarchies are proved to have discrete zero curvature representations associated with a discrete spectral problem, which also shows that the positive and negative hierarchies correspond to positive and negative power expansions of Lax operators with respect to the spectral parameter, respectively. Moreover, the integrable lattice models in the positive hierarchy are of polynomial type, and the integrable lattice models in the negative hierarchy are of rational type. Further, we construct infinite conservation laws about the positive hierarchy.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11971441,11871440,and 11931017)Key Scientific Research Projects of Colleges and Universities in Henan Province,China(Grant No.20A110006).
文摘With the aid of Lenard recursion equations, an integrable hierarchy of nonlinear evolution equations associated with a 2 × 2 matrix spectral problem is proposed, in which the first nontrivial member in the positive flows can be reduced to a new generalization of the Wadati–Konno–Ichikawa(WKI) equation. Further, a new generalization of the Fokas–Lenells(FL) equation is derived from the negative flows. Resorting to these two Lax pairs and Riccati-type equations, the infinite conservation laws of these two corresponding equations are obtained.
文摘The conservation law of J-integral in two-media with a crack paralleling to the interface of the two media was firstly proved by analytical and numerical finite element method. Then a schedule model was established that an interface crack is inserted in four media. According to the J-integral conservation law on multi-media, the energy release ratio of Ⅰ-type crack was considered to be conservation when the middle medium layers are very thin. And the conservation law was also convinced by numerical method. By means of the dimension analysis on the model, the asymptotic results and formula calculating the energy release ratio and complex stress intensity factor are presented.
文摘As the international community place high value on the ocean,protecting the marine environment and promoting the sustainable development of the ocean have become the consensus of all countries around the world.Integrated ocean management,as an important part of the ocean governance mechanism within the framework of the United Nations,have approached to the scene view of human gradually.Taking the United Nations Convention on the Law of the Sea as the framework,the current international integrated ocean management is summarized and the future development trend of ocean management is analyzed,providing reference for China to build a maritime power and a maritime community with shared future.
文摘In modem society, the rule or law is an important pattern of national governance. Realizing human rights is closely related to the rule of law. What the author intends to discuss is how to use integral thinking to understand human rights as well as the relationship between human rights and development and how to comprehend the rule of law and the promotion of human rights through the rule of law in the process of facilitating the cause of human rights.
基金The project supported by the Scientific Research Award Foundation for Outstanding Young and Middle-Aged Scientists of Shandong Province of China
文摘A discrete matrix spectral problem and the associated hierarchy of Lax integrable lattice equations are presented, and it is shown that the resulting Lax integrable lattice equations are all Liouville integrable discrete Hamiltonian systems. A new integrable symplectic map is given by binary Bargmann constraint of the resulting hierarchy. Finally, an infinite set of conservation laws is given for the resulting hierarchy.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11501520 and 11331008)the Outstanding Young Talent Research Fund of Zhengzhou University(Grant No.1521315001)
文摘A novel hierarchy of integrable nonlinear evolution equations related to the combined Ablowitz–Kaup–Newell–Segur(AKNS) and Wadati–Konno–Ichikawa(WKI) spectral problems is proposed,from which the Lax pair for a corresponding negative flow and its infinite many conservation laws are obtained.Furthermore,a reduction of this hierarchy is discussed,by which a generalized sinh-Gordon equation is derived on the basis of its negative flow.
文摘The classical power law relaxation, i.e. relaxation of current with inverse of power of time for a step-voltage excitation to dielectric—as popularly known as Curie-von Schweidler law is empirically derived and is observed in several relaxation experiments on various dielectrics studies since late 19th Century. This relaxation law is also regarded as “universal-law” for dielectric relaxations;and is also termed as power law. This empirical Curie-von Schewidler relaxation law is then used to derive fractional differential equations describing constituent expression for capacitor. In this paper, we give simple mathematical treatment to derive the distribution of relaxation rates of this Curie-von Schweidler law, and show that the relaxation rate follows Zipf’s power law distribution. We also show the method developed here give Zipfian power law distribution for relaxing time constants. Then we will show however mathematically correct this may be, but physical interpretation from the obtained time constants distribution are contradictory to the Zipfian rate relaxation distribution. In this paper, we develop possible explanation that as to why Zipfian distribution of relaxation rates appears for Curie-von Schweidler Law, and relate this law to time variant rate of relaxation. In this paper, we derive appearance of fractional derivative while using Zipfian power law distribution that gives notion of scale dependent relaxation rate function for Curie-von Schweidler relaxation phenomena. This paper gives analytical approach to get insight of a non-Debye relaxation and gives a new treatment to especially much used empirical Curie-von Schweidler (universal) relaxation law.
文摘A new active control method was proposed, in which the analytical control law was deduced by using a step by step integral method to differential equation of motion under the condition of static error being zero. This control law is terse in mathematical expression and convenient for practical use. The simulation results demonstrate that the proposed method can provide much more remarkable peak response reduction of seismically excited structures than the classical LQR method.
基金supported by the National Natural Science Foundation of China(Grant No.11401259)the Fundamental Research Funds for the Central Universities,China(Grant No.JUSRR11407)
文摘In this paper, we propose a variational integrator for nonlinear Schrodinger equations with variable coefficients. It is shown that our variational integrator is naturally multi-symplectic. The discrete multi-symplectic structure of the integrator is presented by a multi-symplectic form formula that can be derived from the discrete Lagrangian boundary function. As two examples of nonlinear Schrodinger equations with variable coefficients, cubic nonlinear Schrodinger equations and Gross-Pitaevskii equations are extensively studied by the proposed integrator. Our numerical simulations demonstrate that the integrator is capable of preserving the mass, momentum, and energy conservation during time evolutions. Convergence tests are presented to verify that our integrator has second-order accuracy both in time and space.