Immersed Boundary Ghost Cell Method

thehexman · January 29, 2010, 05:27

Hello all,

I am currently writing a student research paper on the use of the IBM to model the human heart valves (a not quite new topic, I'm afraid, but in my department this hasn't been done yet). Currently I'm doing literature research and finding out about the various methods of the IB concept. At the moment, I am having trouble understanding the Ghost Cell Method (GCM).

In particular, I do not exactly understand the interpolation. But first let me tell you what I understand so far:

- A forcing term is added to the discrete system of governing equations
- This forcing term is derived from the boundary conditions (e.g. no slip: u=0)
- In general, the derivation of the forcing term is done by estimating the velocity field and and correcting it at the boundary to fit the boundary conditions.
- Ghost cells are defined inside the immersed body, so that each GC has at least one neighbor in the fluid
- If the artificial value for the flow parameters were known in these cells, the computation could stop right there w/o having to solve in the fluid
- However, there is no information in the GC themselves, but nearby there is the boundary condition (as before, a no slip condition would be an example)
- The value of the GC flow parameters can then be extrapolated (using a bilinear or trilinear scheme, depending on the case) from the boundary conditions and the fluid cells close by.

My question is: what is the point of calculating the value of the GC flow parameters if the boundary condition is already known?

Isn't it just more work to do so, because the flow parameters inside the body are purely artificial?

I think this has something to do with the fact that the IB does not necessarily coincide with the mesh points and therefore one cannot impose boundary conditions on a discrete equation if these conditions do not lie on the mesh points. Therefore, another cell has to be added with a value that approximates the boundary conditions between the already existing cells and the new GC. Am I completely off-track with this line of thought?

viraj · February 5, 2010, 06:10

hi

I am working on the same method and implemented it for flow over a cylinder problem. I think you are not off the track. In ghost cell method the boundary doesn't coincide with the Cartesian mesh how to satisfy the boundary condition. Here we take help of ghost points. We interpolate the value at ghost point from neighboring fluid points such that boundary condition is satisfied. We dont solve for the solid cell which are completely inside the solid and doesn't have any fluid point as a neighbor.

I hope some doubts are cleared.

bye

January 29, 2010, 05:27	Immersed Boundary Ghost Cell Method	#1
thehexman New Member Join Date: Jan 2010 Posts: 3 Rep Power: 16	Hello all, I am currently writing a student research paper on the use of the IBM to model the human heart valves (a not quite new topic, I'm afraid, but in my department this hasn't been done yet). Currently I'm doing literature research and finding out about the various methods of the IB concept. At the moment, I am having trouble understanding the Ghost Cell Method (GCM). In particular, I do not exactly understand the interpolation. But first let me tell you what I understand so far: - A forcing term is added to the discrete system of governing equations - This forcing term is derived from the boundary conditions (e.g. no slip: u=0) - In general, the derivation of the forcing term is done by estimating the velocity field and and correcting it at the boundary to fit the boundary conditions. - Ghost cells are defined inside the immersed body, so that each GC has at least one neighbor in the fluid - If the artificial value for the flow parameters were known in these cells, the computation could stop right there w/o having to solve in the fluid - However, there is no information in the GC themselves, but nearby there is the boundary condition (as before, a no slip condition would be an example) - The value of the GC flow parameters can then be extrapolated (using a bilinear or trilinear scheme, depending on the case) from the boundary conditions and the fluid cells close by. My question is: what is the point of calculating the value of the GC flow parameters if the boundary condition is already known? Isn't it just more work to do so, because the flow parameters inside the body are purely artificial? I think this has something to do with the fact that the IB does not necessarily coincide with the mesh points and therefore one cannot impose boundary conditions on a discrete equation if these conditions do not lie on the mesh points. Therefore, another cell has to be added with a value that approximates the boundary conditions between the already existing cells and the new GC. Am I completely off-track with this line of thought?

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Immersed boundary method	sfar	Main CFD Forum	31	June 21, 2017 03:52
ghost fluid method (shock tube problem)	Amir	Main CFD Forum	0	March 1, 2009 06:16
ghost fluid method	ma	Main CFD Forum	1	October 6, 2005 13:26
Spectral Method Boundary condition	Deepak	Main CFD Forum	1	February 25, 2004 09:48
Immersed boundary boundary treatment	AFP	Main CFD Forum	0	January 30, 2004 07:16

February 5, 2010, 06:10		#2
viraj New Member viraj Join Date: Feb 2010 Posts: 2 Rep Power: 0	hi I am working on the same method and implemented it for flow over a cylinder problem. I think you are not off the track. In ghost cell method the boundary doesn't coincide with the Cartesian mesh how to satisfy the boundary condition. Here we take help of ghost points. We interpolate the value at ghost point from neighboring fluid points such that boundary condition is satisfied. We dont solve for the solid cell which are completely inside the solid and doesn't have any fluid point as a neighbor. I hope some doubts are cleared. bye