diff --git a/tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/constant/couplingProperties b/tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/constant/couplingProperties index 01def497..a2033d83 100644 --- a/tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/constant/couplingProperties +++ b/tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/constant/couplingProperties @@ -165,7 +165,7 @@ DiFeliceDragProps KochHillDragProps { - //verbose true; + verbose true; velFieldName "U"; voidfractionFieldName "voidfraction"; interpolation true; diff --git a/tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/system/controlDict b/tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/system/controlDict index 05af8ef1..31108b0e 100644 --- a/tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/system/controlDict +++ b/tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/system/controlDict @@ -93,5 +93,5 @@ functions outputInterval 1; } -}; +} // ************************************************************************* // diff --git a/tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/system/controlDict.foam b/tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/system/controlDict.foam index 7493598b..59cbc811 100644 --- a/tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/system/controlDict.foam +++ b/tutorials/cfdemSolverPiso/ErgunTestMPI/CFD/system/controlDict.foam @@ -52,7 +52,7 @@ maxCo 0.1; libs ("libOpenFOAM.so" "libgroovyBC.so"); functions -( +{ probes { @@ -113,5 +113,5 @@ functions ); factor 1; }*/ -); +} // ************************************************************************* // diff --git a/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/0/N b/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/0/N index 185a7387..5ba28508 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/0/N +++ b/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/0/N @@ -16,7 +16,7 @@ FoamFile dimensions [1 -3 0 0 0 0 0]; -internalField uniform 0; +internalField uniform 1; boundaryField { diff --git a/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/0/U b/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/0/U index aa546d22..55eb783a 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/0/U +++ b/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/0/U @@ -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 diff --git a/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/combustionProperties b/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/combustionProperties index 593fa936..d6fb4041 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/combustionProperties +++ b/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/combustionProperties @@ -15,7 +15,7 @@ FoamFile } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -combustionModel noCombustion;//laminar;////noCombustion;//PaSR; +combustionModel noCombustion; //laminar;////noCombustion;//PaSR; active false; //true; diff --git a/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/couplingProperties b/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/couplingProperties index 4db72597..5f7a73a0 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/couplingProperties +++ b/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/couplingProperties @@ -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; diff --git a/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/foam.dat b/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/foam.dat index 527eb2be..71fcec14 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/foam.dat +++ b/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/foam.dat @@ -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; } -} +} diff --git a/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/thermophysicalProperties b/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/thermophysicalProperties index 84a8434b..17068d88 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/thermophysicalProperties +++ b/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/thermophysicalProperties @@ -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"; diff --git a/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/transportProperties b/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/transportProperties index b8605f77..3851e205 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/transportProperties +++ b/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/transportProperties @@ -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; diff --git a/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/turbulenceProperties b/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/turbulenceProperties index e8fa087c..d1156961 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/turbulenceProperties +++ b/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/constant/turbulenceProperties @@ -17,26 +17,4 @@ FoamFile simulationType laminar; - /*laminar -{ - RASModel laminar; - - turbulence off; - - printCoeffs off; -} */ - -/* - * simulationType RAS; - -RAS -{ - RASModel kOmegaSST; - - turbulence on; - - printCoeffs on; -} */ - - // ************************************************************************* // diff --git a/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/system/controlDict b/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/system/controlDict index cc8daa96..82d6e262 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/system/controlDict +++ b/tutorials/cfdemSolverRhoPimpleChem/test_case/CFD/system/controlDict @@ -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; diff --git a/tutorials/cfdemSolverRhoPimpleChem/test_case/DEM/in.liggghts_init b/tutorials/cfdemSolverRhoPimpleChem/test_case/DEM/in.liggghts_init index 22f12b8a..60562d08 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/test_case/DEM/in.liggghts_init +++ b/tutorials/cfdemSolverRhoPimpleChem/test_case/DEM/in.liggghts_init @@ -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 diff --git a/tutorials/cfdemSolverRhoPimpleChem/test_case/DEM/in.liggghts_run b/tutorials/cfdemSolverRhoPimpleChem/test_case/DEM/in.liggghts_run index 408e706e..2ab52e36 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/test_case/DEM/in.liggghts_run +++ b/tutorials/cfdemSolverRhoPimpleChem/test_case/DEM/in.liggghts_run @@ -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