`Journal of Applied MathematicsVolume 2012, Article ID 465972, 14 pageshttp://dx.doi.org/10.1155/2012/465972`
Research Article

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

School of Materials and Energy, Guangdong University of Technology, Guangzhou 510643, China

Received 15 January 2012; Accepted 3 August 2012

Copyright © 2012 Zhengrong Zhang and Xiangwei Zhang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### 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) . The applications of FDM are affected for unsuitability to complicated structure flow field and inconservation of field variables . The FEM schemes for N-S equations are very complicated, and the computational efficiency of nonstandard FEMs is low . FVM is an effective numerical method for solutions of fluid flow, but exact treatments of the moving boundary are difficult and complicated . 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 . 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 . 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 . 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 .

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 . 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 . 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 aswhere 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 , which meet the requirement of partition of unity , 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: where 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 where , denote the basic series for velocity and pressure cover functions of physical cover ; , are the cover DOF variables, and the DOF numbers are and .

By substituting (2.3) into (2.2), the global cover functions for velocity and pressure in element (e) can be written aswhere 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,

The partial derivatives of element cover functions can be obtained from (2.4a) and (2.4b) as the follows: where 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 that is,where ; and the coefficients for element manifold equations are

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 , 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: 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 where are similar to (2.8a) and (2.8b); , , and are the unknown variables. The equations can be written in the matrix form as that is, 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 where is time step; and are the velocity DOF variables at time step and , respectively. Then as to every manifold element (e), the element manifold equations (2.8b) can be reformulated as the follows: Equations (3.5) uniting (2.8a) can be stated in matrix form as that is, where is the pressure cover DOF variables at time step is the element unknown cover DOF variables at time step is the element flow matrix with respect to the unknown velocity DOF variables at time step 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, where are the coordinates of element (e) center point; are the side lengths.

##### 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 and , the fluid density is , and three dynamic viscosities of fluid are = 0.01, 0.005, 0.0025, and so Re numbers are known as =100, 200, 400, respectively. All the flow field is meshed into 196 elements and 229 nodes by rectangular 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 and rectangular grids, respectively. The NMM results from grids are nearly coincident with the FVM results from 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 , the width is . The size of square cylinder is , 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 is the average flow velocity in the channel, is the channel centerline and , the pressure condition of outlet BC is , and the nonslip boundary conditions are set on other fixed walls, that is, .

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 uniform element grid is used for the domain from to and element grid for the domain from to , the entire flow field is meshed by elements and 901 nodes; the second one is that uniform element grid is used for and element grid for , the flow field is meshed by elements and nodes; the third one is that uniform element grid is used for and element grid for , the flow field is meshed by elements and nodes.

The numerical analyses of steady flow at different Re numbers ( = 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 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 , and the third element grid or finer element grid is more reasonable for flow at and 150.

The streamline patterns and pressure distributions at different Re numbers are shown in Figure 4. At 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. . As Re number increases to Re , Cherry’s results show that two vortexes become unsymmetrical and Karman vortex street is formed behind the square cylinder gradually , but the stable convergent results are obtained from NMM solutions of steady N-S equations at , 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 , 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 (, 75, 100, 150) are also analyzed through NMM solutions of N-S equations in this paper, where numerical solutions for steady flow at are adopted as the initial conditions, the total computational time is 10 s, and the time step is 0.01 s.

When , 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 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].

Figure 7: The streamline patterns for flow field from element grid and  s time step at Re = 150.

The pressure contours for flow field at 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.

Figure 8: The pressure contours for flow field from element grid and s time step at .

#### 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 , 50, 75, 100, 150 and for unsteady flow at , 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).

#### References

1. J. H. Ferziger and M. Perić, Computational Methods for Fluid Dynamics, Springer, New York, NY, USA, 3rd edition, 1996.
2. J. Nordström, J. Gong, E. van der Weide, and M. Svärd, “A stable and conservative high order multi-block method for the compressible Navier-Stokes equations,” Journal of Computational Physics, vol. 228, no. 24, pp. 9020–9035, 2009.
3. L. Zhang and Z. Chen, “A stabilized mixed finite element method for single-phase compressible flow,” Journal of Applied Mathematics, vol. 2011, Article ID 129724, 16 pages, 2011.
4. S. Boivin, F. Cayré, and J. M. Hérard, “A finite volume method to solve the Navier-Stokes equations for incompressible flows on unstructured meshes,” International Journal of Thermal Sciences, vol. 39, no. 8, pp. 806–825, 2000.
5. P. Pořízková, K. Kozel, and J. Horáček, “Numerical simulation of unsteady compressible flow in convergent channel: pressure spectral analysis,” Journal of Applied Mathematics, vol. 2012, Article ID 545120, 9 pages, 2012.
6. G. H. Shi, “Manifold method,” in Proceedings of the 1st International Forum on Discontinuous Deformation Analysis & Simulation of Discontinuous Media, pp. 52–204, Berkeley, Calif, USA, 1996.
7. G. H. Shi, “Manifold method of material analysis,” in Proceedings of the 9th Army Conference on Applied Mathematics and Computing, pp. 51–76, Minneapolis, Minn, USA, 1992.
8. K. Terada and M. Kurumatani, “Performance assessment of generalized elements in the finite cover method,” Finite Elements in Analysis and Design, vol. 41, no. 2, pp. 111–132, 2004.
9. K. Terada and M. Kurumatani, “An integrated procedure for three-dimensional structural analysis with the finite cover method,” International Journal for Numerical Methods in Engineering, vol. 63, no. 15, pp. 2102–2123, 2005.
10. K. Terada, T. Ishii, T. Kyoya, and Y. Kishino, “Finite cover method for progressive failure with cohesive zone fracture in heterogeneous solids and structures,” Computational Mechanics, vol. 39, no. 2, pp. 191–210, 2007.
11. H. Y. Liu, S. R. Lv, and S. Q. Qin, “Simulation of rock shock damage evolution law with numerical manifold method,” Engineering Science, vol. 19, no. 1, pp. 392–396, 2007 (Chinese).
12. G. F. Wei and W. Feng, “Incompatible numerical manifold method based on heat exchange problem,” Chinese Quarterly of Mechanics, vol. 26, no. 3, pp. 451–454, 2005 (Chinese).
13. Z. R. Zhang, X. W. Zhang, and J. H. Yan, “Manifold method coupled velocity and pressure for Navier-Stokes equations and direct numerical solution of unsteady incompressible viscous flow,” Computers & Fluids, vol. 39, no. 8, pp. 1353–1365, 2010.
14. N. J. Cherry, R. Hillier, and M. E. M. P. Latour, “Unsteady measurements in a separated and reattaching flow,” Journal of Fluid Mechanics, vol. 144, pp. 13–46, 1984.
15. A. Mukhopadhyay, G. Biswas, and T. Sundararajan, “Numerical investigation of confined wakes behind a square cylinder in a channel,” International Journal for Numerical Methods in Fluids, vol. 14, no. 12, pp. 1473–1484, 1992.
16. Z. Dang, G. Xi, and S. J. Wang, “Application of the piso method in simulation of unsteady wakes behind a square cylinder,” Journal of Engineering Thermophysics, vol. 20, no. 3, pp. 317–321, 1999 (Chinese).
17. Y. Yoshida and T. Nomura, “A transient solution method for the finite element incompressible Navier-Stokes equations,” International Journal for Numerical Methods in Fluids, vol. 5, no. 10, pp. 873–890, 1985.
18. O. C. Zienkiewicz, S. Qu, R. L. Taylor, and S. Nakazawa, “The patch test for mixed formulations,” International Journal for Numerical Methods in Engineering, vol. 23, no. 10, pp. 1873–1883, 1986.