skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A simple finite element method for the Stokes equations

Abstract

The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in primal velocity-pressure formulation and is so simple such that both velocity and pressure are approximated by piecewise constant functions. Implementation issues as well as error analysis are investigated. A basis for a divergence free subspace of the velocity field is constructed so that the original saddle point problem can be reduced to a symmetric and positive definite system with much fewer unknowns. The numerical experiments indicate that the method is accurate.

Authors:
ORCiD logo [1];  [2]
  1. Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
  2. Univ. of Arkansas, Little Rock, AR (United States)
Publication Date:
Research Org.:
Oak Ridge National Lab. (ORNL), Oak Ridge, TN (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1362244
Grant/Contract Number:
AC05-00OR22725
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Advances in Computational Mathematics
Additional Journal Information:
Journal Volume: 43; Journal Issue: 6; Journal ID: ISSN 1019-7168
Publisher:
Springer
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING

Citation Formats

Mu, Lin, and Ye, Xiu. A simple finite element method for the Stokes equations. United States: N. p., 2017. Web. doi:10.1007/s10444-017-9526-z.
Mu, Lin, & Ye, Xiu. A simple finite element method for the Stokes equations. United States. doi:10.1007/s10444-017-9526-z.
Mu, Lin, and Ye, Xiu. Tue . "A simple finite element method for the Stokes equations". United States. doi:10.1007/s10444-017-9526-z. https://www.osti.gov/servlets/purl/1362244.
@article{osti_1362244,
title = {A simple finite element method for the Stokes equations},
author = {Mu, Lin and Ye, Xiu},
abstractNote = {The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in primal velocity-pressure formulation and is so simple such that both velocity and pressure are approximated by piecewise constant functions. Implementation issues as well as error analysis are investigated. A basis for a divergence free subspace of the velocity field is constructed so that the original saddle point problem can be reduced to a symmetric and positive definite system with much fewer unknowns. The numerical experiments indicate that the method is accurate.},
doi = {10.1007/s10444-017-9526-z},
journal = {Advances in Computational Mathematics},
number = 6,
volume = 43,
place = {United States},
year = {Tue Mar 21 00:00:00 EDT 2017},
month = {Tue Mar 21 00:00:00 EDT 2017}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Save / Share:
  • This paper presents an alternative formulation to the primitive variable form of the Navier-Stokes equations. The approach is based on the use of the velocity and vorticity variables as a logical extension of the stream function-vorticity method which is quite effective for two-dimensional flows. Second-order equations are obtained for the variables and discretization is through a weak-Galerkin finite element method, followed by Newton linearization. The velocity and vorticity components are solved simultaneously by a direct solver. The scheme is demonstrated for two-dimensional and three-dimensional incompressible and subsonic internal enclosed flow problems. 31 refs., 16 figs., 1 tab.
  • This paper deals with a high-order accurate discontinuous finite element method for the numerical solution of the compressible Navier-Stokes equations. We extend a discontinuous finite element discretization originally considered for hyperbolic systems such as the Euler equations to the case of the Navier-Stokes equations by treating the viscous terms with a mixed formulation. The method combines two key ideas which are at the basis of the finite volume and of the finite element method, the physics of wave propagation being accounted for by means of Riemann problems and accuracy being obtained by means of high-order polynomial approximations within elements. Asmore » a consequence the method is ideally suited to compute high-order accurate solution of the Navier-Stokes equations on unstructured grids. The performance of the proposed method is illustrated by computing the compressible viscous flow on a flat plate and around a NACA0012 airfoil for several flow regimes using constant, linear, quadratic, and cubic elements. 23 refs., 24 figs., 3 tabs.« less
  • No abstract prepared.
  • The main purpose of this paper is to develop a finite analytic (FA) numerical solution for unsteady two-dimensional Navier-Stokes equations. The FA method utilizes the analytic solution in a small local element to formulate the algebraic representation of partial differential equations. In this study the combination of linear and exponential functions that satisfy the governing equation is adopted as the boundary function, thereby improving the accuracy of the finite analytic solution. Two flows, one a starting cavity flow and the other a vortex shedding flow behind a rectangular block, are solved by the FA method. The starting square cavity flowmore » is solved for Reynolds numbers of 400, 1000, and 2000 to show the accuracy and stability of the FA solution. The FA solution for flow over a rectangular block (H x H/4) predicts the Strouhal number for Reynolds numbers of 100 and 500 to be 0.156 and 0.125. Details of the flow patterns are given. In addition to streamlines and vorticity distribution, rest-streamlines are given to illustrate the vortex motion downstream of the block.« less