Direct Simulation of Low-Re Flow around a Square Cylinder by Numerical Manifold Method for Navier-Stokes Equations

Abstract

Numerical manifold method was applied to directly solve Navier-Stokes (N-S) equations for incompressible viscous flow in this paper, and numerical manifold schemes for N-S equations coupled velocity and pressure were derived based on Galerkin weighted residuals method as well. Mixed cover with linear polynomial function for velocity and constant function for pressure was employed in finite element cover system. As an application, mixed cover 4-node rectangular manifold element has been used to simulate the incompressible viscous flow around a square cylinder in a channel. Numerical tests illustrate that NMM is an effective and high-order accurate numerical method for incompressible viscous flow N-S equations.

1. Introduction

In computational fluid dynamics (CFD), Navier-Stokes (N-S) equations for incompressible viscous flow can be solved by several numerical methods generally, such as finite difference method (FDM), finite element method (FEM), and finite volume method (FVM) [1]. The applications of FDM are affected for unsuitability to complicated structure flow field and inconservation of field variables [2]. The FEM schemes for N-S equations are very complicated, and the computational efficiency of nonstandard FEMs is low [3]. FVM is an effective numerical method for solutions of fluid flow, but exact treatments of the moving boundary are difficult and complicated [4]. Furthermore, fractional step algorithms are used to solve velocity and pressure variables unintegratedly in all these numerical methods, and the continuity equation and momentums equations are satisfied, respectively [5]. All the defects of these numerical methods will affect the computational efficiency and accuracy.

Numerical manifold method (NMM) also known as manifold method or finite cover method (FCM) is a generalized numerical method proposed by Shi in the early 1990s [6, 7]. The method performs numerical computation with finite element cover system, which is composed of two independent cover grids: mathematical cover grid and physical cover grid. Mathematical covers define the accuracy of approximate solution, and physical covers determine the solution domain. High accurate numerical manifold schemes for mathematical physics equations can be constructed by adopting different cover functions and cover-weighted functions for physical variables in finite element cover system [8]. The method has been successfully applied in some complicated engineering problems, such as numerical simulation of crack initiation and propagation, damage evolvement of fractured rock masses, coupling vibration, and potential flow problem, and the advantages have been demonstrated in these fields [9–12]. As to fluid flow problems, direct-numerical-solution coupled velocity and pressure variables can be implemented by adopting manifold schemes with mixed cover for N-S equations. The authors have applied numerical manifold method into analyses of unsteady incompressible viscous flow, and the results have illustrated the validity of NMM [13].

Flow around a square cylinder is a typical model to validate the performance of numerical methods for solution of incompressible viscous N-S equations. The flow structure has been investigated experimentally and numerically. Experimental investigations have shown that the flow characteristic is different at different Reynolds numbers [14]. FVM, FDM, and others numerical methods have been applied to simulate the flow around a square cylinder, the detail results have been obtained and met the experimental results [15–17]. The experimental and numerical results have been employed to validate the numerical methods in CFD as benchmarks.

In this paper, numerical manifold schemes of direct solutions coupled velocity and pressure for N-S equations were constructed and applied to analyze incompressible viscous flow around a square cylinder in a parallel channel. The validations of numerical schemes to steady and unsteady flow were completed.

2. Numerical Manifold Schemes for N-S Equations of Incompressible Viscous Flow

2.1. Galerkin Integration Expressions of N-S Equations

For numerical solution of incompressible viscous flow, the integration expressions for the continuity equation and N-S equations can be obtained by Galerkin weighted residual method, and the weak solution form can be stated as_Ω𝑣_𝑗𝜕𝜕𝑥_𝑗𝛿𝑝dΩ=_Γ1𝑣_𝑛𝛿𝑝dΓ,_Ω𝜌𝜕𝑣_𝑖𝜕𝑡+𝑣_𝑗𝜕𝑣_𝑖𝜕𝑥_𝑗𝛿𝑣_𝑖+−𝑝𝛿_𝑖𝑗+𝜇𝜕𝑣_𝑖𝜕𝑥_𝑗+𝜕𝑣_𝑗𝜕𝑥_𝑖𝜕𝜕𝑥_𝑗𝛿𝑣_𝑖=dΩ_Γ2𝑝_𝑖𝛿𝑣_𝑖dΓ+_Ω𝜌𝑓_𝑖𝛿𝑣_𝑖dΩ, (2.1a)(2.1b)where 𝑡 denotes the time variable; 𝜌 is the fluid density; 𝜇 is the dynamic viscosity; 𝑣_𝑗 is the velocity in 𝑥_𝑗 direction; 𝑝 is the pressure; 𝑓_𝑖 is the body force per unit volume in 𝑥_𝑖 direction; 𝛿𝑝 is referred to the weight function for continuity equation and 𝛿𝑣_𝑖 for momentum equation; Ω is the domain of solution; 𝑣_𝑛 is the known normal velocity on boundary Γ₁; 𝑝_𝑖 is the known pressure components in 𝑥_𝑖 direction on boundary Γ₂.

2.2. Numerical Manifold Schemes for N-S Equations

2.2.1. Element Velocity and Pressure Cover Functions

When NMM is applied into solution of N-S equations, suggested that there are 𝑀 manifold elements in finite element cover system and each element (e) has 𝑛 physical covers, and the cover weight functions are 𝑁_𝑗(𝑥,𝑦)(𝑗=1,2,…,𝑛), which meet the requirement of partition of unity ∑^𝑛_𝑗=1𝑁_𝑗(𝑥,𝑦)=1(𝑥,𝑦∈𝑒), then the element global velocity and pressure cover functions 𝑣_𝑗(𝑥,𝑦) and 𝑝(𝑥,𝑦) can be given by weighted summation of the cover functions of every physical cover as the follows: 𝑣_𝑖𝑥,𝑦=𝑁_𝑗𝑥,𝑦𝑣_𝑖𝑗,𝑝𝑥,𝑦𝑥,𝑦=𝑁_𝑗𝑥,𝑦𝑝_𝑗𝑥,𝑦, (2.2) where 𝑖=1,2 for two-dimensional (2D) flow problem; 𝑣_𝑖𝑗(𝑥,𝑦) are the velocity component cover functions for physical cover 𝑈_𝑗 of element (𝑒) and 𝑝_𝑗(𝑥,𝑦) is the pressure cover function. If polynomial functions are taken as the cover functions of every physical cover, the cover functions can be stated in the form as 𝑣_𝑖𝑗𝑥,𝑦=𝑓^𝑣_𝑘⋅𝐷^𝑣_𝑖𝑗𝑘,𝑝_𝑗𝑥,𝑦=𝑓^𝑝_𝑙⋅𝐷^𝑝_𝑗𝑙, (2.3) where 𝑓^𝑣_𝑘,𝑓^𝑝_𝑙 denote the basic series for velocity and pressure cover functions of physical cover 𝑈_𝑗; 𝐷^𝑣_𝑖𝑗𝑘,𝐷^𝑝_𝑗𝑙 are the cover DOF variables, and the DOF numbers are 𝑚_𝑣(𝑘=1,2,…,𝑚_𝑣) and 𝑚_𝑝(𝑙=1,2,…,𝑚_𝑝).

By substituting (2.3) into (2.2), the global cover functions for velocity and pressure in element (e) can be written as𝑣_𝑖𝑥,𝑦=𝑁_𝑗𝑥,𝑦𝑣_𝑖𝑗𝑥,𝑦=𝑁_𝑗𝑥,𝑦𝑓^𝑣_𝑘𝐷^𝑣_𝑖𝑗𝑘=𝑇^𝑣_𝑗𝑘𝐷^𝑣_𝑖𝑗𝑘,𝑝𝑥,𝑦=𝑁_𝑗𝑥,𝑦𝑝_𝑗𝑥,𝑦=𝑁_𝑗𝑥,𝑦𝑓^𝑝_𝑙𝐷^𝑝_𝑗𝑙=𝑇^𝑝_𝑗𝑙𝐷^𝑝_𝑗𝑙, (2.4a)(2.4b)where 𝑇^𝑣_𝑗𝑘 is the element velocity cover basic function; 𝑇^𝑝_𝑗𝑙 is the element pressure cover basic function. The element cover basic functions consist of the basic series of cover functions and cover weight functions, and the element DOF variables are composed of all the cover DOF variables.

In Galerkin integration equations (2.1a) and (2.1b), the weight functions for continuity equation and N-S equations can adopt the element basic functions of velocity components and pressure accordingly, that is, 𝛿𝑣_𝑖=𝑇^𝑣_𝑗𝑘;𝛿𝑝=𝑇^𝑝_𝑗𝑙. (2.5)

The partial derivatives of element cover functions can be obtained from (2.4a) and (2.4b) as the follows: 𝑣_𝑖,𝛼=𝑇^𝑣_{𝑗𝑘,𝛼}𝐷^𝑣_𝑖𝑗𝑘;𝑝_,𝛼=𝑇^𝑝_{𝑗𝑙,𝛼}𝐷^𝑝_𝑗𝑙;𝑣_𝑖,𝑡=𝑇^𝑣_𝑗𝑘𝐷^𝑣_{𝑖𝑗𝑘,𝑡}=𝑇^𝑣_𝑗𝑘𝐷̇^𝑣_𝑖𝑗𝑘, (2.6) where 𝛼=1,2 for 2D problem; 𝑇^𝑣_{𝑗𝑘,𝛼}, 𝑇^𝑝_{𝑗𝑙,𝛼} are the partial derivatives of element velocity and pressure cover basic functions with respect to 𝑥_𝛼, respectively; 𝐷̇^𝑣_𝑖𝑗𝑘 are the partial derivatives of element velocity DOF variables with respect to t.

2.2.2. Element Manifold Schemes

For each element (𝑒) in finite element cover system, element manifold equations can be derived by introducing (2.4a), (2.4b), (2.5), and (2.6) into (2.2), and the equations can be stated as _Ω𝑒𝑇^𝑣_𝛽𝐼𝐷^𝑣_𝛼𝛽𝐼𝑇^𝑝_{𝑗𝑙,𝛼}dΩ=_Γ𝑒1𝑣_𝑛𝑇^𝑝_𝑗𝑙dΓ,_Ω𝑒𝜌𝑇^𝑣_𝛽𝐼𝐷̇^𝑣_𝑖𝛽𝐼𝑇^𝑣_𝑗𝑘+𝜌𝑇^𝑣_𝛽𝐼𝐷^𝑣_𝛼𝛽𝐼𝑇^𝑣_{𝛾𝐽,𝛼}𝐷^𝑣_𝑖𝛾𝐽𝑇^𝑣_𝑗𝑘−𝑇^𝑝_𝛽𝐾𝐷^𝑝_𝛽𝐾𝛿_𝑖𝛼𝑇^𝑣_{𝑗𝑘,𝛼}+𝜇𝑇^𝑣_{𝛾𝐽,𝛼}𝐷^𝑣_𝑖𝛾𝐽+𝑇^𝑣_{𝛾𝐽,𝑖}𝐷^𝑣_𝛼𝛾𝐽𝑇^𝑣_{𝑗𝑘,𝛼}=dΩ_Γ𝑒2𝑝_𝑖𝑇^𝑣_𝑗𝑘dΓ+_Ω𝑒𝜌𝑓_𝑖𝑇^𝑣_𝑗𝑘dΩ, (2.7) that is,𝐺^𝑒_{𝑗𝑙𝛼𝛽𝐼}𝐷^𝑣_𝛼𝛽𝐼=𝐻^𝑒_𝑗𝑙,𝐴^𝑒_{𝑗𝑘𝛾𝐽}𝐷̇^𝑣_𝑖𝛾𝐽+𝐵^𝑒_{𝑗𝑘𝛼𝛽𝛾𝐼𝐽}𝐷^𝑣_𝛼𝛽𝐼𝐷^𝑣_𝑖𝛾𝐽+𝐶^𝑒_{𝑖𝑗𝑘𝛽𝐾}𝐷^𝑝_𝛽𝐾+𝐸^𝑒_{𝑗𝑘𝑖𝛼𝛾𝐽}𝐷^𝑣_𝛼𝛾𝐽=𝐹^𝑒_𝑖𝑗𝑘, (2.8a)(2.8b)where 𝛽,𝛾=1,2,…,𝑛;𝐼,𝐽=1,2,…,𝑚_𝑣,𝐾=1,2,…,𝑚_𝑝; and the coefficients for element manifold equations are 𝐴^𝑒_{𝑗𝑘𝛾𝐽}=_Ω𝑒𝜌𝑇^𝑣_𝛾𝐽𝑇^𝑣_𝑗𝑘dΩ,𝐵^𝑒_{𝑗𝑘𝛼𝛽𝛾𝐼𝐽}=_Ω𝑒𝜌𝑇^𝑣_𝛽𝐼𝑇^𝑣_{𝛾𝐽,𝛼}𝑇^𝑣_𝑗𝑘𝐶dΩ,_{𝑒𝑖𝑗𝑘𝛽𝐾}=−_Ω𝑒𝑇^𝑝_𝛽𝐾𝛿_𝑖𝛼𝑇^𝑣_{𝑗𝑘,𝛼}dΩ,𝐸^𝑒_{𝑗𝑘𝑖𝛼𝛾𝐽}=_Ω𝑒𝜇𝑇^𝑣_{𝑗𝑘,𝑠}𝛿_𝑖𝛼𝑇^𝑣_{𝛾𝐽,𝑠}+𝛿_𝑠𝛼𝑇^𝑣_{𝛾𝐽,𝑖}𝐹dΩ,_{𝑒𝑖𝑗𝑘}=_Γ𝑒2𝑝_𝑖𝑇^𝑣_𝑗𝑘dΓ+_Ω𝑒𝜌𝑓_𝑖𝑇^𝑣_𝑗𝑘dΩ,𝐺^𝑒_{𝑗𝑙𝛼𝛽𝐼}=_Ω𝑒𝑇^𝑣_𝛽𝐼𝑇^𝑝_{𝑗𝑙,𝛼}dΩ,𝐻^𝑒_𝑗𝑙=_Γ𝑒1𝑣_𝑛𝑇^𝑝_𝑗𝑙dΓ. (2.9)

Equations (2.8a) and (2.8b) are the nonlinear algebra equations with respect to coupled velocity and pressure DOF variables and first-order linear ordinary differential equations with respect to time.

3. Some Key Techniques

3.1. Mixed Cover Functions for Velocity and Pressure

In FEM, when velocity and pressure field are discretized in the same way, the discrete element can not be ensured to satisfy the Ladyzanskya Babuska Brezzi stability condition (LBB condition), and these kinds of element cannot be applied to directly solve N-S equations coupled velocity and pressure for the spatial oscillation of pressure field. Mixed elements can be built to meet the LBB condition by increasing velocity interpolation node [18], but the elements are inconvenient to apply in practice for the complicated calculation process and difficult element discretization.

In finite element cover system of NMM, the element cover functions are composed of the weight functions and the cover functions of physical covers. Generally, the weight functions will adopt the element shape functions as in FEM, which are usually defined by the element shape and nodes. The cover functions of every physical cover for different physical variables can employ different order functions according to the solving physical equations. When NMM is applied to solve N-S equations for incompressible viscous flow, the velocity and pressure field can be discretized in the same way and the weight functions will employ the same interpolation functions for element velocity and pressure variables. But, the cover functions of every element physical cover can apply different order functions for velocity and pressure variables, so it will form a manifold element with mixed cover functions for velocity and pressure variables, which can meet the requirements of different order approximate functions for velocity and pressure in Galerkin integration expressions (2.1a) and (2.1b) and satisfy the LBB condition for increasing the DOF numbers of velocity variables. The mixed cover manifold elements can be applied in direct solution of incompressible viscous N-S-equations-coupled velocity and pressure variables.

In theory, velocity components and pressure covers can employ very high-order functions, but it will cause a very complicated calculation process, so low-order polynomial functions are favorable. If pressure cover function of every physical cover is defined as constant function, velocity cover functions can adopt linear polynomial functions, and then the basic series in (2.3) can be given by the following localized series for 2D problem: 𝑓^𝑣_𝑘=1𝜂₁𝑥−𝑥_𝑐𝜂₂𝑦−𝑦_𝑐;𝑓^𝑝_𝑙=[1], (3.1) where (𝑥,𝑦) is the coordinate variables; (𝑥_𝑐,𝑦_𝑐) is the center coordinates of physical cover; 𝜂₁, 𝜂₂ are the coefficients of cover functions.

3.2. Solution of Manifold Equations for Steady N-S Equations

As to 2D steady N-S equations, element manifold equations (2.8a) and (2.8b) can be reformatted as 𝐵^𝑒_{𝑗𝑘𝛼𝛽𝛾𝐼𝐽}𝐷^𝑣_𝛼𝛽𝐼𝐷^𝑣_1𝛾𝐽+𝐶^𝑒_{1𝑗𝑘𝛽𝐾}𝐷^𝑝_𝛽𝐾+𝐸^𝑒_{𝑗𝑘1𝛼𝛾𝐽}𝐷^𝑣_𝛼𝛾𝐽=𝐹^𝑒_1𝑗𝑘,𝐵^𝑒_{𝑗𝑘𝛼𝛽𝛾𝐼𝐽}𝐷^𝑣_𝛼𝛽𝐼𝐷^𝑣_2𝛾𝐽+𝐶^𝑒_{2𝑗𝑘𝛽𝐾}𝐷^𝑝_𝛽𝐾+𝐸^𝑒_{𝑗𝑘2𝛼𝛾𝐽}𝐷^𝑣_𝛼𝛾𝐽=𝐹^𝑒_2𝑗𝑘,𝐺^𝑒_{𝑗𝑙𝛼𝛽𝐼}𝐷^𝑣_𝛼𝛽𝐼=𝐻^𝑒_𝑗𝑙, (3.2) where 𝑗,𝑘,𝑙,𝛼,𝛽,𝛾,𝐼,𝐽 are similar to (2.8a) and (2.8b); 𝐷_1𝛾𝐽, 𝐷_2𝛾𝐽, and 𝐷_𝛽𝐾 are the unknown variables. The equations can be written in the matrix form as𝐵^𝑒_{𝑗𝑘1𝛽𝛾𝐼𝐽}𝐷^𝑣_1𝛽𝐼+𝐸^𝑒_{𝑗𝑘11𝛾𝐽}𝐵^𝑒_{𝑗𝑘2𝛽𝛾𝐼𝐽}𝐷^𝑣_1𝛽𝐼+𝐸^𝑒_{𝑗𝑘12𝛾𝐽}𝐶^𝑒_{1𝑗𝑘𝛽𝐾}𝐵^𝑒_{𝑗𝑘1𝛽𝛾𝐼𝐽}𝐷^𝑣_2𝛽𝐼+𝐸^𝑒_{𝑗𝑘21𝛾𝐽}𝐵^𝑒_{𝑗𝑘2𝛽𝛾𝐼𝐽}𝐷^𝑣_2𝛽𝐼+𝐸^𝑒_{𝑗𝑘22𝛾𝐽}𝐶^𝑒_{2𝑗𝑘𝛽𝐾}𝐺^𝑒_{𝑗𝑙1𝛾𝐽}𝐺^𝑒_{𝑗𝑙2𝛾𝐽}0𝐷^𝑣_1𝛾𝐽𝐷^𝑣_2𝛾𝐽𝐷^𝑝_𝛽𝐾=𝐹^𝑒_1𝑗𝑘𝐹^𝑒_2𝑗𝑘𝐻^𝑒_𝑗𝑙, (3.3a) that is, 𝐾^𝑒𝐷^𝑒=𝐹^𝑒, (3.3b)where [𝐾^𝑒] is the element flow matrix with respect to unknown velocity DOF variables; {𝐷^𝑒} is the element DOF variables; {𝐹^𝑒} is the element generalized loads.

The manifold equations (3.3a) and (3.3b) are nonlinear algebra equations and can be solved by two methods of directly linearizing alternate iteration and Newton-Raphson iterative in general.

3.3. Solution of Manifold Equations for Steady N-S Equations

As to 2D unsteady N-S equations, the first-order linear ordinary differential equation (2.8b) with respect to time need be translated into linear algebra equations by time difference replacing the derivative. At any discrete time step 𝑡_𝑖, the nonlinear manifold equations can be solved by implicit iterative methods or explicit methods. If implicit iterative methods are used, the derivative of velocity DOF variables for every physical cover with respect to time in (2.8b) will be discretized by time difference as 𝐷̇^𝑣_𝑖𝑗𝑘=𝐷^𝑣(𝑖+1)_𝑖𝑗𝑘−𝐷^𝑣(𝑖)_𝑖𝑗𝑘,Δ𝑡 (3.4) where Δ𝑡 is time step; 𝐷^𝑣(𝑖)_𝑖𝑗𝑘 and 𝐷^𝑣(𝑖+1)_𝑖𝑗𝑘 are the velocity DOF variables at time step 𝑡_𝑖 and 𝑡_𝑖+1, respectively. Then as to every manifold element (e), the element manifold equations (2.8b) can be reformulated as the follows: 𝐴^𝑒_{𝑗𝑘𝛾𝐽}𝐷^𝑣(𝑖+1)_𝑖𝛾𝐽+Δ𝑡𝐵^𝑒_{𝑗𝑘𝛼𝛽𝛾𝐼𝐽}𝐷^𝑣(𝑖+1)_𝛼𝛽𝐼𝐷^𝑣(𝑖+1)_𝑖𝛾𝐽+Δ𝑡𝐶^𝑒_{𝑖𝑗𝑘𝛽𝐾}𝐷^𝑝(𝑖+1)_𝛽𝐾+Δ𝑡𝐸^𝑒_{𝑗𝑘𝑖𝛼𝛾𝐽}𝐷^𝑣(𝑖+1)_𝛼𝛾𝐽=Δ𝑡𝐹^𝑒_𝑖𝑗𝑘+𝐴^𝑒_{𝑗𝑘𝛽𝐼}𝐷^𝑣(𝑖)_𝑖𝑗𝑘. (3.5) Equations (3.5) uniting (2.8a) can be stated in matrix form as 𝐴^𝑒_{𝑗𝑘𝛾𝐽}+Δ𝑡𝐵^𝑒_{𝑗𝑘1𝛽𝛾𝐼𝐽}𝐷^𝑣(𝑖+1)_1𝛽𝐼+𝐸^𝑒_{𝑗𝑘11𝛾𝐽}Δ𝑡𝐵^𝑒_{𝑗𝑘2𝛽𝛾𝐼𝐽}𝐷^𝑣(𝑖+1)_1𝛽𝐼+𝐸^𝑒_{𝑗𝑘12𝛾𝐽}Δ𝑡𝐶^𝑒_{1𝑗𝑘𝛽𝐾}Δ𝑡𝐵^𝑒_{𝑗𝑘1𝛽𝛾𝐼𝐽}𝐷^𝑣(𝑖+1)_2𝛽𝐼+𝐸^𝑒_{𝑗𝑘21𝛾𝐽}𝐴^𝑒_{𝑗𝑘𝛾𝐽}+Δ𝑡𝐵^𝑒_{𝑗𝑘2𝛽𝛾𝐼𝐽}𝐷^𝑣(𝑖+1)_2𝛽𝐼+𝐸^𝑒_{𝑗𝑘22𝛾𝐽}Δ𝑡𝐶^𝑒_{2𝑗𝑘𝛽𝐾}𝐺^𝑒_{𝑗𝑙1𝛾𝐽}𝐺^𝑒_{𝑗𝑙2𝛾𝐽}0⋅𝐷^𝑣(𝑖+1)_1𝛾𝐽𝐷^𝑣(𝑖+1)_2𝛾𝐽𝐷^𝑝(𝑖+1)_𝛽𝐾=Δ𝑡𝐹^𝑒_1𝑗𝑘+𝐴^𝑒_{𝑗𝑘𝛾𝐽}𝐷^𝑣(𝑖)_1𝛾𝐽Δ𝑡𝐹^𝑒_2𝑗𝑘+𝐴^𝑒_{𝑗𝑘𝛾𝐽}𝐷^𝑣(𝑖)_2𝛾𝐽𝐻^𝑒_𝑗𝑙, (3.6) that is, 𝑘^𝑒𝐷^𝑒(𝑖+1)=𝐹^𝑒, (3.7) where {𝐷^𝑝(𝑖+1)_𝛽𝐾} is the pressure cover DOF variables at time step 𝑡_𝑖+1;{𝐷^𝑒(𝑖+1)} is the element unknown cover DOF variables at time step 𝑡_𝑖+1;[𝑘^𝑒] is the element flow matrix with respect to the unknown velocity DOF variables at time step 𝑡_𝑖+1;{𝐹^𝑒} is the element generalized loads with respect to the boundary conditions and the known velocity DOF variables at time step 𝑡_𝑖.

The manifold equations (3.7) for unsteady flow can be solved by iterative methods. In this paper, Newton-Raphson iterative method is also used here for solution of unsteady flow manifold equations.

4. Numerical Example of Flow around a Square Cylinder in a Channel

To illustrate the validity of the present numerical manifold method, low-Re incompressible viscous flow around a square cylinder in a channel and flow past a step are investigated in details through direct solution of N-S equations in dimensionless form in this paper as two numerical examples. In numerical manifold analysis, standard rectangular finite element cover system with 4-node manifold element as show in Figure 1 is used, the mixed cover functions of physical covers with the basic series equations (3.1) are adopted, and the shape functions of rectangular element are employed as the cover weight functions in (2.2), that is, 𝑁₁=1𝑥,𝑦41−𝑥−𝑥₀𝑎1−𝑦−𝑦₀𝑏,𝑁₂=1𝑥,𝑦41+𝑥−𝑥₀𝑎1−𝑦−𝑦₀𝑏,𝑁₃=1𝑥,𝑦41+𝑥−𝑥₀𝑎1+𝑦−𝑦₀𝑏,𝑁₄=1𝑥,𝑦41−𝑥−𝑥₀𝑎1+𝑦−𝑦₀𝑏, (4.1) where (𝑥₀,𝑦₀) are the coordinates of element (e) center point; 2𝑎×2𝑏 are the side lengths.

Figure 1 

4-node rectangular manifold element.

4.1. Analysis for Steady Flow past a Step

In numerical analysis of flow past a step, finite element cover system and boundary conditions are shown in Figure 2(a), the characteristic velocity of the flow field inlet is 𝑢=1 and 𝑣=0, the fluid density is 𝜌=1, and three dynamic viscosities of fluid are 𝜇 = 0.01, 0.005, 0.0025, and so Re numbers are known as Re=100, 200, 400, respectively. All the flow field is meshed into 196 elements and 229 nodes by 0.1×0.1 rectangular grids.

(a) Flow field configuration and boundary conditions

(b) NNM solutions (0.1×0.1 grids)

(c) FVM solutions (0.1×0.1 grids)

(d) FVM solutions (0.05×0.05 grids)

Figure 2 

The streamline patterns and pressure distributions for flow field past a step.

The streamline and pressure distribution in flow field past a step at Re = 200 are shown in Figure 2(b). The streamline patterns obtained from direct numerical solutions have revealed the vortex structure at the step of flow field, and high accurate pressure distributions have been obtained as well. The numerical solutions from FVM are shown in Figures 2(c) and 2(d) comparatively, which obtained from 0.1×0.1 and 0.05×0.05 rectangular grids, respectively. The NMM results from 0.1×0.1 grids are nearly coincident with the FVM results from 0.05×0.05 grids. The comparative analyses illustrate that numerical manifold schemes with mixed cover are stable numerical schemes for incompressible viscous flow, can implement direct numerical solution of N-S equations coupled velocity and pressure variables and can improve the solution accuracy at the same grids compared with FVM.

4.2. Analysis for Steady Flow around a Square Cylinder in a Channel

The flow field configuration and boundary conditions of flow around a square cylinder in a channel are shown in Figure 3. The length of the channel is 𝐿=4, the width is 𝑊=1. The size of square cylinder is 𝑑×𝑑=0.25×0.25, its horizontal axis coincides with the channel’s centerline, and its vertical axis is perpendicular to the channel’s centreline. The velocity condition of the inlet AD is 𝑢/𝑢_av=1.5[1−(1−𝑦/𝑦_𝑚)²](𝑢_av=1is the average flow velocity in the channel, 𝑦_𝑚=0.5 is the channel centerline) and 𝑣=0, the pressure condition of outlet BC is 𝑝=0, and the nonslip boundary conditions are set on other fixed walls, that is, 𝑢=𝑣=0.

Figure 3 

Flow field configuration and boundary conditions of flow around a square cylinder.

The flow field is meshed by three mixed grids, one is that 32×16 uniform element grid is used for the domain from 𝑥=0 to 𝑥=2 and 20×16 element grid for the domain from 𝑥=2 to 𝑥=4, the entire flow field is meshed by 52×16=832 elements and 53×17= 901 nodes; the second one is that 64×32 uniform element grid is used for 𝑥=0∼2 and 40×32 element grid for 𝑥=2∼4, the flow field is meshed by 104×32=3328 elements and 105×33=3465 nodes; the third one is that 96×48 uniform element grid is used for 𝑥=0∼2 and 64×48 element grid for 𝑥=2∼4, the flow field is meshed by 160×48=7680 elements and 161×49=7889 nodes.

The numerical analyses of steady flow at different Re numbers (Re = 25, 50, 75, 100, 150) are performed by the NMM solutions of steady N-S equations. The streamline patterns and variables distributions of flow field at Re=100 from three different element grids are compared. The results from the first element grid are different to these from two other grids, but the results from the second grid are closely similar to these from the third one, and the average difference of velocity and pressure on the channel centreline is less than 5%, so the second element grid is enough for flow at Re<100, and the third element grid or finer element grid is more reasonable for flow at Re=100 and 150.

The streamline patterns and pressure distributions at different Re numbers are shown in Figure 4. At Re=25 and 50, there are two stable symmetrical vortexes behind the square cylinder, and the flow patterns are similar to the numerical solutions of Mukhopadhyay et al. [15]. As Re number increases to Re ≈[50,200], Cherry’s results show that two vortexes become unsymmetrical and Karman vortex street is formed behind the square cylinder gradually [14], but the stable convergent results are obtained from NMM solutions of steady N-S equations at Re=75, 100, and 150, and there are also two stable symmetrical vortexes behind the square cylinder. The distances from back flow point of the vortexes to square cylinder at Re=25, 50, 75, 100, 150 are 1.0, 1.2, 1.4, 1.6, 1.8, respectively, which are approximately linear with Re numbers. These results illustrate that numerical manifold schemes for steady N-S equations have the strong-stability-preserving property.

Figure 4 

Streamline pattern and pressure distribution for flow field at different Re numbers.

In NMM analysis unlike FVM and FEM analysis, the pressure distributions of flow field are obtained from direct numerical solutions of N-S equations and continuity equation, so it can improve the solution accuracy of pressure field. The pressure distributions of flow field show that high-pressure area is formed in frontage of the square cylinder for flow blockage, low pressure areas are formed near behind two front corner points for flow separating, and large pressure gradient is produced near the corner points. The pressure distributions are symmetrical, and the pressure of flow field will decrease as Re number increasing.

The velocity 𝑢 distributions along horizontal lines through geometric center of flow field are shown in Figure 5. The velocity 𝑢 drops rapidly from inlet to the square cylinder for flow blockage. Negative velocity appears behind the cylinder for back flow, and the velocity will rise gradually with away from the cylinder, and the velocity restoring becomes slow as Re number increases. The pressure distributions along horizontal lines through geometric center of flow field are shown in Figure 6. The pressure rises gradually from inlet to the square cylinder. Positive pressure appears behind the cylinder at low Re number, negative pressure appears at comparatively high Re number, and the pressure will drops or rises to zero at outlet of flow field.

Figure 5 

Comparison of velocity 𝑢 along horizontal lines through geometric center of flow field.

Figure 6 

Comparison of pressure along horizontal lines through geometric center of flow field.

4.3. Analysis for Unsteady Flow around a Square Cylinder in a Channel

The flow field configuration and boundary conditions of unsteady flow around a square cylinder in a channel are the same as these of steady flow, and the flow field is meshed by the third mixed grids. The unsteady flows at different Re numbers (Re=50, 75, 100, 150) are also analyzed through NMM solutions of N-S equations in this paper, where numerical solutions for steady flow at Re=50 are adopted as the initial conditions, the total computational time is 10 s, and the time step is 0.01 s.

When Re≤50, numerical solutions from unsteady flow analyses are closely similar to the results from steady flow analyses, and there are two steady symmetrical vortexes behind the square cylinder. As Re number increasing, two vortexes become unsymmetrical, and Karman Vortex Street is formed behind the square cylinder gradually. The streamline patterns for flow field at Re=150 are shown in Figure 7, and the periodical change process of appearing, developing, moving, separating of the vortexes are clearly displayed in the streamline patterns. The Strouhal number is 0.2438, the difference is 2.2% and 6.5% to [15, 16]. The flow patterns are similar to Mukhopadhyay’s and Yutaka’s results, which is obtained from very fine element grid [15, 17].

(a) T/4

(b) T/2

(c) 3T/4

Figure 7 

The streamline patterns for flow field from 48×160 element grid and Δ𝑡=0.01 s time step at Re = 150.

The pressure contours for flow field at Re=150 are show in Figure 8. The pressure distributions beside and in frontage of the square cylinder are almost identical to the results from steady flow analyses and are different behind the cylinder, which will change with the moving of vortexes.

(a) T/4 

(b) T/2

(c) 3T/4

Figure 8 

The pressure contours for flow field from 48×160 element grid and Δ𝑡=0.01s time step at Re=150.

5. Summary

Numerical manifold method for direct coupled solution of incompressible viscous flow N-S equations has been developed in this paper. Numerical manifold schemes integrated velocity and pressure were derived based on Galerkin-weighted residuals method as well. Mixed cover with linear polynomial function for velocity and constant function for pressure was adopted in finite element cover system. Compared with FVM and FEM, in NMM for incompressible viscous flow, the accuracy of velocity variables approximation can be improved by adopting high-order cover function, direct numerical solution of N-S equations, and continuity equation coupled velocity and pressure variables can be implemented by adopting finite cover system with mixed cover manifold element, so it can improve the solution accuracy of velocity and pressure variables.

As an application, mixed cover 4-node rectangular manifold element has been used to simulate flow around a square cylinder in a channel and past a step at low Re numbers. As to flow around a square cylinder in a channel, accurate numerical results have been presented for steady flow at Re=25, 50, 75, 100, 150 and for unsteady flow at Re=100, 150. The rules on streamline patterns and pressure distributions were analyzed in details. Numerical solutions were in very good agreement with the published data from very fine element meshes. Numerical tests indicate that manifold method is an effective and high-order accurate numerical method for incompressible viscous flow Navier-Stokes equations.

Acknowledgment

Financial support for this project was provided by the National Natural Science Foundation of China (no. 50975050).