Current progress in solving the time-dependent, incompressible Navier-Stokes equations in three-dimensions by (almost) the FEM
We have been pursuing the development of numerical methods which, while originally derived via the conventional Galerkin Finite Element Method (GFEM), have been significantly modified in the interest of cost-effectiveness to the point where we now refer to them as a useful blend of finite element and finite difference methods. Starting with the simplest 3-D isoparametric element and the simplest time integration method (explicit Euler), the principal modifications of GFEM are that the mass matrix is diagonalized via lumping and all Galerkin integrals are evaluated approximately using one-point quadrature. Additional gains in computational speed have been achieved by two modifications of the time integration method: i.e. the introduction of a procedure called subcycling, which permits less frequent updates of the pressure relative to the stability-limited processes of advection and diffusion, and the often expensively small time-step restriction required to assure numerical stability has been eased via the introduction of an additional anisotropic viscosity. The current status of this hybrid method is summarized, and sample numerical results are presented.
- Research Organization:
- Lawrence Livermore National Lab., CA (USA)
- DOE Contract Number:
- W-7405-ENG-48
- OSTI ID:
- 5257180
- Report Number(s):
- UCRL-87445; CONF-820631-1; ON: DE82012259
- Resource Relation:
- Conference: 4. international conference on finite elements in water resources, Hanover, F.R. Germany, 21 Jun 1982; Other Information: Portions of document are illegible
- Country of Publication:
- United States
- Language:
- English
Similar Records
Application of a modified finite element method to the time-dependent thermal convection of a liquid metal
Solution of the multi-dimensional, incompressible Navier-Stokes equations using low-order finite elements and one-point quadrature