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.     

Element-free Galerkin (EFG) method for analysis of the time-fractional partial differential equations

The present paper deals with the numerical solution of time-fractional partial differential equations using the element-free Galerkin (EFG) method, which is based on the moving least-square approximation. Compared with numerical methods based on meshes, the EFG method for time-fractional partial differential equations needs only scattered nodes instead of meshing the domain of the problem. It neither requires element connectivity nor suffers much degradation in accuracy when nodal arrangements are very irregular. In this method, the first-order time derivative is replaced by the Caputo fractional derivative of order α (0<α ≤1). The Galerkin weak form is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. Several numerical examples are presented and the results we obtained are in good agreement with the exact solutions.     

Differential Quadrature Analysis of Moving Load Problems

The differential quadrature method (DQM) has been successfully used in a variety of fields. Similar to the conventional point discrete methods such as the collocation method and finite difference method, however, the DQM has some difficulty in dealing with singular functions like the Dirac-delta function. In this paper, two modifications are introduced to overcome the difficulty encountered in solving differential equations with Dirac-delta functions by using the DQM. The moving point load is work-equivalent to loads applied at all grid points and the governing equation is numerically integrated before it is discretized in terms of the differential quadrature. With these modifications, static behavior and forced vibration of beams under a stationary or a moving point load are successfully analyzed by directly using the DQM. It is demonstrated that the modified DQM can yield very accurate solutions. The compactness and computational efficiency of the DQM are retained in solving the partial differential equations with a time dependent Dirac-delta function.     

Application of homotopy perturbation method to the RLW and generalized modified Boussinesq equations

In this Letter, He's homotopy perturbation method (HPM) is implemented for finding the solitary-wave solutions of the regularized long-wave (RLW) and generalized modified Boussinesq (GMB) equations. We obtain numerical solutions of these equations for the initial conditions. We will show that the convergence of the HPM is faster than those obtained by the Adomian decomposition method (ADM). 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.     

The Method of Fundamental Solutions for Solving Convection-Diffusion Equations with Variable Coefficients

A meshless method based on the method of fundamental solutions (MFS) is proposed to solve the time-dependent partial differential equations with variable coefficients. The proposed method combines the time discretization and the one-stage MFS for spatial discretization. In contrast to the traditional two-stage process, the one-stage MFS approach is capable of solving a broad spectrum of partial differential equations. The numerical implementation is simple since both closed-form approximate particular solution and fundamental solution are easier to find than the traditional approach. The numerical results show that the one-stage approach is robust and stable.     

Simulating non-linear coupled oscillators by an iterative differential quadrature method

An iterative method based on differential quadrature rules is proposed as a new unified frame of resolution for non-linear two-degree-of-freedom systems. Dynamical systems with Duffing-type non-linearity have been considered. Differential quadrature rules have been applied with a careful distribution of sampling points to reduce the governing equation of motion to two second-order non-linear, non-autonomous ordinary differential equations and to solve the time-domain problem. The time domain of the problem is discretized by means of time intervals, with the same distribution of sampling points used to discretize the space domain (which can be seen as a single interval). It will be shown that accurate solutions depend not only on the choice of the distribution of sampling points, but also on the length of the time interval one refers to in the computations. The numerical results, utilized to draw Poincaré maps, are successfully compared with those obtained using the Runge-Kutta method.     

结构动力分析的时域配点DQ半解析法

直接从控制微分方程出发,以梁为对象提出了本计算方法。该方法在空间域采用DQ法(DiferentialQuadratureMethod),在时间域取级数,采用时域配点的方式得到求位移场全部待定参数的可解方程组,求解一次线性方程组即可求得全域的响应位移场。     

Conditional quadrature method of moments for kinetic equations

Kinetic equations arise in a wide variety of physical systems and efficient numerical methods are needed for their solution. Moment methods are an important class of approximate models derived from kinetic equations, but require closure to truncate the moment set. In quadrature-based moment methods (QBMM), closure is achieved by inverting a finite set of moments to reconstruct a point distribution from which all unclosed moments (e.g. spatial fluxes) can be related to the finite moment set. In this work, a novel moment-inversion algorithm, based on 1-D adaptive quadrature of conditional velocity moments, is introduced and shown to always yield realizable distribution functions (i.e. non-negative quadrature weights). This conditional quadrature method of moments (CQMOM) can be used to compute exact N-point quadratures for multi-valued solutions (also known as the multi-variate truncated moment problem), and provides optimal approximations of continuous distributions. In order to control numerical errors arising in volume averaging and spatial transport, an adaptive 1-D quadrature algorithm is formulated for use with CQMOM. The use of adaptive CQMOM in the context of QBMM for the solution of kinetic equations is illustrated by applying it to problems involving particle trajectory crossing (i.e. collision-less systems), elastic and inelastic particle–particle collisions, and external forces (i.e. fluid drag).     

Numerical Simulation of MHD Peristaltic Flow with Variable Electrical Conductivity and Joule Dissipation Using Generalized Differential Quadrature Method

In this paper, the MHD peristaltic flow inside wavy walls of an asymmetric channel is investigated, where the walls of the channel are moving with peristaltic wave velocity along the channel length. During this investigation,the electrical conductivity both in Lorentz force and Joule heating is taken to be temperature dependent. Also, the long wavelength and low Reynolds number assumptions are utilized to reduce the governing partial differential equations into a set of coupled nonlinear ordinary differential equations. The new set of obtained equations is then numerically solved using the generalized differential quadrature method(GDQM). This is the first attempt to solve the nonlinear equations arising in the peristaltic flows using this method in combination with the Newton-Raphson technique. Moreover, in order to check the accuracy of the proposed numerical method, our results are compared with the results of built-in Mathematica command NDSolve. Taking Joule heating and viscous dissipation into account, the effects of various parameters appearing in the problem are used to discuss the fluid flow characteristics and heat transfer in the electrically conducting fluids graphically. In presence of variable electrical conductivity, velocity and temperature profiles are highly decreasing in nature when the intensity of the electrical conductivity parameter is strengthened.     

Direct quadrature method of moments solution of the Fokker–Planck equation

The direct quadrature method of moments is presented as an efficient and accurate means of numerically computing solutions of the Fokker–Planck equation corresponding to stochastic nonlinear dynamical systems. The theoretical details of the solution procedure are first presented. The method is then used to solve Fokker–Planck equations for both 1D and 2D (noisy van der Pol oscillator) processes which possess nonlinear stochastic differential equations. Higher-order moments of the stationary solutions are computed and prove to be very accurate when compared to analytic (1D process) and Monte Carlo (2D process) solutions.     

Numerical Solutions of Fractional Boussinesq Equation

Based upon the Adomian decomposition method, a scheme is developed to obtain numerical solutions of a fractional Boussinesq equation with initial condition, which is introduced by replacing some order time and space derivatives by fractional derivatives. The fractional derivatives are described in the Caputo sense. So the traditional Adomian decomposition method for differential equations of integer order is directly extended to derive explicit and numerical solutions of the fractional differential equations. The solutions of our model equation are calculated in the form of convergent series with easily computable components.     

Stability Analysis and Order Improvement for Time Domain Differential Quadrature Method

The differential quadrature method has been widely used in scientific and engineering computation. However, for the basic characteristics of time domain differential quadrature method, such as numerical stability and calculation accuracy or order, it is still lack of systematic analysis conclusions. In this paper, according to the principle of differential quadrature method, it has been derived and proved that the weighting coefficients matrix of differential quadrature method meets the importantV-transformation feature. Through the equivalence of the differential quadrature method and the implicit Runge-Kutta method, it has been proved that the differential quadrature method is A-stable and $s$-stage $s$-order method. On this basis, in order to further improve the accuracy of the time domain differential quadrature method, a class of improved differential quadrature method of $s$-stage $2s$-order has been proposed by using undetermined coefficients method and Padéapproximations. The numerical results show that the proposed differential quadrature method is more precise than the traditional differential quadrature method.     

Navier-Stokes方程的线性化微分求积法

提出一种新的线性化微分求积法(LDQM),将这种目的应用到流函数和涡量形式的Navier-Stokes方程.通过LDQM,非线性方程很容易被解出来,并且容易处理压力的边界条件.为检验本目的,计算了两个数值算例.     

An element-free Galerkin (EFG) method for generalized Fisher equations (GFE)

A generalized Fisher equation(GFE) relates the time derivative of the average of the intrinsic rate of growth to its variance.The exact mathematical result of the GFE has been widely used in population dynamics and genetics,where it originated.Many researchers have studied the numerical solutions of the GFE,up to now.In this paper,we introduce an element-free Galerkin(EFG) method based on the moving least-square approximation to approximate positive solutions of the GFE from population dynamics.Compared with other numerical methods,the EFG method for the GFE needs only scattered nodes instead of meshing the domain of the problem.The Galerkin weak form is used to obtain the discrete equations,and the essential boundary conditions are enforced by the penalty method.In comparison with the traditional method,numerical solutions show that the new method has higher accuracy and better convergence.Several numerical examples are presented to demonstrate the effectiveness of the method.     

The choosing of reproducing kernel particle shape function with mathematic proof

Many mechanical problems can be induced from differential equations with boundary conditions; there exist analytic and numerical methods for solving the differential equations. Usually it is not so easy to obtain analytic solutions. So it is necessary to give numerical solutions. The reproducing kernel particle (RKP) method is based on the Garlerkin Meshless method. According to the Sobolev space and Fourier transform, the RKP shape function is mathematically proved in this paper.     

Three-dimensional free vibration of thick functionally graded annular plates in thermal environment

The free vibration analysis of functionally graded (FG) thick annular plates subjected to thermal environment is studied based on the 3D elasticity theory. The material properties are assumed to be temperature dependent and graded in the thickness direction. Considering the thermal environment effects and using Hamilton's principle, the equations of motion are derived. The effects of the initial thermal stresses are considered accurately by obtaining them from the 3D thermoelastic equilibrium equations. The differential quadrature method (DQM) as an efficient and accurate numerical tool is used to solve both the thermoelastic equilibrium and free vibration equations. Very fast rate of convergence of the method is demonstrated. Also, the formulation is validated by comparing the results with those obtained based on the first-order shear deformation theory and also with those available in the literature for the limit cases, i.e. annular plates without thermal effects. The effects of temperature rise, material and geometrical parameters on the natural frequencies are investigated. The new results can be used as benchmark solutions for future researches.     

Collisions and umbilic catastrophes

Two procedures to evaluate the hyperbolic umbilic canonical diffraction integral, a function that occurs in many applications of semiclassical collision theory and optics, are described. The first is the integration of the various differential equations satisfied by this function. We present schemes to perform the numerical integration and investigate their range of validity. A second algorithm, which evaluates the integral directly by quadrature, is given. Here, attention is focused on the very rapidly oscillating regions of the integrand. The direct quadrature method is shown to be more widely applicable and better suited to numerical computation.     

干斜压大气拉格朗日原始方程组的半解析解法和非线性密度流数值试验

以差分方程代替微分方程给大气原始方程组求解带来了诸多难以解决的问题, 对于(半)拉格朗日模式来说质点轨迹的计算与Helmholtz方程的求解是两大难题. 本文通过对气压变量代换, 并在积分时间步长内将原始方程组线性化, 近似为常微分方程组, 求出方程组的半解析解, 再采用精细积分法求解半解析解. 半解析方法可同时计算风、气压和位移, 无需求解Helmholtz方程, 质点的位移采用积分风的半解析解得到, 相比采用风速外推的计算方法, 半解析方法更科学合理. 非线性密度流试验检验表明: 半解析模式能够清晰地模拟Kelvin-Helmholtz 切变不稳定涡旋的发生和发展过程; 模拟的气压场和风场环流结构与标准解非常相似, 且数值解是收敛的, 同时, 总质量和总能量具有较好的守恒性. 试验初步证明了采用半解析方法求解大气原始方程组是可行的, 为大气数值模式的构建提供了一个新的思路.     

New Complexiton Solutions of (2+1)-Dimensional Nizhnik-Novikov-Veselov Equations

In this paper, we extend the multiple Riccati equations rational expansion method by introducing a new ansatz. Using this method, many complexiton solutions of the (2 1)-dimensional Nizhnik-Novikov-Veselov equations are obtained which include various combination of hyperbolic and trigonometric periodic function solutions, various combination of hyperbolic and rational function solutions, various combination of trigonometric periodic and rational function solutions, etc. The method can be also used to solve other nonlinear partial differential equations.     

New explicit and exact travelling wave solutions for a system of dispersive long wave equations

A method to construct the new exact solutions of nonlinear partial differential equations (NLPDEs) in a unified way is presented, which is named an improved sine-cosine method. This method is more powerful than the sine-cosine method. Systems of dispersive long wave equations in (1+1) and (2+1) dimensions are chosen to illustrate the method and several types of explicit and exact travelling wave solutions are obtained. These solutions contain Wang's results and other types of solitary wave solutions and new solutions. The method presented here is general and can also be applied to solve more systems of nonlinear partial differential equations, such as the coupled KdV equations.