diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/CO2 b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/CO2 index ff10ed20..a481f4b7 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/CO2 +++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/CO2 @@ -38,7 +38,8 @@ boundaryField inlet { - type zeroGradient; + type fixedValue; + value uniform 0.0; } outlet diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/O2 b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/O2 index 6b0f56c8..64811916 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/O2 +++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/O2 @@ -38,7 +38,8 @@ boundaryField inlet { - type zeroGradient; + type fixedValue; + value uniform 1.0; } outlet diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/U b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/U index aa546d22..cae54f4e 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/U +++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/0/U @@ -17,32 +17,37 @@ FoamFile dimensions [0 1 -1 0 0 0 0]; -internalField uniform (0.0 0 0); +internalField uniform (1.0 0 0); boundaryField { // Walls in no slip condition top { - type zeroGradient; + type fixedValue; + value uniform (0.0 0 0); } bottom { - type zeroGradient; + type fixedValue; + value uniform (0.0 0 0); } side-walls { - type zeroGradient; + type fixedValue; + value uniform (0.0 0 0); } inlet { type fixedValue; - value uniform (0.0 0 0); + value uniform (1.0 0 0); } outlet { type zeroGradient; + // type fixedValue; + // value uniform (0.0 0 0); } } diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/couplingProperties b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/couplingProperties index 382a6874..98c3ee4a 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/couplingProperties +++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/couplingProperties @@ -111,6 +111,7 @@ speciesProps partMoleName "partN"; pressureFieldName "p"; partPName "partP"; + Nevery 10; } SyamlalThermCondProps diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/polyMesh/blockMeshDict b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/polyMesh/blockMeshDict index 45851e7b..9f01da08 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/polyMesh/blockMeshDict +++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/constant/polyMesh/blockMeshDict @@ -31,7 +31,7 @@ vertices blocks ( - hex (0 1 2 3 4 5 6 7)(13 4 4) simpleGrading (1 1 1) + hex (0 1 2 3 4 5 6 7)(29 9 9) simpleGrading (1 1 1) ); diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/controlDict b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/controlDict index 3b00e70b..1f7d0e4a 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/controlDict +++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/controlDict @@ -17,19 +17,19 @@ FoamFile application cfdemSolverRhoPimpleChem; -startFrom latestTime; +startFrom startTime; startTime 0; stopAt endTime; -endTime 0.5; +endTime 20.0; deltaT 0.0001; -writeControl adjustableRunTime; +writeControl timeStep; -writeInterval 0.005; +writeInterval 10000; purgeWrite 0; @@ -53,14 +53,60 @@ maxDeltaT 0.1; // ************************************************************************* // -/* libs ( - "libsimpleFunctionObjects.so" + libs ( + // "libsimpleFunctionObjects.so" "libfieldFunctionObjects.so" ); functions -{ - Average-O2 +{ + probes1 + { + type probes; + + functionObjectLibs ("libsampling.so"); + + #include "probesDict"; + } + + globalMassFrac + { + type volRegion; + libs ("libfieldFunctionObjects.so"); + writeControl timeStep;//outputTime; + writeInterval 1; + log true; + writeFields false; + regionType all; + name c0; + operation weightedVolAverage; + weightField rhoeps; + fields + ( + O2 + CO2 + N2 + ); + } + + globalMass + { + type volRegion; + libs ("libfieldFunctionObjects.so"); + writeControl timeStep;//outputTime; + writeInterval 1; + log true; + writeFields false; + regionType all; + name c1; + operation volIntegrate; + fields + ( + rhoeps + ); + } + +/* Average-O2 { type fieldAverage; functionObjectLibs ( "libfieldFunctionObjects.so" ); @@ -139,10 +185,9 @@ functions ( CO2 ); - } - - -} */ + } +*/ +} // ************************************************************************* // diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/decomposeParDict b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/decomposeParDict index 58b7b05f..c5dcd7c4 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/decomposeParDict +++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/CFD/system/decomposeParDict @@ -15,14 +15,14 @@ FoamFile } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -numberOfSubdomains 4; +numberOfSubdomains 2; //method scotch; method simple; simpleCoeffs { - n (2 2 1); + n (2 1 1); delta 0.001; } diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/DEM/in.liggghts_init b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/DEM/in.liggghts_init index 5d6a8ff3..4a51825c 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/DEM/in.liggghts_init +++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/DEM/in.liggghts_init @@ -9,7 +9,7 @@ boundary f f f newton off units si -processors 2 2 1 +processors 2 1 1 region reg block 0. 15. 0. 5. 0. 5. units box create_box 1 reg @@ -40,7 +40,7 @@ fix zwalls1 all wall/gran model hertz tangential history primitive type 1 zp fix zwalls2 all wall/gran model hertz tangential history primitive type 1 zplane 5.0 # create single particle in a specific spot -create_atoms 1 single 7.5 2 2 units box +create_atoms 1 single 7.5 2.5 2.5 units box set atom 1 diameter 0.2 density 2267 vx 0 vy 0 vz 0 # particle distributions and insertion diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/DEM/in.liggghts_run b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/DEM/in.liggghts_run index 87d97194..b1fc742e 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/DEM/in.liggghts_run +++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/DEM/in.liggghts_run @@ -10,10 +10,11 @@ boundary f f f newton off units si -processors 2 2 1 +processors 2 1 1 # read the restart file read_restart ../DEM/post/restart/liggghts.restart +reset_timestep 0 neighbor 0.0005 bin neigh_modify delay 0 binsize 0.01 @@ -52,7 +53,7 @@ fix cfd3 all couple/cfd/chemistry n_species 3 species_names O2 CO2 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 O2 molMassA 31.99 speciesC CO2 molMassC 44.01 molMassB 12.01 k 2.5e1 rmin 0.005 nevery 1000 # apply nve integration to all particles that are inserted as single particles fix integr all nve/sphere diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/parCFDDEMrun.sh b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/parCFDDEMrun.sh index cd0a4fc5..d0924910 100644 --- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/parCFDDEMrun.sh +++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/parCFDDEMrun.sh @@ -19,7 +19,7 @@ logpath=$casePath headerText="Particle_in_Duct" logfileName="log_$headerText" solverName="cfdemSolverRhoPimpleChem" -nrProcs="4" +nrProcs="2" machineFileName="none" # yourMachinefileName | none debugMode="off" # on | off| strict testHarnessPath="$CFDEM_TEST_HARNESS_PATH" diff --git a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/parDEMrun.sh b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/parDEMrun.sh index d1e3a033..f1a2533e 100755 --- a/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/parDEMrun.sh +++ b/tutorials/cfdemSolverRhoPimpleChem/Shrinking_particle_model/parDEMrun.sh @@ -20,7 +20,7 @@ logpath="$casePath" headerText="run_liggghts_init_DEM" logfileName="log_$headerText" solverName="in.liggghts_init" -nrProcs=4 +nrProcs=2 machineFileName="none" debugMode="off" #--------------------------------------------------------------------------------#