Ultraconvergence for averaging discontinuous finite elements and its applications in Hamiltonian system

This paper discusses the k-degree averaging discontinuous finite element solution for the initial value problem of ordinary differential equations.When k is even,the averaging numerical flux (the average of left and right limits for the discontinuous finite element at nodes) has the optimal-order ultraconvergence 2k + 2.For nonlinear Hamiltonian systems (e.g.,Schro¨dinger equation and Kepler system) with momentum conservation,the discontinuous finite element methods preserve momentum at nodes.These properti...     

Continuous finite element methods for Hamiltonian systems

By applying the continuous finite element methods of ordinary differential equations,the linear element methods are proved having second-order pseudo-symplectic scheme and the quadratic element methods are proved having third-order pseudo- symplectic scheme respectively for general Hamiltonian systems,and they both keep energy conservative.The finite element methods are proved to be symplectic as well as energy conservative for linear Hamiltonian systems.The numerical results are in agree- ment with theory.     

PERTURBATION FINITE ELEMENT METHOD FOR SOLVING GEOMETRICALLY NONLINEAR PROBLEMS OF AXISYMMETRICAL SHELL

In analysing the geometrically nonlinear problem of an axisymmetrical thin-walled shell, the paper combines the perturbation method with the finite element method by introducing the former into the variational equation to obtain a series of linear equations of different orders and then solving the equations with the latter. It is well-known that the finite element method can be used to deal with difficult problems as in the case of structures with complicated shapes or boundary conditions, and the perturbation method can change the nonlinear problems into linear ones. Evidently the combination of the two methods will give an efficient solution to many difficult nonlinear problems and clear away some obstacles resulted from using any of the two methods solely. The paper derives all the formulas concerning an axisym-metric shell of large deformation by means of the perturbation finite element method and gives two numerical examples,the results of which show good convergence characteristics.     

An investigation of an Emden-Fowler equation from thin film flow

A third-order ordinary differential equation (ODE) for thin film flow with both Neumann and Dirichlet boundary conditions is transformed into a second-order nonlinear ODE with Dirichlet boundary conditions.Numerical solutions of the nonlinear second-order ODE are investigated using finite difference schemes.A finite difference formulation to an Emden-Fowler representation of the second-order nonlinear ODE is shown to converge faster than a finite difference formulation of the standard form of the second-order nonlinear ODE.Both finite difference schemes satisfy the von Neumann stability criteria.When mapping the numerical solution of the second-order ODE back to the variables of the original third-order ODE we recover the position of the contact line.A nonlinear relationship between the position of the contact line and physical parameters is obtained.     

Mixed time discontinuous space-time finite element method for convection diffusion equations

A mixed time discontinuous space-time finite element scheme for secondorder convection diffusion problems is constructed and analyzed. Order of the equation is lowered by the mixed finite element method. The low order equation is discretized with a space-time finite element method, continuous in space but discontinuous in time. Stability, existence, uniqueness and convergence of the approximate solutions are proved. Numerical results are presented to illustrate efficiency of the proposed method.     

Numerical solution of nonlinear ordinary differential equation for a turning point problem

By using the method in[3],several useful estimations of the derivatives of the solutionof the boundary value problem for a nonlinear ordinary differential equation with a turningpoint are obtained.With the help of the technique in[4],the uniform convergence on thesmall parameterεfor a difference scheme is proved.At the end of this paper,a numericalexample is given.The numerical result coincides with theoretical analysis.     

TRANSIENT SPHERICAL FLOW OF NON-NEWTONIAN POWER-LAW FLUIDS IN POROUS MEDIA

The transient spherical flow behavior of a slightly compressible, non-Newtonian, power-law fluids in porous media is studied. A nonlinear partial differential equation of parabolic type is derived. The diffusivity equation for spherical flow is a special case of the new equation. We obtain analytical, asymptotic and approximate solutions by using the methods of Laplace transform and weighted mass conservation. The structures of asymptotic and approximate solutions are similar, which enriches the theory of one-dimensional flow of non-Newtonian fluids through porous media.     

A HIGH ACCURACY DIFFERENCE SCHEME FOR THE SINGULAR PERTURBATION PROBLEM OF THE SECOND-ORDER LINEAR ORDINARY DIFFERENTIAL EQUATION IN CONSERVATION FORM

In this paper,combining the idea of difference method and finite element method,weconstruct a difference scheme for a self-adjoint problem in conservation form.Its solutionuniformly converges to that of the original differential equation problem with order h~3.     

Global attractor for Klein-Gordon-Schrodinger lattice system

We considered the longtime behavior of solutions of a coupled lattice dynamical system of Klein-Gordon-Schrodinger equation (KGS lattice system). We first proved the existence of a global attractor for the system considered here by introducing an equivalent norm and using "End Tails" of solutions. Then we estimated the upper bound of the Kolmogorov delta-entropy of the global attractor by applying element decomposition and the covering property of a polyhedron by balls of radii delta in the finite dimensional space. Finally, we presented an approximation to the global attractor by the global attractors of finite-dimensional ordinary differential systems.     

Global attractor for Klein-Gordon-Schroedinger lattice system

We considered the longtime behavior of solutions of a coupled lattice dynamical system of Klein-Gordon-Schroedinger equation （KGS lattice system）. We first proved the existence of a global attractor for the system considered here by introducing an equivalent norm and using ＂End Tails＂ of solutions. Then we estimated the upper bound of the Kolmogorov delta-entropy of the global attractor by applying element decomposition and the covering property of a polyhedron by balls of radii delta in the finite dimensional space. Finally, we presented an approximation to the global attractor by the global attractors of finite-dimensional ordinary differential systems.     

THE FINITE ELEMENT METHOD OF SINGULAR PERTURBATION PROBLEM

In this paper we construct new finite element subspace using polynomials of different degrees and the new finite element scheme is established. The convergence of the scheme and the stability of the reduced difference equation are proved.     

Peristaltic motion of third grade fluid in curved channel

＂ Analysis is performed to study the slip effects on the peristaltic flow of non-Newtonian fluid in a curved channel with wall properties. The resulting nonlinear partial differential equations are transformed to a single ordinary differential equation in a stream function by using the assumptions of long wavelength and low Reynolds number. This differential equation is solved numerically by employing the built-in routine for solving nonlinear boundary value problems （BVPs） through the software Mathematica. In addition, the analytic solutions for small Deborah number are computed with a regular perturbation technique. It is noticed that the symmetry of bolus is destroyed in a curved channel. An intensification in the slip effect results in a larger magnitude of axial velocity. Further, the size and circulation of the trapped boluses increase with an increase in the slip parameter. Different from the case of planar channel, the axial velocity profiles are tilted towards the lower part of the channel. A comparative study between analytic and numerical solutions shows excellent agreement.     

FINITE ELEMENT ANALYSIS OF NON-NEWTONIAN FLUID FLOW IN 2-D BRANCHING CHANNEL

This paper presents finite element analysis of non-Newtonian fluid flow in 2-dbranching channel.The Galerkin method and mixed finite element method are usedHere the fluid is considered as incompressible,non-Newtonian fluid with Oldyorddifferential-type constitutive equation.The non-linear algebraic equation system whichis formulated with finite element method is solved by means of continuous differentialmethod The results show that finite element method is suitable for the analysis of non-Newtonian fluid flow with complex geometry.     

Symmetry solutions of a nonlinear elastic wave equation with third-order anharmonic corrections

Lie symmetry method is applied to analyze a nonlinear elastic wave equation for longitudinal deformations with third-order anharmonic corrections to the elastic energy. Symmetry algebra is found and reductions to second-order ordinary differential equations （ODEs） are obtained through invariance under different symmetries. The reduced ODEs are further analyzed to obtain several exact solutions in an explicit form. It was observed in the literature that anharmonic corrections generally lead to solutions with time-dependent singularities in finite times singularities, we also obtain solutions which Along with solutions with time-dependent do not exhibit time-dependent singularities.     

Self-adaptive one-dimensional nonlinear finite element method based on element energy projection method

The element energy projection （EEP） method for computation of super- convergent resulting in a one-dimensional finite element method （FEM） is successfully used to self-adaptive FEM analysis of various linear problems, based on which this paper presents a substantial extension of the whole set of technology to nonlinear problems. The main idea behind the technology transfer from linear analysis to nonlinear analysis is to use Newton＇s method to linearize nonlinear problems into a series of linear problems so that the EEP formulation and the corresponding adaptive strategy can be directly used without the need for specific super-convergence formulation for nonlinear FEM. As a re- sult, a unified and general self-adaptive algorithm for nonlinear FEM analysis is formed. The proposed algorithm is found to be able to produce satisfactory finite element results with accuracy satisfying the user-preset error tolerances by maximum norm anywhere on the mesh. Taking the nonlinear ordinary differential equation （ODE） of second-order as the model problem, this paper describes the related fundamental idea, the imple- mentation strategy, and the computational algorithm. Representative numerical exam- ples are given to show the efficiency, stability, versatility, and reliability of the proposed approach.     

A CLASS OF VARIATIONAL DIFFERENCE SCHEMES FOR A SINGULAR PERTURBATION PROBLEM

In this paper,a singularly perturbed boundary value problem for second order self-adjoint ordinary differential equation is discussed.A class of variational difference schemesis constructed by the finite element method.Uniform convergence about small parameter isproved under a weaker smooth condition with respect to the coefficients of the equation.The schemes studied in refs.[1],[3],[4]and[5]belong to the class.     

Nonlinear vibration analysis of fractional viscoelastic Euler–Bernoulli nanobeams based on the surface stress theory

The nonlinear vibrations of viscoelastic Euler–Bernoulli nanobeams are studied using the fractional calculus and the Gurtin–Murdoch theory. Employing Hamilton's principle, the governing equation considering surface effects is derived. The fractional integro-partial differential governing equation is first converted into a fractional–ordinary differential equation in the time domain using the Galerkin scheme. Thereafter, the set of nonlinear fractional time-dependent equations expressed in a state-space form is solved using the predictor–corrector method. Finally, the effects of initial displacement, fractional derivative order, viscoelasticity coefficient, surface parameters and thickness-to-length ratio on the nonlinear time response of simply-supported and clamped-free silicon viscoelastic nanobeams are investigated.     

A symplectic finite element method based on Galerkin discretization for solving linear systems

We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems. For the dynamic responses of continuous medium structures, the traditional numerical algorithm is the dissipative algorithm and cannot maintain long-term energy conservation. Thus, a symplectic finite element method with energy conservation is constructed in this paper. A linear elastic system can be discretized into multiple elements, and a Hamiltonian system of each element ...     

Conservation Law and Application of J-Integral in Multi-Materials

The conservation law of J-integral in two-media with a crack paralleling to the interface of the two media was firstly proved by analytical and numerical finite element method. Then a schedule model was established that an interface crack is inserted in four media. According to the J-integral conservation law on multi-media, the energy release ratio of Ⅰ-type crack was considered to be conservation when the middle medium layers are very thin. And the conservation law was also convinced by numerical method. By means of the dimension analysis on the model, the asymptotic results and formula calculating the energy release ratio and complex stress intensity factor are presented.     

ON THE SINGULAR PERTURBATION OF A NONLINEAR ORDINARY DIFFERENTIAL EQUATION WITH TWO PARAMETERS

In this paper,the method of differential inequalities has been applied to study theboundary value problems of nonlinear ordinary differential equation with two parameters.The asymptotic solutions have been found and the remainders have been estimated.