|
[Sponsors] |
September 4, 2008, 21:24 |
Boundary cell looping in posdat.f [v4.06]
|
#1 |
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 |
|
September 5, 2008, 10:22 |
Re: Boundary cell looping in posdat.f [v4.06] *NM*
|
#2 |
Guest
Posts: n/a
|
|
|
September 5, 2008, 10:23 |
Re: Boundary cell looping in posdat.f [v4.06]
|
#3 |
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.
|
|
September 5, 2008, 12:21 |
Re: Boundary cell looping in posdat.f [v4.06]
|
#4 |
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 |
|
September 8, 2008, 03:55 |
Re: Boundary cell looping in posdat.f [v4.06]
|
#5 |
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
|
|
September 14, 2008, 13:42 |
Re: Boundary cell looping in posdat.f [v4.06]
|
#6 |
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 |
|
September 22, 2008, 15:34 |
Re: Boundary cell looping in posdat.f [v4.06]
|
#7 |
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.
|
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |