zhyang-dev/CFDEMcoupling-PFM
1
0
Fork 0
mirror of https://github.com/ParticulateFlow/CFDEMcoupling-PFM.git synced 2025-12-08 06:37:44 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

changed tutorial case for test

Browse Source
This commit is contained in:
ekinaci
2016-06-03 12:38:23 +02:00
parent 9aba825507
commit c9c0e49a0f
5 changed files with 29 additions and 66 deletions
Download Patch File Download Diff File Expand all files Collapse all files

6
tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/CFD/constant/couplingProperties
View File

@ -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

4
tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/CFD/constant/thermophysicalProperties
View File

@ -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;

4
tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/CFD/system/controlDict
View File

@ -23,9 +23,9 @@ startTime 0;
stopAt endTime;
endTime 0.2;
endTime 0.1;
deltaT 0.005;
deltaT 0.001;
writeControl adjustableRunTime;

34
tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/DEM/in.liggghts_init
View File

@ -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

47
tutorials/cfdemSolverRhoPimpleChem/simplifiedSiwek_in_CFDEM/DEM/in.liggghts_run
View File

@ -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