A STABILIZED CRANK-NICOLSON MIXED FINITE VOLUME ELEMENT FORMULATION FOR THE NON-STATIONARY PARABOLIZED NAVIER-STOKES EQUATIONS

A time semi-discrete Crank-Nicolson （CN） formulation with second-order time accuracy for the non-stationary parabolized Navier-Stokes equations is firstly established. And then, a fully discrete stabilized CN mixed finite volume element （SCNMFVE） formu- lation based on two local Gaussian integrals and parameter-free with the second-order time accuracy is established directly from the time semi-discrete CN formulation so that it could avoid the discussion for semi-discrete SCNMFVE formulation with respect to spatial wriables and its theoretical analysis becomes very simple. Finally, the error estimates of SCNMFVE solutions are provided.

A New Elasticity and Finite Element Formulation for Different Young's Modulus When Tension and Compression Loadings

This paper presents a new elasticity and finite element formulation for different Young's modulus when tension and compression loadings in anisotropy media. The case studies, such as anisotropy and isotropy, were investigated. A numerical example was shown to find out the changes of neutral axis at the pure bending beams.

Nodeless variable finite element method for heat transfer analysis by means of flux-based formulation and mesh adaptation

Based on flux-based formulation, a nodeless variable element method is developed to analyze two-dimensional steady-state and transient heat transfer problems. The nodeless variable element employs quadratic interpolation functions to provide higher solution accuracy without necessity to actually generate additional nodes. The flux-based formulation is applied to reduce the complexity in deriving the finite element equations as compared to the conventional finite element method, The solution accuracy is further improved by implementing an adaptive meshing technique to generaie finite element mesh that can adapt and move along corresponding to the solution behavior. The technique generates small elements in the regions of steep solution gradients to provide accurate solution, and meanwhile it generates larger elements in the other regions where the solution gradients are slight to reduce the computational time and the computer memory. The effectiveness of the combined procedure is demonstrated by heat transfer problems that have exact solutions. These problems tire： （a） a steady-state heat conduction analysis in a square plate subjected to a highly localized surface heating, and （b） a transient heat conduction analysis in a long plate subjected to moving heat source.

Discrete formulation of mixed finite element methods for vapor deposition chemical reaction equations

The vapor deposition chemical reaction processes, which are of extremely extensive applications, can be classified as a mathematical model by the following governing nonlinear partial differential equations containing velocity vector, temperature field, pressure field, and gas mass field. The mixed finite element （MFE） method is employed to study the system of equations for the vapor deposition chemical reaction processes. The semidiscrete and fully discrete MFE formulations are derived. And the existence and convergence （error estimate） of the semidiscrete and fully discrete MFE solutions are demonstrated. By employing MFE method to treat the system of equations for the vapor deposition chemical reaction processes, the numerical solutions of the velocity vector, the temperature field, the pressure field, and the gas mass field can be found out simultaneously. Thus, these researches are not only of important theoretical means, but also of extremely extensive applied vistas.

A STABILIZED MIXED FINITE ELEMENT FORMULATION FOR THE NON-STATIONARY INCOMPRESSIBLE BOUSSINESQ EQUATIONS

In this study, we employ mixed finite element （MFE） method, two local Gauss integrals, and parameter-free to establish a stabilized MFE formulation for the non-stationary incompressible Boussinesq equations. We also provide the theoretical analysis of the existence, uniqueness, stability, and convergence of the stabilized MFE solutions for the stabilized MFE formulation.

A BICUBIC B-SPLINE FINITE ELEMENT METHOD FOR FOURTH-ORDER SEMILINEAR PARABOLIC OPTIMAL CONTROL PROBLEMS

A bicubic B-spline finite element method is proposed to solve optimal control problems governed by fourth-order semilinear parabolic partial differential equations.Its key feature is the selection of bicubic B-splines as trial functions to approximate the state and costate variables in two space dimensions.A Crank-Nicolson difference scheme is constructed for time discretization.The resulting numerical solutions belong to C2in space,and the order of the coefficient matrix is low.Moreover,the Bogner-Fox-Schmit element is considered for comparison.Two numerical experiments demonstrate the feasibility and effectiveness of the proposed method.

Multiscale Hybrid-Mixed Finite Element Method for Flow Simulation in Fractured Porous Media

The multiscale hybrid-mixed(MHM)method is applied to the numerical approximation of two-dimensional matrix fluid flow in porous media with fractures.The two-dimensional fluid flow in the reservoir and the one-dimensional flow in the discrete fractures are approximated using mixed finite elements.The coupling of the two-dimensional matrix flow with the one-dimensional fracture flow is enforced using the pressure of the one-dimensional flow as a Lagrange multiplier to express the conservation of fluid transfer between the fracture flow and the divergence of the one-dimensional fracture flux.A zero-dimensional pressure(point element)is used to express conservation of mass where fractures intersect.The issuing simulation is then reduced using the MHM method leading to accurate results with a very reduced number of global equations.A general system was developed where fracture geometries and conductivities are specified in an input file and meshes are generated using the public domain mesh generator GMsh.Several test cases illustrate the effectiveness of the proposed approach by comparing the multiscale results with direct simulations.

A co-rotational updated Lagrangian formulation for a 2D beam element with consideration of the deformed curvature

A tensor-based updated Lagrangian （UL） formulation for the geometrically nonlinear analysis of 2D beam-column structures is developed by using curvilinear coordinates, which has considered the effects of the deformed curvature. Between the known configuration C1 and the desired configuration C2, a configuration C2^＊ derived by rigid-body motion of C1 is introduced to eliminate the element-end transverse displacements between C2^＊ and C2. A stiffness matrix is obtained in C2^＊; and then by a transformation defined by the element-end displacements, the stiffness matrix in C2^＊ is transformed into that in CI. Comparing the stiffness matrix with that in the conventional UL formulation for a 2D beam element, the initial displacement stiffness matrix emerges, which results from the deformed curvature within the element. Numerical examples have verified the accuracy and efficiency of the present formulation, and the results show that the deformed curvatures have significant effects when deformations are large.

ANALYSIS OF AUGMENTED THREE-FIELD MACRO-HYBRID MIXED FINITE ELEMENT SCHEMES

On the basis of composition duality principles, augmented three-field macrohybrid mixed variational problems and finite element schemes are analyzed. The compatibility condition adopted here, for compositional dualization, is the coupling operator surjectivity, property that expresses in a general operator sense the Ladysenskaja-Babulka-Brezzi inf-sup condition. Variational macro-hybridization is performed under the assumption of decomposable primal and dual spaces relative to nonoverlapping domain decompositions. Then, through compositional dualization macro-hybrid mixed problems are obtained, with internal boundary dual traces as Lagrange multipliers. Also, ＂mass＂ preconditioned aug- mentation of three-field formulations are derived, stabilizing macro-hybrid mixed finite element schemes and rendering possible speed up of rates of convergence. Dual mixed incompressible Darcy flow problems illustrate the theory throughout the paper.

A FINITE ELEMENT SOLVER FOR NAVIER-STOKES EQUATIONS VIA VORTICITY AND VELOCITY

The incompressible Navier Stokes equations are solved via variables of vorticity and velocity. Firstly, a rigorous variational framework with the equivalence between the velocity pressure and the vorticity velocity formulations is presented in a Lipschitz domain. Next, a class of Galerkin finite element approximations of the corresponding variational form is introduced, and a convergence analysis is given for the Stokes problem. Finally, an iterative finite element solver for the Navier Stokes problem is proposed.

GEOMETRICALLY NONLINEAR FE FORMULATIONS FOR THE MACRO-ELEMENT UNIPLET OF FOLDABLE STRUCTURES

Geometrically nonlinear stiffness matrix due to large displacement small strain was firstly formulated explicitly for the basic components of pantographic foldable structures,namely, the uniplet, derived from a three node beam element.The formulation of the uniplet stiffness matrix is based on the precise nonlinear finite element theory and the displacement harmonized and internal force constraints are applied directly to the deformation modes of the three node beam element. The formulations were derived in general form, and can be simplified for particular foldable structures, such as flat, cylindrical and spherical structures.Finally, two examples were presented to illustrate the applications of the stiffness matrix evolved.

A FINITE ELEMENT ANALYSIS OF THE LARGE PLASTIC DEFORMATION BEHAVIOR OF AMORPHOUS GLASSY CIRCULAR POLYMERIC BARS

The BPA eight-chain molecular network model is introduced into the finite element formulation of elastic-plastic large deformation. And then, the tensile deformation localization development of the amorphous glassy circular polymeric bars (such as polycarbonates) is numerically simulated. The simulated results are compared with experimental ones, and very good consistence between numerical simulation and experiment is obtained, which shows the efficiency of the finite element analysis. Finally, the influences of the microstructure parameter S-ss on tensile neck-propagation and triaxial stress effect are studied.

Reference Coordinate System of Nonlinear Beam Element Based on the Geometrically Exact Formulation under Large Spatial Rotations and Deformations

Analysis of slender beam structures in a three-dimensional space is widely applicable in mechanical and civil engineering. This paper presents a new procedure to determine the reference coordinate system of a beam element under large rotation and elastic deformation based on a newly introduced physical concept: the zero twist sectional condition, which means that a non-twisted section between two nodes always exists and this section can reasonably be regarded as a reference coordinate system to calculate the internal element forces. This method can avoid the disagreement of the reference coordinates which might occur under large spatial rotations and deformations. Numerical examples given in the paper prove that this procedure guarantees the numerical exactness of the inherent formulation and improves the numerical efficiency, especially under large spatial rotations.

A REDUCED MFE FORMULATION BASED ON POD FOR THE NON-STATIONARY CONDUCTION-CONVECTION PROBLEMS

In this article, a reduced mixed finite element （MFE） formulation based on proper orthogonal decomposition （POD） for the non-stationary conduction-convection problems is presented. Also the error estimates between the reduced MFE solutions based on POD and usual MFE solutions are derived. It is shown by numerical examples that the results of numerical computation are consistent with theoretical conclusions. Moreover, it is shown that the reduced MFE formulation based on POD is feasible and efficient in finding numerical solutions for the non-stationary conduction-convection problems.

A REDUCED FE FORMULATION BASED ON POD METHOD FOR HYPERBOLIC EQUATIONS

A proper orthogonal decomposition（POD） method was successfully used in the reduced-order modeling of complex systems.In this paper,we extend the applications of POD method,namely,apply POD method to a classical finite element（FE） formulation for second-order hyperbolic equations with real practical applied background,establish a reduced FE formulation with lower dimensions and high enough accuracy,and provide the error estimates between the reduced FE solutions and the classical FE solutions and the implementation of algorithm for solving reduced FE formulation so as to provide scientific theoretic basis for service applications.Some numerical examples illustrate the fact that the results of numerical computation are consistent with theoretical conclusions.Moreover,it is shown that the reduced FE formulation based on POD method is feasible and efficient for solving FE formulation for second-order hyperbolic equations.

A REDUCED-ORDER MFE FORMULATION BASED ON POD METHOD FOR PARABOLIC EQUATIONS

In this paper, we extend the applications of proper orthogonal decomposition （POD） method, i.e., apply POD method to a mixed finite element （MFE） formulation naturally satisfied Brezz-Babu^ka for parabolic equations, establish a reduced-order MFE formulation with lower dimensions and sufficiently high accuracy, and provide the error estimates between the reduced-order POD MFE solutions and the classical MFE solutions and the implementation of algorithm for solving reduced-order MFE formulation. Some numerical examples illustrate the fact that the results of numerical computation are consis- tent with theoretical conclusions. Moreover, it is shown that the new reduced-order MFE formulation based on POD method is feasible and efficient for solving MFE formulation for parabolic equations.

A reduced-order extrapolation algorithm based on CNLSMFE formulation and POD technique for two-dimensional Sobolev equations

A reduced-order extrapolation algorithm based on Crank-Nicolson least-squares mixed finite element （CNLSMFE） formulation and proper orthogonal decomposition （POD） technique for two-dimensional （2D） Sobolev equations is established. The error estimates of the reduced-order CNLSMFE solutions and the implementation for the reduced-order extrapolation algorithm are provided. A numerical example is used to show that the results of numerical computations are consistent with theoretical conclusions. Moreover, it is shown that the reduced-order extrapolation algorithm is feasible and efficient for seeking numerical solutions to 2D Sobolev equations.

Finite element formulation based on proper orthogonal decomposition for parabolic equations

A proper orthogonal decomposition (POD) method is applied to a usual finite element (FE) formulation for parabolic equations so that it is reduced into a POD FE formulation with lower dimensions and enough high accuracy. The errors between the reduced POD FE solution and the usual FE solution are analyzed. It is shown by numerical examples that the results of numerical computations are consistent with theoretical conclusions. Moreover, it is also shown that this validates the feasibility and efficiency of POD method.

INTEGRATION OF NON-UNIFORM RATIONAL B-SPLINES GEOMETRY AND RATIONAL ABSOLUTE NODAL COORDINATES FORMULATION FINITE ELEMENT ANALYSIS

This investigation is intended to develop a computer procedure for the integration of NURBS geometry and the rational absolute nodal coordinate formulation （RANCF） finite element analysis. A linear transformation is given that can be used to convert the NURBS curve to RANCF cable element mesh retaining the same geometry and the same degree of continuity, including the discussion of continuity control and mesh refinement. The green strain tensor is used to establish the nonlinear dynamic equations with numerical examples to demonstrate the use of the procedure in the dynamic analysis of flexible bodies.

Finite element corotational formulation for geometric nonlinear analysis of thin shells with large rotation and small strain

Based on the consistent symrnetrizable equilibrated （CSE） corotational formulation, a linear triangular flat thin shell element with 3 nodes and 18~ of freedom, constructed by combination of the optimal membrane element and discrete Kirchhoff trian- gle （DKT） bending plate element, was extended to the geometric nonlinear analysis of thin shells with large rotation and small strain. Through derivation of the consistent tangent stiffness matrix and internal force vector, the corotational nonlinear finite element equations were established. The nonlinear equations were solved by using the Newton-Raphson iteration algorithm combined with an automatic load controlled technology. Three typical case studies, i.e., the slit annular thin plate, top opened hemispherical shell and cylindrical shell, validated the accuracy of the formulation established in this paper.