Naslov
Design and implementation of Immersed Boundary Method with discrete forcing approach for boundary conditions

Autori
Jasak, Hrvoje ; Rigler, Damir ; Tuković, Željko

Vrsta, podvrsta i kategorija rada
Radovi u zbornicima skupova, cjeloviti rad (in extenso), znanstveni

Izvornik
6th European Conference on Computational Fluid Dynamics, ECFD 2014 / Huerta, A. ; Onate, E. ; Oliver, X. - : International Center for Numerical Methods in Engineering (CIMNE), 2014, 5319-5332

ISBN
978-849428447-2

Skup
11th World Congress on Computational Mechanics, WCCM 2014, 5th European Conference on Computational Mechanics, ECCM 2014 and 6th European Conference on Computational Fluid Dynamics, ECFD 2014

Mjesto i datum
Barcelona, Španjolska, 20.07.2014. - 25.07.2014

Vrsta sudjelovanja
Predavanje

Vrsta recenzije
Međunarodna recenzija

Ključne riječi
CFD ; Immersed Boundary Method ; Immersed boundary wall functions ; Turbulent flows

Sažetak
This article describes design, implementation and validation of Immersed Boundary Method (IB) for finite volume based Computational Fluid Dynamics (CFD). In CFD, objects are usually described using body conformal meshes, where boundary conditions are imposed on the boundary faces of a mesh and provide the information necessary to properly describe object geometry and flow parameters. Generation of such meshes still requires some effort, even for simple geometries. If an object changes its position (or form) during time, mesh motion, or even topological changes of the mesh are required. Both require computational effort and can result in an unusable mesh. These obstacles can be avoided with an alternative approach: Immersed Boundary Method (IBM). IB method does not use a body- conformal mesh. Instead, background volume mesh is used for the flow domain together with a separate IB surface mesh for the object outline. Generation of the IB surface mesh is not related to the generation of the volume mesh making it significantly easier. As the object is not described using boundary faces, it is necessary to use a IB surface mesh for boundary condition imposition. The IB surface mesh is then immersed into the domain volume and the calculation matrix is modified to account for the object geometry. This enables very fast solutions for different geometries on the same background mesh. Furthermore, moving the IB surface mesh through the background mesh is significantly easier than any mesh-manipulation method. For this work IB method has been implemented into OpenFOAM using a discrete forcing approach, where discretized equations are only affected in the cells touching the IB (direct imposition of Boundary Conditions). Direct imposition of BC preserves sharpness of the object surface. Wall functions are used to allow high Reynolds number flows. Dirichlet and Neumann BC are implemented using quadratic interpolation, where unknown coefficients are determined using least square method on extended stencil [1]. Implementation into OpenFOAM enables use of polyhedral background meshes and automatic mesh refinement. The IB method implementation is validated with four test cases that cover transient and steady-state, laminar and turbulent, single phase and free surface two-phase flow. Both integral values (such as force on the object) and certain fields (such as pressure, or free surface) are compared with experimental results from literature and show satisfactory agreement.

Izvorni jezik
Engleski

Znanstvena područja
Strojarstvo

POVEZANOST RADA