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

update tutorial case

Browse Source
This commit is contained in:
ekinaci
2017-06-21 11:45:12 +02:00
parent 347b83ee2a
commit e488ef863f
14 changed files with 134 additions and 81 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/constant/couplingProperties
View File

@ -165,7 +165,7 @@ DiFeliceDragProps
KochHillDragProps
{
//verbose true;
verbose true;
velFieldName "U";
voidfractionFieldName "voidfraction";
interpolation true;

2
tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/system/controlDict
View File

@ -93,5 +93,5 @@ functions
outputInterval 1;
}
};
}
// ************************************************************************* //

4
tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/system/controlDict.foam
View File

@ -52,7 +52,7 @@ maxCo 0.1;
libs ("libOpenFOAM.so" "libgroovyBC.so");
functions
(
{
probes
{
@ -113,5 +113,5 @@ functions
);
factor 1;
}*/
);
}
// ************************************************************************* //

2
tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/0/N
View File

@ -16,7 +16,7 @@ FoamFile
dimensions [1 -3 0 0 0 0 0];
internalField uniform 0;
internalField uniform 1;
boundaryField
{

4
tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/0/U
View File

@ -17,7 +17,7 @@ FoamFile
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0.0 0 0);
internalField uniform (0.01 0 0); //uniform (6.64e-7 0 0);
boundaryField
{
@ -37,7 +37,7 @@ boundaryField
inlet
{
type fixedValue;
value uniform (0.0 0 0);
value uniform (0.01 0 0); // uniform (6.64e-7 0 0);
}
outlet

2
tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/combustionProperties
View File

@ -15,7 +15,7 @@ FoamFile
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
combustionModel noCombustion<psiThermoCombustion>;//laminar<psiChemistryCombustion>;////noCombustion<psiThermoCombustion>;//PaSR<psiChemistryCombustion>;
combustionModel noCombustion<psiThermoCombustion>; //laminar<psiChemistryCombustion>;////noCombustion<psiThermoCombustion>;//PaSR<psiChemistryCombustion>;
active false; //true;

18
tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/couplingProperties
View File

@ -60,7 +60,7 @@ forceModels
//DiFeliceDrag
gradPForce
viscForce
//KochHillDrag
// KochHillDrag
//DEMbasedDrag
//RongDrag
//Archimedes
@ -114,11 +114,13 @@ speciesProps
partMoleName "partN";
pressureFieldName "p";
partPName "partP";
// verbose true;
}
diffusionCoefficientsProps
{
ChemistryFile "$casePath/CFD/constant/foam.inp";
verbose true;
ChemistryFile "$casePath/CFD/constant/foam.inp";
tempFieldName "T";
totalMoleFieldName "N";
pressureFieldName "p";
@ -130,13 +132,13 @@ diffusionCoefficientsProps
massTransferCoeffProps
{
velFieldName "U";
voidfractionFieldName "voidfraction";
partNuName "nu";
partReynolds "Rep";
velFieldName "U";
voidfractionFieldName "voidfraction";
partNu "partNu";
partReynolds "partRe";
verbose true;
}
SyamlalThermCondProps
{
voidfractionFieldName "voidfraction";
@ -243,7 +245,7 @@ DiFeliceDragProps
KochHillDragProps
{
//verbose true;
verbose true;
velFieldName "U";
voidfractionFieldName "voidfraction";
interpolation true;

77
tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/foam.dat
View File

@ -15,7 +15,80 @@ FoamFile
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
// Sutherland coefficients not valid must be checked before real cases
// for hConst
/*CO2
{
specie
{
nMoles 1;
molWeight 44.01;
}
thermodynamics
{
Cp 1290;
Hf 0;
}
transport
{
mu 1.37e-05;
Pr 0.72;
}
}
O2
{
specie
{
nMoles 1;
molWeight 31.9988;
}
thermodynamics
{
Cp 1120;
Hf 0;
}
transport
{
mu 1.95e-05;
Pr 0.7;
}
}
N2
{
specie
{
nMoles 1;
molWeight 28.0134;
}
thermodynamics
{
Cp 1212;
Hf 0;
}
transport
{
mu 1.66e-05;
Pr 0.7;
}
}
CO
{
specie
{
nMoles 1;
molWeight 28.01;
}
thermodynamics
{
Cp 1227;
Hf 0;
}
transport
{
mu 1.66e-05;
Pr 0.7;
}
}*/
CO2
{
specie
@ -103,7 +176,7 @@ CO
As 1.67212e-06;
Ts 170.672;
}
}
}

18
tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/thermophysicalProperties
View File

@ -17,18 +17,20 @@ FoamFile
thermoType
{
type hePsiThermo; //covers compressible and lagrangian solvers
mixture reactingMixture; //species and reactions are listed in chemistry file
transport sutherland; //calcualtes viscosity as a function of temperature from Stuherland coefficient As and Ts
thermo janaf;
energy sensibleEnthalpy; //uses enthaly in the solution
equationOfState perfectGas;
specie specie;
type hePsiThermo; //covers compressible and lagrangian solvers
mixture reactingMixture; //species and reactions are listed in chemistry file
//transport const; // assumes constant viscosity and Prandtl number
transport sutherland; //calcualtes viscosity as a function of temperature from Stuherland coefficient As and Ts
thermo janaf; //hConst; //janaf;
energy sensibleEnthalpy; //uses enthaly in the solution
//equationOfState rhoConst; // rho = const.
equationOfState perfectGas; // rho = 1/RT*p
specie specie;
}
chemistryReader foamChemistryReader;
foamChemistryThermoFile "$casePath/CFD/constant/foam.dat";
foamChemistryThermoFile "$casePath/CFD/constant/foam.dat" ;// "$casePath/CFD/constant/foam.dat";
foamChemistryFile "$casePath/CFD/constant/foam.inp";

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

@ -17,7 +17,7 @@ FoamFile
transportModel Newtonian;
nu nu [ 0 2 -1 0 0 0 0 ] 1.328e-05; // dnyamic viscosity of CO : 1.66*10^-5 [kg/ms] -- density: 1.25 kg/m3
nu nu [ 0 2 -1 0 0 0 0 ] 1.328e-05; // dnyamic viscosity of CO : 1.66*10^-5 [kg/ms] -- density: 1.25 kg/m3
// dynamic viscosity of O2 = 2,04-10^-5[kg/ms];
// density of O2 = 1.4290 [kg/m^3]
// kinematic viscosity of air (o2+n2) = 1.48e-05;

22
tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/turbulenceProperties
View File

@ -17,26 +17,4 @@ FoamFile
simulationType laminar;
/*laminar
{
RASModel laminar;
turbulence off;
printCoeffs off;
} */
/*
* simulationType RAS;
RAS
{
RASModel kOmegaSST;
turbulence on;
printCoeffs on;
} */
// ************************************************************************* //

6
tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/system/controlDict
View File

@ -23,11 +23,11 @@ startTime 0;
stopAt endTime;
endTime 3;
endTime 5;
deltaT 0.0001;
deltaT 0.1;
writeControl adjustableRunTime; // timeStep; // //runTime;//
writeControl adjustableRunTime; //runTime;//// timeStep; //
writeInterval 0.5;

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

@ -28,7 +28,7 @@ pair_style gran model hertz tangential history # Hertzian without cohesion
pair_coeff * *
# timestep, gravity
timestep 0.000001
timestep 0.001 #0.00001 #0.000001
fix gravi all gravity 0.0 vector 0.0 -1.0 0.0
# walls

54
tutorials/cfdemSolverRhoPimpleChem/test_case/DEM/in.liggghts_run
View File

@ -18,7 +18,6 @@ read_restart ../DEM/post/restart/liggghts.restart
neighbor 0.0005 bin
neigh_modify delay 0 binsize 0.01
# Material properties required for granular pair styles
fix m1 all property/global youngsModulus peratomtype 5.e6
@ -31,7 +30,7 @@ pair_style gran model hertz tangential history # Hertzian without cohesion
pair_coeff * *
# timestep, gravity
timestep 0.000001
timestep 0.001 #0.00001 #0.000001
fix gravi all gravity 0.0 vector 0.0 -1.0 0.0
# walls
@ -53,10 +52,9 @@ fix cfd3 all couple/cfd/chemistry n_species 4 species_names CO CO2 O2 N2
# 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 0.2
#fix cfd4 all chem/shrink speciesA O2 molMassA 31.99 speciesC CO2 molMassC 44.01 molMassB 12.01 k 2.5e3 rmin 0.005
#fix cfd4 all chem/shrink speciesA CO molMassA 28.01 speciesC CO2 molMassC 44.01 molMassB 12.01 k 2.5e3 rmin 0.005
# Activate for 3-layer unreacted core shrink model
#fix OreReductionCO all chem/shrink/core speciesA O2 molMassA 31.99 speciesC CO2 molMassC 44.01
fix cfd5 all chem/shrink/core speciesA CO molMassA 28.01 speciesC CO2 molMassC 44.01
# Chemical properties for unreacted shrink core (activate only when chem/shrink/core is active)
@ -66,7 +64,7 @@ fix Ea all property/global Ea_cfd5 vector 69488 73674 113859
# particle porosities adn tortuosity
fix porosity all property/global porosity_ scalar 0.15
fix tortuosity all property/global tortuosity_ scalar 3
fix pore_diameter all property/global pore_diameter_ scalar 3.91e-6
fix pore_diameter all property/global pore_diameter_ scalar 0.00001 #3.91e-9
# Material properties for unreacted chemical shrink core (activate only when chem/shrink/core is active)
fix density all property/global density_all vector 7870 5740 5170 5240
@ -80,38 +78,38 @@ fix integr all nve/sphere
###############################################
## DEBUG ##
compute masschange all property/atom mass
compute massreduce all reduce sum c_masschange
fix rmass all ave/time 100 1 100 c_massreduce file rmass.dat
#compute masschange all property/atom mass
#compute massreduce all reduce sum c_masschange
#fix rmass all ave/time 100 1 100 c_massreduce file rmass.dat
compute layerRad1 all reduce sum f_layerRelRad[1]
fix redRad1 all ave/time 100 1 100 c_layerRad1 file relRad1.dat
#compute layerRad1 all reduce sum f_layerRelRad[1]
#fix redRad1 all ave/time 100 1 100 c_layerRad1 file relRad1.dat
compute layerRad2 all reduce sum f_layerRelRad[2]
fix redRad2 all ave/time 100 1 100 c_layerRad2 file relRad2.dat
#compute layerRad2 all reduce sum f_layerRelRad[2]
#fix redRad2 all ave/time 100 1 100 c_layerRad2 file relRad2.dat
compute layerRad3 all reduce sum f_layerRelRad[3]
fix redRad3 all ave/time 100 1 100 c_layerRad3 file relRad3.dat
#compute layerRad3 all reduce sum f_layerRelRad[3]
#fix redRad3 all ave/time 100 1 100 c_layerRad3 file relRad3.dat
compute layerRad4 all reduce sum f_layerRelRad[4]
fix redRad4 all ave/time 100 1 100 c_layerRad4 file relRad4.dat
#compute layerRad4 all reduce sum f_layerRelRad[4]
#fix redRad4 all ave/time 100 1 100 c_layerRad4 file relRad4.dat
variable rad_fe atom "f_layerRelRad[1]*radius"
compute r_fe all reduce sum v_rad_fe
fix r_fe1 all ave/time 100 1 100 c_r_fe file rfe.dat
#variable rad_fe atom "f_layerRelRad[1]*radius"
#compute r_fe all reduce sum v_rad_fe
#fix r_fe1 all ave/time 100 1 100 c_r_fe file rfe.dat
variable rad_w atom "f_layerRelRad[2]*radius"
compute r_w all reduce sum v_rad_w
fix r_w1 all ave/time 100 1 100 c_r_w file rw.dat
#variable rad_w atom "f_layerRelRad[2]*radius"
#compute r_w all reduce sum v_rad_w
#fix r_w1 all ave/time 100 1 100 c_r_w file rw.dat
variable rad_m atom "f_layerRelRad[3]*radius"
compute r_m all reduce sum v_rad_m
fix r_m1 all ave/time 100 1 100 c_r_m file rm.dat
#variable rad_m atom "f_layerRelRad[3]*radius"
#compute r_m all reduce sum v_rad_m
#fix r_m1 all ave/time 100 1 100 c_r_m file rm.dat
variable rad_h atom "f_layerRelRad[4]*radius"
compute r_h all reduce sum v_rad_h
fix r_h1 all ave/time 100 1 100 c_r_h file rh.dat
#variable rad_h atom "f_layerRelRad[4]*radius"
#compute r_h all reduce sum v_rad_h
#fix r_h1 all ave/time 100 1 100 c_r_h file rh.dat
###############################################
# screen output