CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   Phoenics (http://www.cfd-online.com/Forums/phoenics/)
-   -   How can I get the area of a cell? (http://www.cfd-online.com/Forums/phoenics/51497-how-can-i-get-area-cell.html)

 cloudz October 21, 2002 05:22

How can I get the area of a cell?

I want to calculate the airflow between zones. So the velocity and current cell area are need to know,how can I get it? thanks.

 Eugene October 21, 2002 10:45

Re: How can I get the area of a cell?

Try http://htex.mpei.ac.ru/~krinitsky/index1.html software. Load PHI file, enter to analyze mode and click mouse button on cell and you see velocity value and X,Y,Z cell size.

 cloudz October 22, 2002 09:09

Re: How can I get the area of a cell?

i mean,i want to calculate the airflow in ground.for. so i wonder whether there are any symbol used to represent the area of current cell(or the x,y,z extent of current cell)?

 Rami October 22, 2002 10:58

Re: How can I get the area of a cell?

Hi cloudz,

If you wish to access the EAST face (say) of all cells -

get the 0-location index by

l0a = L0B(ASURFE)

loop on the range of indices (e.g., ix1<= ix <=ix2, etc.), and retrieve the area fro the F array by:

do izz = iz1,iz2

jz = Nx*Ny*(izz-1)

do ix = ix1,ix2

do iy = iy1,iy2

icell = iy + Ny*(ix-1)

area = F(l0a + icell + jz)

enddo ! iy

enddo ! ix

enddo ! izz

Note, that part of the above variables (e.g., Nx, Ny etc.) are passed in commons included in some of d_includ/ folder.

I hope this helps. Rami

 cloudz October 23, 2002 09:55

Re: How can I get the area of a cell?

Rami,thank you very much. Can you explain more detailedly,what's the function L0B(ASURFE) meaning? I've looked for the polis,but just find the L0F(variable name) function subroutine.

 Rami October 23, 2002 10:12

Re: How can I get the area of a cell?

See explanation in POLIS encyclopaedia, entry F-ARRAY of EARTH, (k) Block-location Indices (LB's); BFC Geometrical Quantities, file:/phoenics/d_polis/d_enc/enc_fa11.htm

 All times are GMT -4. The time now is 13:32.