CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > Siemens

Boundary cell looping in posdat.f [v4.06]

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 4, 2008, 21:24
Default Boundary cell looping in posdat.f [v4.06]
  #1
hirokata
Guest
 
Posts: n/a
Hello all-

I need to identify neighboring cell pair, front and back, at a (baffle) boundary, in order to calculate flux across it. I am using the FsetBoun looping macro in posdat.f to get the cells attached to the boundary as below:

C************************************************* **********

if (level.eq.1) then

do ir = 1, regi_no

if (regi(ir)%type.eq.BAFF) then

call FsetBoun(fs,0,ND_ALL,NSD_ALL,ir,ISIDE_ALL)

do nset=1,fs%no

do j=fs%ns(nset), fs%ne(nset)

i1 = lfc(1,j)

i2 = lfc(2,j)

! Baffle 1

if (ir.eq.7) then

! some actions

end if

end do

end do

end if

end do

end if

C************************************************* **********

With the above code I can identify all the cells at the boundary. But what I really want to do is to identify the "pairing" cells across the boundary. I was hoping lfc(,) will accomplish this task, but it seems that, at a boundary, lfc(1,j) will return the neighboring cell number and lfc(2,j) will give me negative of the face number at the boundary. The version I am using is V4.06.

As I understand acquiring of pairing cells aross a boundary was possible in Version 3. I am hoping there is a simple way I can accomplish that that I'm not aware of. If any of you know a solution to this problem please let me know since I have been at this for months. Thank you very much in advance for your response.

hirokata
  Reply With Quote

Old   September 5, 2008, 10:22
Default Re: Boundary cell looping in posdat.f [v4.06] *NM*
  #2
Stephanie
Guest
 
Posts: n/a
  Reply With Quote

Old   September 5, 2008, 10:23
Default Re: Boundary cell looping in posdat.f [v4.06]
  #3
Stephanie
Guest
 
Posts: n/a
Are the cells on both side of the baffle belonging to different cell types. I would then use the function call cset and work with the subdomains.
  Reply With Quote

Old   September 5, 2008, 12:21
Default Re: Boundary cell looping in posdat.f [v4.06]
  #4
hirokata
Guest
 
Posts: n/a
Hi Stephanie,

Yes, each side of the baffle belongs to a different cell types. However, what I'm trying to do is to identify the pairing cells during the iterations.

hirokata
  Reply With Quote

Old   September 8, 2008, 03:55
Default Re: Boundary cell looping in posdat.f [v4.06]
  #5
Stephanie
Guest
 
Posts: n/a
well. I would find the pairs at the beginnning of the simulation, only once. I would save the pairs in table variables. If your pair cells belong to different types you can work with subdomains with call cset
  Reply With Quote

Old   September 14, 2008, 13:42
Default Re: Boundary cell looping in posdat.f [v4.06]
  #6
hirokata
Guest
 
Posts: n/a
Thank you Stephanie. I used the call cset and was able to find the pairing cells. I checked the table feature but it only seems to be capable of handling simplest functions. Do you happen to know the way to access 'cell diffusion coefficient' while looping the cells in the domain? I found "af(j)" but I believe this is the face diffusion coefficient. I appreciate for sharing that info if you happen to know it. Thanks!

hirokata
  Reply With Quote

Old   September 22, 2008, 15:34
Default Re: Boundary cell looping in posdat.f [v4.06]
  #7
Stephanie
Guest
 
Posts: n/a
sorry. I don't know how to get the cell diffusion coefficient. A call to Adapco support would certainly help. Or you should open a new forum message.
  Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Looping over nodes of a cell (UDF) error Manoj FLUENT 2 December 1, 2005 01:30
Access cell coordinates and wall data in posdat.f Boban Siemens 5 July 22, 2005 05:25
How can I know a cell is boundary in the posdat.f wy Siemens 1 August 9, 2002 03:53
Looping over faces in a cell Senthil FLUENT 3 April 11, 2002 12:41
May I use global coordinate of cell centers in subroutine Posdat H. Shen Siemens 4 March 29, 2000 10:20


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