
[Sponsors] 
August 27, 2012, 07:48 
newbie question: finite difference, immersed boundary method

#1 
New Member
Peter Maday
Join Date: Aug 2011
Posts: 14
Rep Power: 7 
Dear forum members,
I would like to solve a PDE arising from a variational problem related to a mass conservation problem plus regularization terms. The preferred solution method is finite differences due to the simplicity of implementation. The application area essentially is the extraction of flow information from image sequences depicting fluid motion in curved tubular structures (vessels). The question I would like to address is how to handle the domain boundaries within a FD formulation appropriately where the boundary is given by a 2D curve that is not necessarily aligned to the uniform grid. So far I came across the Immersed Boundary method that as far as I understand introduces additional force terms to the PDE dependent on a set of arbitrarily placed points that describe the boundary. My question is essentially if the IB method is the proper approach to this problem, or are there other possibilities that work with FD and my be an alternative? Additionally if anyone knew about an example C/MATLAB code for IB, that would certainly be of help. In case, to make the problem simpler, and a boundary with edge aligned segments is considered is there a recommended way to handle the BCs? For rectangular boundaries the addition of "ghost points" (virtual domain points for which setting fixed values at each iteration depending on the values of neighboring cells at the previous step enable the handling of common BC types in an elegant manner) seems as a good solution, however for domains containing "step" boundary segments the approach cannot be directly applied, as a ghost cell may have multiple neighbours within the domain and thus its value cannot be set to satisfy multiple BCs at once. Is there an established way of using such ideas with non rectangular domains? Thanks, and looking forward to your feedback, Peter 

August 27, 2012, 08:01 

#2  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,797
Rep Power: 33 
Quote:
Hi, the IB technique is only one possibility among others... but you can use finite volume method that is very simple and can be generalized properly over volumes of general shape without using complicate mapping between physical and computational domain. As suggestion, you can read the book of Peric and Ferziger. 

August 27, 2012, 09:34 

#3 
Member
Join Date: Mar 2009
Posts: 72
Rep Power: 9 
If you have simple and smooth geometry a transfinite grid interpolation along with coordinate mapping works well with FD methods.
The book by Andersson describes this well and it is quite easy to implement. If the geometry becomes more complicated then there is a high risk that you won't be able to use this method however. Here unstructured FV methods might be better, but they are, imho, more difficult to code. gl hf
__________________
"Trying is the first step to failure."  Homer Simpson 

August 27, 2012, 16:00 

#4 
Senior Member
Join Date: Aug 2011
Posts: 251
Rep Power: 8 
The purpose of Immersed Boundary methods as its name indicates it (immersed) is to handle obtsacle within the domain and not really the boundaries of the domain,even if it is always possible. But I think in the foundation paper of Peskin he designed it to handle moving artery wall for blood flow.
But anyway according to your problem as it has been suggested by Filippo and Fordprefect, herearethe different options: 1) If you really want to use FD you have to use coordinate transformation (mapping) to switch from physical domaintocomputational domain 2) use finite volume designed for non orthogonal.the book of Peric is a must and you will find some fortran code to deal with such configuration. 3) you can use also Finite Element method,but muchmorecomplicated to implement,so I would advise you to avoid this option. As Filippo stated it above 2) is not necessarly more complicated than 1) but it is really more powerfull. So for me 2) wins!! Good luck 

August 28, 2012, 05:02 
FD vs FV

#5 
New Member
Peter Maday
Join Date: Aug 2011
Posts: 14
Rep Power: 7 
Thanks a lot for the replies from all of you!
Actually the question has arisen in the context of implementing an image processing algorithm that is based on a variational PDE approach where the original paper does not address the question with nonrectangular domains, however for practical applications with vessels it has to be considered while the paper only focuses on cases with straight tube arrangements. The preference for FD comes from the fact that the paper uses a finite difference discretization. The domain transformation approach mentioned in Chapter 5. of Andersson seems feasible to me, however to have a complete view of the possibilities I would like to ask if any of you knew a book / paper in which the handling of generic curved boundaries is addressed on the original FD domain by using more sophisticated handling of BCs. In 8.1.1 of the Peric book a method is briefly described where the boundary is approximated with a "stepwise" one that is aligned to the edges of the FD grid. Unfortunately no details are disclosed (only a reference is given to "Largeeddy simulation of turbulent boundary layer flow over a hemisphere" by M Manhart, H Wengle that I cannot access). Do you have any suggestion for such methods? Thanks again for your attention. 

August 28, 2012, 05:24 

#6  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,797
Rep Power: 33 
Quote:
As the IB method is concerned, I suggest these three papers: 1) M.D. de Tullio, A. Cristallo, E. Balaras, R. Verzicco, Direct numerical simulation of the pulsatile flow through an aortic bileaflet mechanical heart valve, Journal of Fluid Mechanics, Vol. 662, 2009, pp. 259290 (ISSN 00221120). 2) M.D. de Tullio, P. De Palma, G. Iaccarino, G. Pascazio, M. Napolitano, An immersed boundary method for compressible flows using local grid refinement, Journal of Computational Physics, Vol. 225, 2007, pp. 20982117 (ISSN 00219991). 3) http://www.annualreviews.org/doi/abs....061903.175743 and if you want reading a method mixing IB and FV on unstructured grid: http://onlinelibrary.wiley.com/doi/1...d.278/abstract Good work 

August 28, 2012, 05:27 

#7  
Member
Join Date: Mar 2009
Posts: 72
Rep Power: 9 
Quote:
Unfortunately I have my Ferziger/Peric book at home and can not check what you are referring to. Perhaps I misread you but do you want to implement stepwise boundaries? This is usually not preferred since your numerical surface normal is not aligned with your real surface normal. I guess it depends on what type of equations you are trying to solve however.
__________________
"Trying is the first step to failure."  Homer Simpson 

August 28, 2012, 05:31 

#8  
Senior Member
Join Date: Aug 2011
Posts: 251
Rep Power: 8 
Quote:
First you have to discretize your equations on the domain. 1) you choose to map your physical domain into a computational one where you will apply your FD method 2) you will manage your boundary conditions excatly how you use to do it. or 1) you choose FV method and then you grid your physical domain with non orthogonal control volumes and you discretize your equations on these non orthogonal volumes (see Peric) 2) you manage your boundary conditions how it should be with finite volume (see Peric). Avoid the use of stepwise boundaries. It is not how we should do CFD in 2013 

August 28, 2012, 07:18 

#9  
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 566
Rep Power: 13 
The issue is very complicated.
The reason as to why there is not much available on it is that there is no perfect way of dealing with general boundaries in immersed boundary codes. Having written immersed boundary solver for both cartesian grid and in general unstructured grids, I can tell you that you will not go far if you just follow what is mentioned in the books. It took me almost 8 months of effort to get drag and lifts to match experimental data for our calculations. In this I ended up doing things the way NOT mentioned in any literature. I won't be able to share what I exactly did to make it work because my last company is in process of patenting it. So until they are done with it, I can not tell you what it was. But here are few pointers for you: (1) If you have cartesian type mesh then try to find out extactly where in X,Y and Z axis the solid boundary lies. (2) Write the various descretization operators while accounting for position of boundaries. (3) if you can not do above two things then spend your time with Ghost cell methods and the idea behind it, that would give you the best results specially when things are moving and rotating. (4) You can apply the forcing in Matrices too by just modifying the source terms. I did that compared to calculating forces and distributing them etc etc. Read (5) as to why. (5) Think that the solid regions are the regions of the mesh where velocity field is known. Now how you force them is upto you. Ultimately what you are doing is that calculating the source terms for momentum equation such that when you solve them in matrix form the velocities that matrix solver would give is same as solid velocities in solid regions. (6) The major headache and time consuming part is solid marking in immersed boundaries. In this if things are moving then you will have to mark the solid and fluid regions of the mesh after every movement. This could be very problematic in non uniform cartesian meshes and in unstructured meshes. This typically involves solving nearest neighbor problems which are time consuming. The situation is grave in explicit codes where time step is small and courant dependent. (7) If used properly and immersed boundaries are very powerful tool. Quote:
Quote:


August 28, 2012, 08:55 

#10 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 2,797
Rep Power: 33 
I agree, no books or papers will give you the tricks to make an IB code properly working...


August 28, 2012, 10:21 

#11 
Senior Member
Join Date: Aug 2011
Posts: 251
Rep Power: 8 
We can not advise someone to use IB just to deal with his boundaries which are curved !
The only justification of using IB compared to body fitted grids is if we have moving boundaries or moving obstacles within the domain. In other cases using FV on arbitrary shape control volumes is much more powerful and accurate. Or at least use transform coordinate, mappe your domain and use FD if you want. But using setpwise boundaries combined with IB is definitely not the best option to deal with a boundary. By the way don't you know that most of IB methods are just first order? After this everyone can do what he wants in life 

March 18, 2014, 08:37 
Finite difference, Immersed boundary method

#12 
New Member
yanbin
Join Date: Apr 2013
Posts: 16
Rep Power: 5 
I am combining the Immersed Boundary Method and finite difference method. However, for the geometry corner, it is difficult to determine which velocity is used to interpolated, because u or v forcing are near the both two boundary line. Do you have the same experience about this problems?


March 18, 2014, 08:54 

#13  
New Member
yanbin
Join Date: Apr 2013
Posts: 16
Rep Power: 5 
I have met the similar problem. How do you solved that?
Quote:


March 18, 2014, 10:57 

#14 
Senior Member

The way to handle such situations is basically based on using multiple values. In practice, the cell geometry is fixed but for a given variable you will need as many internal values as many fluid cells will need them.
Consider for example the corner cells for a square cylinder with different temperatures on different sides. You can also have a fully structured grid with the cylinder surface aligned with the grid but, in the end, each corner solid cell has to store information for each of the neighboring cells. 

March 20, 2014, 03:58 
Ibm

#15  
New Member
yanbin
Join Date: Apr 2013
Posts: 16
Rep Power: 5 
But how do we use so much stored information? In other words, which measure do you take to handle the problem?
Quote:


March 20, 2014, 06:27 

#16 
Senior Member

From the practical point of view the answer is: I don't know, as i never programmed it.
In general, if i would like to make it completely transparent to the solver i would add an integer index for each multivalued ghost cell, each index pointing to the values for the given neighbor fluid cell. then i would treat such multiple values as different unknowns in the system. Then you would just need an if clause (based on a multivalue flag and the integer index) to fill the system matrix. Actually i don't know if the ghost cells are usually treated implicitly or not. If not then the previous procedure is even more trivial (i guess). However, i'm not an expert on immersed boundary. An article treating such multivalue method is the following: http://www.sciencedirect.com/science...45793007000394 which, by the way, might be different from my guess 

March 21, 2014, 03:59 

#17 
Senior Member

Consider also the following paper:
http://www.sciencedirect.com/science...21999108000235 I think its description of the overall immersed boundary method is made crystal clear, and it also describes the case of thin surfaces. I might not appreciate their somehow arbitrary treatment of sharpedges but, still, it seems to work. 

March 21, 2014, 04:37 
Ibm

#18  
New Member
yanbin
Join Date: Apr 2013
Posts: 16
Rep Power: 5 
Thanks a lot. I will try to read these paper
Quote:


May 12, 2014, 08:18 
Problems

#19 
New Member
yanbin
Join Date: Apr 2013
Posts: 16
Rep Power: 5 

May 12, 2014, 08:20 
IBM problems

#20  
New Member
yanbin
Join Date: Apr 2013
Posts: 16
Rep Power: 5 
Hi,
For instance, rectangular with bottom line connected to the bottom boundary, the forcing in horizontal line is very small and that in vertical line is much larger. Do you have similar experience ? Quote:


Tags 
ghost point, immersed boundary method 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Immersed boundary method  sfar  Main CFD Forum  29  June 20, 2015 13:06 
Immersed boundary method in OpenFoam  daihu000  OpenFOAM  10  December 22, 2014 15:50 
Question about immersed solid method in CFX  spwater  CFX  6  May 24, 2012 13:05 
finite difference method for navierstokes problem  dallybird  Main CFD Forum  5  February 17, 2003 23:00 
immersed boundary method with unstructured mesh  x.zhang  Main CFD Forum  0  August 25, 2002 21:57 