scholarly journals Multi-Symplectic Method for the Logarithmic-KdV Equation

Symmetry ◽  
2020 ◽  
Vol 12 (4) ◽  
pp. 545
Author(s):  
Yu Zhang ◽  
Shaohua Li

The multi-symplectic integrator, as a numerical integration approach with symmetry, is known to have the characteristic of preserving the qualitative features and geometric properties of certain systems. Using the multi-symplectic integrator, the numerical simulation of the Gaussian solitary wave propagation of the logarithmic Korteweg–de Vries (logarithmic-KdV) equation was investigated. The multi-symplectic formulation of the logarithmic-KdV equation was explored by introducing some intermediate variables. A fully implicit version of the centered box scheme was used to discretize the multi-symplectic equations. In addition, numerical experiments were carried out to show the conservative properties of the proposed scheme.

2013 ◽  
Vol 14 (2) ◽  
pp. 393-411 ◽  
Author(s):  
Shanshan Jiang ◽  
Lijin Wang ◽  
Jialin Hong

AbstractIn this paper we propose stochastic multi-symplectic conservation law for stochastic Hamiltonian partial differential equations, and develop a stochastic multi-symplectic method for numerically solving a kind of stochastic nonlinear Schrödinger equations. It is shown that the stochastic multi-symplectic method preserves the multi-symplectic structure, the discrete charge conservation law, and deduces the recurrence relation of the discrete energy. Numerical experiments are performed to verify the good behaviors of the stochastic multi-symplectic method in cases of both solitary wave and collision.


Materials ◽  
2021 ◽  
Vol 14 (13) ◽  
pp. 3553
Author(s):  
Dengwang Wang ◽  
Yong Gao ◽  
Sheng Wang ◽  
Jie Wang ◽  
Haipeng Li

Carbon/Phenolic (C/P), a typical anisotropic material, is an important component of aerospace and often used to protect the thermodynamic effects of strong X-ray radiation. In this paper, we establish the anisotropic elastic-plastic constitutive model, which is embedded in the in-house code “RAMA” to simulate a two-dimensional thermal shock wave induced by X-ray. Then, we compare the numerical simulation results with the thermal shock wave stress generated by the same strong current electron beam via experiment to verify the correctness of the numerical simulation. Subsequently, we discuss and analyze the rules of thermal shock wave propagation in C/P material by further numerical simulation. The results reveal that the thermal shock wave represents different shapes and mechanisms by the radiation of 1 keV and 3 keV X-rays. The vaporization recoil phenomenon appears as a compression wave under 1 keV X-ray irradiation, and X-ray penetration is caused by thermal deformation under 3 keV X-ray irradiation. The thermal shock wave propagation exhibits two-dimensional characteristics, the energy deposition of 1 keV and 3 keV both decays exponentially, the energy deposition of 1 keV-peak soft X-ray is high, and the deposition depth is shallow, while the energy deposition of 3 keV-peak hard X-ray is low, and the deposition depth is deep. RAMA can successfully realize two-dimensional orthotropic elastoplastic constitutive relation, the corresponding program was designed and checked, and the calculation results for inspection are consistent with the theory. This study has great significance in the evaluation of anisotropic material protection under the radiation of intense X-rays.


Author(s):  
Z. Y. Song ◽  
C. Cheng ◽  
F. M. Xu ◽  
J. Kong

Based on the analytical solution of one-dimensional simplified equation of damping tidal wave and Heuristic stability analysis, the precision of numerical solution, computational time and the relationship between the numerical dissipation and the friction dissipation are discussed with different numerical schemes in this paper. The results show that (1) when Courant number is less than unity, the explicit solution of tidal wave propagation has higher precision and requires less computational time than the implicit one; (2) large time step is allowed in the implicit scheme in order to reduce the computational time, but the precision of the solution also reduce and the calculation precision should be guaranteed by reducing the friction factor: (3) the friction factor in the implicit solution is related to Courant number, presented as the determined friction factor is smaller than the natural value when Courant number is larger than unity, and their relationship formula is given from the theoretical analysis and the numerical experiments. These results have important application value for the numerical simulation of the tidal wave.


Sign in / Sign up

Export Citation Format

Share Document