Numerical Complexiton Solutions of Complex KdV Equation

In this paper, we directly extend the applications of the Adomian decomposition method to investigate the complex KdV equation. By choosing different forms of wave functions as the initial values, three new types of realistic numerical solutions： numerical positon, negaton solution, and particularly the numerical analytical complexiton solution are obtained, which can rapidly converge to the exact ones obtained by Lou et al. Numerical simulation figures are used to illustrate the efficiency and accuracy of the proposed method.     

A Higher-Dimensional Hirota Condition and Its Judging Method

When a one-dimensional nonlinear evolution equation could be transformed into a bilinear differential form as F（Dt, Dx）f . f = O, Hirota proposed a condition for the above evolution equation to have arbitrary N-soliton solutions, we call it the 1-dimensional Hirota condition. As far as higher-dimensional nonlinear evolution equations go, a similar condition is established in this paper, also we call it a higher-dimensional Hirota condition, a corresponding judging theory is given. As its applications, a few two-dimensional KdV-type equations possessing arbitrary N-soliton solutions are obtained.     

Painleve Property and Complexiton Solutions of a Special Coupled KdV Equation

A special coupled KdV equation is proved to be the Painleve property by the Kruskal＇s simplification of WTC method. In order to search new exact solutions of the coupled KdV equation, Hirota＇s bilinear direct method and the conjugate complex number method of exponential functions are applied to this system. As a result, new analytical eomplexiton and soliton solutions are obtained synchronously in a physical field. Then their structures, time evolution and interaction properties are further discussed graphically.     

Analytic Solutions to Forced KdV Equation

The Korteweg-de Vries equation with a forcing term is established by recent studies as a simple mathematical model of describing the physics of a shallow layer of fluid subject to external forcing. In the present paper, we study the analytic solutions to the KdV equation with forcing term by using Hirota＇s direct method. Several exact solutions are given as examples, from which one can see that the same type soliton solutions can be excited by different forced term.     

Wronskian Form of N-Solitonic Solution for a Variable-Coefficient Korteweg-de Vries Equation with Nonuniformities

By the symbolic computation and Hirota method, the bilinear form and an auto-Backlund transformation for a variable-coemcient Korteweg-de Vries equation with nonuniformities are given. Then, the N-solitonic solution in terms of Wronskian form is obtained and verified. In addition, it is shown that the （N - 1）- and N-solitonic solutions satisfy the auto-Backlund transformation through the Wronskian technique.     

Similarity Reductions and Exact Solutions of a Coupled Korteweg de Vries System

For a special coupled Korteweg de Vries （KdV） system, its similarity solutions and reduction equations are obtained by the Clarkson and Kruskal＇s direct method. In addition, its new explicit soliton solutions and traveling wave solutions are found by the deformation and mapping method.     

New Exact Periodic Solitary-Wave Solutions for New （2＋1）-Dimensional KdV Equation

The extended homoclinic test function method is a kind of classic, efficient and well-developed method to solve nonlinear evolution equations. In this paper, with the help of this approach, we obtain new exact solutions （including kinky periodic solitary-wave solutions, periodic soliton solutions, and cross kink-wave solutions） for the new （2＋1）-dimensional KdV equation. These results enrich the variety of the dynamics of higher-dimensionai nonlinear wave field.     

Effects of Adiabatic Dust Charge Fluctuation and Particles Collisions on Dust-Acoustic Solitary Waves in Three-Dimensional Magnetized Dusty Plasmas

Taking into account the combined effects of the external magnetic field, adiabatic dust charge fluctuation and collisions occurring between the charged dust grains and neutral gas particles （dust-neutral collisions）, the dust-acoustic solitary waves in three-dimensional uniform dusty plasmas are investigated analytically. By using the reductive perturbation method, the Korteweg-de Vries （KdV） equation governing the dnst-aconstic solitary waves is obtained. The present analytical results show that only rarefactive solitary waves exist in this system. It is also found that the effects of the wave vector along the z-direction, dust charge variation, collisional frequency, the plasma density, and temperature ratio can significantly influence the characteristics of low-frequency wave modes. Moreover, for the collisional dusty plasmas, there is a certain critical value μc of the plasma density ratio μ, if μ 〈 μc, the width of the waves increases with μ, otherwise the width of waves decreases with μ.     

Similarity Reductions of Nearly Concentric KdV Equation

Basing on the direct method developed by Clarkson and Kruskal, the nearly concentric Korteweg-de Vries （ncKdV） equation can be reduced to three types of （1＋1）-dimensional variable coefficients partial differential equations （PDEs） and three types of variable coefficients ordinary differential equation. Furthermore, three types of （1＋1）-dimensional variable coefficients PDEs are all reduced to constant coefficients PDEs by some transformations.     

Solitary and Periodic Waves in Coupled KdV Equations with Different Linear Dispersion Relations

Based on the invariant expansion method, some reasonable approximate solutions of coupled Korteweg-de Vries （KdV） equations with different linear dispersion relations have been obtained. These solutions contain not only bell type soliton solutions but also periodic wave solutions that expressed by Jacob/elliptic functions. The results also show that if the arbitrary constants are selected suitably, the approximate solutions may become the exact ones.     

Algebraic-Geometric Solution to （2＋1）-Dimensional Sawada-Kotera Equation

Special solution of the （2＋1）-dimensional Sawada Kotera equation is decomposed into three （0＋1）- dimensional Bargmann flows. They are straightened out on the Jacobi variety of the associated hyperelliptic curve. Explicit algebraic-geometric solution is obtained on the basis of a deeper understanding of the KdV hierarchy.     

（G′/G）-Expansion Method Equivalent to Extended Tanh Function Method

In a recent article [Physics Letters A 372 （2008） 417], Wang et al. proposed a method, which is called the （G′/G）-expansion method, to look for travelling wave solutions of nonlinear evolution equations. The travelling wave solutions involving parameters of the KdV equation, the mKdV equation, the variant Boussinesq equations, and the Hirota-Satsuma equations are obtained by using this method. They think the （G′/G）-expansion method is a new method and more travelling wave solutions of many nonlinear evolution equations can be obtained. In this paper, we will show that the （G′/G）-expansion method is equivalent to the extended tanh function method.     

A New （2＋1）-Dimensional KdV Equation and Its Localized Structures

A new （2＋1）-dimensional KdV equation is constructed by using Lax pair generating technique. Exact solutions of the new equation are studied by means of the singular manifold method. Bgcklund transformation in terms of the singular manifold is obtained. And localized structures are also investigated.     

Impact of symmetrized and Burt—Foreman Hamiltonians on spurious solutions and energy levels of InAs/GaAs quantum dots

<正>We present a systematic investigation of calculating quantum dots(QDs) energy levels using the finite element method in the frame of the eight-band k·p method.Numerical results including piezoelectricity,electron and hole levels,as well as wave functions are achieved.In the calculation of energy levels,we do observe spurious solutions(SSs) no matter Burt-Foreman or symmetrized Hamiltonians are used.Different theories are used to analyse the SSs,we find that the ellipticity theory can give a better explanation for the origin of SSs and symmetrized Hamiltonian is easier to lead to SSs.The energy levels simulated with the two Hamiltonians are compared to each other after eliminating SSs,different Hamiltonians cause a larger difference on electron energy levels than that on hole energy levels and this difference decreases with the increase of QD size.     

Multi-physics analysis of the RFQ for Injector Scheme II of C-ADS driver linac

A 162.5 MHz, 2.1 MeV radio frequency quadruples （RFQ） structure is being designed for the Injector Scheme Ⅱ of the China Accelerator Driven Sub-critical System （C-ADS） linac. The RFQ will operate in continuous wave （CW） mode as required. For the CW normal conducting machine, the heat management will be one of the most important issues, since the temperature fluctuation may cause cavity deformation and lead to the resonant frequency shift. Therefore a detailed multi-physics analysis is necessary to ensure that the cavity can stably work at the required power level. The multi-physics analysis process includes RF electromagnetic analysis, thermal analysis, mechanical analysis, and this process will be iterated for several cycles until a satisfactory solution can be found. As one of the widely accepted measures, the cooling water system is used for frequency fine tunning, so the tunning capability of the cooling water system is also studied under different conditions. The results indicate that with the cooling water system, both the temperature rise and the frequency shift can be controlled at an acceptable level.     

Dynamic investigation of the finite dissolution of silicon particles in aluminum melt with a lower dissolution limit

The finite dissolution model of silicon particles in the aluminum melt is built and calculated by the finite difference method, and the lower dissolution limit of silicon particles in the aluminum melt is proposed and verified by experiments, which could be the origin of microinhomogeneity in aluminum-silicon melts. When the effects of curvature and interface reaction on dissolution are not considered; the dissolution rate first decreases and later increases with time. When the effects of curvature and interface reaction on dissolution are considered, the dissolution rate first decreases and later increases when the interface reaction coefficient （k） is larger than 10 1, and the dissolution rate first decreases and later tends to be constant when k is smaller than 10-3. The dissolution is controlled by both diffusion and interface reaction when k is larger than 10-3, while the dissolution is controlled only by the interface reaction when k is smaller than 10-4.     

Critical Behaviors in a Stochastic Local Limited One-Dimensional Rice-Pile Model

A stochastic local limited one-dimensional rice-pile model is numerically investigated. The distributions for avalanche sizes have a clear power-law behavior and it displays a simple finite size scaling. We obtain the avalanche exponents Ts= 1.54±0.10,βs = 2.17±0.10 and TT = 1.80±0.10, βT =1.46 ± 0.10. This self-organized critical model belongs to the same universality class with the Oslo rice-pile model studied by K. Christensen et al. [Phys. Rev. Lett. 77 （1996） 107], a rice-pile model studied by L.A.N. Amaral et al. [Phys. Rev. E 54 （1996） 4512], and a simple deterministic self-organized critical model studied by M.S. Vieira [Phys. Rev. E 61 （2000） 6056].     

Deforming the Travelling Wave Solutions of the KdV Equation to Those of the Generalized Fifth Order KdV Equation

Using some limiting procedures, the solutions of the fifth order KdV equation u_t + (μu²+ υu_xx + αuu_xx + βu_x² + γu³ + δu_xxxx)_x = 0 would degenerate into the solutions of a simple equation, say KdV equation. In this letter, we analyze the possibility of the inverse procedure of the limiting process mentioned above for the travelling wave solutions. The results show that the procedure for deforming a travelling wave solution of the KdV equation to that of the generalized fifth order KdV equation can be accomplished by some pure algebraic tricks. Moreover, this inverse procedure is not unique in general.     

Design on a highly birefringent and nonlinear photonic crystal fiber in the C waveband

In this paper, a novel photonic crystal fiber （PCF） with high birefringence and nonlinearity is designed. The charac- teristics of birefringence, dispersion and nonlinearity are studied by using the full-vector finite element method （FVFEM）. The numerical results show that the phase birefringence and nonlinear coefficient of PCF can be up to 4.51× 10-3 and 32.8972 w-l.km-1 at 1.55 μm, respectively. The proposed PCF could be found to have important applications in the polarization-dependent nonlinear optics such as the pulse compress and reshaping in the C waveband.     

Finite size specimens with cracks of icosahedral AI-Pd-Mn quasicrystals

Icosahedral quasicrystals are the most important and thermodynamically stable in all about 200 kinds of quasicrystals currently observed. Beyond the scope of classical elasticity, apart from a phonon displacement field, there is a phason displacement field in the elasticity of the quasicrystal, which induces an important effect on the mechanical properties of the material and makes an analytical solution difficult to obtain. In this paper, a finite element algorithm for the static elasticity of icosahedral quasicrystals is developed by transforming the elastic boundary value problem of the icosahedral quasicrystals into an equivalent variational problem. Analytical and numerical solutions for an icosahedral A1-Pd-Mn quasicrystal cuboid subjected to a uniaxial tension with different phonon-phason coupling parameters are given to verify the validity of the numerical approach. A comparison between the analytical and numerical solutions of the specimen demonstrates the accuracy and efficiency of the present algorithm. Finally, in order to reveal the fracture behavior of the icosahedral A1-Pd-Mn quasicrystal, a cracked specimen with a finite size of matter is investigated, both with and without phonon-phason coupling. Meanwhile, the geometry factors are calculated, including the stress intensity factor and the crack opening displacement for the finite-size specimen. Computational results reveal the importance of pbonon-phason coupling effect on the icosahedral A1-Pd-Mn quasicrystal. Furthermore, the finite element procedure can be used to solve more complicated boundary value problems.