Incompressible computational fluid dynamics and the continuity constraint method for the three-dimensional Navier-Stokes equations
- Oak Ridge National Lab., TN (United States). Computational Physics and Engineering Div.
- Univ. of Tennessee, Knoxville, TN (United States). Dept. of Engineering Science and Mechanics
As computational fluid dynamics (CFD) continues to mature, algorithms are required to exploit the most recent advances in approximation theory, numerical mathematics, computing architectures, and hardware. Meeting this requirement is particularly challenging in incompressible fluid mechanics. One key to accomplishing this goal is recognition of the dual role assumed by the pressure, i.e., a mechanism for instantaneously enforcing conservation of mass and a force in the mechanical balance law for conservation of momentum. Proving this assertion has motivated the development of a new, primitive-variable, incompressible, CFD algorithm called the continuity constraint method (CCM). The theoretical basis for the CCM consists of a finite-element spatial semidiscretization of a Galerkin weak statement, equal-order interpolation for all state variables, a {theta}-implicit time-integration scheme, and a quasi-Newton iterative procedure extended by a Taylor weak statement formulation for dispersion error control. This paper presents: the formulation of the unsteady evolution of the divergence error; an investigation of the role of nonsmoothness in the discretized continuity-constraint function; the development of a uniformly H{sup 1} Galerkin weak statement for the Reynolds-averaged Navier-Stokes pressure Poisson equation; and a derivation of physically and numerically well-posed boundary conditions. In contrast to the general family of pressure-relaxation incompressible CFD algorithms, in the CCM, the mathematically smooth and physically motivated genuine pressure is an underlying replacement for the nonsmooth continuity-constraint function to control inherent dispersive-error mechanisms. The genuine pressure is calculated by the diagnostic pressure Poisson equation, evaluated using the verified solenoidal velocity field. This produces a clear view of the distinct boundary conditions required for the continuity-constraint function and genuine pressure.
- Sponsoring Organization:
- USDOE
- OSTI ID:
- 242677
- Journal Information:
- Numerical Heat Transfer. Part B, Fundamentals, Vol. 29, Issue 2; Other Information: PBD: Mar 1996
- Country of Publication:
- United States
- Language:
- English
Similar Records
NSFnets (Navier-Stokes flow nets): Physics-informed neural networks for the incompressible Navier-Stokes equations
Accuracy, convergence and stability of TWS finite element algorithms for the incompressible Navier-Stokes equations. [Taylor Weak Statement]