Fictitious Domain and Immersed Boundary methods in OpenFOAM: Application to Complex Geometries

Federico Municchi, Stefan Radl

Research output: Contribution to conference(Old data) Lecture or PresentationResearch

Abstract

Building high quality meshes can be a difficult task when dealing with complex geometries. Furthermore, many industrial applications require topological optimization procedures in order to, for example, minimize the drag force exerted on an immersed body. In other situations, it may be necessary to simulate a complex moving object, like a stirrer. The most straightforward (but sub-optimal) approach consists in defining moving patches in OpenFOAM and to perform a re-meshing procedure that results in a drastic change of the mesh topology. Furthermore, a subsequent mapping of previous fields into the new mesh is required. Both procedures are time consuming, since the meshing algorithm has to ensure a minimum mesh quality, and the overall quality of the resulting solution tends to decrease.

In our present contribution we show how to address these problems using two alternative strategies: a Fictitious Domain (FD) and Immersed Boundary (IB) methodology. Both methods rely on a suitable representation of immersed objects by means of additional terms in the governing equations rather than imposing complex boundary conditions. This allows us (i) to use structured, fixed grids, and (ii) to avoid remeshing operations. Clearly, this results in major savings in computational time, improved stability of the calculation, and increased accuracy due to the fixed, high-quality meshes, and the lack of the field mapping step.

We will show how the method can be implemented in OpenFOAM in a conservative and numerically stable way. We do this in the frame if a projection method for pressure-based solvers. Finally, we will show examples of applications, as well as present results of a verification study.
Original languageEnglish
Publication statusPublished - 10 Nov 2016
EventPFAU XIII: OpenFOAM user meeting - Vienna, Austria
Duration: 10 Nov 201610 Nov 2016

Workshop

WorkshopPFAU XIII: OpenFOAM user meeting
Abbreviated titlePFAUXIII
CountryAustria
CityVienna
Period10/11/1610/11/16

Fingerprint

Geometry
Industrial applications
Drag
Topology
Boundary conditions

Keywords

  • Immersed Boundary, Fictitious Domain, OpenFOAM, gas-particle

Cite this

Municchi, F., & Radl, S. (2016). Fictitious Domain and Immersed Boundary methods in OpenFOAM: Application to Complex Geometries. PFAU XIII: OpenFOAM user meeting, Vienna, Austria.

Fictitious Domain and Immersed Boundary methods in OpenFOAM: Application to Complex Geometries. / Municchi, Federico; Radl, Stefan.

2016. PFAU XIII: OpenFOAM user meeting, Vienna, Austria.

Research output: Contribution to conference(Old data) Lecture or PresentationResearch

Municchi, F & Radl, S 2016, 'Fictitious Domain and Immersed Boundary methods in OpenFOAM: Application to Complex Geometries' PFAU XIII: OpenFOAM user meeting, Vienna, Austria, 10/11/16 - 10/11/16, .
Municchi F, Radl S. Fictitious Domain and Immersed Boundary methods in OpenFOAM: Application to Complex Geometries. 2016. PFAU XIII: OpenFOAM user meeting, Vienna, Austria.
Municchi, Federico ; Radl, Stefan. / Fictitious Domain and Immersed Boundary methods in OpenFOAM: Application to Complex Geometries. PFAU XIII: OpenFOAM user meeting, Vienna, Austria.
@conference{0b057a2a9138469a81779e38b4f7d046,
title = "Fictitious Domain and Immersed Boundary methods in OpenFOAM: Application to Complex Geometries",
abstract = "Building high quality meshes can be a difficult task when dealing with complex geometries. Furthermore, many industrial applications require topological optimization procedures in order to, for example, minimize the drag force exerted on an immersed body. In other situations, it may be necessary to simulate a complex moving object, like a stirrer. The most straightforward (but sub-optimal) approach consists in defining moving patches in OpenFOAM and to perform a re-meshing procedure that results in a drastic change of the mesh topology. Furthermore, a subsequent mapping of previous fields into the new mesh is required. Both procedures are time consuming, since the meshing algorithm has to ensure a minimum mesh quality, and the overall quality of the resulting solution tends to decrease.In our present contribution we show how to address these problems using two alternative strategies: a Fictitious Domain (FD) and Immersed Boundary (IB) methodology. Both methods rely on a suitable representation of immersed objects by means of additional terms in the governing equations rather than imposing complex boundary conditions. This allows us (i) to use structured, fixed grids, and (ii) to avoid remeshing operations. Clearly, this results in major savings in computational time, improved stability of the calculation, and increased accuracy due to the fixed, high-quality meshes, and the lack of the field mapping step. We will show how the method can be implemented in OpenFOAM in a conservative and numerically stable way. We do this in the frame if a projection method for pressure-based solvers. Finally, we will show examples of applications, as well as present results of a verification study.",
keywords = "Immersed Boundary, Fictitious Domain, OpenFOAM, gas-particle",
author = "Federico Municchi and Stefan Radl",
year = "2016",
month = "11",
day = "10",
language = "English",
note = "PFAU XIII: OpenFOAM user meeting, PFAUXIII ; Conference date: 10-11-2016 Through 10-11-2016",

}

TY - CONF

T1 - Fictitious Domain and Immersed Boundary methods in OpenFOAM: Application to Complex Geometries

AU - Municchi, Federico

AU - Radl, Stefan

PY - 2016/11/10

Y1 - 2016/11/10

N2 - Building high quality meshes can be a difficult task when dealing with complex geometries. Furthermore, many industrial applications require topological optimization procedures in order to, for example, minimize the drag force exerted on an immersed body. In other situations, it may be necessary to simulate a complex moving object, like a stirrer. The most straightforward (but sub-optimal) approach consists in defining moving patches in OpenFOAM and to perform a re-meshing procedure that results in a drastic change of the mesh topology. Furthermore, a subsequent mapping of previous fields into the new mesh is required. Both procedures are time consuming, since the meshing algorithm has to ensure a minimum mesh quality, and the overall quality of the resulting solution tends to decrease.In our present contribution we show how to address these problems using two alternative strategies: a Fictitious Domain (FD) and Immersed Boundary (IB) methodology. Both methods rely on a suitable representation of immersed objects by means of additional terms in the governing equations rather than imposing complex boundary conditions. This allows us (i) to use structured, fixed grids, and (ii) to avoid remeshing operations. Clearly, this results in major savings in computational time, improved stability of the calculation, and increased accuracy due to the fixed, high-quality meshes, and the lack of the field mapping step. We will show how the method can be implemented in OpenFOAM in a conservative and numerically stable way. We do this in the frame if a projection method for pressure-based solvers. Finally, we will show examples of applications, as well as present results of a verification study.

AB - Building high quality meshes can be a difficult task when dealing with complex geometries. Furthermore, many industrial applications require topological optimization procedures in order to, for example, minimize the drag force exerted on an immersed body. In other situations, it may be necessary to simulate a complex moving object, like a stirrer. The most straightforward (but sub-optimal) approach consists in defining moving patches in OpenFOAM and to perform a re-meshing procedure that results in a drastic change of the mesh topology. Furthermore, a subsequent mapping of previous fields into the new mesh is required. Both procedures are time consuming, since the meshing algorithm has to ensure a minimum mesh quality, and the overall quality of the resulting solution tends to decrease.In our present contribution we show how to address these problems using two alternative strategies: a Fictitious Domain (FD) and Immersed Boundary (IB) methodology. Both methods rely on a suitable representation of immersed objects by means of additional terms in the governing equations rather than imposing complex boundary conditions. This allows us (i) to use structured, fixed grids, and (ii) to avoid remeshing operations. Clearly, this results in major savings in computational time, improved stability of the calculation, and increased accuracy due to the fixed, high-quality meshes, and the lack of the field mapping step. We will show how the method can be implemented in OpenFOAM in a conservative and numerically stable way. We do this in the frame if a projection method for pressure-based solvers. Finally, we will show examples of applications, as well as present results of a verification study.

KW - Immersed Boundary, Fictitious Domain, OpenFOAM, gas-particle

M3 - (Old data) Lecture or Presentation

ER -