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

latest working tutorial case

Browse Source
This commit is contained in:
ekinaci
2016-09-12 13:49:52 +02:00
parent 8e5c2df15d
commit 9439089c05
24 changed files with 75 additions and 72 deletions
Download Patch File Download Diff File Expand all files Collapse all files

1
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/Allclean
View File

@ -10,6 +10,7 @@ casePath="$(dirname "$(readlink -f ${BASH_SOURCE[0]})")"
cd $casePath/CFD
cleanCase
rm $casePath/CFD/shrink*.dat
rm $casePath/CFD/rmass*.dat
rm $casePath/log*
rm $casePath/run_error.log
#cp -R 0/Org/p 0/p

4
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/Allrun.sh
View File

@ -27,5 +27,5 @@ else
fi
#- run parallel CFD-DEM in new terminal
#gnome-terminal --title='cfdemSolverPiso ErgunTestMPI CFD' -e "bash $casePath/parCFDDEMrun.sh"
bash $casePath/parCFDDEMrun.sh
gnome-terminal -e "bash $casePath/parCFDDEMrun.sh"
#bash $casePath/parCFDDEMrun.sh

6
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/CO
View File

@ -23,17 +23,17 @@ boundaryField
{
top
{
type slip;
type zeroGradient;
}
bottom
{
type slip;
type zeroGradient;
}
side-walls
{
type slip;
type zeroGradient;
}
inlet

6
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/CO2
View File

@ -23,17 +23,17 @@ boundaryField
{
top
{
type slip;
type zeroGradient;
}
bottom
{
type slip;
type zeroGradient;
}
side-walls
{
type slip;
type zeroGradient;
}
inlet

6
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/Ksl
View File

@ -23,17 +23,17 @@ boundaryField
{
top
{
type slip;
type zeroGradient;
}
bottom
{
type slip;
type zeroGradient;
}
side-walls
{
type slip;
type zeroGradient;
}
inlet

6
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/O2
View File

@ -23,17 +23,17 @@ boundaryField
{
top
{
type slip;
type zeroGradient;
}
bottom
{
type slip;
type zeroGradient;
}
side-walls
{
type slip;
type zeroGradient;
}
inlet

6
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/Qsource
View File

@ -23,17 +23,17 @@ boundaryField
{
top
{
type slip;
type zeroGradient;
}
bottom
{
type slip;
type zeroGradient;
}
side-walls
{
type slip;
type zeroGradient;
}
inlet

8
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/T
View File

@ -17,23 +17,23 @@ FoamFile
dimensions [0 0 0 1 0 0 0];
internalField uniform 1073; //Taken from paper of Smith, I.W. Nineteenth Symposium (International) on Combustion, 1982, p. 1045
internalField uniform 293.15;
boundaryField
{
top
{
type slip;
type zeroGradient;
}
bottom
{
type slip;
type zeroGradient;
}
side-walls
{
type slip;
type zeroGradient;
}
inlet

10
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/U
View File

@ -17,27 +17,27 @@ FoamFile
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0.0218 0 0);
internalField uniform (0.0204 0 0);
boundaryField
{
// Walls in no slip condition
top
{
type slip;
type zeroGradient;
}
bottom
{
type slip;
type zeroGradient;
}
side-walls
{
type slip;
type zeroGradient;
}
inlet
{
type fixedValue;
value uniform (0.0218 0 0);
value uniform (0.0204 0 0);
}
outlet

6
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/Us
View File

@ -23,15 +23,15 @@ boundaryField
{
top
{
type slip;
type zeroGradient;
}
bottom
{
type slip;
type zeroGradient;
}
side-walls
{
type slip;
type zeroGradient;
}
inlet

6
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/alphat
View File

@ -23,15 +23,15 @@ boundaryField
{
top
{
type slip;
type zeroGradient;
}
bottom
{
type slip;
type zeroGradient;
}
side-walls
{
type slip;
type zeroGradient;
}
inlet
{

6
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/epsilon
View File

@ -23,15 +23,15 @@ boundaryField
{
top
{
type slip;
type zeroGradient;
}
bottom
{
type slip;
type zeroGradient;
}
side-walls
{
type slip;
type zeroGradient;
}
inlet
{

6
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/k
View File

@ -23,15 +23,15 @@ boundaryField
{
top
{
type slip;
type zeroGradient;
}
bottom
{
type slip;
type zeroGradient;
}
side-walls
{
type slip;
type zeroGradient;
}
inlet
{

6
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/mut
View File

@ -23,15 +23,15 @@ boundaryField
{
top
{
type slip;
type zeroGradient;
}
bottom
{
type slip;
type zeroGradient;
}
side-walls
{
type slip;
type zeroGradient;
}
inlet
{

6
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/nut
View File

@ -23,17 +23,17 @@ boundaryField
{
top
{
type slip;
type zeroGradient;
}
bottom
{
type slip;
type zeroGradient;
}
side-walls
{
type slip;
type zeroGradient;
}
inlet

9
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/p
View File

@ -23,15 +23,15 @@ boundaryField
{
top
{
type slip;
type zeroGradient;
}
bottom
{
type slip;
type zeroGradient;
}
side-walls
{
type slip;
type zeroGradient;
}
inlet
{
@ -40,7 +40,8 @@ boundaryField
outlet
{
type zeroGradient;
type fixedValue;
value uniform 101325;
}
}

11
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/p.org
View File

@ -23,27 +23,28 @@ boundaryField
{
top
{
type slip;
type zeroGradient;
}
bottom
{
type slip;
type zeroGradient;
}
side-walls
{
type slip;
type zeroGradient;
}
inlet
{
type slip;
type zeroGradient;
}
outlet
{
type slip;
type fixedValue;
value uniform 101325;
}
}

6
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/rho
View File

@ -22,17 +22,17 @@ boundaryField
{
top
{
type slip;
type zeroGradient;
}
bottom
{
type slip;
type zeroGradient;
}
side-walls
{
type slip;
type zeroGradient;
}
inlet

6
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/thDiff
View File

@ -24,17 +24,17 @@ boundaryField
{
top
{
type slip;
type zeroGradient;
}
bottom
{
type slip;
type zeroGradient;
}
side-walls
{
type slip;
type zeroGradient;
}
inlet

6
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/0/voidfraction
View File

@ -23,17 +23,17 @@ boundaryField
{
top
{
type slip;
type zeroGradient;
}
bottom
{
type slip;
type zeroGradient;
}
side-walls
{
type slip;
type zeroGradient;
}
inlet

2
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/CFD/constant/transportProperties
View File

@ -17,7 +17,7 @@ FoamFile
transportModel Newtonian;
nu nu [ 0 2 -1 0 0 0 0 ] 1.5327e-05; // ḱinematic viscosity of O2 = 1.525e-05; // kinematic viscosity of air (o2+n2) = 1.48e-05;
nu nu [ 0 2 -1 0 0 0 0 ] 1.4276e-05; // ḱinematic viscosity of O2 = 1.525e-05; // kinematic viscosity of air (o2+n2) = 1.48e-05;
kf kf [ 1 1 -3 -1 0 0 0 ] 0; // comes from energy model

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

@ -23,13 +23,13 @@ startTime 0;
stopAt endTime;
endTime 10;
endTime 20;
deltaT 0.0001;
writeControl adjustableRunTime;
writeInterval 0.25;
writeInterval 0.50;
purgeWrite 0;

2
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_init
View File

@ -52,7 +52,7 @@ region bc block 0. 1.5 0. 0.5 0. 0.5 units box
fix pts1 all particletemplate/sphere 1 atom_type 1 density constant 2267 radius constant 0.035
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 once overlapcheck yes all_in yes particles_in_region 100 region bc
fix ins all insert/pack seed 100001 distributiontemplate pdd1 vel constant 0. 0. 0. insert_every once overlapcheck yes all_in yes particles_in_region 1 region bc
# apply nve integration to all particles that are inserted as single particles
#fix integr all nve/sphere

12
tutorials/cfdemSolverRhoPimpleChem/Particle_in_Duct/DEM/in.liggghts_run
View File

@ -51,7 +51,7 @@ fix cfd3 all couple/cfd/chemistry n_species 2 species_names O2 CO2
# this should shrink the particle
#fix cfd4 all chem/shrink speciesA O2 molMassA 31.99 speciesC CO2 molMassC 44.01 molMassB 12.01 k 2.5e3 rdef hertzpct 0.2
fix cfd4 all chem/shrink speciesA O2 molMassA 31.99 speciesC CO2 molMassC 44.01 molMassB 12.01 k 25 rmin 0.005
fix cfd4 all chem/shrink speciesA O2 molMassA 31.99 speciesC CO2 molMassC 44.01 molMassB 12.01 k 2.5e3 rmin 0.005
# apply nve integration to all particles that are inserted as single particles
fix integr all nve/sphere
@ -59,15 +59,15 @@ fix integr all nve/sphere
# print radius change according to timestep
compute radChange all property/atom radius
compute radreduce all reduce sum c_radChange
fix shrink all ave/time 50 1 50 c_radreduce file shrink.dat
fix shrink all ave/time 25 1 25 c_radreduce file shrink.dat
compute radChange1 all property/atom radius
compute radreduce1 all reduce max c_radChange1
fix shrink1 all ave/time 50 1 50 c_radreduce1 file shrink-1.dat
fix shrink1 all ave/time 25 1 25 c_radreduce1 file shrink-1.dat
compute radChange2 all property/atom radius
compute radreduce2 all reduce min c_radChange2
fix shrink2 all ave/time 50 1 50 c_radreduce2 file shrink-2.dat
compute masschange all property/atom mass
compute massreduce all reduce sum c_masschange
fix rmass all ave/time 25 1 25 c_massreduce file rmass.dat
# screen output
compute rke all erotate/sphere