A Conservative Multi-block Algorithm for Two-dimensional Numerical Model

Authors: Yaoxin Zhang, Yafei Jia, Sam S.Y. Wang

Abstract:

A multi-block algorithm and its implementation in two-dimensional finite element numerical model CCHE2D are presented. In addition to a conventional Lagrangian Interpolation Method (LIM), a novel interpolation method, called Consistent Interpolation Method (CIM), is proposed for more accurate information transfer across the interfaces. The consistent interpolation solves the governing equations over the auxiliary elements constructed around the interpolation nodes using the same numerical scheme used for the internal computational nodes. With the CIM, the momentum conservation can be maintained as well as the mass conservation. An imbalance correction scheme is used to enforce the conservation laws (mass and momentum) across the interfaces. Comparisons of the LIM and the CIM are made using several flow simulation examples. It is shown that the proposed CIM is physically more accurate and produces satisfactory results efficiently.

Keywords: Multi-block algorithm, conservation, interpolation, numerical model, flow simulation.

Digital Object Identifier (DOI): doi.org/10.5281/zenodo.1055851

References:

[1] Benek, J.A., Steger, J.L., and Dougherby, F.C. (1983), ''A Flexiable Grid Embedding Technique with Application to the Euler Equations'', AIAA Paper 83-1944, 1983.
[2] Berger, M. (1987), ''On Conservation at Grid interface'', SIAM J. Numer. Anal, 24(5).
[3] Brakkee, E., Wesseling, P., and Kassels, C.G. (2000), ''Schwarz domain decomposition for the incompressible Navier-Stokes equations in general coordinates'', Int. J. Numer. Methods Fluids, 32:141-173.
[4] Caruso, S.C.(1985), ''Adaptive grid techniques for elliptic fluid-flow problems'', Ph.D. Dissertation, Department of Mechanical Engineering, Stanford University, Stanford, CA.
[5] Chen, H.C., and Chen, M.(1998), ''Chimera RANS simulation of a berthing DDG-51 ship in translational and rotational motions'', Int. J. Offshore and Polar Engrg., 8(3): 182-191.
[6] Chen, H.C., and Liu, T.(1999), ''Turbulent flow Induced by Full-scale Ship in Harbor'', ASCE, J. Engrg. Mech.,125(7): 827-835.
[7] Chesshire, G., and Henshaw, W.D.(1994), ''A scheme for conservative interpolation on overlapping grids'', SIAM J. of Scientific Computing, 15(4): 819-845.
[8] Francisco, J.S., and Sabadell, F.J.(2000), ''The multi-block method---a new strategy based on domain decomposition for the solution of wave propagation problems'', Int. J. Numer. Methods Engrg., 49: 1353-1376.
[9] Furukawa, M., Yamasaki, M., and Inoue, M.(1989), ''A zonal approach for solving the Navier-Stokes Equations using a TVD Finite Volume Method'', International Symposium on CFD, Nagoya, Nagoya Trade & Industry Center, Japan, August 28-31, 1989.
[10] Gresho, P.M. and Sani, R.L.(1987), ''On pressure boundary conditions for the incompressible Navier-Stokes equations'', Int. J. Numer. Methods Fluids, 7: 1111-1145.
[11] Henshaw, W.D. (1985), ''Part I: The numerical solution of hyperbolic systems of conservation laws. Part II: Composite overlapping grid techniques'', Ph.D. Dissertation, California Institute of Technology, Pasadena, CA, 1985.
[12] Henshaw, W.D. (1994), ''A fourth-order accurate method for the incompressible Navier-Stokes equations on overlapping grids'', J. Comp. Phys, 113(1):13-25.
[13] Hessinius, K., and Rai, M. (1986), ''Three Dimensional Conservative, Euler Computations using Patched Grid Systems and Explicit Methods'', AIAA Paper 86-1088, 1986.
[14] Jia, Y.F., and Wang, S. S-Y.(1999), ''Numerical Model for channel flow and morphological changes studies'', ASCE J. Hydraulic Engrg, 125(9): 924-933.
[15] Jia, Y.F., and Wang, Sam S.Y. (2001a), ''CCHE2D: Two-dimensional Hydrodynamic and Sedimenst Transport Model for Unsteady Open Channel Flows Over Loose Bed'', Technical Report, NCCHE-TR-2001-1, Aug. 2001.
[16] Jia, Y.F., and Wang, Sam S.Y. (2001b), ''CCHE2D Verification and Validation Tests Documentation'', Technical Report, NCCHE-TR-2001-2, Aug. 2001.
[17] Jia, Y.F., Wang, S.S.Y., and Xu, Y.C (2002), ''Validation and Application of a 2D Model to Channels with Complex Geometry'', Int. J. Comp. Engr. Sci., 3(1): 57-71.
[18] Klopfer, G.H., and Molvik, G.A. (1991), ''Conservative multi-zonal interface algorithm for 3D Navier-Stokes equations'', AIAA paper 91-1601-CP, AIAA 10th CFD Conference, Honolulu, Hawaii.
[19] Meakin, R.L. (1986), ''Application of boundary conforming coordinate and domain decomposition principles to environmental flows'', Ph.D. Thesis, Department of Civil Engineering, Stanford University.
[20] Meakin, R.L., and Street, R.L. (1988), ''Simulation of environmental flow problems in geometrically complex domains, part 2: A domain-splitting method'', Compu. Methods Appl. Mech. Eng., 68: 311-331.
[21] Moon, Y. J. and Liou, M.-S. (1989), ''Conservative treatment of boundary interfaces for overlaid grids and multi-level grid adaptations'', AIAA Paper 89-1980-CP, AIAA 9th CFD Conference, Buffalo, NY.
[22] Perng, C.Y., and Street, R.L. (1991), ''A coupled multi-domain-splitting technique for simulating incompressible flows in geometrically complex domains'', Int. J. Numer. Methods Fluids, 13: 269-286.
[23] Rai, M. M. (1986a), ''A conservative treatment of zonal boundaries for Euler equation calculations'', J. Comput. Phys., 62: 472-503.
[24] Rai, M.M. (1986b), An Implicit Conservative Zonal Boundary Scheme for Euler Equation Calculations, Comp. Fluids, 14: 295-319.
[25] Rajaratnam, N., and Nwachukwu, B.A.(1983), ''Flow near groid-like structures'', J. of Hydraulic Engineering, 109(3): 463-481
[26] Shetter, A.S., and Murphy, K.K. (1996), ''A numerical study of division flow in open channels'', Journal of Hydraulic Research, IAHR, 34(5): 651-675.
[27] Sheng, C., Taylor, L.K., and Whitefield, D.L. (1995), ''Multiblock multigrid solution of three-dimensional incompressible turbulent flows about appended submarine configurations'', AIAA-95-0203, 33rd Aerospace sciences meeting and exhibit, January 9-12, 1995, Reno, NV.
[28] Sinha, K. S., Sotiropoulos, F., and Odgaard, A. J.(1998), ''Three-dimensional numerical model for flow through natural rivers'', J. Hydr. Engrg., ASCE, 124(1):13-24.
[29] Steger, J.L. (1991), ''Thoughts on the Chimera Method of Simulation of Three-Dimensional Viscous Flow'', in Proceedings, CFD symposium on Aeropropulsion, Cleveland, Ohio, NASA CP-3078, 1991, 1-10.
[30] Tu, J.Y., and Fuchs, L.(1992), ''Overlapping grids and multi-grid methods for 3D unsteady flow calculation in I.C engine'', Int. J. Numer. Methods Fluids, 15: 693-714.
[31] Van Leer, B. (1982), ''Flux-Vector Splitting for the Euler Equations'', Lecture Notes in Physics, Vol. 170: 507-512.
[32] Wang, S.S.Y., and Hu, K.K. (1997), ''Improved methodology for formulating finite-element hydrodynamic models'', Finite element in fluids, T.J.Chung, ed., Vol. 8, Hemisphere Publishing Corp., 457-478.
[33] Wang, Z.J. (1994), ''Conservative Chimera for 3D Euler Equations on Structured/Structured, Structured/Unstructured Grids'', Proceedings of the 2nd Overset Composite Grid and Solution Technology Symposium, Ft Walton Beach, Florida, October 1994.
[34] Wright, J., and Shyy, W.(1993), ''A pressure-based composite grid method for the Navier-Stokes equations'', J. Comp. Phy., 107: 225-238.
[35] Xie, B.L.(1996), ''Experiment on flow in a sudden-expanded channel'', Technical Report, Wuhan University of Hydraulic and Electric Engineering, China.
[36] Zang, Y., and Street, R.L. (1995), ''A composite multi-grid method for calculating unsteady incompressible flows in geometrically complex domains'', Int. J. Numer. Methods Fluids, 20: 341-361.