Generator2p has been developped by A. Albrecht to initialize a number of particles in partial equilibrium with Homogeneous Isotropic Turbulence.

Input :

• turbulent field at NTMIX format (from THI3D), already partitioned in parralel case
• characteristics data file

Output :

• particle data file at NTMIX format.
• progress README output file

Internal

• You have to verify and modify the sources for Albrecht 2001 add on to particles restart file (sources terms, here set to 0)
• Npartmax, nprocmax, nsizema,x overlapmax are internal parameters for maximum memory allocation : to be extended if needed

• NEC SX5, Monoprocessor

• Homogeneous and isotropic turbulence
• Spatially randomly distribution of particles
• Possible polydispersion (NOT VALIDATED)
• Particles velocity  are initialized according to [P. Fevrier, 2000] as the sum of two parts:

Correlated part :

with the Stokes number :

Random part :


 
 

Characteristics data file :

Line number/data

1 ttime (has to be the same as in continuous field file)
2 NX (line 2->14 : same as TEST_IN of NTMIX)
3 NY
4 NZ
5 nprocx
6 nprocy
7 nprocz
8 overlap
9 xmin
10 xmax
11 ymin
12 ymax
13 zmin
14 zmax
15 ndroplet (number of droplets)
16 nd_tirage (has to be > to ndroplet*nprocx*nprocy*nprocz)
17 polydispersion ? (.true. or .false.)
18 variance (rd' **2 )
19 rd_init (radius)
20 td_init (temperature)
21 equi_drop (.true. <=> partial equilibrium initalization)
22 tr_d_init  (response time of droplet for partial equilibrium)
23 u_drop (fixed velocity, equi_drop=.false.)
24 v_drop (fixed velocity)
25 w_drop (fixed velocity)
26 nparticle  (number of particles)
27 np_tirage  (has to be > to nparticle*nprocx*nprocy*nprocz)
28 rp_init (radius)
29 tp_init (temperature)
30 equi_part (.true. <=> partial equilibrium initalization)
31 tr_p_init (response time of particle for partial equilibrium)
32 u_part (fixed velocity, equi_part=.false.)
33 v_part  (fixed velocity)
34 w_part  (fixed velocity)
35 tracers (.true. or .false.)
36 homo (0: aleatoire, 1:homogene)
37 gas restart file name (char*12, without proc number)
38 q11_fluid (gas turbulent kinetic energy)
39 t11_fluid (response time of fluid)

Note author : J. Hélie, IMFT, 2002

Back to script page

Back to main page