HIGH ACCURACY FINITE VOLUME ELEMENT METHOD FOR TWO-POINT BOUNDARY VALUE PROBLEM OF SECOND ORDER ORDINARY DIFFERENTIAL EQUATIONS

In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.

AN ADI GALERKIN METHOD WITH MOVING FINITE ELEMENT SPACES FOR A CLASS OF SECOND-ORDER HYPERBOLIC EQUATIONS

An alternating direction implicit (ADI) Galerkin method with moving finite element spaces is formulated for a class of second order hyperbolic equations in two space variables. A priori H 1 error estimate is derived.

Numerical Experiments Using MATLAB: Superconvergence of Conforming Finite, Element Approximation for Second Order, Elliptic Problems

The superconvergence in the finite element method is a phenomenon in which the finite element approximation converges to the exact solution at a rate higher than the optimal order error estimate. Wang proposed and analyzed superconvergence of the conforming finite element method by L2-projections. The goal of this paper is to perform numerical experiments using MATLAB to support and to verify the theoretical results in Wang for the superconvergence of the conforming finite element method (CFEM) for the second order elliptic problems by L2-projection methods. MATLAB codes are published at https://github.com/annaleeharris/Superconvergence-CFEM for anyone to use and to study.

Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems

The superconvergence in the finite element method is a phenomenon in which the fi-nite element approximation converges to the exact solution at a rate higher than the optimal order error estimate. Wang proposed and analyzed superconvergence of the conforming finite element method by L2-projections. However, since the conforming finite element method (CFEM) requires a strong continuity, it is not easy to construct such finite elements for the complex partial differential equations. Thus, the nonconforming finite element method (NCFEM) is more appealing computationally due to better stability and flexibility properties compared to CFEM. The objective of this paper is to establish a general superconvergence result for the nonconforming finite element approximations for second-order elliptic problems by L2-projection methods by applying the idea presented in Wang. MATLAB codes are published at https://github.com/annaleeharris/Superconvergence-NCFEM for anyone to use and to study. The results of numerical experiments show great promise for the robustness, reliability, flexibility and accuracy of superconvergence in NCFEM by L2- projections.

SQUEEZE FLOW OF A SECOND-ORDER FLUID BETWEEN TWO PARALLEL DISKS OR TWO SPHERES

The normal viscous force of squeeze flow between two arbitrary rigid spheres with an interstitial second-order fluid was studied for modeling wet granular materials using the discrete element method. Based on the Reynolds' lubrication theory, the small parameter method was introduced to approximately analyze velocity field and stress distribution between the two disks. Then a similar procedure was carried out for analyzing the normal interaction between two nearly touching, arbitrary rigid spheres to obtain the pressure distribution and the resulting squeeze force. It has been proved that the solutions can be reduced to the case of a Newtonian fluid when the non-Newtonian terms are neglected.

A Second Order Characteristic Mixed Finite Element Method for Convection Diffusion Reaction Equations

A combined approximate scheme is defined for convection-diffusion-reaction equations. This scheme is constructed by two methods. Standard mixed finite element method is used for diffusion term. A second order characteristic finite element method is presented to handle the material derivative term, that is, the time derivative term plus the convection term. The stability is proved and the L2-norm error estimates are derived for both the scalar unknown variable and its flux. The scheme is of second order accuracy in time increment, symmetric, and unconditionally stable.

Dynamic Analysis of Kineto-Elastic Beam System with Second-order Effect

Dynamic equations of motional flexible beam elements were derived considering second-order effect. Non-linear finite element method and three-node Euler-Bernoulli beam elements were used. Because accuracy is higher in non-linear structural analysis,three-node beam elements are used to deduce shape functions and stiffness matrices in dynamic equations of flexible elements. Static condensation method was used to obtain the finial dynamic equations of three-node beam elements. According to geometrical relations of nodal displacements in concomitant and global coordinate system,dynamic equations of elements can be transformed to global coordinate system by concomitant coordinate method in order to build the global dynamic equations. Analyzed amplitude condition of flexible arm support of a port crane,the results show that second-order effect should be considered in kinetic-elastic analysis for heavy load machinery of big flexibility.

Finite element analysis of second order wave radiation by a group of cylinders in the time domain

A finite element based numerical method is employed to analyze the wave radiation by multiple or a group of cylinders in the time domain. The nonlinear free surface and body surface boundary conditions are satisfied based on the perturbation method up to the second order. The first- and second-order velocity potential problems at each time step are solved through a Finite Element Method （FEM）. The matrix equation of the FEM is solved through iteration and the initial solution is obtained from the result at the previous time step. The three-dimensional （3-D） mesh required is generated based on a two-dimensional （2-D） hybrid mesh on a horizontal plane and its extension in the vertical direction. The hybrid mesh is generated by combining an unstructured grid away from cylinders and two structured grids near the cylinder and the artificial boundary. The fluid velocity on the free surface and the cylinder surface are calculated by using a differential method. Results for various configurations including the cases of two cylinders and four cylinders and a group of eighteen cylinders are obtained to show the joint influences of cylinders on the first- and second- order waves and forces, including the effects of spacing ratios and wave frequency on the second order waves and the mean force, in particular.

Non-Split PML Boundary Condition for Finite Element Time-Domain Modeling of Ground Penetrating Radar

As a highly efficient absorbing boundary condition, Perfectly Matched Layer (PML) has been widely used in Finite Difference Time Domain (FDTD) simulation of Ground Penetrating Radar (GPR) based on the first order electromagnetic wave equation. However, the PML boundary condition is difficult to apply in GPR Finite Element Time Domain (FETD) simulation based on the second order electromagnetic wave equation. This paper developed a non-split perfectly matched layer (NPML) boundary condition for GPR FETD simulation based on the second order electromagnetic wave equation. Taking two-dimensional TM wave equation as an example, the second order frequency domain equation of GPR was derived according to the definition of complex extending coordinate transformation. Then it transformed into time domain by means of auxiliary differential equation method, and its FETD equation is derived based on Galerkin method. On this basis, a GPR FETD forward program based on NPML boundary condition is developed. The merits of NPML boundary condition are certified by compared with wave field snapshots, signal and reflection errors of homogeneous medium model with split and non-split PML boundary conditions. The comparison demonstrated that the NPML algorithm can reduce memory occupation and improve calculation efficiency. Furthermore, numerical simulation of a complex model verifies the good absorption effects of the NPML boundary condition in complex structures.

Convergence analysis of the rectangular Morley element scheme for second order problem in arbitrary dimensions

We present the convergence analysis of the rectangular Morley element scheme utilised on the second order problem in arbitrary dimensions. Specifically, we prove that the convergence of the scheme is of O(h) order in energy norm and of O(h^2) order in L^2 norm on general d-rectangular triangulations. Moreover, when the triangulation is uniform, the convergence rate can be of O(h^2) order in energy norm, and the convergence rate in L^2 norm is still of O(h^2) order, which cannot be improved. Numerical examples are presented to demonstrate our theoretical results.

CONVERGENCE ANALYSIS OF MIXED VOLUME ELEMENT-CHARACTERISTIC MIXED VOLUME ELEMENT FOR THREE-DIMENSIONAL CHEMICAL OIL-RECOVERY SEEPAGE COUPLED PROBLEM

The physical model is described by a seepage coupled system for simulating numerically three-dimensional chemical oil recovery, whose mathematical description includes three equations to interpret main concepts. The pressure equation is a nonlinear parabolic equation, the concentration is defined by a convection-diffusion equation and the saturations of different components are stated by nonlinear convection-diffusion equations. The transport pressure appears in the concentration equation and saturation equations in the form of Darcy velocity, and controls their processes. The flow equation is solved by the conservative mixed volume element and the accuracy is improved one order for approximating Darcy velocity. The method of characteristic mixed volume element is applied to solve the concentration, where the diffusion is discretized by a mixed volume element method and the convection is treated by the method of characteristics. The characteristics can confirm strong computational stability at sharp fronts and it can avoid numerical dispersion and nonphysical oscillation. The scheme can adopt a large step while its numerical results have small time-truncation error and high order of accuracy. The mixed volume element method has the law of conservation on every element for the diffusion and it can obtain numerical solutions of the concentration and adjoint vectors. It is most important in numerical simulation to ensure the physical conservative nature. The saturation different components are obtained by the method of characteristic fractional step difference. The computational work is shortened greatly by decomposing a three-dimensional problem into three successive one-dimensional problems and it is completed easily by using the algorithm of speedup. Using the theory and technique of a priori estimates of differential equations, we derive an optimal second order estimates in 12 norm. Numerical examples are given to show the effectiveness and practicability and the method is testified as a powerful tool to solve the important problems.

A SECOND ORDER UNCONDITIONALLY CONVERGENT FINITE ELEMENT METHOD FOR THE THERMAL EQUATION WITH JOULE HEATING PROBLEM

In this paper,we study the finite element approximation for nonlinear thermal equation.Because the nonlinearity of the equation,our theoretical analysis is based on the error of temporal and spatial discretization.We consider a fully discrete second order backward difference formula based on a finite element method to approximate the temperature and electric potential,and establish optimal L^(2) error estimates for the fully discrete finite element solution without any restriction on the time-step size.The discrete solution is bounded in infinite norm.Finally,several numerical examples are presented to demonstrate the accuracy and efficiency of the proposed method.

A 3D pyramid spline element

In this paper,a 13-node pyramid spline element is derived by using the tetrahedron volume coordinates and the B-net method,which achieves the second order completeness in Cartesian coordinates.Some appropriate examples were employed to evaluate the performance of the proposed element.The numerical results show that the spline element has much better performance compared with the isoparametric serendipity element Q20 and its degenerate pyramid element P13 especially when mesh is distorted,and it is comparable to the Lagrange element Q27.It has been demonstrated that the spline finite element method is an efficient tool for developing high accuracy elements.

Finite Element Analysis to Two-Dimensional Nonlinear Sloshing Problems

A two-dimensional nonlinear sloshing problem is analyzed by means of the fully nonlinear theory and time domain second order theory of water waves. Liquid sloshing in a rectangular container Subjected to a horizontal excitation is simulated by the finite element method. Comparisons between the two theories are made based on their numerical results. It is found that good agreement is obtained for the case of small amplitude oscillation and obvious differences occur for large amplitude excitation. Even though, the second order solution can still exhibit typical nonlinear features of nonlinear wave and can be used instead of the fully nonlinear theory.

METHOD OF NONCONFORMING MIXED FINITE ELEMENTFOR SECOND ORDER ELLIPTIC PROBLEMSMETHOD OF NONCONFORMING MIXED FINITE ELEMENTFOR SECOND ORDER ELLIPTIC PROBLEM

In this paper, the method of non-conforming mixed finite element for second order elliptic problems is discussed and a format of real optimal order for the lowest order error estimate.

A Method for Calculation of Low-Frequency Slow Drift Motions Based on NURBS for Floating Bodies

Through a higher-order boundary element method based on NURBS （Non-uniform Rational B-splines）, the calculation of second-order low-frequency forces and slow drift motions is conducted for floating bodies. In the floating body＇s inner domain, an auxiliary equation is obtained by applying a Green function which satisfies the solid surface condition. Then, the auxiliary equation and the velocity potential equation are combined in the fluid domain to remove the solid angle coefficient and the singularity of the double layer potentials in the integral equation. Thus, a new velocity potential integral equation is obtained. The new equation is extended to the inner domain to reheve the irregular frequency effects; on the basis of the order analysis, the comparison is made about the contribution of all integral terms with the result in the second-order tow-frequency problem; the higher-order boundary element method based on NURBS is apphed to calculate the geometric position and velocity potentials; the slow drift motions are calculated by the spectrum analysis method. Removing the solid angle coefficient can apply NURBS technology to the hydrodynamic calculation of floating bodies with complex surfaces, and the extended boundary integral method can reduce the irregular frequency effects. Order analysis shows that free surface integral can be neglected, and the numerical results can also prove the correctness of order analysis. The results of second-order low-frequency forces and slow drift motions and the comparison with the results from references show that the application of the NURBS technology to the second-order low-frequency problem is of high efficiency and credible results.

NUMERICAL METHOD OF MIXED FINITE VOLUME-MODIFIED UPWIND FRACTIONAL STEP DIFFERENCE FOR THREE-DIMENSIONAL SEMICONDUCTOR DEVICE TRANSIENT BEHAVIOR PROBLEMS

Transient behavior of three-dimensional semiconductor device with heat conduc- tion is described by a coupled mathematical system of four quasi-linear partial differential equations with initial-boundary value conditions. The electric potential is defined by an ellip- tic equation and it appears in the following three equations via the electric field intensity. The electron concentration and the hole concentration are determined by convection-dominated diffusion equations and the temperature is interpreted by a heat conduction equation. A mixed finite volume element approximation, keeping physical conservation law, is used to get numerical values of the electric potential and the accuracy is improved one order. Two con- centrations and the heat conduction are computed by a fractional step method combined with second-order upwind differences. This method can overcome numerical oscillation, dispersion and decreases computational complexity. Then a three-dimensional problem is solved by computing three successive one-dimensional problems where the method of speedup is used and the computational work is greatly shortened. An optimal second-order error estimate in L2 norm is derived by using prior estimate theory and other special techniques of partial differential equations. This type of mass-conservative parallel method is important and is most valuable in numerical analysis and application of semiconductor device.

A Mixed-finite Volume Element Coupled with the Method of Characteristic Fractional Step Difference for Simulating Transient Behavior of Semiconductor Device of Heat Conductor And Its Numerical Analysis

The mathematical system is formulated by four partial differential equations combined with initial- boundary value conditions to describe transient behavior of three-dimensional semiconductor device with heat conduction. The first equation of an elliptic type is defined with respect to the electric potential, the successive two equations of convection dominated diffusion type are given to define the electron concentration and the hole concentration, and the fourth equation of heat conductor is for the temperature. The electric potential appears in the equations of electron concentration, hole concentration and the temperature in the formation of the intensity. A mass conservative numerical approximation of the electric potential is presented by using the mixed finite volume element, and the accuracy of computation of the electric intensity is improved one order. The method of characteristic fractional step difference is applied to discretize the other three equations, where the hyperbolic terms are approximated by a difference quotient in the characteristics and the diffusion terms are discretized by the method of fractional step difference. The computation of three-dimensional problem works efficiently by dividing it into three one-dimensional subproblems and every subproblem is solved by the method of speedup in parallel. Using a pair of different grids （coarse partition and refined partition）, piecewise threefold quadratic interpolation, variation theory, multiplicative commutation rule of differential operators, mathematical induction and priori estimates theory and special technique of differential equations, we derive an optimal second order estimate in L2-norm. This numerical method is valuable in the simulation of semiconductor device theoretically and actually, and gives a powerful tool to solve the international problem presented by J. Douglas, Jr.

Numerical Simulation of Two-dimensional Nonlinear Sloshing Problems

Numerical simulation of a two-dimensional nonlinear sloshing problem is preceded by the finite element method. Two theories are used. One is fully nonlinear theory; the other is time domain second order theory. A liquid sloshing in a rectangular container subjected to a horizontal excitation is simulated using these two theories. Numerical results are obtained and comparisons are made. It is found that a good agreement is obtained for the case of small amplitude oscillation. For the situation of large amplitude excitation, although the differences between using the two theories are obvious the second order solution can still exhibit typical nonlinear features of nonlinear wave.

塔式起重机二阶侧移计算

为计算塔式起重机塔身二阶效应引起的二阶水平侧移,结合工程实际情况,用理论推导方式得到近似计算公式,对QTP125型塔式起重机进行现场测验和仿真模拟,验证计算公式的合理性。结果表明:计算公式可以用于快速计算实际工程中塔身静止状态下的二阶侧移,能为塔式起重机结构安全设计提供参考。