CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   CD-adapco (http://www.cfd-online.com/Forums/cd-adapco/)
-   -   Extracting information on a particle track (http://www.cfd-online.com/Forums/cd-adapco/53552-extracting-information-particle-track.html)

Bo Jensen November 17, 2003 03:18

Extracting information on a particle track
 
Hi!

I would like to get the temperature in each cell that a particle track passes through. Is that possible?

Regards

Bo

Michiel November 19, 2003 08:20

Re: Extracting information on a particle track
 
Hi Bo,

Yes it is. The variable idc(idr) gives you the cell number a particle is in. You can print these data from the cell number. If idc.eq.0 the particle is stuck in the domain, If idc.lt.0 it has left the domain.

Michiel

This was sent to mee as an example in POSDAT.f by starcd

COMMON /C33001/ IDC

C --- DROPLET DATA

OPEN (UNIT=100,FILE='domain.out')

OPEN (UNIT=101,FILE='evap.out')

OPEN (UNIT=102,FILE='stuck.out')

REWIND (100)

REWIND (101)

REWIND (102)

MASSWO = 0.0

HEATWO = 0.0

DO IDR=1,NDR C --- DROPLET (IDR) IS LOCATED IN CELL (IC)

IC=IDC(IDR)

MASFRA = DRMAIC(IDR,1)/DRM(IDR)

IF (IC.GT.0) THEN

WRITE (100,911) IDR,IDT(IDR),DRD(IDR),DRM(IDR),DRT(IDR),MASFRA

ELSE IF (IC.EQ.0) THEN

WRITE (101,911) IDR,IDT(IDR),DRD(IDR),DRM(IDR),DRT(IDR),MASFRA

ELSE

WRITE (102,911) IDR,IDT(IDR),DRD(IDR),DRM(IDR),DRT(IDR),MASFRA

ENDIF

MASSWO = DRMAIC(IDR,1) + MASSWO

HEATWO = DRT(IDR)*

&(4181.72*DRMAIC(IDR,1)+1750*(DRM(IDR)-DRMAIC(IDR,1)))+HEATWO

ENDDO

CLOSE (100)

CLOSE (101)

CLOSE (102)


All times are GMT -4. The time now is 16:51.