hy! is there any solution to identify the particles diameter for the ones that had disappeared from the domain, for the case of a poly-disperse Lagrangian injection?! thankx in advance! regards, lelia

For STAR-CCM+ use the incident mass field function and divide by particle mass.

For STAR-CD there is a small Fortran program which used to be on the adapco online web site that reads the .trk file and then you can do any type of analysis:

C----------------------------------------------------------------------C

Program PTread C----------------------------------------------------------------------C C C Read STAR Trackfile, and write new trackfiles with only specific C Tracks for post processing visualization. C C Note: Only works on Big-Endian Platforms (i.e. IBM RS6000, NOT Linux) C C----------------------------------------------------------------------C C

CHARACTER FID*4,IDENT*4,BTYPE*4,FTYPE*8,ANS*4,FILE33*80

LOGICAL ITRACK,LSUMM,LPRIN

INTEGER ZERO(11) C C----------------------------------------------------------------------C C C --- Open Files and write out header information

OPEN(UNIT=33,FILE='vat-1.trk',STATUS='OLD',

& FORM='UNFORMATTED') C C --- Read in and write out Header

READ(33) FID,JVERS,IDENT,ZERO

WRITE(*,*) FID,JVERS,IDENT,ZERO C----------------------------------------------------------------------C C

ii=1 C --- Read in track 5 READ(33,ERR=5,END=800)IPART,X,Y,Z,U,V,W,ITRACK,TIM E,ICEL,

& TEMP,DIAM,DMASS,COUNT

print *,ii,TIME

if(TIME.gt.2.79)then

print *,'Inside the time range'

write(*,*)

WRITE(36,*)IPART,TIME,X,Y,Z,ICEL

WRITE(37,*)IPART,TIME,TEMP,DIAM,DMASS,COUNT,COUNT

endif C----------------------------------------------------------------------C C

GO TO 5 800 CONTINUE

CLOSE(33)

CLOSE(36)

CLOSE(37) C C----------------------------------------------------------------------C C 111 stop

END

Dear Ping,

I tried to use the ccm+ and incident mass field function, but apparently I am getting 20 times more mass coming out then injected (tested on a single straight canal). have you encounter such a problem? I cannot get a balance of what I inject and what the ccm count at the outlet.

any hint welcome.

matej