亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        NUMERICAL METHOD FOR MULTI-BODY FLUID INTERACTION BASED ON IMMERSED BOUNDARY METHOD*

        2011-06-27 05:54:02MINGPingjianZHANGWenping

        MING Ping-jian, ZHANG Wen-ping

        School of Power and Energy Engineering, Harbin Engineering University, Harbin 150001, China, E-mail: pingjianming@hrbeu.edu.cn

        NUMERICAL METHOD FOR MULTI-BODY FLUID INTERACTION BASED ON IMMERSED BOUNDARY METHOD*

        MING Ping-jian, ZHANG Wen-ping

        School of Power and Energy Engineering, Harbin Engineering University, Harbin 150001, China, E-mail: pingjianming@hrbeu.edu.cn

        A Cartesian grid based on Immersed Boundary Method (IBM), proposed by the present authors, is extended to unstructured grids. The advantages of IBM and Body Fitted Grid (BFG) are taken to enhance the computation efficiency of the fluid structure interaction in a complex domain. There are many methods to generate the BFG, among which the unstructured grid method is the most popular. The concept of Volume Of Solid (VOS) is used to deal with the multi rigid body and fluid interaction. Each body surface is represented by a set of points which can be traced in an anti-clockwise order with the solid area on the left side of surface. An efficient Lagrange point tracking algorithm on the fixed grid is applied to search the moving boundary grid points. This method is verified by low Reynolds number flows in the range from =100 Re to 1 000 in the cavity with a moving lid. The results are in a good agreement with experimental data in literature. Finally, the flow past two moving cylinders is simulated to test the capability of the method.

        fluid-structure interaction, Immersed Boundary Method (IBM), Volume Of Solid (VOS), unstructured grids

        Introduction

        The multi-body and fluid interaction can be widely found in marine and ocean engineering, and is a very important issue in Computation Fluid Dynamics (CFD). Numerical methods were developed as a complementary to experimental and analytical studies, including the Lagrangian moving mesh method[1]and the Chimera grid method[2,3]. The Lagrangian moving mesh is very inefficient or even impossible to simulate large amplitude motions, for which a remeshing procedure[4]is needed to deal with problems related with the low mesh quality. The Chimera grid method can deal with large amplitude moving body problems, but it is very expensive for the interpolation of variable values between moving grids and background grids. Recently, an immersed boundary method was proposed to simulate the fluid flow and rigid body interaction on fixed meshes. The immersed boundary method was originally proposed by Peskin in 1972 to analyze the blood flow in a human heart model. In his simulation, only a simple Cartesian mesh system wasused, with grids not necessarily conforming to the geometry of the heart[5]. The immersed boundary method is put into a mathematical formulation and the effects of solid boundaries are reflected in forces in the Navier-Stokes equations. There comes the socalled “Immersed Boundary (IB)” method, which gradually becomes a basic concept that the body is simply ‘immersed’ in the background mesh without necessarily fitting the underlying mesh points. Many modified IBMs were proposed, which are quite different from the most conventional CFD approaches[6]. More details of IBM history can be found in literature[7].

        The Immersed Boundary Method (IBM) is originally implemented in the Cartesian grid, with an obvious drawback that the grid number is increased drastically with the increase of the Reynolds number[8]. In some situations, a great number of extra outside grids are required in the Cartesian grid, for example, in the lobe pump as shown in Fig.1(a). In order to reduce the extra grid number, IBM is extended to curvilinear grids to generate a body fitted grid through a coordinate transform[9]. In this article, the immersed boundary method is further improved, and then is used for unstructured grids as shown in Fig.1(b). The scopeof the immersed boundary method is thus extended. Combined with the concept of Volume Of Solid (VOS), this article focuses on the multi-body and fluid interaction problem.

        Fig.1 Examples of computational grids

        The immersed boundary methods can be roughly divided into two classes[10]. In the first class, the immersed solid surface is regarded as a “diffuse”interface of finite thickness to represent its effects on the surrounding fluid. The thickness is usually taken as the local mesh spacing. In this class includes most methods that employ forces and/or mass sources as well as those that use the volume-of-fluid and phasefield approaches. In the second class, the boundary is tracked as a sharp interface and the communication between the moving boundary and the flow solver is usually achieved directly through modifying the computational stencil near the immersed boundary.

        In this article, the IBM method of the first class[11]is extended to deal with the multi-body and fluid interaction with unstructured grids, because it is easy to be used with only a few modifications of the exiting CFD codes. This method is implemented in an inhouse code to solve multi-body fluid interaction problems. The VOS can be further simplified by assuming that the body interior is occupied by the same incompressible fluid as that outside with a prescribed divergence-free velocity field. In this view, a fluidbody interface is similar to a fluid-fluid interface in the Volume Of Fluid (VOF) method for the two-fluid flow problems. The body can thus be identified by the VOS function as the VOF function[12,13]. Since the body is assumed to be incompressible, the pressure field inside the body obeys the same governing equation as that outside. It is the basic idea which is used to develop this immersed boundary method.

        1. Computation scheme

        1.1Governing equations

        The governing equations for viscous flow of incompressible fluid can be written in the integral form as:

        Continuity

        whereΩis the integral domain and ?Ωis the boundary, dS, dΩare the elemental area vector and the elemental domain, respectively,Uis the velocity vector (u,v),ρis the fluid density,μis the fluid dynamic viscosity,pis the static pressure andFis the body force vector (fx,fy).

        1.2Numerical scheme

        1.2.1 Momentum prediction

        The momentum equation can be solved for the cell center velocity vector with a given face flux and pressure. They can be either the initial values or intermediate values.

        The Navier-Stokes equation is discretized in space with the collocated finite volume method and in time with the Euler implicit method. The linear equation set can be obtained as:

        The VOS concept is used to consider the immersed boundary effect. The following formulations

        are applied to calculate the body forces, whereu′′andv′′ represent the second correction values which can be linearly interpolated from the fluid and the immersed rigid body velocities as done by NG[13]:

        whereαis the total solid fraction within a cell. The body force can be rewritten as follows

        From the Newton’s third law, we can obtain the forces on the rigid body.

        1.2.2 Pressure correction

        The velocity does not always satisfy the continuity equation after the momentum prediction and it is necessary to correct it based on the mass conservation equation. A pressure correction equation can be derived as in literature[14].

        The pressure can be corrected by solving the above equation. The face flux and the velocity can also be updated in this manner.

        1.2.3 Numerical method of VOS

        A link data structure is used for the immersed boundary containing a number of points. An efficient Lagrange point tracking algorithm on a fixed grid is applied to search the moving boundary grid points to improve the code performance[15]. The link direction is decided by the solid areas always on the left side of the boundary. In order to calculate the solid volume fraction in each cell, an array of pointer data structure is defined and a local polygon link is formed in the shadowed area. After the local link formation, the transformation from link to graph is completed. The solid volume fraction is calculated based on these graphs[14]. There are some special cases which require careful treatments, the details of which can be found in Ref.[16].

        Fig.2 Sketch of immersed boundary fronts

        After calculating the solid volume fraction of all cells which cross the immersed boundary, it is necessary to find the cells which are totally contained in the solid area. In this article, a method based on the concept of the advancing front is proposed to identify those cells quickly. As shown in Fig.2, one seeks the cell faces that are in the solid area but not cross the immersed boundary, which would form the front link, and then one advances the fronts until they disappear. The unstructured gird around an ellipse and the partially enlarged diagram are shown in Fig.2. The bold solid line represents the edge link, the thin line represents the immersed boundary and the dotted line represents the grid mesh meeting with the immersed boundary.

        In order to deal with the multi-body case, the procedure is repeated for each body surface. As shown in Fig.3, there are two cylinders, so the VOS calculation procedure will be carried out twice.

        Fig.3 Sketch of multi-body case

        2. Presentation of results

        2.1Lid driven cavity flow

        In order to verify the present immersed boundary method, a benchmark problem is considered. This method is also used to simulate the flow past a stationary cylinder as in literature[17-20].

        Fig.4 The sketch of the computational domain and grids

        Fig.5 VOS distribution

        Fig.7 Results forRe=400

        Fig.6 Results forRe=100

        Fig.8 Results forRe=1000

        Fig.9 Computational domain

        Fig.10 Contours of pressure around two cylinders moving to each other atRe=40 at three different stages. The first graph is obtained with a coarse grid and the second one is obtained with a refined grid

        The lid driven cavity flow is a classic benchmark case, which is usually used to evaluate numerical methods and codes[21,22]The computational domain is shown in Fig.4, which is a 1.1 m×1.1 m cavity square. In order to test the present numerical method and codes, a boundary layer of 0.1m in width is assumed to be the solid zone on bottom and right sides. The VOS distribution is shown in Fig.5.

        Fig.11 Contours of vorticity around two cylinders moving to each other atRe=40 at three different stages. The first graph is obtained with a coarse grid and the second one is obtained with a refined grid

        Ghia’s results are a set of the most popular experimental data[23], which are extensively used, andhere they are regarded as the standard. In this article, simulations with different Reynolds numbers ranging from 100 to 1 000 are carried out and the results are shown in Figs.6 through 8. The Reynolds number is expressed as

        whereUmaxis the lid moving velocity anddis the width of the cavity.

        As shown in Fig.6, the velocity distribution agrees well with the experimental data. There are two small eddies and one dominant eddy. The two small eddies grow with the increase of Reynolds number and the center of the dominant eddy moves downwards. The difference of the velocity distribution between numerical results and experimental data increases with the increase of Reynolds number.

        2.2Two cylinders moving to each other

        This case was firstly used by Wang in a series articles for verification. Two cylinders move to each other with a velocity of 1 m/s and the diameter of both cylinders is 1 m as shown in Fig.9.

        Two sets of meshes with 464 052 and 629 268 triangle cells with refinement near the cylinder motion zone are used to simulate the multi-body and fluid flow interaction. The contours of the pressure and the vorticity at three different instants are shown in Fig.10 and Fig.12, respectively, and the top one is the result obtained with a coarse grid and the bottom one is obtained with a refined grid. As shown in Fig.10 and Fig.11, the pressure and velocity contours are similar at different instants.

        Fig.12 Contours of vorticity

        The comparison of the vorticity att=24 with that in literature[2]as shown in Fig.12 shows that the present method is reasonable and accurate.

        3. Conclusions

        A numerical method on the fluid-structure interaction based on the immersed boundary method and the volume of solid method is extended to unstructured grids. This method can be applied to any existing pressure based finite volume CFD solver without complex procedures.

        To verify the present method, the flows past a stationary and an oscillating cylinder are simulated and the results agree well with those in literature. The results for two cylinders moving to each other are reasonable and consistent with Wang’s data, which indicates that this method is accurate.

        Acknowledgement

        This work was supported by the Fundamental Research Fund of Harbin Engineering University (Grant No. HEUFT 09005).

        [1] FARAHAT C., GEUZAINE P. and GRANDMONT C. The discrete geometric conservation law and nonlinear stability of ALE schemes for solution of flow problems on moving grids[J].Journal of Computational Physics,2001, 174(2): 669-694.

        [2] WANG Z. J., PARTHASARATHY V. A fully automated Chimera methodology for multiple moving body problems[J].International Journal for Numerical Methods in Fluids,2000, 33(7): 919-938.

        [3] CHEN H. C., LIU T. Time domain simulation of large amplitude ship roll motions by a Chimera RANS method[J].International Journal of Offshore and Polar Engineering,2002, 12(3): 206-212.

        [4] GARCIA J., ONATE E. An unstructured finite element solver for ship hydrodynamics problems[J].Journal of Applied Mechanics,2003, 70(1): 18-27.

        [5] PESKIN C. S. Flow patterns around heart valves: A numerical method[J].Journal of Computational Physics,1972, 10(2): 252-271.

        [6] DUTSCH H., DURST F. and BECKER S. Low Reynolds number flow around an oscillating circular cylinder at low Keulegan-Carpenter numbers[J].Journal of Fluid Mechanics,1998, 360: 249-271.

        [7] MITTAL R., IACCARINO G. Immersed boundary methods[J].Annual Review Fluid Mechanics,2005, 37: 239-61.

        [8] GHIAS R., MITTAL R. and DONG H. A sharp interface immersed boundary method for compressible viscous flows[J].Journal of Computational Physics,2007, 225(1): 528-553.

        [9] ROMAN F., NAPOLI E. and MILICI B. et al. An improved immersed boundary method for curvilinear grids[J].Computers and Fluids,2009, 38(8): 1510-1527.

        [10] UDAYKUMAR H. S., MITTAL R. and RAMPUNGGOON P. et al. A sharp interface cartesian grid method for simulating flows with complex movingboundaries[J].Journal of Computational Physics,2001, 174(1): 345-380.

        [11] MING Ping-jian, ZHANG Wen-ping. Numerical simulation of low Reynolds number fluid-structure interaction with immersed boundary method[J].Chinese Journal of Aeronautics,2009, 22(5): 480-485.

        [12] PAN D. An immersed boundary method for incompressible flows using volume of body function[J].International Journal for Numerical Methods in Fluids,2006, 50(6): 733-750.

        [13] NG K. C. A collocated finite volume embedding method for simulation of flow past stationary and moving body[J].Computers and Fluids,2009, 38(2): 347-357.

        [14] MING P. J., SUN Y. Z. and DUAN W. Y. et al. Unstructured grid immersed boundary method for numerical simulation of fluid structure interaction[J].Journal of Marine Science and Application,2010, 9(2): 181-186.

        [15] LIU Yong-feng, MING Ping-jian and ZHANG Wenping et al. An efficient Lagrange point tracking algorithm for fixed grids[J].Chinese Journal of Computational Physics,2010, 27(4): 527-532(in Chinese).

        [16] MING Ping-jian, ZHANG Wen-ping and LU Xi-qun et al. Numerical simulation of fluid structure interaction based on combination of immersed boundary method and unstructured grids[J].Chinese Journal of Hydrodynamics,2010, 25(3): 323-331(in Chinese).

        [17] KIM J., KIM D. and CHOI H. An immersed boundary finite volume method for simulations of flow in complex geometries[J].Journal of Computational Physics,2001, 171(1): 132-150.

        [18] CALHOUN D. A Cartesian grid method for solving the two-dimensional stream function-vorticity equations in irregular region[J].Journal of Computational Physics,2002, 176(2): 231-275

        [19] RUSSELL D., WANG Z. J. A Cartesian grid method for modeling multiple moving objects in 2D incompressible viscous flow[J].Journal of Computational Physics,2003, 191(1): 177-205.

        [20] JUNG-Il C., ROSHAN C. O. and JACK R. E. et al. An immersed boundary method for complex incompressible flows[J].Journal of Computational Physics,2007, 224(2): 757-784.

        [21] SCHREIBER R., KELLER H. B. Driven cavity flows by efficient numerical techniques[J].Journal of Computational Physics,1983, 49(2): 310-333.

        [22] JUN Z. Numerical simulation of 2D square driven cavity using fourth-order compact finite difference schemes[J].Computers and Mathematics with Applications,2003, 45(1-3): 43-52.

        [23] GHIA U., GHIA K. N. and SHIN C. T. High-resolutions for incompressible flow using the Navier-Stokes equations and a multigrid method[J].Journal of Computational Physics,1982, 48(3): 387-411.

        February 12, 2011, Revised April 7, 2011)

        * Biography: MING Ping-jian (1980-), Male, Ph. D., Lecturer

        2011,23(4):476-482

        10.1016/S1001-6058(10)60138-2

        精品国产一区二区三区不卡在线| 日本香蕉久久一区二区视频| 毛片av中文字幕一区二区| 日本一区二区三区熟女俱乐部| 色综合久久中文娱乐网| 国产肉丝袜在线观看| 日本精品视频一视频高清| 亚洲性日韩一区二区三区| 久久五月精品中文字幕| 在线观看视频日本一区二区| 日本丰满熟妇videossex8k| 国产在线手机视频| 色婷婷久久99综合精品jk白丝 | 国产一区二区三区视频网 | 女人被躁到高潮嗷嗷叫免费软| 人妖啪啪综合av一区| 亚洲视频在线观看| 日日摸夜夜添夜夜添无码免费视频| 国产主播一区二区在线观看 | 91偷自国产一区二区三区| 亚洲中文字幕无码中文字| 久久精品国产一区二区电影| 天堂女人av一区二区| 亚洲第一幕一区二区三区在线观看| 国产精品毛片一区二区| 国模精品二区| 日本中文字幕官网亚洲| 亚洲熟女综合色一区二区三区| 亚洲男同帅gay片在线观看| 亚洲精品天堂在线观看| 丝袜美腿国产一区二区| 狠狠色婷婷久久综合频道日韩| 亚洲久无码中文字幕热| 亚洲人妻av综合久久| 久久久久九九精品影院| 97人人超碰国产精品最新o| 性色av成人精品久久| 国语对白福利在线观看 | 无套内内射视频网站| 丰满爆乳无码一区二区三区| 国产精品一区二区av白丝在线|