A finite element variational multiscale method for incompressible flows based on two local gauss integrations

In this article, we present a finite element variational multiscale (VMS) method for incompressible flows based on two local Gauss integrations, and compare it with common VMS method which is defined by a low order finite element space L_h$L_h$ on the same grid as X_h$X_h$ for the velocity deformation tensor and a stabilization parameter α$\alpha$. The best algorithmic feature of our method is using two local Gauss integrations to replace projection operator. We theoretically discuss the relationship between our method and common VMS method for the Taylor–Hood elements, and show that the nonlinear system derived from our method by finite element discretization is much smaller than that of common VMS method computationally.     

Fictitious Domain approach with hp-finite element approximation for incompressible fluid flow

We consider the application of Fictitious Domain approach combined with least squares spectral elements for the numerical solution of fluid dynamic incompressible equations. Fictitious Domain methods allow problems formulated on a complicated shaped domain Ω$\Omega$ to be solved on a simpler domain Π$\Pi$ containing Ω$\Omega$. Least Squares Spectral Element Method has been used to develop the discrete model, as this scheme combines the generality of finite element methods with the accuracy of spectral methods. Moreover the least squares methods have theoretical and computational advantages in the algorithmic design and implementation. This paper presents the formulation and validation of the Fictitious Domain Least Squares Spectral Element approach for the steady incompressible Navier–Stokes equations. The convergence of the approximated solution is verified solving two-dimensional benchmark problems, demonstrating the predictive capability of the proposed formulation.     

A massively parallel fractional step solver for incompressible flows

This paper presents a parallel implementation of fractional solvers for the incompressible Navier–Stokes equations using an algebraic approach. Under this framework, predictor–corrector and incremental projection schemes are seen as sub-classes of the same class, making apparent its differences and similarities. An additional advantage of this approach is to set a common basis for a parallelization strategy, which can be extended to other split techniques or to compressible flows.     

Numerical Analysis and Comparison of Four Stabilized Finite Element Methods for the Steady Micropolar Equations

In this paper, four stabilized methods based on the lowest equal-order finite element pair for the steady micropolar Navier–Stokes equations (MNSE) are presented, which are penalty, regular, multiscale enrichment, and local Gauss integration methods. A priori properties, existence, uniqueness, stability, and error estimation based on Fem approximation of all the methods are proven for the physical variables. Finally, some numerical examples are displayed to show the numerical characteristics of these methods.     

Explicit one-step time discretizations for discontinuous Galerkin and finite volume schemes based on local predictors

We consider a family of explicit one-step time discretizations for finite volume and discontinuous Galerkin schemes, which is based on a predictor-corrector formulation. The predictor remains local taking into account the time evolution of the data only within the grid cell. Based on a space–time Taylor expansion, this idea is already inherent in the MUSCL finite volume scheme to get second order accuracy in time and was generalized in the context of higher order ENO finite volume schemes. We interpret the space–time Taylor expansion used in this approach as a local predictor and conclude that other space–time approximate solutions of the local Cauchy problem in the grid cell may be applied. Three possibilities are considered in this paper: (1) the classical space–time Taylor expansion, in which time derivatives are obtained from known space-derivatives by the Cauchy–Kovalewsky procedure; (2) a local continuous extension Runge–Kutta scheme and (3) a local space–time Galerkin predictor with a version suitable for stiff source terms. The advantage of the predictor–corrector formulation is that the time evolution is done in one step which establishes optimal locality during the whole time step. This time discretization scheme can be used within all schemes which are based on a piecewise continuous approximation as finite volume schemes, discontinuous Galerkin schemes or the recently proposed reconstructed discontinuous Galerkin or P_NP_M schemes. The implementation of these approaches is described, advantages and disadvantages of different predictors are discussed and numerical results are shown.     

Multigrid method based on the transformation-free HOC scheme on nonuniform grids for 2D convection diffusion problems

In this paper, a multigrid method based on the high order compact (HOC) difference scheme on nonuniform grids, which has been proposed by Kalita et al. [J.C. Kalita, A.K. Dass, D.C. Dalal, A transformation-free HOC scheme for steady convection–diffusion on non-uniform grids, Int. J. Numer. Methods Fluids 44 (2004) 33–53], is proposed to solve the two-dimensional (2D) convection diffusion equation. The HOC scheme is not involved in any grid transformation to map the nonuniform grids to uniform grids, consequently, the multigrid method is brand-new for solving the discrete system arising from the difference equation on nonuniform grids. The corresponding multigrid projection and interpolation operators are constructed by the area ratio. Some boundary layer and local singularity problems are used to demonstrate the superiority of the present method. Numerical results show that the multigrid method with the HOC scheme on nonuniform grids almost gets as equally efficient convergence rate as on uniform grids and the computed solution on nonuniform grids retains fourth order accuracy while on uniform grids just gets very poor solution for very steep boundary layer or high local singularity problems. The present method is also applied to solve the 2D incompressible Navier–Stokes equations using the stream function–vorticity formulation and the numerical solutions of the lid-driven cavity flow problem are obtained and compared with solutions available in the literature.