Mixture of a New Integral Transform and Homotopy Perturbation Method for Solving Nonlinear Partial Differential Equations

In this paper, we present a new method, a mixture of homotopy perturbation method and a new integral transform to solve some nonlinear partial differential equations. The proposed method introduces also He’s polynomials [1]. The analytical results of examples are calculated in terms of convergent series with easily computed components [2].

Exact Solutions of Two Nonlinear Partial Differential Equations by the First Integral Method

In recent years, many methods have been used to find the exact solutions of nonlinear partial differential equations. One of them is called the first integral method, which is based on the ring theory of commutative algebra. In this paper, exact travelling wave solutions of the Non-Boussinesq wavepacket model and the (2 + 1)-dimensional Zoomeron equation are studied by using the first integral method. From the solving process and results, the first integral method has the characteristics of simplicity, directness and effectiveness about solving the exact travelling wave solutions of nonlinear partial differential equations. In other words, tedious calculations can be avoided by Maple software;the solutions of more accurate and richer travelling wave solutions are obtained. Therefore, this method is an effective method for solving exact solutions of nonlinear partial differential equations.

Dynamic stability analysis of porous functionally graded beams under hygro-thermal loading using nonlocal strain gradient integral model

We present a study on the dynamic stability of porous functionally graded(PFG)beams under hygro-thermal loading.The variations of the properties of the beams across the beam thicknesses are described by the power-law model.Unlike most studies on this topic,we consider both the bending deformation of the beams and the hygro-thermal load as size-dependent,simultaneously,by adopting the equivalent differential forms of the well-posed nonlocal strain gradient integral theory(NSGIT)which are strictly equipped with a set of constitutive boundary conditions(CBCs),and through which both the stiffness-hardening and stiffness-softening effects of the structures can be observed with the length-scale parameters changed.All the variables presented in the differential problem formulation are discretized.The numerical solution of the dynamic instability region(DIR)of various bounded beams is then developed via the generalized differential quadrature method(GDQM).After verifying the present formulation and results,we examine the effects of different parameters such as the nonlocal/gradient length-scale parameters,the static force factor,the functionally graded(FG)parameter,and the porosity parameter on the DIR.Furthermore,the influence of considering the size-dependent hygro-thermal load is also presented.

An Integral Collocation Approach Based on Legendre Polynomials for Solving Riccati, Logistic and Delay Differential Equations

In this paper, we propose and analyze some schemes of the integral collocation formulation based on Legendre polynomials. We implement these formulae to solve numerically Riccati, Logistic and delay differential equations with variable coefficients. The properties of the Legendre polynomials are used to reduce the proposed problems to the solution of non-linear system of algebraic equations using Newton iteration method. We give numerical results to satisfy the accuracy and the applicability of the proposed schemes.

Euler’s First-Order Explicit Method–Peridynamic Differential Operator for Solving Population Balance Equations of the Crystallization Process

Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.

Well-posedness of two-phase local/nonlocal integral polar models for consistent axisymmetric bending of circular microplates

Previous studies have shown that Eringen’s differential nonlocal model would lead to the ill-posed mathematical formulation for axisymmetric bending of circular microplates.Based on the nonlocal integral models along the radial and circumferential directions,we propose nonlocal integral polar models in this work.The proposed strainand stress-driven two-phase nonlocal integral polar models are applied to model the axisymmetric bending of circular microplates.The governing differential equations and boundary conditions(BCs)as well as constitutive constraints are deduced.It is found that the purely strain-driven nonlocal integral polar model turns to a traditional nonlocal differential polar model if the constitutive constraints are neglected.Meanwhile,the purely strain-and stress-driven nonlocal integral polar models are ill-posed,because the total number of the differential orders of the governing equations is less than that of the BCs plus constitutive constraints.Several nominal variables are introduced to simplify the mathematical expression,and the general differential quadrature method(GDQM)is applied to obtain the numerical solutions.The results from the current models(CMs)are compared with the data in the literature.It is clearly established that the consistent softening and toughening effects can be obtained for the strain-and stress-driven local/nonlocal integral polar models,respectively.The proposed two-phase local/nonlocal integral polar models(TPNIPMs)may provide an efficient method to design and optimize the plate-like structures for microelectro-mechanical systems.

On the consistency of two-phase local/nonlocal piezoelectric integral model

In this paper,we propose general strain-and stress-driven two-phase local/nonlocal piezoelectric integral models,which can distinguish the difference of nonlocal effects on the elastic and piezoelectric behaviors of nanostructures.The nonlocal piezoelectric model is transformed from integral to an equivalent differential form with four constitutive boundary conditions due to the difficulty in solving intergro-differential equations directly.The nonlocal piezoelectric integral models are used to model the static bending of the Euler-Bernoulli piezoelectric beam on the assumption that the nonlocal elastic and piezoelectric parameters are coincident with each other.The governing differential equations as well as constitutive and standard boundary conditions are deduced.It is found that purely strain-and stress-driven nonlocal piezoelectric integral models are ill-posed,because the total number of differential orders for governing equations is less than that of boundary conditions.Meanwhile,the traditional nonlocal piezoelectric differential model would lead to inconsistent bending response for Euler-Bernoulli piezoelectric beam under different boundary and loading conditions.Several nominal variables are introduced to normalize the governing equations and boundary conditions,and the general differential quadrature method(GDQM)is used to obtain the numerical solutions.The results from current models are validated against results in the literature.It is clearly established that a consistent softening and toughening effects can be obtained for static bending of the Euler-Bernoulli beam based on the general strain-and stress-driven local/nonlocal piezoelectric integral models,respectively.

ITERATIVE ALGORITHM FOR AXIALLY ACCELERATING STRINGS WITH INTEGRAL CONSTITUTIVE LAW

A numerical method is proposed to simulate the transverse vibrations of a viscoelastic moving string constituted by an integral law. In the numerical computation, the Galerkin method based on the Hermite functions is applied to discretize the state variables, and the Runge- Kutta method is applied to solve the resulting differential-integral equation system. A linear iterative process is designed to compute the integral terms at each time step, which makes the numerical method more efficient and accurate. As examples, nonlinear parametric vibrations of an axially moving viscoelastic string are analyzed.

CERTAIN FAMILY OF INTEGRAL OPERATORS PRESERVING SUBORDINATION AND SUPERORDINATION

We obtain subordination, superordination and sandwich-preserving new theorems for certain integral operators defined on the space of normalized analytic functions in the open unit disk. The sandwich-type theorem for these integral operators is also derived, and the results generalize some recently ones.

A Contour Integral Method for Linear Differential Equations in Complex Plane

This paper presents the contour integral method for solving the linear constant coefficient ordinary differential equations in complex plane,and obtains the uniform expressions of the general solutions.Firstly,by using Residue Theorem,the general form of the contour integral representation for the homogeneous complex differential equation is obtained,which can be degenerated to classical results in real line.As for inhomogeneous complex differential equations with constant coefficients,we construct the integral expression of the particular solution for any continuous forcing term,and give rigorous proof via Residue Theorem.Thus the general solutions of inhomogeneous complex differential equations are also given.The main purpose of this paper is to give a foundation for a complete theory of linear complex differential equations with constant coefficients by a contour integral method.The results can not only solve the inhomogeneous complex differential equation well,but also explain the forms that are difficult to be understood in the classical solutions.

Finite Difference-Peridynamic Differential Operator for Solving Transient Heat Conduction Problems

Transient heat conduction problems widely exist in engineering.In previous work on the peridynamic differential operator(PDDO)method for solving such problems,both time and spatial derivatives were discretized using the PDDO method,resulting in increased complexity and programming difficulty.In this work,the forward difference formula,the backward difference formula,and the centered difference formula are used to discretize the time derivative,while the PDDO method is used to discretize the spatial derivative.Three new schemes for solving transient heat conduction equations have been developed,namely,the forward-in-time and PDDO in space(FT-PDDO)scheme,the backward-in-time and PDDO in space(BT-PDDO)scheme,and the central-in-time and PDDO in space(CT-PDDO)scheme.The stability and convergence of these schemes are analyzed using the Fourier method and Taylor’s theorem.Results show that the FT-PDDO scheme is conditionally stable,whereas the BT-PDDO and CT-PDDO schemes are unconditionally stable.The stability conditions for the FT-PDDO scheme are less stringent than those of the explicit finite element method and explicit finite difference method.The convergence rate in space for these three methods is two.These constructed schemes are applied to solve one-dimensional and two-dimensional transient heat conduction problems.The accuracy and validity of the schemes are verified by comparison with analytical solutions.

Weighted variation inequalities for differential operators and singular integrals in higher dimensions

We prove weighted q-variation inequalities with 2<q<∞for sharp truncations of singular integral operators in higher dimensions.The vector-valued extensions of these inequalities are also given.Parallel results are proven for differential operators.

A New Representation for Second Order Stochastic Integral-differential Operators and Its Applications

In this paper, we＇ll prove new representation theorems for a kind of second order stochastic integral- differential operator by stochastic differential equations （SDEs） and backward stochastic differential equations （BSDEs） with jumps, and give some applications.

LOCAL GAUSSIAN-COLLOCATION SCHEME TO APPROXIMATE THE SOLUTION OF NONLINEAR FRACTIONAL DIFFERENTIAL EQUATIONS USING VOLTERRA INTEGRAL EQUATIONS

This work describes an accurate and effective method for numerically solving a class of nonlinear fractional differential equations.To start the method,we equivalently convert these types of differential equations to nonlinear fractional Volterra integral equations of the second kind by integrating from both sides of them.Afterward,the solution of the mentioned Volterra integral equations can be estimated using the collocation method based on locally supported Gaussian functions.The local Gaussian-collocation scheme estimates the unknown function utilizing a small set of data instead of all points in the solution domain,so the proposed method uses much less computer memory and volume computing in comparison with global cases.We apply the composite non-uniform Gauss-Legendre quadrature formula to estimate singular-fractional integrals in the method.Because of the fact that the proposed scheme requires no cell structures on the domain,it is a meshless method.Furthermore,we obtain the error analysis of the proposed method and demon-strate that the convergence rate of the approach is arbitrarily high.Illustrative examples clearly show the reliability and efficiency of the new technique and confirm the theoretical error estimates.

A Kernel-Free Boundary Integral Method for Variable Coefficients Elliptic PDEs

This work proposes a generalized boundary integral method for variable coefficients elliptic partial differential equations(PDEs),including both boundary value and interface problems.The method is kernel-free in the sense that there is no need to know analytical expressions for kernels of the boundary and volume integrals in the solution of boundary integral equations.Evaluation of a boundary or volume integral is replaced with interpolation of a Cartesian grid based solution,which satisfies an equivalent discrete interface problem,while the interface problem is solved by a fast solver in the Cartesian grid.The computational work involved with the generalized boundary integral method is essentially linearly proportional to the number of grid nodes in the domain.This paper gives implementation details for a secondorder version of the kernel-free boundary integral method in two space dimensions and presents numerical experiments to demonstrate the efficiency and accuracy of the method for both boundary value and interface problems.The interface problems demonstrated include those with piecewise constant and large-ratio coefficients and the heterogeneous interface problem,where the elliptic PDEs on two sides of the interface are of different types.

DIFFERENTIAL GEOMETRICAL METHOD IN ELASTICCOMPOSITE WITH IMPERFECT INTERFACES

A differential geometrical method is for the first time used to calculate the effective moduli of a two-phase elastic composite materials with imperfect interface which the inclusions are assumed to be ellipsoidal of revolutions. All of the interface integral items participating in forming the potential and complementary energy functionals of the composite materials are expressed in terms of intrinsic quantities of the ellipsoidal of revolutions. Based on this, the upper and the lower bound for the effective elastic moduli of the composite materials with inclusions described above have been derived. Under three limiting conditions of sphere, disk and needle shaped inclusions, the results of this paper will return to the bounds obtained by Hashin([6]) (1992).

Hamilton-Jacobi method for solving ordinary differential equations

The Hamilton-Jacobi method for solving ordinary differential equations is presented in this paper. A system of ordinary differential equations of first order or second order can be expressed as a Hamilton system under certain conditions. Then the Hamilton-Jacobi method is used in the integration of the Hamilton system and the solution of the original ordinary differential equations can be found. Finally, an example is given to illustrate the application of the result.

Exponential Stability of Traveling Pulse Solutions of a Singularly Perturbed System of Integral Differential Equations Arising From Excitatory Neuronal Networks

We establish the exponential stability of fast traveling pulse solutions to nonlinear singularly per-turbed systems of integral differential equations arising from neuronal networks.It has been proved that expo-nential stability of these orbits is equivalent to linear stability.Let (?) be the linear differential operator obtainedby linearizing the nonlinear system about its fast pulse,and let σ((?)) be the spectrum of (?).The linearizedstability criterion says that if max{Reλ:λ∈σ((?)),λ≠0}(?)-D,for some positive constant D,and λ=0 is asimple eigenvalue of (?)(ε),then the stability follows immediately (see [13] and [37]).Therefore,to establish theexponential stability of the fast pulse,it suffices to investigate the spectrum of the operator (?).It is relativelyeasy to find the continuous spectrum,but it is very difficult to find the isolated spectrum.The real part ofthe continuous spectrum has a uniformly negative upper bound,hence it causes no threat to the stability.Itremains to see if the isolated spectrum is safe.Eigenvalue functions (see [14] and [35,36]) have been a powerful tool to study the isolated spectrum of the as-sociated linear differential operators because the zeros of the eigenvalue functions coincide with the eigenvaluesof the operators.There have been some known methods to define eigenvalue functions for nonlinear systems ofreaction diffusion equations and for nonlinear dispersive wave equations.But for integral differential equations,we have to use different ideas to construct eigenvalue functions.We will use the method of variation of param-eters to construct the eigenvalue functions in the complex plane C.By analyzing the eigenvalue functions,wefind that there are no nonzero eigenvalues of (?) in {λ∈C:Reλ(?)-D} for the fast traveling pulse.Moreoverλ=0 is simple.This implies that the exponential stability of the fast orbits is true.

Improved precise integration method for differential Riccati equation

An improved precise integration method （IPIM） for solving the differential Riccati equation （DRE） is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise integration method （PIM） for solving the DRE is connected with the scaling and squaring method for computing the exponential of a matrix. The error analysis of the scaling and squaring method for the exponential of a matrix is applied to the PIM of the DRE. Based ,on the error analysis, the criterion for choosing two parameters of the PIM is given. Three kinds of IPIMs for solving the DRE are proposed. The numerical examples machine accuracy solutions. show that the IPIM is stable and gives the

Hermite Matrix Polynomial Collocation Method for Linear Complex Differential Equations and Some Comparisons

In this paper, we introduce a Hermite operational matrix collocation method for solving higher-order linear complex differential equations in rectangular or elliptic domains. We show that based on a linear algebra theorem, the use of different polynomials such as Hermite, Bessel and Taylor in polynomial collocation methods for solving differential equations leads to an equal solution, and the difference in the numerical results arises from the difference in the coefficient matrix of final linear systems of equations. Some numerical examples will also be given.