An embedding finite element method for viscous incompressible flows with complex immersed boundaries on Cartesian grids. Issue 4 (2014)

• Record Type: Journal Article
• Authors: D.C. Lo; Chih‐Min Hsieh; D.L. Young
• Is Part Of: Engineering computations. Volume 31:Issue 4(2014)
• Page Start: 656
• Page End: 680
• Publication Date: 2014
• Subjects: Computer-aided engineering -- Periodicals; Computer graphics -- Periodicals; 620.00285
• DOI: 10.1108/EC-04-2012-0090
• Languages: English
• ISSNs: 0264-4401

Abstract

<abstract> <title> <x xml:space="preserve"> Abstract </x> </title> <p> <bold>Purpose</bold> – The main advantage of the proposed method is that the computations can be performed on a Cartesian grid with complex immersed boundaries (IBs). The purpose of this paper is to device a numerical scheme based on an embedding finite element method for the solution of two‐dimensional (2D) Navier‐Stokes equations. <bold>Design/methodology/approach</bold> – Geometries featuring the stationary solid obstacles in the flow are embedded in the Cartesian grid with special discretizations near the embedded boundary to ensure the accuracy of the solution in the cut cells. To comprehend the complexities of the viscous flows with IBs, the paper adopts a compact interpolation scheme near the IBs that allows to satisfy the second‐order accuracy and the conservation property of the solver. The interpolation scheme is designed by virtue of the shape function in the finite element scheme. <bold>Findings</bold> – Three numerical examples are selected to demonstrate the accuracy and flexibility of the proposed methodology. Simulation of flow past a circular cylinder for a range of Re=20‐200 shows excellent agreements with other results using different numerical schemes. Flows around a pair of tandem cylinders and several bodies are particularly investigated. The paper simulates the time‐based variation of the flow phenomena for uniform flow past a pair of cylinders with various streamwise gaps between two cylinders. The results in terms of drag coefficient and Strouhal number show excellent agreements with the results available in the literature. <bold>Originality/value</bold> – Details of the flow characteristics, such as velocity distribution, pressure and vorticity fields are presented. It is concluded the combined embedding boundary method and FE discretizations are robust and accurate for solving 2D fluid flows with complex IBs.</p> <ack> <title> <x xml:space="preserve"> Acknowledgements </x> </title> <p>The support by the National Science Council, Taiwan, under the grant No. 100‐2221‐E‐022‐015 for the work is gratefully acknowledged.</p> </ack> </abstract>