Two types of first-order decoupled conservative schemes are firstly proposed for the strongly coupled nonlinear Schrodinger system by using pseudospectral method in space and coordinate increment discrete gradient met...
详细信息
Two types of first-order decoupled conservative schemes are firstly proposed for the strongly coupled nonlinear Schrodinger system by using pseudospectral method in space and coordinate increment discrete gradient method in time. And then, in order to improve the solution accuracy in time, the composition methods are employed to construct secondand fourth-order schemes. The proposed schemes are efficient for the system in d >= 2 dimensions and also easy to code because of their decoupled feature. A fast solver is proposed to speed up the computation. Ample numerical examples including the motion of single soliton and interaction of multiple solitary waves are carried out to exhibit the performance of the schemes. (C) 2020 IMACS. Published by Elsevier B.V. All rights reserved.
In this paper, we use the Galerkin weak form to construct a structure-preserving scheme for Klein-Gordon-Schrodinger equation and analyze its conservative and convergent properties. We first discretize the underlying ...
详细信息
In this paper, we use the Galerkin weak form to construct a structure-preserving scheme for Klein-Gordon-Schrodinger equation and analyze its conservative and convergent properties. We first discretize the underlying equation in space direction via a selected finite element method, and the Hamiltonian partial differential equation can be casted into Hamiltonian ordinary differential equations based on the weak form of the system afterwards. Then, the resulted ordinary differential equations are solved by the symmetric discrete gradient method, which yields a charge-preserving and energy-preserving scheme. Moreover, the numerical solution of the proposed scheme is proved to be bounded in the discrete L infinity norm and convergent with the convergence order of O(h2+tau 2) in the discrete L2 norm without any grid ratio restrictions, where h and tau are space and time step, respectively. Numerical experiments conducted last to verify the theoretical analysis.
For more than half a century, most of the plasma scientists have encountered a violation of the conservation laws of charge, momentum, and energy whenever they have numerically solved the first-principle equations of ...
详细信息
For more than half a century, most of the plasma scientists have encountered a violation of the conservation laws of charge, momentum, and energy whenever they have numerically solved the first-principle equations of kinetic plasmas, such as the relativistic Vlasov-Maxwell system. This fatal problem is brought by the fact that both the Vlasov and Maxwell equations are indirectly associated with the conservation laws by means of some mathematical manipulations. Here we propose a quadratic conservative scheme, which can strictly maintain the conservation laws by discretizing the relativistic Vlasov-Maxwell system. A discrete product rule and summation-by-parts are the key players in the construction of the quadratic conservative scheme. Numerical experiments of the relativistic two-stream instability and relativistic Weibel instability prove the validity of our computational theory, and the proposed strategy will open the doors to the first-principle studies of mesoscopic and macroscopic plasma physics. (C) 2018 The Author(s). Published by Elsevier Inc.
structure-preserving algorithms for solving conservative PDEs with added linear dissipation are generalized to systems with time dependent damping/driving terms. This study is motivated by several PDE models of physic...
详细信息
structure-preserving algorithms for solving conservative PDEs with added linear dissipation are generalized to systems with time dependent damping/driving terms. This study is motivated by several PDE models of physical phenomena, such as Korteweg-de Vries, Klein-Gordon, Schrodinger, and Camassa-Holm equations, all with damping/driving terms and time-dependent coefficients. Since key features of the PDEs under consideration are described by local conservation laws, which are independent of the boundary conditions, the proposed (second-order in time) discretizations are developed with the intent of preserving those local conservation laws. The methods are respectively applied to a damped-driven nonlinear Schrodinger equation and a damped Camassa-Holm equation. Numerical experiments illustrate the structure-preserving properties of the methods, as well as favorable results over other competitive schemes. (C) 2018 Elsevier B.V. All rights reserved.
Symplectic schemes are powerful methods for numerically integrating Hamiltonian systems, and their long-term accuracy and fidelity have been proved both theoretically and numerically. However direct applications of st...
详细信息
Symplectic schemes are powerful methods for numerically integrating Hamiltonian systems, and their long-term accuracy and fidelity have been proved both theoretically and numerically. However direct applications of standard symplectic schemes to relativistic charged particle dynamics result in implicit and electromagnetic gauge-dependent algorithms. In the present study, we develop explicit high-order gauge-independent noncanonical symplectic algorithms for relativistic charged particle dynamics using a Hamiltonian splitting method in the 8D phase space. It is also shown that the developed algorithms can be derived as variational integrators by appropriately discretizing the action of the dynamics. Numerical examples are presented to verify the excellent long-term behavior of the algorithms. (C) 2019 Elsevier B.V. All rights reserved.
We give a systematic method for discretizing Hamiltonian partial differential equations. An application of the method to the Hamiltonian form of the coupled Schrodinger-KdV equations yields a temporal first-order cons...
详细信息
We give a systematic method for discretizing Hamiltonian partial differential equations. An application of the method to the Hamiltonian form of the coupled Schrodinger-KdV equations yields a temporal first-order conservative scheme. Temporal second- and fourth-order schemes are developed by employing the composition method. All the schemes are decoupled and exactly conserve three invariants simultaneously. Numerical results show good performance of the schemes and verify the theoretical results. (C) 2018 Elsevier Ltd. All rights reserved.
Based on the multi-symplectic Hamiltonian formula of the generalized Rosenau-type equation, a multi-symplectic scheme and an energy-preserving scheme are proposed. To improve the accuracy of the solution, we apply the...
详细信息
Based on the multi-symplectic Hamiltonian formula of the generalized Rosenau-type equation, a multi-symplectic scheme and an energy-preserving scheme are proposed. To improve the accuracy of the solution, we apply the composition technique to the obtained schemes to develop high-order schemes which are also multi-symplectic and energy-preserving respectively. Discrete fast Fourier transform makes a significant improvement to the computational efficiency of schemes. Numerical results verify that all the proposed schemes have satisfactory performance in providing accurate solution and preserving the discrete mass and energy invariants. Numerical results also show that although each basic time step is divided into several composition steps, the computational efficiency of the composition schemes is much higher than that of the non-composite schemes. (C) 2017 Elsevier B.V. All rights reserved.
We develop two local energy-preserving integrators and a global energy-preserving integrator for the general multisymplectic Hamiltonian system. When applied to the 1D and multi-dimensional N-coupled nonlinear Schrodi...
详细信息
We develop two local energy-preserving integrators and a global energy-preserving integrator for the general multisymplectic Hamiltonian system. When applied to the 1D and multi-dimensional N-coupled nonlinear Schrodinger equations, the given schemes have the exact preservation of the local/global conservation law and are decoupled in the components psi(n,) n = 1, 2, ... , N, i.e., each of the components can be solved independently. The decoupled feature is significant and helpful for overcoming the computational difficulty of the N-coupled (N >= 3) nonlinear Schrodinger equations, especially of the multi-dimensional case. The composition method is employed to improve the accuracy of the schemes in time and the discrete fast Fourier transform is used to reduce the computational complexity. Several numerical experiments are carried out to exhibit the behaviors of the wave solutions. Numerical results confirm the theoretical results. (C) 2018 Elsevier Inc. All rights reserved.
The multisymplectic schemes have been used in numerical simulations for the RLW-type equation successfully. They well preserve the local geometric property, but not other local conservation laws. In this article, we p...
详细信息
The multisymplectic schemes have been used in numerical simulations for the RLW-type equation successfully. They well preserve the local geometric property, but not other local conservation laws. In this article, we propose three novel efficient local structure-preserving schemes for the RLW-type equation, which preserve the local energy exactly on any time-space region and can produce richer information of the original problem. The schemes will be mass- and energy-preserving as the equation is imposed on appropriate boundary conditions. Numerical experiments are presented to verify the efficiency and invariant-preserving property of the schemes. Comparisons with the existing nonconservative schemes are made to show the behavior of the energy affects the behavior of the solution.(c) 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1678-1691, 2017
In this article, we obtain local energy and momentum conservation laws for the Klein-Gordon-Schrodinger equations, which are independent of the boundary condition and more essential than the global conservation laws. ...
详细信息
In this article, we obtain local energy and momentum conservation laws for the Klein-Gordon-Schrodinger equations, which are independent of the boundary condition and more essential than the global conservation laws. Based on the rule that the numerical methods should preserve the intrinsic properties as much as possible, we propose local energy- and momentum-preserving schemes for the equations. The merit of the proposed schemes is that the local energy/momentum conservation law is conserved exactly in any time-space region. With suitable boundary conditions, the schemes will be charge- and energy-/momentum-preserving. Nonlinear analysis shows LEP schemes are unconditionally stable and the numerical solutions converge to the exact solutions with order O(2+h2). The theoretical properties are verified by numerical experiments. (c) 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1329-1351, 2017
暂无评论