Multi-symplectic method for generalized Boussinesq equation

The generalized Boussinesq equation that represents a group of important nonlinear equations possesses many interesting properties. Multi-symplectic formulations of the generalized Boussinesq equation in the Hamilton space are introduced in this paper. And then an implicit multi-symplectic scheme equivalent to the multi-symplectic Box scheme is constructed to solve the partial differential equations （PDEs） derived from the generalized Boussinesq equation. Finally, the numerical experiments on the soliton solutions of the generalized Boussinesq equation are reported. The results show that the multi-symplectic method is an efficient algorithm with excellent long-time numerical behaviors for nonlinear partial differential equations.     

New explicit multi-symplectic scheme for nonlinear wave equation简

Based on splitting multi-symplectic structures, a new multi-symplectic scheme is proposed and applied to a nonlinear wave equation. The explicit multi-symplectic scheme of the nonlinear wave equation is obtained, and the corresponding multi-symplectic conservation property is proved. The backward error analysis shows that the explicit multi-symplectic scheme has good accuracy. The sine-Gordon equation and the Klein-Gordon equation are simulated by an explicit multi-symplectic scheme. The numerical results show that the new explicit multi-symplectic scheme can well simulate the solitary wave behaviors of the nonlinear wave equation and approximately preserve the relative energy error of the equation.     

Multi-symplectic method for the generalized (2+1)-dimensional KdV-mKdV equation

In the present paper,a general solution involving three arbitrary functions for the generalized(2+1)dimensional KdV-mKdV equation,which is derived fromthe generalized(1+1)-dimensional KdV-mKdV equation,is first introduced by means of the Wiess,Tabor,Carnevale(WTC) truncation method.And then multisymplectic formulations with several conservation lawstaken into account are presented for the generalized(2+1)dimensional KdV-mKdV equation based on the multisymplectic theory of Bridges.Subsequently,in order tosimulate the periodic wave solutions in terms of rationalfunctions of the Jacobi elliptic functions derived from thegeneral solution,a semi-implicit multi-symplectic schemeis constructed that is equivalent to the Preissmann scheme.From the results of the numerical experiments,we can conclude that the multi-symplectic schemes can accurately simulate the periodic wave solutions of the generalized(2+1)dimensional KdV-mKdV equation while preserve approximately the conservation laws.     

Analyzing dynamic response of non-homogeneous string fixed at both ends

To survey the local conservation properties of complex dynamic problems, a structure-preserving numerical method, named as generalized multi-symplectic method, is proposed to analyze the dynamic response of a non-homogeneous string fixed at both ends. Firstly, based on the multi-symplectic idea, a generalized multi-symplectic form derived from the vibration equation of a non-homogeneous string is presented. Secondly, several conservation laws are deduced from the generalized multi-symplectic form to illustrate the local properties of the system. Thirdly, a centered box difference scheme satisfying the discrete local momentum conservation law exactly, named as a generalized multi-symplectic scheme, is constructed to analyze the dynamic response of the non-homogeneous string fixed at both ends. Finally, numerical experiments on the generalized multi-symplectic scheme are reported. The results illustrate the high accuracy, the good local conservation properties as well as the excellent long-time numerical behavior of the generalized multi-symplectic scheme well.     

非线性弹性杆中纵波传播过程的数值模拟

基于Hamilton空间体系的多辛理论研究了包含材料非线性效应和几何弥散效应的非线性弹性杆中纵波传播问题。导出了其Bridges意义下的多辛形式及其多种守恒律,并构造了等价于Box多辛格式的隐式多辛格式对不考虑材料非线性效应和几何弥散效应、只考虑材料非线性效应、只考虑几何弥散效应、同时考虑材料非线性效应和几何弥散效应四种情况下不同截面参数的圆杆中的纵波传播过程进行数值模拟,模拟结果不仅全面地反映了非线性效应和几何弥散效应对纵波传播的影响,而且也反映了多辛方法的两大优点:精确的保持多种守恒律和良好的长时间数值行为。     

THE MULTI-SYMPLECTIC ALGORITHM FOR "GOOD" BOUSSINESQ EQUATION

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｒｅｃｅｎｔｙｅａｒｓａｒｅｍａｒｋａｂｌｅｄｅｖｅｌｏｐｍｅｎｔｈａｓｔａｋｅｎｐｌａｃｅｉｎｔｈｅｓｔｕｄｙｏｆｎｏｎｌｉｎｅａｒｅｖｏｌｕｔｉｏｎａｒｙｐａｒｔｉａｌｄｉｆｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｓ.Ａｎｅｘａｍｐｌｅｉｓｔｈｅ“Ｇｏｏｄ”Ｂｏｕｓｓｉｎｅｓｑ (Ｇ .Ｂ .)ｅｑｕａｔｉｏｎｕｔｔ =-ｕｘｘｘｘ+ｕｘｘ+ (ｕ2 ) ｘｘ ( 1 )ｗｈｉｃｈｄｅｓｃｒｉｂｅｓｓｈａｌｌｏｗｗａｔｅｒｗａｖｅｓｐｒｏｐａｇａｔｉｎｇｉｎｂｏｔｈｄｉｒｅｃｔｉｏｎｓ.Ｔｈｅａ…     

The evolution of nonlinear waves in active media

An evolution equation that describes the behaviour of single waves in active media is derived. The basic mathematical model corresponds to pulse transmission in nerve fibers according to the hyperbolic telegraph equation. A numerical experiment is carried out in which the evolution equation is solved by the pseudospectral method and the corresponding stationary wave equation by the standard Runge-Kutta method. The evolution equation has a stationary solution in the form of an unsymmetric solitary wave with a refractive tail. The numerical simulation of the process gives physically admissible results, namely, the suitable wave profile and the existence of the threshold and asymptotic values. The situation analyzed here in an active medium with energy influx is in a certain sense similar to the formation of solitary waves in a conservative medium.     

Multi-symplectic methods for membrane free vibration equation

In this paper,the multi-symplectic formulations of the membrane free vi- bration equation with periodic boundary conditions in Hamilton space are considered. The complex method is introduced and a semi-implicit twenty-seven-points scheme with certain discrete conservation laws—a multi-symplectic conservation law(CLS),a local energy conservation law(ECL)as well as a local momentum conservation law(MCL)—is constructed to discrete the PDEs that are derived from the membrane free vibra- tion equation.The results of the numerical experiments show that the multi-symplectic scheme has excellent long-time numerical behavior.     

THE MULTI-SYMPLECTIC ALGORITHM FOR “GOOD” BOUSSINESQ EQUATION

The multi-symplectic formulations of the “Good” Boussinesq equation were considered. For the multi-symplectic formulation, a new fifteen-point difference scheme which is equivalent to the multi-symplectic Preissman integrator was derived. The numerical experiments show that the multi-symplectic scheme have excellent long-time numerical behavior. Foundation items: the Foundation for Key Laboratory of Scientific/Engineering Computing Institute of Computational Mathematics and Scientific/Engineering Computing, Chinese Academy of Sciences; the Natural Science Foundation of Huaqiao University. Biography: ZENG Wen-ping (1940-), Professor (E-mail: qmz@1sec.cc.ac.cn)     

Particle trajectories under interactions between solitary waves and a linear shear current

This paper is concerned with particle trajectories beneath solitary waves when a linear shear current exists. The fluid is assumed to be incompressible and inviscid, lying on a flat bed. Classical asymptotic expansion is used to obtain a Korteweg-de Vries(Kd V) equation, then a forth-order Runge-Kutta method is applied to get the approximate particle trajectories. On the other hand, our particular attention is paid to the direct numerical simulation(DNS) to the original Euler equations. A conformal map is used to solve the nonlinear boundary value problem. Highaccuracy numerical solutions are then obtained through the fast Fourier transform(FFT) and compared with the asymptotic solutions, which shows a good agreement when wave amplitude is small. Further, it also yields that there are different types of particle trajectories. Most surprisingly,periodic motion of particles could exist under solitary waves, which is due to the wave-current interaction.     

Dynamic symmetry breaking and structure-preserving analysis on the longitudinal wave in an elastic rod with a variable cross-section

The longitudinal wave propagating in an elastic rod with a variable cross-section owns wide engineering background, in which the longitudinal wave dissipation determines some important performances of the slender structure. To reproduce the longitudinal wave dissipation effects on an elastic rod with a variable cross-section, a structure-preserving approach is developed based on the dynamic symmetry breaking theory. For the dynamic model controlling the longitudinal wave propagating in the elast...     

轴向运动复合材料圆柱壳的非线性振动研究

采用Runge–Kutta法和多尺度法对轴向运动分层复合材料薄壁圆柱壳的非线性振动特性进行了研究。首先根据层合壳理论建立轴向运动分层复合材料薄壁圆柱壳的波动方程,利用Galerkin法对方程进行离散,得到相互耦合模态方程组。然后应用Runge –Kutta法分析了不同参数条件下的幅频特性曲线,得到了系统由于固有频率接近所导致的内共振现象,以及系统呈现软特性等非线性特性。最后采用多尺度法进行了系统1:1内共振时的近似解析分析,对系统在不同参数下的振动研究表明,激振力幅值、阻尼、速度等参数对位移响应幅值、共振区间、模态间的耦合度及系统软特性程度均有影响,其结论与数值计算结果一致,并同时对解的稳定性进行了研究。     

Numerical simulation of standing wave with 3D predictor-corrector finite difference method for potential flow equations

A three-dimensional （3D） predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.     

Hamilton体系下介电弹性体圆形薄膜的动力学建模与辛求解

采用辛算法研究了Hamilton体系下介电弹性体圆形薄膜的动力学响应。首先,将该问题引入Hamilton对偶变量体系,借助Legendre变换,给出系统的广义动量和Hamilton函数,通过对Hamilton函数作用量的变分,得到Hamilton体系下的正则方程。其次,对于得到的正则方程给出了辛Runge-Kutta的计算格式。最后,采用二级四阶辛Runge-Kutta算法对动力学系统进行了数值求解,和四级四阶经典Runge-Kutta算法进行对比,结果表明,二级四阶辛Runge-Kutta算法具有保能量以及长时间数值稳定的优势,同时说明四级四阶经典Runge-Kutta算法对于步长依赖的局限性。     

A WAVE EQUATION MODEL TO SOLVE THE MULTIDIMENSIONAL TRANSPORT EQUATION

The wave equation model, originally developed to solve the advection–diffusion equation, is extended to the multidimensional transport equation in which the advection velocities vary in space and time. The size of the advection term with respect to the diffusion term is arbitrary. An operator-splitting method is adopted to solve the transport equation. The advection and diffusion equations are solved separate ly at each time step. During the advection phase the advection equation is solved using the wave equation model. Consistency of the first-order advection equation and the second-order wave equation is established. A finite element method with mass lumping is employed to calculate the three-dimensional advection of both a Gaussian cylinder and sphere in both translational and rotational flow fields. The numerical solutions are accurate in comparison with the exact solutions. The numerical results indicate that (i) the wave equation model introduces minimal numerical oscillation, (ii) mass lumping reduces the computational costs and does not significantly degrade the numerical solutions and (iii) the solution accuracy is relatively independent of the Courant number provided that a stability constraint is satisfied. © 1997 by John Wiley & Sons, Ltd.     

外激励作用下亚音速二维壁板的复杂响应研究

研究了亚音速流中二维壁板在外激励作用下的复杂响应问题。采用迦辽金方法将非线性运动控制方程离散为常微分方程组,采用数值方法进行计算,研究了壁板系统的复杂响应。应用最大李亚普诺夫指数和庞加莱截面方法对系统的运动性质进行了判定。结果表明,系统随着参数的变化呈现出复杂的响应,系统的周期运动与混沌运动会相间出现;系统由周期运动进...     

基于S-A湍流模型的Runge-Kutta有限元算法

采用一方程S-A模型(Spalart-Allmaras模型)封闭雷诺时均N-S方程(RANS方程)进行湍流数值计算,可以减少方程求解数量,节约计算时间。本文对其进行了有限元数值算法研究,首先通过沿流线坐标变换,得到无对流项RANS方程,并引入三阶Runge-Kutta法对其进行时间离散;然后利用沿流线的Taylor展开解决坐标变换带来的网格更新的困难;最后采用Galerkin法进行空间离散,得到湍流模型的有限元算法。基于方柱绕流和覆冰输电线绕流模型,与试验结果进行对比,验证了该算法的有效性,与一阶数值算法相比,该算法在精度和收敛性方面更具优势。     

NUMERICAL SIMULATION OF SHALLOW WATER WAVE PROPAGATION USING A BOUNDARY ELEMENT WAVE EQUATION MODEL

The present paper makes use of a wave equation formulation of the primitive shallow water equations to simulate one-dimensional free surface flow. A numerical formulation of the boundary element method is then developed to solve the wave continuity equation using a time-dependent fundamental solution, while an explicit finite difference scheme is used to derive velocities from the primitive momentum equation. One-dimensional free surface flows in open channels are treated and the results compared with analytical and numerical solutions. © 1997 John Wiley & Sons, Ltd.