CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   Siemens (https://www.cfd-online.com/Forums/siemens/)
-   -   poly-disperse lagrangian injection! (https://www.cfd-online.com/Forums/siemens/56579-poly-disperse-lagrangian-injection.html)

lelia April 16, 2008 05:25

poly-disperse lagrangian injection!
 
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

ping April 17, 2008 01:25

Re: poly-disperse lagrangian injection!
 
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

matej May 2, 2008 04:48

Re: poly-disperse lagrangian injection!
 
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


All times are GMT -4. The time now is 17:49.