SPACE-TIME FINITE ELEMENT METHOD FOR SCHRDINGER EQUATION AND ITS CONSERVATION

Energy conservation of nonlinear Schrodinger ordinary differential equation was proved through using continuous finite element methods of ordinary differential equation; Energy integration conservation was proved through using space-time continuous fully discrete finite element methods and the electron nearly conservation with higher order error was obtained through using time discontinuous only space continuous finite element methods of nonlinear Schrodinger partial equation. The numerical results are in accordance with the theory.

THE SPACE-TIME FINITE ELEMENT METHOD FOR PARABOLIC PROBLEMS

Adaptive space-time finite element method, continuous in space but discontinuous in time for semi-linear parabolic problems is discussed. The approach is based on a combination of finite element and finite difference techniques. The existence and uniqueness of the weak solution are proved without any assumptions on choice of the spacetime meshes. Basic error estimates in L-infinity (L-2) norm, that is maximum-norm in time, L-2-norm in space are obtained. The numerical results are given in the last part and the analysis between theoretic and experimental results are 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.

THE ELASTIC SOLUTION OF CONCENTRATED FORCE ACTING IN ORTHOGONAL ANISOTROPIC HALF-PLANE AND CONSTANT ELEMENT FUNDAMENTAI FORMULAE OF BOUNDARY ELEMENT METHOD

In this paper, the elastic solutions of concentrated force acting in orthogonal anisotropic half-plane are derived by imaginal method and the formulae of coefficient matrix for constant element are put forward. To solve half-plane problems numerically by BEM, this paper provides the necessary formulae. Because the expressions of fundamental solutions are very simple, the. object functions could be obtained for every integral of constant element and higher order element of indirect BEM. Thus, the procedure of integration could be avoided in calculation program

SPECTRAL/HP ELEMENT METHOD WITH HIERARCHICAL RECONSTRUCTION FOR SOLVING NONLINEAR HYPERBOLIC CONSERVATION LAWS

The hierarchical reconstruction （HR） [Liu, Shu, Tadmor and Zhang, SINUM ＇07] has been successfully applied to prevent oscillations in solutions computed by finite volume, Runge-Kutta discontinuous Galerkin, spectral volume schemes for solving hyperbolic conservation laws. In this paper, we demonstrate that HR can also be combined with spectral/hp element method for solving hyperbolic conservation laws. An orthogonal spectral basis written in terms of Jacobi polynomials is applied. High computational efficiency is obtained due to such matrix-free algorithm. The formulation is conservative, and essential nomoscillation is enforced by the HR limiter. We show that HR preserves the order of accuracy of the spectral/hp element method for smooth solution problems and generate essentially non-oscillatory solutions profiles for capturing discontinuous solutions without local characteristic decomposition. In addition, we introduce a postprocessing technique to improve HR for limiting high degree numerical solutions.

THE ANALYTICAL SOLUTIONS BASED ON THE CONCEPT OF FINITE ELEMENT METHODS

On the basis of the concept of finite element methods, the rigorous analytical solutions of structural response in terms of the design variables are researched in this paper. The spatial trusses are taken as an example for the solution of the analytical expressions of the explicit displacements which are proved mathematically; then some conclusions are reached that are useful to structural sensitivity analysis and optimization. In the third part of the paper, a generalized geometric programming method is sugguested for the optimal model with the explicit displacement. Finally, the analytical solutions of the displacements of three trusses are given as examples.

Spectral Element Viscosity Methods for Nonlinear Conservation Laws on the Semi-Infinite Interval

In this paper we propose a spectral element: vanishing viscosity (SEW) method for the conservation laws on the semi-infinite interval. By using a suitable mapping, the problem is first transformed into a modified conservation law in a bounded interval, then the well-known spectral vanishing viscosity technique is generalized to the multi-domain case in order to approximate this trarsformed equation more efficiently. The construction details and convergence analysis are presented. Under a usual assumption of boundedness of the approximation solutions, it is proven that the solution of the SEW approximation converges to the uniciue entropy solution of the conservation laws. A number of numerical tests is carried out to confirm the theoretical results.

FAST SOLUTION FOR LARGE SCALE LINEAR ALGEBRAIC EQUATIONS IN FINITE ELEMENT ANALYSIS

The computational efficiency of numerical solution of linearalgebraic equations in finite elements can be improved in two ways.One is to decrease the fill-in numbers, which are new non-ze- ronumbers in the matrix of global stiffness generated during theprocess of elimination. The other is to reduce the computationaloperation of multiplying a real number by zero. Based on the factthat the order of elimination can determine how many fill-in numbersshould be generated, we present a new method for optimization ofnumbering nodes. This method is quite different from bandwidthoptimiza- tion. Fill-in numbers can be decreased in a large scale bythe use of this method. The bi-factorization method is adopted toavoid multiplying real numbers by zero. For large scale finiteelement analysis, the method presented in this paper is moreefficient than the traditional LDLT method.

A modified discrete element method for concave granular materials based on energy-conserving contact model

The development of a general discrete element method for irregularly shaped particles is the core issue of the simulation of the dynamic behavior of granular materials.The general energy-conserving contact theory is used to establish a universal discrete element method suitable for particle contact of arbitrary shape.In this study,three dimentional(3D)modeling and scanning techniques are used to obtain a triangular mesh representation of the true particles containing typical concave particles.The contact volumebased energy-conserving model is used to realize the contact detection between irregularly shaped particles,and the contact force model is refined and modified to describe the contact under real conditions.The inelastic collision processes between the particles and boundaries are simulated to verify the robustness of the modified contact force model and its applicability to the multi-point contact mode.In addition,the packing process and the flow process of a large number of irregular particles are simulated with the modified discrete element method(DEM)to illustrate the applicability of the method of complex problems.

Solution to Equivalent Power Element of a Six Degrees of Freedom Simulator

The dynamical equation of a six degrees of freedom(DOF) Stewart-type simulator with load and impact force is derived. By associating the direct solution to acceleration with the inverse solution to force, an equivalent power element model of the simulator is further presented, which offers the basis for analysis and design of the control system of the simulator.

Numerical simulation of oxide nanoparticle growth characteristics under the gas detonation chemical reaction by space-time conservation element-solution element method

Under harsh conditions （such as high temperature, high pressure, and millisecond lifetime chemical reaction）, a long-standing challenge remains to accurately predict the growth characteristics of nanosize spherical particles and to determine the rapid chemical reaction flow field characteristics, The growth characteristics of similar spherical oxide nanoparticles are further studied by successfully introducing the space-time conservation element-solution element （CE/SE） algorithm with the monodisperse Kruis model. This approach overcomes the nanosize particle rapid growth limit set and successfully captures the characteristics of the rapid gaseous chemical reaction process. The results show that this approach quantitatively captures the characteristics of the rapid chemical reaction, nanosize particle growth and size distribution. To reveal the growth mechanism for numerous types of oxide nanoparticles, it is very important to choose a rational numerical method and particle physics model.

Self-adaptive strategy for one-dimensional finite element method based on EEP method with optimal super-convergence order

Based on the newly-developed element energy projection （EEP） method with optimal super-convergence order for computation of super-convergent results, an improved self-adaptive strategy for one-dimensional finite element method （FEM） is proposed. In the strategy, a posteriori errors are estimated by comparing FEM solutions to EEP super-convergent solutions with optimal order of super-convergence, meshes are refined by using the error-averaging method. Quasi-FEM solutions are used to replace the true FEM solutions in the adaptive process. This strategy has been found to be simple, clear, efficient and reliable. For most problems, only one adaptive step is needed to produce the required FEM solutions which pointwise satisfy the user specified error tolerances in the max-norm. Taking the elliptical ordinary differential equation of the second order as the model problem, this paper describes the fundamental idea, implementation strategy and computational algorithm and representative numerical examples are given to show the effectiveness and reliability of the proposed approach.

SELF-ADAPTIVE STRATEGY FOR ONE-DIMENSIONAL FINITE ELEMENT METHOD BASED ON ELEMENT ENERGY PROJECTION METHOD

Based on the newly-developed element energy projection （EEP） method for computation of super-convergent results in one-dimensional finite element method （FEM）, the task of self-adaptive FEM analysis was converted into the task of adaptive piecewise polynomial interpolation. As a result, a satisfactory FEM mesh can be obtained, and further FEM analysis on this mesh would immediately produce an FEM solution which usually satisfies the user specified error tolerance. Even though the error tolerance was not completely satisfied, one or two steps of further local refinements would be sufficient. This strategy was found to be very simple, rapid, cheap and efficient. Taking the elliptical ordinary differential equation of second order as the model problem, the fundamental idea, implementation strategy and detailed algorithm are described. Representative numerical examples are given to show the effectiveness and reliability of the proposed approach.

The space time CE/SE method for solving one-dimensional batch crystallization model with fines dissolution

This article is concerned with the numerical investigation of one-dimensional population balance models for batch crystallization process with fines dissolution.In batch crystallization,dissolution of smaller unwanted nuclei below some critical size is of vital importance as it improves the quality of product.The crystal growth rates for both size-independent and size-dependent cases are considered.A delay in recycle pipe is also included in the model.The space–time conservation element and solution element method,originally derived for non-reacting flows,is used to solve the model.This scheme has already been applied to a range of PDEs,mainly in the area of fluid mechanics.The numerical results are compared with those obtained from the Koren scheme,showing that the proposed scheme is more efficient.

Spline Fictitious Boundary Element Alternating Method for Edge Crack Problems with Mixed Boundary Conditions

The alternating method based on the fundamental solutions of the infinite domain containing a crack,namely Muskhelishvili’s solutions,divides the complex structure with a crack into a simple model without crack which can be solved by traditional numerical methods and an infinite domain with a crack which can be solved by Muskhelishvili’s solutions.However,this alternating method cannot be directly applied to the edge crack problems since partial crack surface of Muskhelishvili’s solutions is located outside the computational domain.In this paper,an improved alternating method,the spline fictitious boundary element alternating method(SFBEAM),based on infinite domain with the combination of spline fictitious boundary element method(SFBEM)and Muskhelishvili’s solutions is proposed to solve the edge crack problems.Since the SFBEM and Muskhelishvili’s solutions are obtained in the framework of infinite domain,no special treatment is needed for solving the problem of edge cracks.Different mixed boundary conditions edge crack problems with varies of computational parameters are given to certify the high precision,efficiency and applicability of the proposed method compared with other alternating methods and extend finite element method.

Axisymmetric Finite Element Method for Analysis of Plate on Layered Soil

To obtain the fundamental solution of soil has become the key problem for the semi-analytical and semi-numerical (SASN) method in analyzing plate on layered soil. By applying axisymmetric finite element method (FEM),an expression relating the surface settlement and the reaction of the layered soil can be obtained. Such a reaction can be treated as load acting on the applied external load. Having the plate modelled by four-node elements,the governing equation of the plate can be formed and solved. In this case, the fundamental solution can be introduced into the global soil stiffness matrix and five-node or nine-node element soil stiffness matrix.The existing commercial FEM software can be used to solve the fundamental solution of soil, which can bypass the complicated formula derivation and boasts high computational efficiency as well.

UNIQUENESS OF SOLUTION OF FIELD POINT OF SINGULAR SOURCE OUTSIDE-REGION-DISTRIBUTION METHOD

The uniqueness of solution of field point, inside a convex region due to singular source(s) with kernel function decreasing with distance increasing, outside-region-distribution(s) such that the boundary condition expressed by the response of the source(s) is satisfied, is proved by using the condition of kernel function decreasing with distance increasing anal an integral inequality. Examples of part of these singular sources such as Kelvin's point force, Point-Ring-Couple (PRC) etc. are given. The proof of uniqueness of solution of field point in a twisted shaft of revolution due to PRC distribution is given as an example of application.

AN ITERATIVE PARALLEL ALGORITHM OF FINITE ELEMENT METHOD

In this paper, a parallel algorithm with iterative form for solving finite element equation is presented. Based on the iterative solution of linear algebra equations, the parallel computational steps are introduced in this method. Also by using the weighted residual method and choosing the appropriate weighting functions, the finite element basic form of parallel algorithm is deduced. The program of this algorithm has been realized on the ELXSI-6400 parallel computer of Xi'an Jiaotong University. The computational results show the operational speed will be raised and the CPU time will be cut down effectively. So this method is one kind of effective parallel algorithm for solving the finite element equations of large-scale structures.

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.

A Conceptual Numerical Model of the Wave Equation Using the Complex Variable Boundary Element Method

In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier series. The technique described in this work is suitable for modeling initial-boundary value problems governed by the wave equation on a rectangular domain with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The new numerical scheme is based on the standard approach of decomposing the global initial-boundary value problem into a steady-state component and a time-dependent component. The steady-state component is governed by the Laplace PDE and is modeled with the CVBEM. The time-dependent component is governed by the wave PDE and is modeled using a generalized Fourier series. The approximate global solution is the sum of the CVBEM and generalized Fourier series approximations. The boundary conditions of the steady-state component are specified as the boundary conditions from the global BVP. The boundary conditions of the time-dependent component are specified to be identically zero. The initial condition of the time-dependent component is calculated as the difference between the global initial condition and the CVBEM approximation of the steady-state solution. Additionally, the generalized Fourier series approximation of the time-dependent component is fitted so as to approximately satisfy the derivative of the initial condition. It is shown that the strong formulation of the wave PDE is satisfied by the superposed approximate solutions of the time-dependent and steady-state components.