CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

How to deal with "dummy" cells by FVM

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   November 12, 2002, 19:35
Default How to deal with "dummy" cells by FVM
  #1
Smith Gearheart
Guest
 
Posts: n/a
hi, there

I am using structured grid finite volume method to compute for a case. My problem is that there will be some cells that is occupied with solid. I think these cells definitely should not be included in the computation, so I called them "dummy" cells.(please correct me if the name is not right). However, since I am using an Euler solver, it seems there will be flux into these cells. For example, my domain is something like below(sorry for the bad drawing, maybe here should let us upload some pictures)

------------------------------- |_|_|_|_|_|_|_|_|_|_|_|_|_|_|_| |_|_|_|_|_|_|_|_|_|_|_|_|_|_|_| |_|_|_|_|_|_|_|_|_|_|_|_|*|*|*|<--solid cells |_|_|_|_|_|_|_|_|_|_|_|_|*|*|*|<--solid cells |_|_|_|_|_|_|_|_|_|_|_|_|_|_|_| | | | | | | | | | | | | | | | | -------------------------------

Thanks for advice.
  Reply With Quote

Old   November 12, 2002, 19:38
Default Re: How to deal with "dummy" cells by FVM
  #2
Smith Gearheart
Guest
 
Posts: n/a
sorry for the drawing,I should preview them. it is something like: for example: 4x4 grid, with (2,3)(2,4)(3,3)(3,4) are solids.
  Reply With Quote

Old   November 13, 2002, 07:49
Default Re: How to deal with "dummy" cells by FVM
  #3
onno
Guest
 
Posts: n/a
There are several techniques. One general method would be a technique very similar to that used multi-phase flows by introducing a solid/fluid fraction say alpha. Alpha for the solid (solid fraction) indicates which part of the control volume is occupied by the solid: thus a value of one indicate completely occupied by a solid, a value of zero that the control volume is completely occupied by the fluid. Values in between indicates that percentage occupied by the solid. Thus the (fluid fraction) = 1 (solid fraction) The conservation equations used is very much similar to the multi-phase flow equations, all the terms contains the fluid fraction. When discretised with the finite volume method, the effect of Alpha*Volume indicates that part of the control volume that is available for the fluid. Thus a cell containing a solid will be zero * Volume. A subsequent issue is that the fluid fraction at the control volume faces should be defined as well. Thus, a solid covering 30% of a control volume face makes only 70% flow region available. According to your drawing (if I understand your "bad drawing" correctly) the faces values of your control volumes will either be one or zero. Once this solid fraction starts to move you have several additional problems that could be addressed quite efficiently.

Another very simple technique is to hard code the fluxes at solid control volume faces to be zero and for the purpose of the solver to make the matrix central coefficient of the control volumes containing the solids very large, say E+10. The inverse E-10 is very small and the result of the solver is an unchanged value for the control volumes containing the solids.

Hope this help
  Reply With Quote

Old   November 13, 2002, 10:39
Default Re: How to deal with "dummy" cells by FVM
  #4
Jim Park
Guest
 
Posts: n/a
Don't believe that there's an universally accepted terminology for these cells. I've seen "solid". The Los Alamos folks used "obstacle" for years.

The ring of cells outside your active region that can be used to set boundary condition has - in my experience - been called "ghost", "imaginary", "dummy", "image" and probably some others I've forgotten. So "dummy" could lead to some confusion if it's published.

The condition you need to impose is of course that there's no flow through the face of the obstacle cell. I think a previous post has given you a good solution for that.
  Reply With Quote

Old   November 14, 2002, 03:40
Default Re: How to deal with "dummy" cells by FVM
  #5
versi
Guest
 
Posts: n/a
My experience, as the two above, is that you should set the numerical flux at the wall surface to original physical flux. Ignore all the contributions from the ghost cells, else you will introduce artificial flux into the flow domain.
  Reply With Quote

Old   November 14, 2002, 16:05
Default Re: How to deal with "dummy" cells by FVM
  #6
Smith Gearheart
Guest
 
Posts: n/a
to me, it seems to set zero flux by incorperate huge coefficient is a simple way. Thank you very much. But can you tell me how the boundary conditions on the solid surface are applied? Since I am solving Euler equation, a zero velocity at surface is not correct. thanks again.
  Reply With Quote

Old   November 14, 2002, 16:09
Default thanks, Jim and versi
  #7
Smith Gearheart
Guest
 
Posts: n/a
  Reply With Quote

Old   November 15, 2002, 01:35
Default Re: How to deal with "dummy" cells by FVM
  #8
onno
Guest
 
Posts: n/a
Apart from a possible pressure extrapolation problem where special care needs to be taken, should boundary conditions be catered for in a natural way. Unless you want to apply wall functions for turbulence modelling at these solid cell interfaces. Here some research needs to be done, especially if the cells are partially occupied by a solid.

Thus, if you need wall functions at these solid cell interfaces, you are in for some research and testing (maybe a PhD or two). Here, a faster remedy might be to adjust your coding in such a manner that you can solve for block structured meshes. Alternatively, you should use a very fine mesh near the solid cells.

GOODLUCK
  Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
SnappyHexMesh for internal Flow vishwa OpenFOAM Native Meshers: snappyHexMesh and Others 23 August 6, 2014 03:50
Import netgen mesh to OpenFOAM hsieh Open Source Meshers: Gmsh, Netgen, CGNS, ... 32 September 13, 2011 05:50
snappyHexMesh won't work - zeros everywhere! sc298 OpenFOAM Native Meshers: snappyHexMesh and Others 2 March 27, 2011 21:11
snappyHexMesh aborting Tobi OpenFOAM Native Meshers: snappyHexMesh and Others 0 November 10, 2010 04:23
physical boundary error!! kris CD-adapco 2 August 3, 2005 00:32


All times are GMT -4. The time now is 17:15.