Multiscale phenomena: Green's functions, the Dirichlet-to-Neumann formulation, subgrid scale models, bubbles and the origins of stabilized methods

Finite element formulations based on stabilized bilinear and linear equal-order-interpolation velocity-pressure elements are presented for computation of steady and unsteady incompressible flows. The stabilization procedure involves a slightly modified Galerkin/least-squares formulation of the steady-state equations. The pressure field is interpolated by continuous functions for both the quadrilateral and triangular elements used. These elements are employed in conjunction with the one-step and multi-step time integration of the Navier-Stokes equations. The three test cases chosen for the performance evaluation of these formulations are the standing vortex problem, the lid-driven cavity flow at Reynolds number 400, and flow past a cylinder at Reynolds number 100.

Incompressible flow computations with stabilized bilinear and linear equal-order-interpolation velocity-pressure elements

A new strategy based on the stabilized space-time finite element formulation is proposed for computations involving moving boundaries and interfaces. In the deforming-spatial-domain/space-time (DSD/ST) procedure the variational formulation of a problem is written over its space-time domain, and therefore the deformation of the spatial domain with respect to time is taken into account automatically. Because the space-time mesh is generated over the space-time domain of the problem, within each time step, the boundary (or interface) nodes move with the boundary (or interface). Whether the motion of the boundary is specified or not, the strategy is nearly the same. If the motion of the boundary is unknown, then the boundary nodes move as defined by the other unknowns at the boundary (such as the velocity or the displacement). At the end of each time step a new spatial mesh covers the new spatial domain. For computational feasibility, the finite element interpolation functions are chosen to be discontinuous in time, and the fully discretized equations are solved one space-time slab at a time.

A new strategy for finite element computations involving moving boundaries and interfaces—the deforming-spatial-domain/space-time procedure. I: The concept and the preliminary numerical tests

Publisher Summary This chapter discusses stabilized finite element formulations for incompressible flow computations. Finite element computation of incompressible flows involve two main sources of potential numerical instabilities associated with the Galerkin formulation of a problem. The stabilization techniques that are reviewed more extensively than others are the Galerkin/ least-squares (GLS), streamline-upwind/ Petrov–Galerkin (SUPG), and pressure-stabilizing/Petrov–Galerkin (PSPG) formulations. The SUPG stabilization for incompressible flows is achieved by adding to the Galerkin formulation a series of terms, each in the form of an integral over a different element. These integrals involve the product of the residual of the momentum equation and the advective operator acting on the test function. The natural boundary conditions are the conditions on the stress components, and these are the conditions assumed to be imposed at the remaining part of the boundary. The interpolation functions used for velocity and pressure are piecewise bilinear in space and piecewise linear in time. These computations involve no global coefficient matrices, and therefore need substantially less computer memory and time compared to noniterative solution of the fully discrete equations. It is suggested that for two-liquid flows, the solution and variational function spaces for pressure should include the functions that are discontinuous across the interface.

Gregory M. Hulbert

Papers

A new finite element formulation for computational fluid dynamics: VIII. The galerkin/least-squares method for advective-diffusive equations