弹性力学的重构核粒子边界无单元法

将重构核粒子法(RKPM)和边界积分方程方法结合，提出了一种新的边界积分方程无网格方法——重构核粒子边界无单元法(RKP-BEFM).对弹性力学问题，推导了其重构核粒子边界无单元法的公式，研究其数值积分方案，建立了重构核粒子边界无单元法离散化边界积分方程，并推导了重构核粒子边界无单元法的内点位移和应力积分公式.重构核粒子法形成的形函数具有重构核函数的光滑性，且能再现多项式在插值点的精确值，所以本方法具有更高的精度.最后给出了数值算例，验证了本方法的有效性和正确性. 关键词： 重构核粒子法 弹性力学 边界无单元法     

A comparison of differential and integral equations of radiative transfer

The equations of radiative transfer and of statistical equilibrium of a two-level atom are solved by means of differential and integral equations for a one-dimensional medium. The numerical solutions are compared to the analytic solution. It is found that the integral equation for piecewise quadratic source functions gives more accurate results than does the differential equation.     

New exact solutions to the space–time fractional nonlinear wave equation obtained by the ansatz and functional variable methods

In this paper, the ansatz method and the functional variable method are employed to find new analytic solutions for the space–time nonlinear fractional wave equation, the space–time fractional Kadomtsev–Petviashvili–Benjamin–Bona–Mahony equation and the space–time fractional modified Korteweg–de Vries–Zakharov–Kuznetsov equation. As a result, some exact solutions are obtained in terms of hyperbolic and periodic functions. It is shown that the proposed methods provide a more powerful mathematical tool for constructing exact solutions for many other nonlinear fractional differential equations occurring in nonlinear physical phenomena. We have also presented the numerical simulations for these equations by means of three dimensional plots.     

Analytical solution of transverse oscillation in cyclotron using LP method

We have carried out an approximate analytical solution to precisely consider the influence of magnetic field on the transverse oscillation of particles in a cyclotron.The differential equations of transverse oscillation are solved from the Lindstedt-Poincare method.After careful deduction,accurate first-order analytic solutions are obtained.The analytical solutions are applied to the magnetic field from an isochronous cyclotron with four spiral sectors.The accuracy of these analytical solutions is verified and confirmed from comparison with a numerical method.Finally,we discussed the transverse oscillation at v_0=N/2,using the same analytical solution.     

A second-order accurate numerical method for the two-dimensional fractional diffusion equation

Spatially fractional order diffusion equations are generalizations of classical diffusion equations which are used in modeling practical superdiffusive problems in fluid flow, finance and others. In this paper, we present an accurate and efficient numerical method to solve a fractional superdiffusive differential equation. This numerical method combines the alternating directions implicit (ADI) approach with a Crank–Nicolson discretization and a Richardson extrapolation to obtain an unconditionally stable second-order accurate finite difference method. The stability and the consistency of the method are established. Numerical solutions for an example super-diffusion equation with a known analytic solution are obtained and the behavior of the errors are analyzed to demonstrate the order of convergence of the method.     

A reliable treatment of the homotopy analysis method for viscous flow over a non-linearly stretching sheet in presence of a chemical reaction and under influence of a magnetic field

The similarity solution for the steady two-dimensional flow of an incompressible viscous and electrically conducting fluid over a non-linearly semi-infinite stretching sheet in the presence of a chemical reaction and under the influence of a magnetic field gives a system of non-linear ordinary differential equations. These non-linear differential equations are analytically solved by applying a newly developed method, namely the Homotopy Analysis Method (HAM). The analytic solutions of the system of non-linear differential equations are constructed in the series form. The convergence of the obtained series solutions is carefully analyzed. Graphical results are presented to investigate the influence of the Schmidt number, magnetic parameter and chemical reaction parameter on the velocity and concentration fields. It is noted that the behavior of the HAM solution for concentration profiles is in good agreement with the numerical solution given in reference [A. Raptis, C. Perdikis, Int. J. Nonlinear Mech. 41, 527 (2006)].      

Application of the Generalized Differential Quadrature Method in Solving Burgers'Equations

The aim of this paper is to obtain numerical solutions of the one-dimensional, two-dimensional and coupled Burgers' equations through the generalized differential quadrature method (GDQM). The polynomial-based differential quadrature (PDQ) method is employed and the obtained system of ordinary differential equations is solved via the total variation diminishing Runge-Kutta (TVD-RK) method. The numerical solutions are satisfactorily coincident with
the exact solutions. The method can compete against the methods applied in the literature.     

分数阶Oldroyd-B黏弹性Poiseuille流的Laplace数值反演分析

采用Laplace数值反演的Stehfest算法研究了分数阶Oldroyd-B粘弹性流体在两平板间非定常的Poiseuille流动问题.首先,通过数值解与近似解析解的比较验证了Stehfest算法的有效性.其次,运用Stehfest算法对平板Poiseuille流动进行了研究,揭示了分数阶黏弹性平板流的速度过冲和应力过冲现象,指出这些现象对分数导数的阶数存在明显的依赖性.同时,数值结果表明,整数阶本构方程仅仅是分数阶本构方程的特例,分数阶本构方程较整数阶本构方程具有更广泛的适用性。     

Analytic solutions of the optical bistability equations for a standing wave cavity

The problem of optical bistability in a standing wave cavity in the steady state leads to a pair of coupled, nonlinear, ordinary differential equations for the forward and backward waves. Only numerical solutions have so far been presented for these equations. We give their exact analytic solutions and find good agreement with the numerical results. The exact solutions are shown to reduce to the mean field equation for the input and output fields in the double limits T → 0 and αL → 0 for the mirror transmission and the linear absorption absorption, respectively.     

The finite element method for the coupled Schrödinger-KdV equations

In this Letter, we employ finite element method to study a periodic initial value problem for the coupled Schrödinger-KdV equations. For the case of one dimension, this problem is reduced to a system of ordinary differential equations by using a semi-discrete scheme. The conservation properties of this scheme, the existence and uniqueness of the discrete solutions, and error estimates are presented. In numerical experiments, the resulting system of ordinary differential equations are solved by Runge-Kutta method at each time level. The superior accuracy of this scheme is shown by comparing the numerical solutions with the exact solutions.     

Gauge invariance and formal integrability of the Yang-Mills-Higgs equations

In the framework of the formal theory of overdetermined systems of partial differential equations, it is shown that the Yang-Mills-Higgs equations are an involutive, and hence formally integrable, system. To this end a key role is played by the gauge invariance of the theory and the resulting differential identities involving the field equations themselves. By applying a theorem of Malgrange, an existence theorem for the solutions of the Yang-Mills-Higgs field equations in the analytic context is thus obtained. The approach is within differential geometry.     

Explicit Solutions for Generalized (2+1)-Dimensional NonlinearZakharov-Kuznetsov Equation

The exact solutions of the generalized (2+1)-dimensional nonlinearZakharov-Kuznetsov (Z-K) equation are explored by the method of the improved generalized auxiliary differential equation. Many explicit analytic solutions of the Z-K equation are obtained. The methods used to solve the Z-K equation can be employed in further work to establish new solutions for other nonlinear partial differential equations.     

求解对流扩散方程的一种有限分析方法

本文提出了一种求解线性与拟线性对流扩散方程的有限分析格式,证明了这种格式的解的存在唯一性及绝对稳定性与广义弱稳定性等,大量算例表明它具有较高的精度与很好的稳定性。     

Application of He's homotopy perturbation method to boundary layer flow and convection heat transfer over a flat plate

In this Letter, the problem of forced convection over a horizontal flat plate is presented and the homotopy perturbation method (HPM) is employed to compute an approximation to the solution of the system of nonlinear differential equations governing on the problem. It has been attempted to show the capabilities and wide-range applications of the homotopy perturbation method in comparison with the previous ones in solving heat transfer problems. The obtained solutions, in comparison with the exact solutions admit a remarkable accuracy. A clear conclusion can be drawn from the numerical results that the HPM provides highly accurate numerical solutions for nonlinear differential equations.     

Connections defining representations of zero curvature and the solitons of sine-Gordon and Korteweg-de Vries equations

It is claimed that solutions of travelling-wave type (and, in particular, soliton solutions) of partial differential equations can be created by using connections defining representations of zero curvature. In this paper, we construct solitons of the sine-Gordon and Korteweg-de Vries equations. By previous results of the author, the connections defining representations of zero curvature for a given differential equation generate Bäcklund transformations for this equation. It can be shown that the well-known Lax system (the so-called Lax pair) for the Korteweg-de Vries equation is a special case of a Bäcklund system (i.e., the system of partial differential equations defining a Bäcklund transformation). Note that the creation of solitons by means of the inverse scattering method is in fact a creation of solitons by means of the Lax system (without using connections defining the representations of zero curvature from the very beginning). Moreover, the inverse scattering method is essentially more labor-consuming than the method suggested in the present paper. Further, it is not required to involve any physical notions when using the suggested method. In the final section of the paper, we consider the so-called 2-soliton solutions of sine-Gordon and Korteweg — de Vries equations. Here we systematically use the invariant analytic method developed by G. F. Laptev, which is well-known in differential geometry under the title of Cartan-Laptev method.     

Numerical integration of an inhomogeneous boundary value problem

A numerical method is developped for the integration of stiff inhomogeneous coupled ordinary differential equations. It is shown how to compute integrals with integrands containing the solution of the differential equations. The method is stable and avoids storing of large amounts of intermediate results. The present method can also be applied to problems involving solutions of homogeneous coupled equations.     

Simplified analytic solutions and a novel fast algorithm for Yb-doped double-clad fiber lasers

In this paper, continuous wave Yb³⁺-doped double-clad fiber lasers (DCFLs) with linear-cavity are investigated theoretically and numerically using the rate equations. Under the steady state conditions, the simplified analytic solutions of Yb³⁺-doped DCFLs under considering the scattering loss are deduced in the strongly pump condition. Compared with the known analytic solutions in published literatures, our analytic solutions are more accurate, especially, at higher reflectivity of output mirror. In addition, a fast and stable algorithm based on the Newton-Raphson method is proposed to simulate numerically Yb³⁺-doped DCFLs. The results by simplified analytic solutions are in good agreement with those by the numerical simulation. Moreover, we have performed the optimization of an Yb³⁺-doped DCFL using the simplified analytic solutions and the numerical simulations, respectively.     

光学Maxwell-Bloch方程的数值算法研究及其应用

建立了高准确度快速求解均匀展宽二能级体系光学Maxwell-Bloch耦合方程的数值算法.通过与特定条件得到的解析解的比较,验证了算法所具有的高收敛性和稳定性,并可保持算法的误差阶数,因此算法是可靠并实用的.应用该算法数值求解了一般条件下的MB方程,并由计算结果分析了失谐量、弛豫时间、初始光强对光脉冲在介质中的传播及对Bloch矢量演化的影响.所建立的数值算法对MB方程以及修正的这类偏微分方程组具有普适性.     

Variational iteration method for solving compressible Euler equations

This paper applies the variational iteration method to obtain approximate analytic solutions of compressible Euler equations in gas dynamics.This method is based on the use of Lagrange multiplier for identification of optimal values of parameters in a functional.Using this method,a rapid convergent sequence is produced which converges to the exact solutions of the problem.Numerical results and comparison with other two numerical solutions verify that this method is very convenient and efficient.     

二阶常微分方程组边值问题的一种数值解法

本文用压缩变挟方法,求解了函数解随自变量增加而迅速增长的一类常微分方程的边值问题。结果表明,这种方法,可加快收敛速度。在奇点W=0附近,采用局部解析近似处理,得到了比较满意的数值结果。