#Contact model atom_style granular atom_modify map array boundary f f f newton off communicate single vel yes processors 1 4 1 units si region reg block -5.000 3.200 -3.500 6.0000 -3.000 3.000 units box create_box 17 reg neighbor 0.02 bin neigh_modify delay 0 #Material properties required for new pair styles fix m1 all property/global youngsModulus peratomtype 6000000 6000000 6000000 6000000 6000000 6000000 6000000 6000000 6000000 6000000 6000000 6000000 6000000 6000000 6000000 6000000 999000000 fix m2 all property/global poissonsRatio peratomtype .2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 0.2000 fix m3 all property/global coefficientRestitution peratomtypepair 17 .4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 0.4000 fix m4 all property/global coefficientFriction peratomtypepair 17 .5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000 fix m5 all property/global coefficientRollingFriction peratomtypepair 17 .0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 0.0200 fix 11 all property/global minSeparationDistanceRatio scalar 1.01 1.01 1.01 1.01 1.01 1.01 1.01 1.01 1.01 1.01 1.01 1.01 1.01 1.01 1.01 1.01 #(value=value for the minimum separation distance, recommended as 1.01) fix 12 all property/global maxSeparationDistanceRatio scalar 1.1 1.1 1.1 1.1 1.1 1.1 1.1 1.1 1.1 1.1 1.1 1.1 1.1 1.1 1.1 1.1 #(value=value for the maximum separation distance, recommended as 1.1) fix 21 all property/global surfaceLiquidContentInitial scalar 22 20 25 5 10 38 4 15 30 40 30 60 28 32 44 55 #(value=value for the initial surface liquid volume in % of the solid volume) fix 22 all property/global surfaceTension scalar 0.07275 0.07275 0.07275 0.07275 0.07275 0.07275 0.07275 0.07275 0.07275 0.07275 0.07275 0.07275 0.07275 0.07275 0.07275 0.07275 #(value=value for the surface tension of liquid (SI units N/m)) fix 31 all property/global fluidViscosity scalar 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 0.1 #(value=value for the fluidViscosity (SI units Pas)) fix 32 all property/global contactAngle peratomtype 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 15 #(value_i=value for contact angle of atom type i and fluid) #New pair style pair_style gran model hertz tangential history cohesion washino/capillary/viscous #Hertzian with cohesion pair_coeff * * timestep 0.000005 #0.0000005 fix gravi all gravity 9.81 vector -1.0 0.0 0.0 fix CercleE all mesh/surface file Cercle.Entree.stl type 17 scale 0.001 rotate axis 0 0 1 angle 3.57 fix Cylindre all mesh/surface file cylindre5.stl type 17 scale 0.001 rotate axis 0 0 1 angle 3.57 fix Sortie all mesh/surface file Sortie5.stl type 17 scale 0.001 rotate axis 0 0 1 angle 3.57 curvature 1e-5 fix Entree all mesh/surface file EN.stl type 17 scale 0.001 rotate axis 0 0 1 angle 3.57 curvature 1e-5 fix in all mesh/surface file in.stl type 17 scale 0.001 rotate axis 0 0 1 angle 3.57 curvature 1e-5 #fix spiral all mesh/surface file spiral.stl type 17 scale 0.001 rotate axis 0 0 1 angle 3.57 fix granwalls all wall/gran model hertz tangential history mesh n_meshes 5 meshes Entree in CercleE Cylindre Sortie #spiral #Particle specification fix pts1 type1_ particletemplate/sphere 15485863 atom_type 1 density constant 950 radius constant 0.001000000 fix pts2 type2_ particletemplate/sphere 15485867 atom_type 2 density constant 950 radius constant 0.001133333 fix pts3 type3_ particletemplate/sphere 32452843 atom_type 3 density constant 950 radius constant 0.001266667 fix pts4 type4_ particletemplate/sphere 49979687 atom_type 4 density constant 950 radius constant 0.001400000 fix pts5 type5_ particletemplate/sphere 49979693 atom_type 5 density constant 950 radius constant 0.001533333 fix pts6 type6_ particletemplate/sphere 67867967 atom_type 6 density constant 950 radius constant 0.002666667 fix pts7 type7_ particletemplate/sphere 67867979 atom_type 7 density constant 950 radius constant 0.002800000 fix pts8 type8_ particletemplate/sphere 86028121 atom_type 8 density constant 950 radius constant 0.003033333 fix pts9 type9_ particletemplate/sphere 865877 atom_type 9 density constant 950 radius constant 0.0036066667 fix pts10 type10_ particletemplate/sphere 887819 atom_type 10 density constant 950 radius constant 0.0042200000 fix pts11 type11_ particletemplate/sphere 887827 atom_type 11 density constant 950 radius constant 0.004333333 fix pts12 type12_ particletemplate/sphere 887837 atom_type 12 density constant 950 radius constant 0.004466667 fix pts13 type13_ particletemplate/sphere 887839 atom_type 13 density constant 950 radius constant 0.004600000 fix pts14 type14_ particletemplate/sphere 887849 atom_type 14 density constant 950 radius constant 0.005033333 fix pts15 type15_ particletemplate/sphere 887867 atom_type 15 density constant 950 radius constant 0.00566667 fix pts16 type16_ particletemplate/sphere 888011 atom_type 16 density constant 950 radius constant 0.006000000 fix ins_mesh all mesh/surface/planar file Insertion.stl type 17 scale 0.001 rotate axis 0 0 1 angle 3.57 #Particle distribution for insertion fix pdd1 all particledistribution/discrete 89783 16 pts1 0.0625 pts2 0.0625 pts3 0.0625 pts4 0.0625 pts5 0.0625 pts6 0.0625 pts7 0.0625 pts8 0.0625 pts9 0.0625 pts10 0.0625 pts11 0.0625 pts12 0.0625 pts13 0.0625 pts14 0.0625 pts15 0.0625 pts16 0.0625 fix ins all insert/stream seed 86028157 distributiontemplate pdd1 & nparticles 15000 particlerate 2500 insert_every 10000 overlapcheck yes all_in no vel constant -1 0 0 & insertion_face ins_mesh #extrude_length 0.3 #************* Dump ************************************ run 1 dump myDump2 all local 10000 post1/coll*.csv f_surfaceLiquidContent[0] f_liquidFlux[0] f_liquidSource[0] #******** insert particles ****************** run 1000001 upto unfix ins run 1000000