diff --git a/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/CFD/constant/couplingProperties b/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/CFD/constant/couplingProperties index 39c88ddd..3058f8ae 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/CFD/constant/couplingProperties +++ b/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/CFD/constant/couplingProperties @@ -31,7 +31,7 @@ syncMode false; modelType "A"; // A or B -couplingInterval 1000; +couplingInterval 1;//1000; voidFractionModel divided;//centre;// @@ -96,9 +96,9 @@ speciesProps { ChemistryFile "$casePath/CFD/constant/foam.inp"; tempFieldName "T"; - partTempName "Tp"; + partTempName "partTemp"; densityFieldName "rho"; - partRhoName "rhop"; + partRhoName "partRho"; } SyamlalThermCondProps diff --git a/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/CFD/constant/thermophysicalProperties b/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/CFD/constant/thermophysicalProperties index ce42a25e..eac1d19d 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/CFD/constant/thermophysicalProperties +++ b/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/CFD/constant/thermophysicalProperties @@ -28,9 +28,9 @@ thermoType chemistryReader foamChemistryReader; -foamChemistryThermoFile "$casePath/CFD/constant/foam.dat"; +foamChemistryThermoFile "$casePath/CFD/constant/foam.dat"; -foamChemistryFile "$casePath/CFD/constant/foam.inp"; +foamChemistryFile "$casePath/CFD/constant/foam.inp"; inertSpecie N2; diff --git a/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/CFD/system/controlDict b/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/CFD/system/controlDict index d55e964d..40b92a3f 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/CFD/system/controlDict +++ b/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/CFD/system/controlDict @@ -23,9 +23,9 @@ startTime 0; stopAt endTime; -endTime 0.2; +endTime 0.1; -deltaT 0.005; +deltaT 0.001; writeControl adjustableRunTime; diff --git a/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/DEM/in.liggghts_init b/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/DEM/in.liggghts_init index d3930d36..4d6181f1 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/DEM/in.liggghts_init +++ b/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/DEM/in.liggghts_init @@ -1,4 +1,4 @@ -# Pour granular particles into chute container, then induce flow +# Particle insertion into domain atom_style granular atom_modify map array echo both @@ -11,15 +11,13 @@ newton off units si processors 2 2 1 -region reg block 0. 0.1 0. 0.1 0. 0.05 units box +region reg block 0. 0.1 0. 0.1 0. 0.1 units box create_box 1 reg -neighbor 0.002 bin -neigh_modify delay 0 - +neighbor 0.01 bin +neigh_modify delay 0 # Material properties required for granular pair styles - fix m1 all property/global youngsModulus peratomtype 5.e6 fix m2 all property/global poissonsRatio peratomtype 0.45 fix m3 all property/global coefficientRestitution peratomtypepair 1 0.3 @@ -30,45 +28,37 @@ pair_style gran model hertz tangential history # Hertzian without cohesion pair_coeff * * # timestep, gravity -timestep 0.00001 +timestep 0.001 #0.00001 fix gravi all gravity 9.81 vector 0.0 0.0 -1.0 -# heat transfer (taken from packedBedTemp) -fix ftco all property/global thermalConductivity peratomtype 5. -fix ftca all property/global thermalCapacity peratomtype 10. -fix heattransfer all heat/gran initial_temperature 263. - # walls fix xwalls1 all wall/gran model hertz tangential history primitive type 1 xplane 0.0 fix xwalls2 all wall/gran model hertz tangential history primitive type 1 xplane 0.1 fix ywalls1 all wall/gran model hertz tangential history primitive type 1 yplane 0.0 fix ywalls2 all wall/gran model hertz tangential history primitive type 1 yplane 0.1 fix zwalls1 all wall/gran model hertz tangential history primitive type 1 zplane 0.0 -fix zwalls2 all wall/gran model hertz tangential history primitive type 1 zplane 0.05 +fix zwalls2 all wall/gran model hertz tangential history primitive type 1 zplane 0.1 # particle distributions and insertion -region bc block 0. 0.1 0. 0.1 0. 0.05 units box -fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 2000 radius constant 0.001 +region bc block 0. 0.1 0. 0.1 0. 0.1 units box +fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 2000 radius constant 0.01 fix pdd1 all particledistribution/discrete 1 1 pts1 1.0 -fix ins all insert/pack seed 100001 distributiontemplate pdd1 vel constant 0. 0. 0. insert_every 100 overlapcheck yes all_in yes particles_in_region 20 region bc +fix ins all insert/pack seed 1000001 distributiontemplate pdd1 vel constant 0. 0. 0. insert_every 100 overlapcheck yes all_in yes particles_in_region 5 region bc # apply nve integration to all particles that are inserted as single particles fix integr all nve/sphere # screen output compute rke all erotate/sphere -thermo_style custom step atoms ke c_rke vol +thermo_style custom step atoms ke c_rke vol thermo 1000 thermo_modify lost ignore norm no compute_modify thermo_temp dynamic yes # insert the first particles so that dump is not empty -run 1 -dump dmp all custom 50 post/dump.liggghts_init id type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius f_heattransfer[0] f_heatFlux[0] -#force : f_couple_cfd[0] f_couple_cfd[1] f_couple_cfd[2] -#node : f_couple_cfd[6] -#cell id : f_couple_cfd[7] +run 1 +dump dmp all custom 50 post/dump.liggghts_init id type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius run 10000 upto write_restart post/restart/liggghts.restart diff --git a/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/DEM/in.liggghts_run b/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/DEM/in.liggghts_run index 5c4d2561..2044b4c2 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/DEM/in.liggghts_run +++ b/tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/DEM/in.liggghts_run @@ -15,8 +15,9 @@ processors 2 2 1 # read the restart file read_restart ../DEM/post/restart/liggghts.restart -neighbor 0.003 bin -neigh_modify delay 0 binsize 0.01 +neighbor 0.01 bin +neigh_modify delay 0 binsize 0.0 ##0.01 + # Material properties required for granular pair styles @@ -30,7 +31,7 @@ pair_style gran model hertz tangential history # Hertzian without cohesion pair_coeff * * # timestep, gravity -timestep 0.00001 +timestep 0.001 #0.00001 fix gravi all gravity 9.81 vector 0.0 -1.0 0.0 # walls @@ -39,48 +40,20 @@ fix xwalls2 all wall/gran model hertz tangential history primitive type 1 xp fix ywalls1 all wall/gran model hertz tangential history primitive type 1 yplane 0.0 fix ywalls2 all wall/gran model hertz tangential history primitive type 1 yplane 0.1 fix zwalls1 all wall/gran model hertz tangential history primitive type 1 zplane 0.0 -fix zwalls2 all wall/gran model hertz tangential history primitive type 1 zplane 0.05 - -# heat transfer -fix ftco all property/global thermalConductivity peratomtype 5. # lambda in [W/(K*m)] -fix ftca all property/global thermalCapacity peratomtype 0.1 # cp in [J/(kg*K)] -fix heattransfer all heat/gran initial_temperature 600. - -# set particle temperature for the bed -run 0 -region total block INF INF INF INF INF INF units box -set region total property/atom Temp 600. - -# change the particles density -# set group all density 2000 +fix zwalls2 all wall/gran model hertz tangential history primitive type 1 zplane 0.1 # cfd coupling -fix cfd all couple/cfd couple_every 100 mpi +fix cfd all couple/cfd couple_every 1 mpi fix cfd2 all couple/cfd/force -# this one invokes heat transfer calculation, transfers per-particle temperature and adds convective heat flux to particles -fix cfd3 all couple/cfd/convection T0 600 - # this should invoke chemistry -fix cfd4 all couple/cfd/chemistry n_species 6 species_names N2 O2 CH4 H2 CO2 H2O +fix cfd3 all couple/cfd/chemistry notequalto_n_species should give error + +#n_species 6 species_names_ N2 O2 CH4 H2 CO2 H2O # apply nve integration to all particles that are inserted as single particles fix integr all nve/sphere -# cfd coupling -#fix cfd all couple/cfd couple_every 100 mpi -#fix cfd2 all couple/cfd/force/implicit -#fix cfd2 all couple/cfd/force/accumulator RongDrag 10 1.5e-3 -#fix cfd2 all couple/cfd/force/implicit/accumulated #CrankNicolson 0.5 -# apply nve integration to all particles that are inserted as single particles -#fix integr all nve/sphere - -# center of mass -#compute centerOfMass all com - -# compute total dragforce -#compute dragtotal all reduce sum f_dragforce[1] f_dragforce[2] f_dragforce[3] - # screen output compute rke all erotate/sphere thermo_style custom step atoms ke c_rke vol @@ -88,6 +61,6 @@ thermo 1000 thermo_modify lost ignore norm no compute_modify thermo_temp dynamic yes -dump dmp all custom 50 ../DEM/post/dump.liggghts_run id type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius f_Temp[0] f_heatFlux[0] +dump dmp all custom 50 ../DEM/post/dump.liggghts_run id type x y z ix iy iz vx vy vz fx fy fz omegax omegay omegaz radius run 1