The k-Version of Finite Element Method for Self-Adjoint Operators in BVP

In this paper a new mathematical and computational framework for boundary value problems described by self-adjoint differential operators is presented. In this framework, numerically computed solutions, when converged, possess the same degree of global smoothness in terms of differentiability up to any desired order as the theoretical solutions. This is accomplished using spaces Ĥk,p that contain basis functions of degree p and order k - 1 (or the order of the space k). It is shown that the order of space k is an intrinsically important independent parameter in all finite element computational processes in addition to the discretization characteristic length h and the degree of basis functions p when the theoretical solutions are analytic. Thus, in all finite element computations, all quantities of interest (e.g., quadratic functional, error or residual functional, norms and seminorms, error norms, etc.) are dependent on h, p as well as k. Therefore, for fixed h and p, convergence of the finite element process can also be investigated by changing k, hence k-convergence and thus the k-version of finite element method. With h, p, and k as three independent parameters influencing all finite element processes, we now have k, hk, pk, and hpk versions of finite element methods. The issue of minimally conforming finite element spaces is reexamined and it is demonstrated that the definition of currently believed minimally conforming space which permit weak convergence of the highest-order derivatives of the dependent variables appearing in the bilinear form is not justifiable mathematically or from physics view point. A new criterion is proposed for establishing the minimally conforming spaces which is more in agreement with the physics and mathematics of the BVP. Significant features and merits of the proposed mathematical and computational framework are presented, discussed, illustrated, and substantiated mathematically as well as numerically with the Galerkin and least-squares finite element formulations for self-adjoint boundary-value problems.

[1]  I. Norman Katz,et al.  The p-Version of the Finite Element Method for Problems Requiring $C^1$-Continuity , 1985 .

[2]  Graham F. Carey,et al.  Least‐squares finite element method and preconditioned conjugate gradient solution , 1987 .

[3]  B. Rubinsky,et al.  A mixed‐variable continuously deforming finite element method for parabolic evolution problems. Part III: Numerical implementation and computational results , 1989 .

[4]  O. C. Zienkiewicz,et al.  An ‘upwind’ finite element scheme for two‐dimensional convective transport equation , 1977 .

[5]  E. Eason.  A review of least-squares methods for solving partial differential equations , 1976 .

[6]  Alberto Peano,et al.  Hierarchies of conforming finite elements for plane elasticity and plate bending , 1976 .

[7]  B. Jiang.  The Least-Squares Finite Element Method , 1998 .

[8]  C. D. Mote,et al.  Solution of non‐linear boundary value problems by discrete least squares , 1977 .

[9]  L. R. Scott,et al.  A nodal basis for ¹ piecewise polynomials of degree ≥5 , 1975 .

[10]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[11]  Alberto Peano,et al.  Conforming approximations for Kirchhoff plates and shells , 1979 .

[12]  T. Hughes,et al.  A new finite element formulation for computational fluid dynamics. X - The compressible Euler and Navier-Stokes equations , 1991 .

[13]  Max Gunzburger,et al.  On least squares approximations to indefinite problems of the mixed type , 1978 .

[14]  Philip M. Gresho,et al.  Finite element simulations of steady, two-dimensional, viscous incompressible flow over a step , 1981 .

[15]  Karan S. Surana,et al.  p‐version least‐squares finite element formulation for convection‐diffusion problems , 1993 .

[16]  Paul P. Lynn.  Least squares finite element analysis of laminar boundary layer flows , 1974 .

[17]  B. Rubinsky,et al.  A mixed-variable continuously deforming finite element method for parabolic evolution problems. Part I: The variational formulation for a single evolution equation , 1989 .

[18]  K. S. Surana,et al.  Nonweak/Strong Solutions of Linear and Nonlinear Hyperbolic and Parabolic Equations Resulting from a Single Conservation Law , 2000, Int. J. Comput. Eng. Sci..

[19]  O. C. Zienkiewicz,et al.  Quadratic finite element schemes for two-dimensional convective-transport problems , 1977 .

[20]  Graham F. Carey,et al.  Least‐squares finite elements for first‐order hyperbolic systems , 1988 .

[21]  Paul P. Lynn,et al.  Efficient least squares finite elements for two‐dimensional laminar boundary layer analysis , 1976 .

[22]  Thomas J. R. Hughes,et al.  A new finite element formulation for computational fluid dynamics: IV. A discontinuity-capturing operator for multidimensional advective-diffusive systems , 1986 .

[23]  Karan S. Surana,et al.  p‐version least‐squares finite element formulation of Burgers' equation , 1993 .

[24]  Karan S. Surana,et al.  Investigation of diffusion in p-version ‘LSFE’ and ‘STLSFE’ formulations , 1995 .

[25]  Robert L. Lee,et al.  Don''t suppress the wiggles|they''re telling you something! Computers and Fluids , 1981 .

[26]  T. Hughes,et al.  A new finite element formulation for computational fluid dynamics: II. Beyond SUPG , 1986 .

[27]  G. Carey,et al.  A stable least‐squares finite element method for non‐linear hyperbolic problems , 1988 .

[28]  Graham F. Carey,et al.  Adaptive refinement for least‐squares finite elements with element‐by‐element conjugate gradient solution , 1987 .

[29]  K. S. Surana,et al.  Nonweak/Strong Solutions in Gas Dynamics: A C11p-Version STLSFEF in Lagrangian Frame of Reference Using r, U, T Primitive Variables , 2001, Int. J. Comput. Eng. Sci..

[30]  Karan S. Surana,et al.  A SPACE-TIME COUPLED P-VERSION LEAST SQUARES FINITE ELEMENT FORMULATION FOR UNSTEADY TWO-DIMENSIONAL NAVIER-STOKES EQUATIONS , 1996 .

[31]  T. Hughes,et al.  Stabilized finite element methods. I: Application to the advective-diffusive model , 1992 .

[32]  Karan S. Surana,et al.  p‐version least squares finite element formulation for two‐dimensional, incompressible fluid flow , 1994 .

[33]  On quadratic elements in finite element solutions of steady‐state convection—diffusion equation , 1980 .

[34]  B. Jiang,et al.  Least-squares finite element method for fluid dynamics , 1990 .

[35]  J. F. T. Pittman,et al.  Isoparametric Hermite elements , 1994 .

[36]  Karan S. Surana,et al.  Non‐weak/strong solutions in gas dynamics: a C11 p‐version STLSFEF in Lagrangian frame of reference using ρ, u, p primitive variables , 2002 .

[37]  Karan S. Surana,et al.  p-Version least squares finite element formulation for axisymmetric incompressible non-Newtonian fluid flow , 1994 .

[38]  Graham F. Carey,et al.  Nonlinear preconditioned conjugate gradient and least-squares finite elements , 1987 .

[39]  Ching L. Chang,et al.  An error analysis of least-squares finite element method of velocity-pressure-vorticity formulation for Stokes problem , 1990 .

[40]  A new mixed formulation for 2D incompressible flows , 1984 .

[41]  O. C. Zienkiewicz,et al.  Reduced integration technique in general analysis of plates and shells , 1971 .

[42]  Thomas J. R. Hughes,et al.  A new finite element formulation for computational fluid dynamics: III. The generalized streamline operator for multidimensional advective-diffusive systems , 1986 .

[43]  Karan S. Surana,et al.  A space–time coupled p‐version least‐squares finite element formulation for unsteady fluid dynamics problems , 1994 .

[44]  Thomas J. R. Hughes,et al.  A simple scheme for developing ‘upwind’ finite elements , 1978 .

[45]  O. C. Zienkiewicz,et al.  Finite element methods for second order differential equations with significant first derivatives , 1976 .

[46]  C. Fletcher.  A primitive variable finite element formulation for inviscid, compressible flow , 1979 .

[47]  Graham F. Carey,et al.  Convergence studies of least-squares finite elements for first-order systems , 1989 .

[48]  B. C. Bell,et al.  p-version least squares finite element formulation for two-dimensional, incompressible, non-Newtonian isothermal and non-isothermal fluid flow , 1994 .

[49]  Bruce M. Irons,et al.  A conforming quartic triangular element for plate bending , 1969 .

[50]  O. C. Zienkiewicz,et al.  A note on upwinding and anisotropic balancing dissipation in finite element approximations to convective diffusion problems , 1980 .

[51]  P. P. Lynn,et al.  Use of the least squares criterion in the finite element formulation , 1973 .