A boundary condition-enforced immersed boundary method for compressible viscous flows. (10th September 2016)

• Record Type: Journal Article
• Authors: Qiu, Y.L.; Shu, C.; Wu, J.; Sun, Y.; Yang, L.M.; Guo, T.Q.
• Is Part Of: Computers & fluids. Volume 136(2016)
• Page Start: 104
• Page End: 113
• Publication Date: 2016-09-10
• Subjects: Immersed boundary method -- Compressible flows -- Cartesian mesh -- Diffuse interface IBM; Fluid dynamics -- Data processing -- Periodicals; 532.050285
• DOI: 10.1016/j.compfluid.2016.06.004
• Languages: English
• ISSNs: 0045-7930

Abstract

Highlights: It is the first work to develop diffuse interface IBM for compressible flows. A fractional step method is adopted, including predictor step and corrector step. Physical boundary conditions are accurately enforced. Force calculation is simply done from Newton's 2nd law. It is easy for applications. Abstract: A boundary condition-enforced immersed boundary method (IBM) for simulating compressible viscous flows is presented in this paper. It is developed from the implicit velocity correction-based IBM proposed by Wu and Shu (Journal of Computational Physics 228, 1963–1979 (2009)) for simulation of incompressible viscous flows. In their work, the velocity field is corrected by implementing the no-slip boundary condition. However, their method cannot be directly applied to simulate compressible flows. This is because the immersed boundary has a great effect on velocity, temperature, density and pressure fields in compressible flows, but their method has no mechanism to correct the density and pressure fields. To overcome this difficulty, in this work, the momentum field is firstly corrected by converting the no-slip velocity condition to the momentum condition. After that, from the continuity equation, the density correction can be made from the momentum correction. Then, the corrected velocity field can be given from the corrected momentum field and corrected density. Similar to the correction of momentum field, the temperature field can be corrected from the given temperature condition. Finally, the pressure can be obtained from the equation of state by using the corrected density and temperature. In this way, all the flow variables are corrected and the given physical boundary conditions are accurately implemented. To validate the proposed method, the flows over a stationary cylinder and one or two airfoils at different Mach and Reynolds numbers are simulated. Compared to the results in the literature, good agreement is achieved.