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

clean up rcfdmSolverRhoPimpleChem/3D_hot_reacting

Browse Source
This commit is contained in:
danielque
2023-07-04 16:08:02 +02:00
parent b0e4905ffd
commit fee08f3e08
113 changed files with 1163 additions and 1369 deletions
Download Patch File Download Diff File Expand all files Collapse all files

13
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/constant/RASProperties
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -15,11 +15,10 @@ FoamFile
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
RASModel laminar;//kEpsilon;
RASModel laminar;
turbulence off;
printCoeffs off;
// ************************************************************************* //

32
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/constant/couplingProperties
View File

@ -1,29 +1,19 @@
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.4 |
| \\ / A nd | Web: http://www.openfoam.org |
| \\/ M anipulation | |
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
root "";
case "";
instance "";
local "";
class dictionary;
object couplingProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
//===========================================================================//
// sub-models & settings
modelType "A";
@ -78,7 +68,8 @@ turbulenceModelType "turbulenceProperties";
getParticleTypes true;
getParticleDensities false;
getParticleEffVolFactors true;
//===========================================================================//
// sub-model properties
implicitCoupleProps
@ -102,7 +93,6 @@ viscForceProps
densityFieldName "rho";
}
BeetstraDragProps
{
velFieldName "U";
@ -148,12 +138,6 @@ SyamlalThermCondProps
{
}
centreProps
{
alphaMin 0.10;
}
engineProps
{
treeSearch true;

11
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/constant/g
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -18,5 +18,4 @@ FoamFile
dimensions [0 1 -2 0 0 0 0];
value ( 0 0 -9.81 );
// ************************************************************************* //

41
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/constant/liggghtsCommands
View File

@ -1,26 +1,17 @@
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.4 |
| \\ / A nd | Web: http://www.openfoam.org |
| \\/ M anipulation | |
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
root "";
case "";
instance "";
local "";
class dictionary;
object liggghtsCommands;
version 2.0;
format ascii;
class dictionary;
object liggghtsCommands;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
liggghtsCommandModels
@ -28,15 +19,5 @@ liggghtsCommandModels
runLiggghts
writeLiggghts
);
// ************************************************************************* //
/*runLiggghtsProps
{
preNo false;
}
writeLiggghtsProps
{
writeLast off;
overwrite off;
}*/

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/constant/symmetrizerPropertiesPost_1
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/constant/symmetrizerPropertiesPost_2
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/constant/symmetrizerPropertiesPost_3
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/constant/symmetrizerPropertiesPost_4
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/constant/symmetrizerPropertiesPre_1
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/constant/symmetrizerPropertiesPre_2
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/constant/symmetrizerPropertiesPre_3
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/constant/symmetrizerPropertiesPre_4
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/constant/thermophysicalProperties
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

14
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/constant/transportProperties
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -19,9 +19,9 @@ transportModel Newtonian;
nu nu [ 0 2 -1 0 0 0 0 ] 1e-05;
kf kf [ 1 1 -3 -1 0 0 0 ] 0.026;
kf kf [ 1 1 -3 -1 0 0 0 ] 0.026;
Cp Cp [ 0 2 -2 -1 0 0 0 ] 1005;
Cp Cp [ 0 2 -2 -1 0 0 0 ] 1005;
CrossPowerLawCoeffs
{

11
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/constant/turbulenceProperties
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -17,5 +17,4 @@ FoamFile
simulationType laminar;
// ************************************************************************* //

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/orig.0/Ksl
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

16
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/orig.0/T
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -28,19 +28,19 @@ boundaryField
bottom
{
type fixedValue;
value uniform 293;
value uniform 293;
}
walls
{
type fixedValue;
value uniform 293;
value uniform 293;
}
inlet
{
type fixedValue;
value uniform 2600;
value uniform 2600;
}
}

16
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/orig.0/U
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -39,9 +39,9 @@ boundaryField
inlet
{
// type uniformFixedValue;
// uniformValue constant (0 0 0);
// actual inflow of 112.6 would require small time step, use smaller inflow instead
//type uniformFixedValue;
//uniformValue constant (0 0 0);
// actual inflow of 112.6 would require small time step, use smaller inflow instead
type flowRateInletVelocity;
massFlowRate constant 10;
rhoInlet 4.0; // fallback value

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/orig.0/Us
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

14
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/orig.0/addSource
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -14,9 +14,9 @@ FoamFile
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -3 0 0 0 0];
dimensions [1 -1 -3 0 0 0 0];
internalField uniform 0;
internalField uniform 0;
boundaryField
{

12
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/orig.0/p
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -24,7 +24,7 @@ boundaryField
{
type zeroGradient;
}
top
{
type fixedValue;

12
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/orig.0/partEffThermCondField
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -34,7 +34,7 @@ boundaryField
{
// fixed heat loss, put all in gas phase
type fixedValue;
value uniform 0.0;
value uniform 0.0;
}
inlet

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/orig.0/partTemp
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/orig.0/particleForces
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/orig.0/voidfraction
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/orig.0/wallHeatLoss
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

346
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/system/blockMeshDict
View File

@ -1,194 +1,190 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 2.1.x |
| \ / A nd | Web: www.OpenFOAM.org |
| \/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
convertToMeters 1.0;
vertices
(
( 0.0 0.0 -1.12 ) // mmm0 = 0
( 2.7577164466275357 2.7577164466275352 -1.12 ) // c0h0 = 1
( 4.242640687119286 4.242640687119285 -1.12 ) // o0h0 = 2
( 2.1410771502600355e-16 3.4966443414554087 -1.12 ) // c1h0 = 3
( 3.6739403974420594e-16 6.0 -1.12 ) // o1h0 = 4
( -2.7577164466275352 2.7577164466275357 -1.12 ) // c2h0 = 5
( -4.242640687119285 4.242640687119286 -1.12 ) // o2h0 = 6
( -3.4966443414554087 4.282154300520071e-16 -1.12 ) // c3h0 = 7
( -6.0 7.347880794884119e-16 -1.12 ) // o3h0 = 8
( -2.757716446627536 -2.7577164466275352 -1.12 ) // c4h0 = 9
( -4.242640687119286 -4.242640687119285 -1.12 ) // o4h0 = 10
( -6.423231450780106e-16 -3.4966443414554087 -1.12 ) // c5h0 = 11
( -1.102182119232618e-15 -6.0 -1.12 ) // o5h0 = 12
( 2.757716446627535 -2.757716446627536 -1.12 ) // c6h0 = 13
( 4.242640687119284 -4.242640687119286 -1.12 ) // o6h0 = 14
( 3.4966443414554087 -8.564308601040142e-16 -1.12 ) // c7h0 = 15
( 6.0 -1.4695761589768238e-15 -1.12 ) // o7h0 = 16
( 0.0 0.0 0.746 ) // mmm1 = 17
( 2.7577164466275357 2.7577164466275352 0.746 ) // c0h1 = 18
( 4.242640687119286 4.242640687119285 0.746 ) // o0h1 = 19
( 2.1410771502600355e-16 3.4966443414554087 0.746 ) // c1h1 = 20
( 3.6739403974420594e-16 6.0 0.746 ) // o1h1 = 21
( -2.7577164466275352 2.7577164466275357 0.746 ) // c2h1 = 22
( -4.242640687119285 4.242640687119286 0.746 ) // o2h1 = 23
( -3.4966443414554087 4.282154300520071e-16 0.746 ) // c3h1 = 24
( -6.0 7.347880794884119e-16 0.746 ) // o3h1 = 25
( -2.757716446627536 -2.7577164466275352 0.746 ) // c4h1 = 26
( -4.242640687119286 -4.242640687119285 0.746 ) // o4h1 = 27
( -6.423231450780106e-16 -3.4966443414554087 0.746 ) // c5h1 = 28
( -1.102182119232618e-15 -6.0 0.746 ) // o5h1 = 29
( 2.757716446627535 -2.757716446627536 0.746 ) // c6h1 = 30
( 4.242640687119284 -4.242640687119286 0.746 ) // o6h1 = 31
( 3.4966443414554087 -8.564308601040142e-16 0.746 ) // c7h1 = 32
( 6.0 -1.4695761589768238e-15 0.746 ) // o7h1 = 33
( 0.0 0.0 3.913 ) // mmm2 = 34
( 3.3919912293518686 3.391991229351868 3.913 ) // c0h2 = 35
( 5.2184480451567214 5.218448045156721 3.913 ) // o0h2 = 36
( 2.6335248948198435e-16 4.300872539990152 3.913 ) // c1h2 = 37
( 4.518946688853733e-16 7.38 3.913 ) // o1h2 = 38
( -3.391991229351868 3.3919912293518686 3.913 ) // c2h2 = 39
( -5.218448045156721 5.2184480451567214 3.913 ) // o2h2 = 40
( -4.300872539990152 5.267049789639687e-16 3.913 ) // c3h2 = 41
( -7.38 9.037893377707467e-16 3.913 ) // o3h2 = 42
( -3.391991229351869 -3.391991229351868 3.913 ) // c4h2 = 43
( -5.2184480451567214 -5.218448045156721 3.913 ) // o4h2 = 44
( -7.900574684459529e-16 -4.300872539990152 3.913 ) // c5h2 = 45
( -1.3556840066561198e-15 -7.38 3.913 ) // o5h2 = 46
( 3.3919912293518673 -3.391991229351869 3.913 ) // c6h2 = 47
( 5.21844804515672 -5.2184480451567214 3.913 ) // o6h2 = 48
( 4.300872539990152 -1.0534099579279374e-15 3.913 ) // c7h2 = 49
( 7.38 -1.8075786755414934e-15 3.913 ) // o7h2 = 50
( 0.0 0.0 6.408 ) // mmm3 = 51
( 3.3919912293518686 3.391991229351868 6.408 ) // c0h3 = 52
( 5.2184480451567214 5.218448045156721 6.408 ) // o0h3 = 53
( 2.6335248948198435e-16 4.300872539990152 6.408 ) // c1h3 = 54
( 4.518946688853733e-16 7.38 6.408 ) // o1h3 = 55
( -3.391991229351868 3.3919912293518686 6.408 ) // c2h3 = 56
( -5.218448045156721 5.2184480451567214 6.408 ) // o2h3 = 57
( -4.300872539990152 5.267049789639687e-16 6.408 ) // c3h3 = 58
( -7.38 9.037893377707467e-16 6.408 ) // o3h3 = 59
( -3.391991229351869 -3.391991229351868 6.408 ) // c4h3 = 60
( -5.2184480451567214 -5.218448045156721 6.408 ) // o4h3 = 61
( -7.900574684459529e-16 -4.300872539990152 6.408 ) // c5h3 = 62
( -1.3556840066561198e-15 -7.38 6.408 ) // o5h3 = 63
( 3.3919912293518673 -3.391991229351869 6.408 ) // c6h3 = 64
( 5.21844804515672 -5.2184480451567214 6.408 ) // o6h3 = 65
( 4.300872539990152 -1.0534099579279374e-15 6.408 ) // c7h3 = 66
( 7.38 -1.8075786755414934e-15 6.408 ) // o7h3 = 67
( 0.0 0.0 24.035 ) // mmm4 = 68
( 2.1579131194860466 2.157913119486046 24.035 ) // c0h4 = 69
( 3.3198663376708413 3.3198663376708404 24.035 ) // o0h4 = 70
( 1.6753928700784779e-16 2.7361241971888575 24.035 ) // c1h4 = 71
( 2.874858360998412e-16 4.695 24.035 ) // o1h4 = 72
( -2.157913119486046 2.1579131194860466 24.035 ) // c2h4 = 73
( -3.3198663376708404 3.3198663376708413 24.035 ) // o2h4 = 74
( -2.7361241971888575 3.3507857401569558e-16 24.035 ) // c3h4 = 75
( -4.695 5.749716721996824e-16 24.035 ) // o3h4 = 76
( -2.157913119486047 -2.157913119486046 24.035 ) // c4h4 = 77
( -3.3198663376708417 -3.3198663376708404 24.035 ) // o4h4 = 78
( -5.026178610235433e-16 -2.7361241971888575 24.035 ) // c5h4 = 79
( -8.624575082995235e-16 -4.695 24.035 ) // o5h4 = 80
( 2.157913119486046 -2.157913119486047 24.035 ) // c6h4 = 81
( 3.31986633767084 -3.3198663376708417 24.035 ) // o6h4 = 82
( 2.7361241971888575 -6.7015714803139115e-16 24.035 ) // c7h4 = 83
( 4.695 -1.1499433443993648e-15 24.035 ) // o7h4 = 84
( 0.0 0.0 26.035 ) // mmm5 = 85
( 2.1579131194860466 2.157913119486046 26.035 ) // c0h5 = 86
( 3.3198663376708413 3.3198663376708404 26.035 ) // o0h5 = 87
( 1.6753928700784779e-16 2.7361241971888575 26.035 ) // c1h5 = 88
( 2.874858360998412e-16 4.695 26.035 ) // o1h5 = 89
( -2.157913119486046 2.1579131194860466 26.035 ) // c2h5 = 90
( -3.3198663376708404 3.3198663376708413 26.035 ) // o2h5 = 91
( -2.7361241971888575 3.3507857401569558e-16 26.035 ) // c3h5 = 92
( -4.695 5.749716721996824e-16 26.035 ) // o3h5 = 93
( -2.157913119486047 -2.157913119486046 26.035 ) // c4h5 = 94
( -3.3198663376708417 -3.3198663376708404 26.035 ) // o4h5 = 95
( -5.026178610235433e-16 -2.7361241971888575 26.035 ) // c5h5 = 96
( -8.624575082995235e-16 -4.695 26.035 ) // o5h5 = 97
( 2.157913119486046 -2.157913119486047 26.035 ) // c6h5 = 98
( 3.31986633767084 -3.3198663376708417 26.035 ) // o6h5 = 99
( 2.7361241971888575 -6.7015714803139115e-16 26.035 ) // c7h5 = 100
( 4.695 -1.1499433443993648e-15 26.035 ) // o7h5 = 101
( 0.0 0.0 -1.12 ) // mmm0 = 0
( 2.7577164466275357 2.7577164466275352 -1.12 ) // c0h0 = 1
( 4.242640687119286 4.242640687119285 -1.12 ) // o0h0 = 2
( 2.1410771502600355e-16 3.4966443414554087 -1.12 ) // c1h0 = 3
( 3.6739403974420594e-16 6.0 -1.12 ) // o1h0 = 4
( -2.7577164466275352 2.7577164466275357 -1.12 ) // c2h0 = 5
( -4.242640687119285 4.242640687119286 -1.12 ) // o2h0 = 6
( -3.4966443414554087 4.282154300520071e-16 -1.12 ) // c3h0 = 7
( -6.0 7.347880794884119e-16 -1.12 ) // o3h0 = 8
( -2.757716446627536 -2.7577164466275352 -1.12 ) // c4h0 = 9
( -4.242640687119286 -4.242640687119285 -1.12 ) // o4h0 = 10
( -6.423231450780106e-16 -3.4966443414554087 -1.12 ) // c5h0 = 11
( -1.102182119232618e-15 -6.0 -1.12 ) // o5h0 = 12
( 2.757716446627535 -2.757716446627536 -1.12 ) // c6h0 = 13
( 4.242640687119284 -4.242640687119286 -1.12 ) // o6h0 = 14
( 3.4966443414554087 -8.564308601040142e-16 -1.12 ) // c7h0 = 15
( 6.0 -1.4695761589768238e-15 -1.12 ) // o7h0 = 16
( 0.0 0.0 0.746 ) // mmm1 = 17
( 2.7577164466275357 2.7577164466275352 0.746 ) // c0h1 = 18
( 4.242640687119286 4.242640687119285 0.746 ) // o0h1 = 19
( 2.1410771502600355e-16 3.4966443414554087 0.746 ) // c1h1 = 20
( 3.6739403974420594e-16 6.0 0.746 ) // o1h1 = 21
( -2.7577164466275352 2.7577164466275357 0.746 ) // c2h1 = 22
( -4.242640687119285 4.242640687119286 0.746 ) // o2h1 = 23
( -3.4966443414554087 4.282154300520071e-16 0.746 ) // c3h1 = 24
( -6.0 7.347880794884119e-16 0.746 ) // o3h1 = 25
( -2.757716446627536 -2.7577164466275352 0.746 ) // c4h1 = 26
( -4.242640687119286 -4.242640687119285 0.746 ) // o4h1 = 27
( -6.423231450780106e-16 -3.4966443414554087 0.746 ) // c5h1 = 28
( -1.102182119232618e-15 -6.0 0.746 ) // o5h1 = 29
( 2.757716446627535 -2.757716446627536 0.746 ) // c6h1 = 30
( 4.242640687119284 -4.242640687119286 0.746 ) // o6h1 = 31
( 3.4966443414554087 -8.564308601040142e-16 0.746 ) // c7h1 = 32
( 6.0 -1.4695761589768238e-15 0.746 ) // o7h1 = 33
( 0.0 0.0 3.913 ) // mmm2 = 34
( 3.3919912293518686 3.391991229351868 3.913 ) // c0h2 = 35
( 5.2184480451567214 5.218448045156721 3.913 ) // o0h2 = 36
( 2.6335248948198435e-16 4.300872539990152 3.913 ) // c1h2 = 37
( 4.518946688853733e-16 7.38 3.913 ) // o1h2 = 38
( -3.391991229351868 3.3919912293518686 3.913 ) // c2h2 = 39
( -5.218448045156721 5.2184480451567214 3.913 ) // o2h2 = 40
( -4.300872539990152 5.267049789639687e-16 3.913 ) // c3h2 = 41
( -7.38 9.037893377707467e-16 3.913 ) // o3h2 = 42
( -3.391991229351869 -3.391991229351868 3.913 ) // c4h2 = 43
( -5.2184480451567214 -5.218448045156721 3.913 ) // o4h2 = 44
( -7.900574684459529e-16 -4.300872539990152 3.913 ) // c5h2 = 45
( -1.3556840066561198e-15 -7.38 3.913 ) // o5h2 = 46
( 3.3919912293518673 -3.391991229351869 3.913 ) // c6h2 = 47
( 5.21844804515672 -5.2184480451567214 3.913 ) // o6h2 = 48
( 4.300872539990152 -1.0534099579279374e-15 3.913 ) // c7h2 = 49
( 7.38 -1.8075786755414934e-15 3.913 ) // o7h2 = 50
( 0.0 0.0 6.408 ) // mmm3 = 51
( 3.3919912293518686 3.391991229351868 6.408 ) // c0h3 = 52
( 5.2184480451567214 5.218448045156721 6.408 ) // o0h3 = 53
( 2.6335248948198435e-16 4.300872539990152 6.408 ) // c1h3 = 54
( 4.518946688853733e-16 7.38 6.408 ) // o1h3 = 55
( -3.391991229351868 3.3919912293518686 6.408 ) // c2h3 = 56
( -5.218448045156721 5.2184480451567214 6.408 ) // o2h3 = 57
( -4.300872539990152 5.267049789639687e-16 6.408 ) // c3h3 = 58
( -7.38 9.037893377707467e-16 6.408 ) // o3h3 = 59
( -3.391991229351869 -3.391991229351868 6.408 ) // c4h3 = 60
( -5.2184480451567214 -5.218448045156721 6.408 ) // o4h3 = 61
( -7.900574684459529e-16 -4.300872539990152 6.408 ) // c5h3 = 62
( -1.3556840066561198e-15 -7.38 6.408 ) // o5h3 = 63
( 3.3919912293518673 -3.391991229351869 6.408 ) // c6h3 = 64
( 5.21844804515672 -5.2184480451567214 6.408 ) // o6h3 = 65
( 4.300872539990152 -1.0534099579279374e-15 6.408 ) // c7h3 = 66
( 7.38 -1.8075786755414934e-15 6.408 ) // o7h3 = 67
( 0.0 0.0 24.035 ) // mmm4 = 68
( 2.1579131194860466 2.157913119486046 24.035 ) // c0h4 = 69
( 3.3198663376708413 3.3198663376708404 24.035 ) // o0h4 = 70
( 1.6753928700784779e-16 2.7361241971888575 24.035 ) // c1h4 = 71
( 2.874858360998412e-16 4.695 24.035 ) // o1h4 = 72
( -2.157913119486046 2.1579131194860466 24.035 ) // c2h4 = 73
( -3.3198663376708404 3.3198663376708413 24.035 ) // o2h4 = 74
( -2.7361241971888575 3.3507857401569558e-16 24.035 ) // c3h4 = 75
( -4.695 5.749716721996824e-16 24.035 ) // o3h4 = 76
( -2.157913119486047 -2.157913119486046 24.035 ) // c4h4 = 77
( -3.3198663376708417 -3.3198663376708404 24.035 ) // o4h4 = 78
( -5.026178610235433e-16 -2.7361241971888575 24.035 ) // c5h4 = 79
( -8.624575082995235e-16 -4.695 24.035 ) // o5h4 = 80
( 2.157913119486046 -2.157913119486047 24.035 ) // c6h4 = 81
( 3.31986633767084 -3.3198663376708417 24.035 ) // o6h4 = 82
( 2.7361241971888575 -6.7015714803139115e-16 24.035 ) // c7h4 = 83
( 4.695 -1.1499433443993648e-15 24.035 ) // o7h4 = 84
( 0.0 0.0 26.035 ) // mmm5 = 85
( 2.1579131194860466 2.157913119486046 26.035 ) // c0h5 = 86
( 3.3198663376708413 3.3198663376708404 26.035 ) // o0h5 = 87
( 1.6753928700784779e-16 2.7361241971888575 26.035 ) // c1h5 = 88
( 2.874858360998412e-16 4.695 26.035 ) // o1h5 = 89
( -2.157913119486046 2.1579131194860466 26.035 ) // c2h5 = 90
( -3.3198663376708404 3.3198663376708413 26.035 ) // o2h5 = 91
( -2.7361241971888575 3.3507857401569558e-16 26.035 ) // c3h5 = 92
( -4.695 5.749716721996824e-16 26.035 ) // o3h5 = 93
( -2.157913119486047 -2.157913119486046 26.035 ) // c4h5 = 94
( -3.3198663376708417 -3.3198663376708404 26.035 ) // o4h5 = 95
( -5.026178610235433e-16 -2.7361241971888575 26.035 ) // c5h5 = 96
( -8.624575082995235e-16 -4.695 26.035 ) // o5h5 = 97
( 2.157913119486046 -2.157913119486047 26.035 ) // c6h5 = 98
( 3.31986633767084 -3.3198663376708417 26.035 ) // o6h5 = 99
( 2.7361241971888575 -6.7015714803139115e-16 26.035 ) // c7h5 = 100
( 4.695 -1.1499433443993648e-15 26.035 ) // o7h5 = 101
);
blocks
(
hex ( 0 15 1 3 17 32 18 20 ) ( 12 12 5 ) simpleGrading ( 1 1 1 ) // 0
hex ( 1 2 4 3 18 19 21 20 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 1
hex ( 3 4 6 5 20 21 23 22 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 2
hex ( 0 3 5 7 17 20 22 24 ) ( 12 12 5 ) simpleGrading ( 1 1 1 ) // 3
hex ( 5 6 8 7 22 23 25 24 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 4
hex ( 7 8 10 9 24 25 27 26 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 5
hex ( 0 7 9 11 17 24 26 28 ) ( 12 12 5 ) simpleGrading ( 1 1 1 ) // 6
hex ( 9 10 12 11 26 27 29 28 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 7
hex ( 11 12 14 13 28 29 31 30 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 8
hex ( 0 11 13 15 17 28 30 32 ) ( 12 12 5 ) simpleGrading ( 1 1 1 ) // 9
hex ( 13 14 16 15 30 31 33 32 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 10
hex ( 15 16 2 1 32 33 19 18 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 11
hex ( 17 32 18 20 34 49 35 37 ) ( 12 12 9 ) simpleGrading ( 1 1 1 ) // 12
hex ( 18 19 21 20 35 36 38 37 ) ( 11 12 9 ) simpleGrading ( 1 1 1 ) // 13
hex ( 20 21 23 22 37 38 40 39 ) ( 11 12 9 ) simpleGrading ( 1 1 1 ) // 14
hex ( 17 20 22 24 34 37 39 41 ) ( 12 12 9 ) simpleGrading ( 1 1 1 ) // 15
hex ( 22 23 25 24 39 40 42 41 ) ( 11 12 9 ) simpleGrading ( 1 1 1 ) // 16
hex ( 24 25 27 26 41 42 44 43 ) ( 11 12 9 ) simpleGrading ( 1 1 1 ) // 17
hex ( 17 24 26 28 34 41 43 45 ) ( 12 12 9 ) simpleGrading ( 1 1 1 ) // 18
hex ( 26 27 29 28 43 44 46 45 ) ( 11 12 9 ) simpleGrading ( 1 1 1 ) // 19
hex ( 28 29 31 30 45 46 48 47 ) ( 11 12 9 ) simpleGrading ( 1 1 1 ) // 20
hex ( 17 28 30 32 34 45 47 49 ) ( 12 12 9 ) simpleGrading ( 1 1 1 ) // 21
hex ( 30 31 33 32 47 48 50 49 ) ( 11 12 9 ) simpleGrading ( 1 1 1 ) // 22
hex ( 32 33 19 18 49 50 36 35 ) ( 11 12 9 ) simpleGrading ( 1 1 1 ) // 23
hex ( 34 49 35 37 51 66 52 54 ) ( 12 12 7 ) simpleGrading ( 1 1 1 ) // 24
hex ( 35 36 38 37 52 53 55 54 ) ( 11 12 7 ) simpleGrading ( 1 1 1 ) // 25
hex ( 37 38 40 39 54 55 57 56 ) ( 11 12 7 ) simpleGrading ( 1 1 1 ) // 26
hex ( 34 37 39 41 51 54 56 58 ) ( 12 12 7 ) simpleGrading ( 1 1 1 ) // 27
hex ( 39 40 42 41 56 57 59 58 ) ( 11 12 7 ) simpleGrading ( 1 1 1 ) // 28
hex ( 41 42 44 43 58 59 61 60 ) ( 11 12 7 ) simpleGrading ( 1 1 1 ) // 29
hex ( 34 41 43 45 51 58 60 62 ) ( 12 12 7 ) simpleGrading ( 1 1 1 ) // 30
hex ( 43 44 46 45 60 61 63 62 ) ( 11 12 7 ) simpleGrading ( 1 1 1 ) // 31
hex ( 45 46 48 47 62 63 65 64 ) ( 11 12 7 ) simpleGrading ( 1 1 1 ) // 32
hex ( 34 45 47 49 51 62 64 66 ) ( 12 12 7 ) simpleGrading ( 1 1 1 ) // 33
hex ( 47 48 50 49 64 65 67 66 ) ( 11 12 7 ) simpleGrading ( 1 1 1 ) // 34
hex ( 49 50 36 35 66 67 53 52 ) ( 11 12 7 ) simpleGrading ( 1 1 1 ) // 35
hex ( 51 66 52 54 68 83 69 71 ) ( 12 12 50 ) simpleGrading ( 1 1 1 ) // 36
hex ( 52 53 55 54 69 70 72 71 ) ( 11 12 50 ) simpleGrading ( 1 1 1 ) // 37
hex ( 54 55 57 56 71 72 74 73 ) ( 11 12 50 ) simpleGrading ( 1 1 1 ) // 38
hex ( 51 54 56 58 68 71 73 75 ) ( 12 12 50 ) simpleGrading ( 1 1 1 ) // 39
hex ( 56 57 59 58 73 74 76 75 ) ( 11 12 50 ) simpleGrading ( 1 1 1 ) // 40
hex ( 58 59 61 60 75 76 78 77 ) ( 11 12 50 ) simpleGrading ( 1 1 1 ) // 41
hex ( 51 58 60 62 68 75 77 79 ) ( 12 12 50 ) simpleGrading ( 1 1 1 ) // 42
hex ( 60 61 63 62 77 78 80 79 ) ( 11 12 50 ) simpleGrading ( 1 1 1 ) // 43
hex ( 62 63 65 64 79 80 82 81 ) ( 11 12 50 ) simpleGrading ( 1 1 1 ) // 44
hex ( 51 62 64 66 68 79 81 83 ) ( 12 12 50 ) simpleGrading ( 1 1 1 ) // 45
hex ( 64 65 67 66 81 82 84 83 ) ( 11 12 50 ) simpleGrading ( 1 1 1 ) // 46
hex ( 66 67 53 52 83 84 70 69 ) ( 11 12 50 ) simpleGrading ( 1 1 1 ) // 47
hex ( 68 83 69 71 85 100 86 88 ) ( 12 12 5 ) simpleGrading ( 1 1 1 ) // 48
hex ( 69 70 72 71 86 87 89 88 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 49
hex ( 71 72 74 73 88 89 91 90 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 50
hex ( 68 71 73 75 85 88 90 92 ) ( 12 12 5 ) simpleGrading ( 1 1 1 ) // 51
hex ( 73 74 76 75 90 91 93 92 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 52
hex ( 75 76 78 77 92 93 95 94 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 53
hex ( 68 75 77 79 85 92 94 96 ) ( 12 12 5 ) simpleGrading ( 1 1 1 ) // 54
hex ( 77 78 80 79 94 95 97 96 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 55
hex ( 79 80 82 81 96 97 99 98 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 56
hex ( 68 79 81 83 85 96 98 100 ) ( 12 12 5 ) simpleGrading ( 1 1 1 ) // 57
hex ( 81 82 84 83 98 99 101 100 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 58
hex ( 83 84 70 69 100 101 87 86 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 59
hex ( 0 15 1 3 17 32 18 20 ) ( 12 12 5 ) simpleGrading ( 1 1 1 ) // 0
hex ( 1 2 4 3 18 19 21 20 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 1
hex ( 3 4 6 5 20 21 23 22 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 2
hex ( 0 3 5 7 17 20 22 24 ) ( 12 12 5 ) simpleGrading ( 1 1 1 ) // 3
hex ( 5 6 8 7 22 23 25 24 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 4
hex ( 7 8 10 9 24 25 27 26 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 5
hex ( 0 7 9 11 17 24 26 28 ) ( 12 12 5 ) simpleGrading ( 1 1 1 ) // 6
hex ( 9 10 12 11 26 27 29 28 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 7
hex ( 11 12 14 13 28 29 31 30 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 8
hex ( 0 11 13 15 17 28 30 32 ) ( 12 12 5 ) simpleGrading ( 1 1 1 ) // 9
hex ( 13 14 16 15 30 31 33 32 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 10
hex ( 15 16 2 1 32 33 19 18 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 11
hex ( 17 32 18 20 34 49 35 37 ) ( 12 12 9 ) simpleGrading ( 1 1 1 ) // 12
hex ( 18 19 21 20 35 36 38 37 ) ( 11 12 9 ) simpleGrading ( 1 1 1 ) // 13
hex ( 20 21 23 22 37 38 40 39 ) ( 11 12 9 ) simpleGrading ( 1 1 1 ) // 14
hex ( 17 20 22 24 34 37 39 41 ) ( 12 12 9 ) simpleGrading ( 1 1 1 ) // 15
hex ( 22 23 25 24 39 40 42 41 ) ( 11 12 9 ) simpleGrading ( 1 1 1 ) // 16
hex ( 24 25 27 26 41 42 44 43 ) ( 11 12 9 ) simpleGrading ( 1 1 1 ) // 17
hex ( 17 24 26 28 34 41 43 45 ) ( 12 12 9 ) simpleGrading ( 1 1 1 ) // 18
hex ( 26 27 29 28 43 44 46 45 ) ( 11 12 9 ) simpleGrading ( 1 1 1 ) // 19
hex ( 28 29 31 30 45 46 48 47 ) ( 11 12 9 ) simpleGrading ( 1 1 1 ) // 20
hex ( 17 28 30 32 34 45 47 49 ) ( 12 12 9 ) simpleGrading ( 1 1 1 ) // 21
hex ( 30 31 33 32 47 48 50 49 ) ( 11 12 9 ) simpleGrading ( 1 1 1 ) // 22
hex ( 32 33 19 18 49 50 36 35 ) ( 11 12 9 ) simpleGrading ( 1 1 1 ) // 23
hex ( 34 49 35 37 51 66 52 54 ) ( 12 12 7 ) simpleGrading ( 1 1 1 ) // 24
hex ( 35 36 38 37 52 53 55 54 ) ( 11 12 7 ) simpleGrading ( 1 1 1 ) // 25
hex ( 37 38 40 39 54 55 57 56 ) ( 11 12 7 ) simpleGrading ( 1 1 1 ) // 26
hex ( 34 37 39 41 51 54 56 58 ) ( 12 12 7 ) simpleGrading ( 1 1 1 ) // 27
hex ( 39 40 42 41 56 57 59 58 ) ( 11 12 7 ) simpleGrading ( 1 1 1 ) // 28
hex ( 41 42 44 43 58 59 61 60 ) ( 11 12 7 ) simpleGrading ( 1 1 1 ) // 29
hex ( 34 41 43 45 51 58 60 62 ) ( 12 12 7 ) simpleGrading ( 1 1 1 ) // 30
hex ( 43 44 46 45 60 61 63 62 ) ( 11 12 7 ) simpleGrading ( 1 1 1 ) // 31
hex ( 45 46 48 47 62 63 65 64 ) ( 11 12 7 ) simpleGrading ( 1 1 1 ) // 32
hex ( 34 45 47 49 51 62 64 66 ) ( 12 12 7 ) simpleGrading ( 1 1 1 ) // 33
hex ( 47 48 50 49 64 65 67 66 ) ( 11 12 7 ) simpleGrading ( 1 1 1 ) // 34
hex ( 49 50 36 35 66 67 53 52 ) ( 11 12 7 ) simpleGrading ( 1 1 1 ) // 35
hex ( 51 66 52 54 68 83 69 71 ) ( 12 12 50 ) simpleGrading ( 1 1 1 ) // 36
hex ( 52 53 55 54 69 70 72 71 ) ( 11 12 50 ) simpleGrading ( 1 1 1 ) // 37
hex ( 54 55 57 56 71 72 74 73 ) ( 11 12 50 ) simpleGrading ( 1 1 1 ) // 38
hex ( 51 54 56 58 68 71 73 75 ) ( 12 12 50 ) simpleGrading ( 1 1 1 ) // 39
hex ( 56 57 59 58 73 74 76 75 ) ( 11 12 50 ) simpleGrading ( 1 1 1 ) // 40
hex ( 58 59 61 60 75 76 78 77 ) ( 11 12 50 ) simpleGrading ( 1 1 1 ) // 41
hex ( 51 58 60 62 68 75 77 79 ) ( 12 12 50 ) simpleGrading ( 1 1 1 ) // 42
hex ( 60 61 63 62 77 78 80 79 ) ( 11 12 50 ) simpleGrading ( 1 1 1 ) // 43
hex ( 62 63 65 64 79 80 82 81 ) ( 11 12 50 ) simpleGrading ( 1 1 1 ) // 44
hex ( 51 62 64 66 68 79 81 83 ) ( 12 12 50 ) simpleGrading ( 1 1 1 ) // 45
hex ( 64 65 67 66 81 82 84 83 ) ( 11 12 50 ) simpleGrading ( 1 1 1 ) // 46
hex ( 66 67 53 52 83 84 70 69 ) ( 11 12 50 ) simpleGrading ( 1 1 1 ) // 47
hex ( 68 83 69 71 85 100 86 88 ) ( 12 12 5 ) simpleGrading ( 1 1 1 ) // 48
hex ( 69 70 72 71 86 87 89 88 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 49
hex ( 71 72 74 73 88 89 91 90 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 50
hex ( 68 71 73 75 85 88 90 92 ) ( 12 12 5 ) simpleGrading ( 1 1 1 ) // 51
hex ( 73 74 76 75 90 91 93 92 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 52
hex ( 75 76 78 77 92 93 95 94 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 53
hex ( 68 75 77 79 85 92 94 96 ) ( 12 12 5 ) simpleGrading ( 1 1 1 ) // 54
hex ( 77 78 80 79 94 95 97 96 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 55
hex ( 79 80 82 81 96 97 99 98 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 56
hex ( 68 79 81 83 85 96 98 100 ) ( 12 12 5 ) simpleGrading ( 1 1 1 ) // 57
hex ( 81 82 84 83 98 99 101 100 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 58
hex ( 83 84 70 69 100 101 87 86 ) ( 11 12 5 ) simpleGrading ( 1 1 1 ) // 59
);

31
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/system/controlDict
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -54,7 +54,7 @@ functions
type fieldAverage;
libs ("libfieldFunctionObjects.so");
writeControl writeTime;
timeStart 18.0;
timeStart 18.0;
fields
(
@ -65,14 +65,14 @@ functions
base time;
window 5.0001;
}
Us
Us
{
mean on;
prime2Mean off;
base time;
window 5.0001;
}
voidfraction
voidfraction
{
mean on;
prime2Mean off;
@ -95,7 +95,7 @@ functions
window 2.0001;
windowName w1;
}
voidfraction
voidfraction
{
mean on;
prime2Mean off;
@ -110,9 +110,9 @@ functions
{
type writeObjects;
functionObjectLibs ( "libutilityFunctionObjects.so" );
exclusiveWriting true;
exclusiveWriting true;
objects ("Us" "UsMean" "voidfraction" "voidfractionMean" "partType" "dSauter" "dSauterMean");
writeControl timeStep;
writeControl timeStep;
writeInterval 10000;
}
@ -120,9 +120,9 @@ functions
{
type probes;
functionObjectLibs ("libsampling.so");
functionObjectLibs ("libsampling.so");
#include "probesDict";
#include "probesDict";
}
cellObj1
@ -148,10 +148,9 @@ functions
{
type abort;
libs ("libutilityFunctionObjects.so");
action writeNow;
file "ABORT";
action writeNow;
file "ABORT";
}
}
// ************************************************************************* //

1
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/system/createPatchDict
View File

@ -12,7 +12,6 @@ FoamFile
class dictionary;
object createPatchDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
pointSync false;

16
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/system/decomposeParDict
View File

@ -1,10 +1,10 @@
/*-------------------------------*- C++ -*---------------------------------*\
| ========= |
| \\ / OpenFOAM |
| \\ / |
| \\ / The Open Source CFD Toolbox |
| \\/ http://www.OpenFOAM.org |
\*-------------------------------------------------------------------------*/
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
@ -14,7 +14,6 @@ FoamFile
location "system";
object decomposeParDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
numberOfSubdomains 32;
@ -28,4 +27,3 @@ simpleCoeffs
}
// ************************************************************************* //

12
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/system/fvOptions
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -17,7 +17,7 @@ FoamFile
source1
{
type limitTemperature;
type limitTemperature;
active yes;
limitTemperatureCoeffs
{

18
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/system/fvSchemes
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -24,15 +24,15 @@ gradSchemes
{
default Gauss linear;
grad(p) cellLimited Gauss linear 1;
// grad(U) Gauss linear;
//grad(U) Gauss linear;
grad(U) cellLimited Gauss linear 1;
}
divSchemes
{
default Gauss linear;
div(phi,U) Gauss upwind;
// div(phi,U) bounded Gauss linearUpwindV grad(U);
div(phi,U) Gauss upwind;
//div(phi,U) bounded Gauss linearUpwindV grad(U);
div(phi,k) Gauss limitedLinear 1;
div(phi,epsilon) Gauss limitedLinear 1;
div(phi,R) Gauss limitedLinear 1;
@ -68,7 +68,7 @@ interpolationSchemes
snGradSchemes
{
// default corrected;
//default corrected;
default limited 0.333;
}

14
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/system/fvSolution
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -57,7 +57,7 @@ solvers
{
//solver smoothSolver;
//smoother GaussSeidel;
solver PBiCG;
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0;
@ -67,7 +67,7 @@ solvers
{
//solver smoothSolver;
//smoother GaussSeidel;
solver PBiCG;
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0;

61
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/system/probesDict
View File

@ -1,21 +1,19 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 2.1.x |
| \ / A nd | Web: www.OpenFOAM.org |
| \/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object probesDict;
version 2.0;
format ascii;
class dictionary;
object probesDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
fields
(
voidfraction
@ -36,27 +34,26 @@ writeInterval 40;
// Locations to be probed.
probeLocations
(
(0.001 0.001 -1.0)
(0.001 0.001 0.0)
(0.001 0.001 1.0)
(0.001 0.001 2.0)
(0.001 0.001 3.0)
(0.001 0.001 4.0)
(0.001 0.001 5.0)
(0.001 0.001 6.0)
(0.001 0.001 7.0)
(0.001 0.001 8.0)
(0.001 0.001 9.0)
(0.001 0.001 10.0)
(0.001 0.001 12.0)
(0.001 0.001 14.0)
(0.001 0.001 16.0)
(0.001 0.001 18.0)
(0.001 0.001 20.0)
(0.001 0.001 22.0)
(0.001 0.001 24.0)
(0.001 0.001 26.0)
(0.001 0.001 -1.0)
(0.001 0.001 0.0)
(0.001 0.001 1.0)
(0.001 0.001 2.0)
(0.001 0.001 3.0)
(0.001 0.001 4.0)
(0.001 0.001 5.0)
(0.001 0.001 6.0)
(0.001 0.001 7.0)
(0.001 0.001 8.0)
(0.001 0.001 9.0)
(0.001 0.001 10.0)
(0.001 0.001 12.0)
(0.001 0.001 14.0)
(0.001 0.001 16.0)
(0.001 0.001 18.0)
(0.001 0.001 20.0)
(0.001 0.001 22.0)
(0.001 0.001 24.0)
(0.001 0.001 26.0)
);
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

12
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/system/topoSetDict_OtherRegs
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -12,8 +12,8 @@ FoamFile
class dictionary;
object topoSetDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
actions
(
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/system/topoSetDict_RWFaces
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2012 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/CFD/system/topoSetDict_RWReg
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2012 |
| \\ / A nd | Website: www.openfoam.com |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

70
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/DEM/DEMparameters
View File

@ -24,26 +24,26 @@ variable roll_C equal 0.5 # coefficient of rolling friction
# mean density coke = 950
# bulk density = 500
variable dens_C equal 950 # density
variable rad_C equal 0.088 # c.g. radius (cg = 4)
variable poly_C equal 0.81 # polydisp factor
variable dens_C equal 950 # density
variable rad_C equal 0.088 # c.g. radius (cg = 4)
variable poly_C equal 0.81 # polydisp factor
# DEM-parameters FINE COKE
variable yM_fC equal 5.e7 # Young's modulus
variable pR_fC equal 0.25 # Poisson ratio
variable rest_fC equal 0.5 # coefficient of restitution
variable fric_fC equal 0.9 # coefficient of friction
variable roll_fC equal 0.5 # coefficient of rolling friction
variable yM_fC equal 5.e7 # Young's modulus
variable pR_fC equal 0.25 # Poisson ratio
variable rest_fC equal 0.5 # coefficient of restitution
variable fric_fC equal 0.9 # coefficient of friction
variable roll_fC equal 0.5 # coefficient of rolling friction
# mean diameter fine coke = 0.02 m
# mean density fine coke = 950
# bulk density fine coke = 570
variable dens_fC equal 950 # density
variable rad_fC equal 0.08 # c.g. radius (cg = 8)
#variable poly_fC equal 0.92 # polydisp factor
variable poly_fC equal 1.2 # polydisp factor accounting for near-wall inaccuracy
variable dens_fC equal 950 # density
variable rad_fC equal 0.08 # c.g. radius (cg = 8)
#variable poly_fC equal 0.92 # polydisp factor
variable poly_fC equal 1.2 # polydisp factor accounting for near-wall inaccuracy
@ -52,23 +52,23 @@ variable yM_O equal 5.e7 # Young's modulus
variable pR_O equal 0.25 # Poisson ratio
variable rest_O equal 0.5 # coefficient of restitution
variable fric_O equal 0.8 # coefficient of friction
variable roll_O equal 0.4 # coefficient of rolling friction
variable roll_O equal 0.4 # coefficient of rolling friction
# values according to Modern Blast Furnace Ironmaking
# mean diameter ore mixture = 0.015 m
# mean density ore mixture = 3300
# bulk density = 2200
variable dens_O equal 3300 # density
variable rad_O equal 0.0825 # c.g. radius (cg = 11)
variable poly_O equal 1.03 # polydisp factor
variable dens_O equal 3300 # density
variable rad_O equal 0.0825 # c.g. radius (cg = 11)
variable poly_O equal 1.03 # polydisp factor
# DEM-parameters FINE ORE (represents the mixture of sinter, pellets and actual ore)
variable yM_fO equal 5.e7 # Young's modulus
variable pR_fO equal 0.25 # Poisson ratio
variable rest_fO equal 0.5 # coefficient of restitution
variable fric_fO equal 0.8 # coefficient of friction
variable yM_fO equal 5.e7 # Young's modulus
variable pR_fO equal 0.25 # Poisson ratio
variable rest_fO equal 0.5 # coefficient of restitution
variable fric_fO equal 0.8 # coefficient of friction
variable roll_fO equal 0.4 # coefficient of rolling friction
# values according to G. Holzinger's literature study
@ -76,10 +76,10 @@ variable roll_fO equal 0.4 # coefficient of rolling friction
# mean density fine ore mixture = 3300
# bulk density = 2200
variable dens_fO equal 3300 # density
variable rad_fO equal 0.072 # c.g. radius (cg = 24)
#variable poly_fO equal 1.03 # polydisp factor
variable poly_fO equal 1.2 # polydisp factor accounting for near-wall inaccuracy
variable dens_fO equal 3300 # density
variable rad_fO equal 0.072 # c.g. radius (cg = 24)
#variable poly_fO equal 1.03 # polydisp factor
variable poly_fO equal 1.2 # polydisp factor accounting for near-wall inaccuracy
@ -88,19 +88,19 @@ variable poly_fO equal 1.2 # polydisp factor accounting for near-wall inac
variable yM_W equal 5.e6 # Young's modulus
variable pR_W equal 0.25 # Poisson ratio
variable rest_W equal 0.5 # coefficient of restitution
variable fric_W equal 0.85 # coefficient of friction
variable fric_W equal 0.85 # coefficient of friction
variable roll_W equal 0.45 # coefficient of rolling friction
# DEM-parameters mixed
variable rest_CO equal 0.2 # coefficient of restitution
variable rest_WO equal 0.2 # coefficient of restitution
variable rest_CW equal 0.5 # coefficient of restitution
variable fric_CO equal 0.8 # coefficient of friction
variable fric_WO equal 0.2 # coefficient of friction
variable fric_CW equal 0.5 # coefficient of friction
variable roll_CO equal 0.02 # for particles and walls with friction
variable roll_WO equal 0.02 # for particles and walls with friction
variable roll_CW equal 0.02 # for particles and walls with friction
variable rest_CO equal 0.2 # coefficient of restitution
variable rest_WO equal 0.2 # coefficient of restitution
variable rest_CW equal 0.5 # coefficient of restitution
variable fric_CO equal 0.8 # coefficient of friction
variable fric_WO equal 0.2 # coefficient of friction
variable fric_CW equal 0.5 # coefficient of friction
variable roll_CO equal 0.02 # for particles and walls with friction
variable roll_WO equal 0.02 # for particles and walls with friction
variable roll_CW equal 0.02 # for particles and walls with friction
variable cV equal 0.1 # characteristic velocity
variable cV equal 0.1 # characteristic velocity

20
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/DEM/in.coupling_integration
View File

@ -7,18 +7,18 @@ variable Tin equal 293
group fullyDeformed empty
fix cfd1 all couple/cfd couple_every ${Ncouple} mpi
fix cfd2 all couple/cfd/force/implicit transfer_density no transfer_type yes transfer_property name effvolfactor type scalar-atom
fix cfd3 all couple/cfd/convection T0 293
fix cfd4 all couple/cfd/parttempfield T0 293
fix cfd5 all couple/cfd/deform rmin 0.01
fix cfd1 all couple/cfd couple_every ${Ncouple} mpi
fix cfd2 all couple/cfd/force/implicit transfer_density no transfer_type yes transfer_property name effvolfactor type scalar-atom
fix cfd3 all couple/cfd/convection T0 293
fix cfd4 all couple/cfd/parttempfield T0 293
fix cfd5 all couple/cfd/deform rmin 0.01
# make sure particle outside CFD domain are initialized with reasonable temperature
# this is executed at the beginning of each run
# this is executed at the beginning of each run
region outsideCFD block INF INF INF INF 26.0 INF units box
variable NFixTemp equal ${Ncouple}+2
fix eeTemp all execute ${NFixTemp} " set region outsideCFD property/atom Temp ${Tin} " once
fix eeTemp all execute ${NFixTemp} " set region outsideCFD property/atom Temp ${Tin} " once
fix integr all nve/sphere
fix limV1 all limit/vel 1.0
fix limV2 RW limit/vel 0.01
fix integr all nve/sphere
fix limV1 all limit/vel 1.0
fix limV2 RW limit/vel 0.01

24
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/DEM/in.liggghts_CFDDEM
View File

@ -3,23 +3,23 @@ include ../DEM/DEMparameters
###############################################
atom_style granular
atom_modify map array
communicate single vel yes
atom_style granular
atom_modify map array
communicate single vel yes
boundary f f f
newton off
units si
boundary f f f
newton off
units si
processors 2 2 8
read_restart ../DEM/liggghts.restart
reset_timestep 0
read_restart ../DEM/liggghts.restart
reset_timestep 0
#read_restart ../DEM/restart.B
#reset_timestep 400000
#read_restart ../DEM/restart.B
#reset_timestep 400000
neighbor ${skin} bin
neigh_modify delay 0
neighbor ${skin} bin
neigh_modify delay 0
timestep ${dt}

42
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/DEM/in.material_settings
View File

@ -2,28 +2,28 @@
# material properties required for new pair styles
##################################################
fix m1 all property/global youngsModulus peratomtype ${yM_C} ${yM_C} ${yM_O} ${yM_O} ${yM_W}
fix m2 all property/global poissonsRatio peratomtype ${pR_C} ${pR_C} ${pR_O} ${pR_O} ${pR_W}
fix m3 all property/global coefficientRestitution peratomtypepair 5 ${rest_C} ${rest_C} ${rest_CO} ${rest_CO} ${rest_CW} ${rest_C} ${rest_C} ${rest_CO} ${rest_CO} ${rest_CW} ${rest_CO} ${rest_CO} ${rest_O} ${rest_O} ${rest_WO} ${rest_CO} ${rest_CO} ${rest_O} ${rest_O} ${rest_WO} ${rest_CW} ${rest_CW} ${rest_WO} ${rest_WO} ${rest_W}
fix m4 all property/global coefficientFriction peratomtypepair 5 ${fric_C} ${fric_C} ${fric_CO} ${fric_CO} ${fric_CW} ${fric_C} ${fric_C} ${fric_CO} ${fric_CO} ${fric_CW} ${fric_CO} ${fric_CO} ${fric_O} ${fric_O} ${fric_WO} ${fric_CO} ${fric_CO} ${fric_O} ${fric_O} ${fric_WO} ${fric_CW} ${fric_CW} ${fric_WO} ${fric_WO} ${fric_W}
fix m5 all property/global coefficientRollingFriction peratomtypepair 5 ${roll_C} ${roll_C} ${roll_CO} ${roll_CO} ${roll_CW} ${roll_C} ${roll_C} ${roll_CO} ${roll_CO} ${roll_CW} ${roll_CO} ${roll_CO} ${roll_O} ${roll_O} ${roll_WO} ${roll_CO} ${roll_CO} ${roll_O} ${roll_O} ${roll_WO} ${roll_CW} ${roll_CW} ${roll_WO} ${roll_WO} ${roll_W}
fix m6 all property/global characteristicVelocity scalar ${cV}
fix m1 all property/global youngsModulus peratomtype ${yM_C} ${yM_C} ${yM_O} ${yM_O} ${yM_W}
fix m2 all property/global poissonsRatio peratomtype ${pR_C} ${pR_C} ${pR_O} ${pR_O} ${pR_W}
fix m3 all property/global coefficientRestitution peratomtypepair 5 ${rest_C} ${rest_C} ${rest_CO} ${rest_CO} ${rest_CW} ${rest_C} ${rest_C} ${rest_CO} ${rest_CO} ${rest_CW} ${rest_CO} ${rest_CO} ${rest_O} ${rest_O} ${rest_WO} ${rest_CO} ${rest_CO} ${rest_O} ${rest_O} ${rest_WO} ${rest_CW} ${rest_CW} ${rest_WO} ${rest_WO} ${rest_W}
fix m4 all property/global coefficientFriction peratomtypepair 5 ${fric_C} ${fric_C} ${fric_CO} ${fric_CO} ${fric_CW} ${fric_C} ${fric_C} ${fric_CO} ${fric_CO} ${fric_CW} ${fric_CO} ${fric_CO} ${fric_O} ${fric_O} ${fric_WO} ${fric_CO} ${fric_CO} ${fric_O} ${fric_O} ${fric_WO} ${fric_CW} ${fric_CW} ${fric_WO} ${fric_WO} ${fric_W}
fix m5 all property/global coefficientRollingFriction peratomtypepair 5 ${roll_C} ${roll_C} ${roll_CO} ${roll_CO} ${roll_CW} ${roll_C} ${roll_C} ${roll_CO} ${roll_CO} ${roll_CW} ${roll_CO} ${roll_CO} ${roll_O} ${roll_O} ${roll_WO} ${roll_CO} ${roll_CO} ${roll_O} ${roll_O} ${roll_WO} ${roll_CW} ${roll_CW} ${roll_WO} ${roll_WO} ${roll_W}
fix m6 all property/global characteristicVelocity scalar ${cV}
fix ftco all property/global thermalConductivity peratomtype 1.7 1.7 80 80 1.0
fix ftca all property/global thermalCapacity peratomtype 850 850 600 600 1.0
fix heattransfer all heat/gran initial_temperature 293 contact_area projection
fix ftco all property/global thermalConductivity peratomtype 1.7 1.7 80 80 1.0
fix ftca all property/global thermalCapacity peratomtype 850 850 600 600 1.0
fix heattransfer all heat/gran initial_temperature 293 contact_area projection
# pair style
pair_style gran model hertz tangential history rolling_friction cdt
pair_coeff * *
# gravity
fix gravi all gravity 9.81 vector 0.0 0.0 -1.0
fix gravi all gravity 9.81 vector 0.0 0.0 -1.0
fix bottom all mesh/surface file ../../CFDDEM/DEM/stl_files/surface_bottom.stl type 5
fix wall all mesh/surface file ../../CFDDEM/DEM/stl_files/surface_walls.stl type 5
fix bottom all mesh/surface file ../../CFDDEM/DEM/stl_files/surface_bottom.stl type 5
fix wall all mesh/surface file ../../CFDDEM/DEM/stl_files/surface_walls.stl type 5
fix granwalls all wall/gran model hertz tangential history mesh n_meshes 2 meshes bottom wall
fix granwalls all wall/gran model hertz tangential history mesh n_meshes 2 meshes bottom wall
#################################
@ -35,14 +35,14 @@ group finecoke type 2
group ore type 3
group fineore type 4
fix pts1 all particletemplate/sphere 2 atom_type 1 density constant ${dens_C} radius constant ${rad_C}
fix pdd1 all particledistribution/discrete 3 1 pts1 1.0
fix pts1 all particletemplate/sphere 2 atom_type 1 density constant ${dens_C} radius constant ${rad_C}
fix pdd1 all particledistribution/discrete 3 1 pts1 1.0
fix pts2 all particletemplate/sphere 5 atom_type 3 density constant ${dens_O} radius constant ${rad_O}
fix pdd2 all particledistribution/discrete 7 1 pts2 1.0
fix pts2 all particletemplate/sphere 5 atom_type 3 density constant ${dens_O} radius constant ${rad_O}
fix pdd2 all particledistribution/discrete 7 1 pts2 1.0
fix pts3a all particletemplate/sphere 11 atom_type 2 density constant ${dens_fC} radius constant ${rad_fC}
fix pts3b all particletemplate/sphere 13 atom_type 4 density constant ${dens_fO} radius constant ${rad_fO}
fix pdd3 all particledistribution/discrete 17 2 pts3a 0.2 pts3b 0.8
fix pts3a all particletemplate/sphere 11 atom_type 2 density constant ${dens_fC} radius constant ${rad_fC}
fix pts3b all particletemplate/sphere 13 atom_type 4 density constant ${dens_fO} radius constant ${rad_fO}
fix pdd3 all particledistribution/discrete 17 2 pts3a 0.2 pts3b 0.8
fix fp all property/atom/polydispparcel effvolfactor ${poly_C} ${poly_fC} ${poly_O} ${poly_fO} 1.0
fix fp all property/atom/polydispparcel effvolfactor ${poly_C} ${poly_fC} ${poly_O} ${poly_fO} 1.0

12
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/DEM/in.output
View File

@ -3,10 +3,10 @@
########
#screen output
compute 1 all erotate/sphere
thermo_style custom step atoms ke c_1 vol
thermo 10000
thermo_modify lost ignore norm no
compute_modify thermo_temp dynamic yes
compute 1 all erotate/sphere
thermo_style custom step atoms ke c_1 vol
thermo 10000
thermo_modify lost ignore norm no
compute_modify thermo_temp dynamic yes
dump dmp1 all custom/vtk ${Ndump} ../DEM/post/dump*.CFDDEM.vtk id type radius x y z vx vy vz fx fy fz mass f_Temp[0] f_partDeformations[0]
dump dmp1 all custom/vtk ${Ndump} ../DEM/post/dump*.CFDDEM.vtk id type radius x y z vx vy vz fx fy fz mass f_Temp[0] f_partDeformations[0]

1
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/DEM/in.regions
View File

@ -14,4 +14,3 @@ group DMcoke region DM
variable DMheight05 equal ${DMheight}+0.5
region hearth block -7 7 -7 7 -1.0 ${DMheight05}

0
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/DEM/monitor/.gitkeep → tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/CFDDEM/DEM/monitor/.gitignore vendored
View File

12
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/constant/RASProperties
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -15,7 +15,7 @@ FoamFile
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
RASModel laminar;//kEpsilon;
RASModel laminar;
turbulence off;

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/constant/chemistryProperties
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/constant/combustionProperties
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

47
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/constant/couplingProperties
View File

@ -1,29 +1,19 @@
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.4 |
| \\ / A nd | Web: http://www.openfoam.org |
| \\/ M anipulation | |
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
root "";
case "";
instance "";
local "";
class dictionary;
object couplingProperties;
version 2.0;
format ascii;
class dictionary;
object couplingProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
//===========================================================================//
// sub-models & settings
modelType "A";
@ -87,7 +77,6 @@ chemistryModels
turbulenceModelType "turbulenceProperties";
getParticleTypes true;
getParticleDensities true;
getParticleEffVolFactors true;
@ -97,7 +86,8 @@ coupleRecForce false;
setForceField true;
nEveryFlow 1;
//===========================================================================//
// sub-model properties
implicitCoupleProps
@ -121,7 +111,6 @@ viscForceProps
densityFieldName "rho";
}
BeetstraDragProps
{
velFieldName "U";
@ -169,13 +158,10 @@ directedDiffusiveRelaxationProps
particleDeformationProps
{
refFieldName "partTempRef";
// partTypes (3 4);
partTypes (1 2 3 4);
defaultDeformCellsName "noOreReg";
// same as Huo et al. (2017)
// lowerBounds (1475 1475);
// upperBounds (1675 1675);
lowerBounds (1475 1475 1475 1475);
lowerBounds (1475 1475 1475 1475);
upperBounds (1675 1675 1675 1675);
// backgroundVoidage 0.05; // ATTENTION: if coke is also deformed, full hearth gets background voidage
voidfractionFieldName "voidfractionRec";
@ -233,7 +219,7 @@ diffusionCoefficientsProps
ChemistryFile "$FOAM_CASE/constant/foam.inp";
diffusantGasNames (
CO CO2 O2
);
);
}
massTransferCoeffProps
@ -267,11 +253,6 @@ constDiffSmoothingProps
smoothingLength 0.5;
}
centreProps
{
alphaMin 0.10;
}
engineProps
{
treeSearch true;

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/constant/foam_hconst.dat
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

11
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/constant/g
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -18,5 +18,4 @@ FoamFile
dimensions [0 1 -2 0 0 0 0];
value ( 0 0 -9.81 );
// ************************************************************************* //

29
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/constant/liggghtsCommands
View File

@ -1,26 +1,17 @@
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.4 |
| \\ / A nd | Web: http://www.openfoam.org |
| \\/ M anipulation | |
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
root "";
case "";
instance "";
local "";
class dictionary;
object liggghtsCommands;
version 2.0;
format ascii;
class dictionary;
object liggghtsCommands;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
liggghtsCommandModels

36
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/constant/recProperties
View File

@ -1,29 +1,19 @@
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.4 |
| \\ / A nd | Web: http://www.openfoam.org |
| \\/ M anipulation | |
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
root "";
case "";
instance "";
local "";
class dictionary;
object recProperties;
version 2.0;
format ascii;
class dictionary;
object recProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
//===========================================================================//
// sub-models & settings
@ -50,14 +40,9 @@ surfaceScalarFields
// phi
);
//verbose true;
//couplingSubStep 1;
initialRecSteps 20000; //10;
//===========================================================================//
// sub-model properties
standardRecModelProps
@ -70,7 +55,6 @@ standardRecModelProps
EuclideanNormProps
{
fieldType "volScalarField";
// fieldName "partType";
fieldName "voidfraction";
}

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/constant/symmetrizerProperties_1
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/constant/symmetrizerProperties_2
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/constant/symmetrizerProperties_3
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/constant/symmetrizerProperties_4
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

22
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/constant/thermophysicalProperties
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -14,18 +14,6 @@ FoamFile
object thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
/*
thermoType
{
type heRhoThermo; //hePsiThermo;
mixture reactingMixture; //species and reactions are listed in chemistry file
transport const; //sutherland; //calcualtes viscosity as a function of temperature from Stuherland coefficient As and Ts
thermo hConst; //janaf;
equationOfState perfectGas;
energy sensibleInternalEnergy; //sensibleEnthalpy;
specie specie;
}
*/
thermoType
{

14
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/constant/transportProperties
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -19,9 +19,9 @@ transportModel Newtonian;
nu nu [ 0 2 -1 0 0 0 0 ] 1e-05;
kf kf [ 1 1 -3 -1 0 0 0 ] 0.026;
kf kf [ 1 1 -3 -1 0 0 0 ] 0.026;
Cp Cp [ 0 2 -2 -1 0 0 0 ] 1005;
Cp Cp [ 0 2 -2 -1 0 0 0 ] 1005;
Cv Cv [ 0 2 -2 -1 0 0 0 ] 718;

11
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/constant/turbulenceProperties
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -17,5 +17,4 @@ FoamFile
simulationType laminar;
// ************************************************************************* //

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/orig.0/CO
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/orig.0/CO2
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/orig.0/Ksl
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/orig.0/N2
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/orig.0/O2
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

16
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/orig.0/T
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -28,19 +28,19 @@ boundaryField
bottom
{
type fixedValue;
value uniform 293;
value uniform 293;
}
walls
{
type fixedValue;
value uniform 293;
value uniform 293;
}
inlet
{
type fixedValue;
value uniform 1600;
value uniform 1600;
}
}

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/orig.0/U
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/orig.0/Us
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

14
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/orig.0/addSource
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -14,9 +14,9 @@ FoamFile
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -3 0 0 0 0];
dimensions [1 -1 -3 0 0 0 0];
internalField uniform 0;
internalField uniform 0;
boundaryField
{

12
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/orig.0/p
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -24,7 +24,7 @@ boundaryField
{
type zeroGradient;
}
top
{
type fixedValue;

12
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/orig.0/partEffThermCondField
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -33,7 +33,7 @@ boundaryField
walls
{
type fixedValue;
value uniform 2.0;
value uniform 2.0;
}
inlet

2
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/orig.0/partTemp
View File

@ -17,7 +17,7 @@ FoamFile
dimensions [0 0 0 1 0 0 0];
internalField nonuniform List<scalar>
internalField nonuniform List<scalar>
144528
(
2700

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/orig.0/particleForces
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/orig.0/thCond
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/orig.0/voidfraction
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

10
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/orig.0/wallHeatLoss
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{

136
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/system/controlDict
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -79,14 +79,14 @@ functions
base time;
window 60.0001;
}
Us
Us
{
mean on;
prime2Mean off;
base time;
window 60.0001;
}
voidfraction
voidfraction
{
mean on;
prime2Mean off;
@ -117,74 +117,74 @@ functions
);
}
probes1
{
type probes;
probes1
{
type probes;
functionObjectLibs ("libsampling.so");
functionObjectLibs ("libsampling.so");
#include "probesDict1";
}
probes2
{
type probes;
#include "probesDict1";
}
functionObjectLibs ("libsampling.so");
probes2
{
type probes;
#include "probesDict2";
}
functionObjectLibs ("libsampling.so");
massInflow
{
type surfaceFieldValue;
libs ("libfieldFunctionObjects.so");
writeControl timeStep;
writeInterval 20;
log true;
writeFields false;
regionType patch;
name inlet;
operation sum;
fields
(
phi
);
}
#include "probesDict2";
}
massOutflow
{
$massInflow;
name top;
}
massInflow
{
type surfaceFieldValue;
libs ("libfieldFunctionObjects.so");
writeControl timeStep;
writeInterval 20;
log true;
writeFields false;
regionType patch;
name inlet;
operation sum;
fields
(
phi
);
}
energyAndSpeciesInflow
{
type surfaceFieldValue;
libs ("libfieldFunctionObjects.so");
writeControl timeStep;
writeInterval 20;
log true;
writeFields false;
regionType patch;
name inlet;
operation weightedSum;
weightField phi;
fields
(
e
CO
CO2
N2
O2
);
}
massOutflow
{
$massInflow;
name top;
}
energyAndSpeciesOutflow
{
$energyAndSpeciesInflow;
name top;
}
energyAndSpeciesInflow
{
type surfaceFieldValue;
libs ("libfieldFunctionObjects.so");
writeControl timeStep;
writeInterval 20;
log true;
writeFields false;
regionType patch;
name inlet;
operation weightedSum;
weightField phi;
fields
(
e
CO
CO2
N2
O2
);
}
energyAndSpeciesOutflow
{
$energyAndSpeciesInflow;
name top;
}
wallheatFluxes
{
@ -204,7 +204,7 @@ functions
QPartPart
);
}
totReactionHeatField
{
type volFieldValue;

14
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/system/createPatchDict
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -48,7 +48,7 @@ pointSync false;
// Patches to create.
patches
(
{
{
name inletLeft;
patchInfo
{
@ -58,7 +58,7 @@ patches
set inletLeft;
}
{
{
name inletRight;
patchInfo
{

16
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/system/decomposeParDict
View File

@ -1,10 +1,10 @@
/*-------------------------------*- C++ -*---------------------------------*\
| ========= |
| \\ / OpenFOAM |
| \\ / |
| \\ / The Open Source CFD Toolbox |
| \\/ http://www.OpenFOAM.org |
\*-------------------------------------------------------------------------*/
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
@ -22,7 +22,7 @@ numberOfSubdomains 32;
method scotch; //simple;
simpleCoeffs
{
{
n (2 2 8);
delta 0.001;
}

66
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/system/fvOptions
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -17,7 +17,7 @@ FoamFile
source1
{
type limitTemperature;
type limitTemperature;
active yes;
limitTemperatureCoeffs
{
@ -28,57 +28,3 @@ source1
}
}
/*
source2
{
type fixedTemperatureConstraint;
selectionMode cellSet;
cellSet rw_left;
mode uniform;
temperature 3450;
}
source3
{
type fixedTemperatureConstraint;
selectionMode cellSet;
cellSet rw_right;
mode uniform;
temperature 3450;
}
*/
/*
heatSource1
{
type scalarSemiImplicitSource;
selectionMode cellSet;
cellSet rw_left;
volumeMode absolute; //specific;
injectionRateSuSp
{
// dm/dt = 1.65 kg / s / RW, C = 850, T = 2100
e (3e6 0);
}
}
heatSource2
{
type scalarSemiImplicitSource;
selectionMode cellSet;
cellSet rw_right;
volumeMode absolute; //specific;
injectionRateSuSp
{
// dm/dt = 1.65 kg / s / RW, C = 850, T = 2100
e (3e6 0);
}
}
*/

14
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/system/fvSchemes
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -24,7 +24,7 @@ gradSchemes
{
default Gauss linear;
grad(p) cellLimited Gauss linear 1;
// grad(U) Gauss linear;
//grad(U) Gauss linear;
grad(U) cellLimited Gauss linear 1;
}
@ -79,7 +79,7 @@ interpolationSchemes
snGradSchemes
{
// default corrected;
//default corrected;
default limited 0.333;
}

36
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/system/fvSolution
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -70,7 +70,7 @@ solvers
relTol 0;
}
"(voidfraction|Ksl|UsNext|voidfractionNext|sSmoothField|vSmoothField|fSmooth)"
{
@ -111,10 +111,10 @@ relaxationFactors
{
fields
{
U 0.3;
p 0.3;
UFinal 0.3;
pFinal 0.3;
U 0.3;
p 0.3;
UFinal 0.3;
pFinal 0.3;
CO 0.3;
COFinal 0.3;
CO2 0.3;
@ -123,10 +123,10 @@ relaxationFactors
O2Final 0.3;
N2 0.3;
N2Final 0.3;
e 0.3;
he 0.3;
eFinal 0.3;
heFinal 0.3;
e 0.3;
he 0.3;
eFinal 0.3;
heFinal 0.3;
}
equations
{
@ -138,10 +138,10 @@ relaxationFactors
// O2Final 0.3;
// N2 0.3;
// N2Final 0.3;
// e 0.3;
// he 0.3;
// eFinal 0.3;
// heFinal 0.3;
// e 0.3;
// he 0.3;
// eFinal 0.3;
// heFinal 0.3;
".*" 0.3;
}
}

6
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/system/monitorPressureWork
View File

@ -36,14 +36,14 @@
forAll(patches, patchi)
{
const fvPatch& curPatch = patches[patchi];
word patchName = mesh().boundary()[patchi].name();
if (patchName == "inlet")
{
forAll(curPatch, facei)
{
label faceCelli = curPatch.faceCells()[facei];
inpU += phi.boundaryField()[patchi][facei] * p.boundaryField()[patchi][facei] / rho.boundaryField()[patchi][facei];
inpU += phi.boundaryField()[patchi][facei] * p.boundaryField()[patchi][facei] / rho.boundaryField()[patchi][facei];
}
}
else if (patchName == "top")
@ -51,7 +51,7 @@
forAll(curPatch, facei)
{
label faceCelli = curPatch.faceCells()[facei];
outpU += phi.boundaryField()[patchi][facei] * p.boundaryField()[patchi][facei] / rho.boundaryField()[patchi][facei];
outpU += phi.boundaryField()[patchi][facei] * p.boundaryField()[patchi][facei] / rho.boundaryField()[patchi][facei];
}
}
}

66
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/system/probesDict1
View File

@ -1,21 +1,19 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 2.1.x |
| \ / A nd | Web: www.OpenFOAM.org |
| \/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object probesDict;
version 2.0;
format ascii;
class dictionary;
object probesDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
fields
(
voidfraction
@ -43,30 +41,28 @@ writeInterval 10;
// Locations to be probed.
probeLocations
(
(0.05 0.05 -1.0)
(0.05 0.05 0.0)
(0.05 0.05 1.0)
(0.05 0.05 2.0)
(0.05 0.05 3.0)
(0.05 0.05 4.0)
(0.05 0.05 5.0)
(0.05 0.05 6.0)
(0.05 0.05 7.0)
(0.05 0.05 8.0)
(0.05 0.05 9.0)
(0.05 0.05 10.0)
(0.05 0.05 11.0)
(0.05 0.05 12.0)
(0.05 0.05 13.0)
(0.05 0.05 14.0)
(0.05 0.05 15.0)
(0.05 0.05 16.0)
(0.05 0.05 17.0)
(0.05 0.05 18.0)
(0.05 0.05 19.0)
(0.05 0.05 20.0)
(0.05 0.05 -1.0)
(0.05 0.05 0.0)
(0.05 0.05 1.0)
(0.05 0.05 2.0)
(0.05 0.05 3.0)
(0.05 0.05 4.0)
(0.05 0.05 5.0)
(0.05 0.05 6.0)
(0.05 0.05 7.0)
(0.05 0.05 8.0)
(0.05 0.05 9.0)
(0.05 0.05 10.0)
(0.05 0.05 11.0)
(0.05 0.05 12.0)
(0.05 0.05 13.0)
(0.05 0.05 14.0)
(0.05 0.05 15.0)
(0.05 0.05 16.0)
(0.05 0.05 17.0)
(0.05 0.05 18.0)
(0.05 0.05 19.0)
(0.05 0.05 20.0)
);
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

69
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/system/probesDict2
View File

@ -1,21 +1,19 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 2.1.x |
| \ / A nd | Web: www.OpenFOAM.org |
| \/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object probesDict;
version 2.0;
format ascii;
class dictionary;
object probesDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
fields
(
voidfraction
@ -42,31 +40,30 @@ writeInterval 10;
// Locations to be probed.
probeLocations
(
(0.05 -0.05 -1.0)
(0.05 -0.05 0.0)
(0.05 -0.05 1.0)
(0.05 -0.05 2.0)
(0.05 -0.05 3.0)
(0.05 -0.05 4.0)
(0.05 -0.05 5.0)
(0.05 -0.05 6.0)
(0.05 -0.05 7.0)
(0.05 -0.05 8.0)
(0.05 -0.05 9.0)
(0.05 -0.05 10.0)
(0.05 -0.05 10.0)
(0.05 -0.05 11.0)
(0.05 -0.05 12.0)
(0.05 -0.05 13.0)
(0.05 -0.05 14.0)
(0.05 -0.05 15.0)
(0.05 -0.05 16.0)
(0.05 -0.05 17.0)
(0.05 -0.05 18.0)
(0.05 -0.05 19.0)
(0.05 -0.05 20.0)
(
(0.05 -0.05 -1.0)
(0.05 -0.05 0.0)
(0.05 -0.05 1.0)
(0.05 -0.05 2.0)
(0.05 -0.05 3.0)
(0.05 -0.05 4.0)
(0.05 -0.05 5.0)
(0.05 -0.05 6.0)
(0.05 -0.05 7.0)
(0.05 -0.05 8.0)
(0.05 -0.05 9.0)
(0.05 -0.05 10.0)
(0.05 -0.05 10.0)
(0.05 -0.05 11.0)
(0.05 -0.05 12.0)
(0.05 -0.05 13.0)
(0.05 -0.05 14.0)
(0.05 -0.05 15.0)
(0.05 -0.05 16.0)
(0.05 -0.05 17.0)
(0.05 -0.05 18.0)
(0.05 -0.05 19.0)
(0.05 -0.05 20.0)
);
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

14
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/system/topoSetDictCZ
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -12,11 +12,11 @@ FoamFile
class dictionary;
object topoSetDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
actions
(
{
{
name aboveCZ;
type cellSet;
action new;

30
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/CFD/system/topoSetDictMonitor
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -12,11 +12,11 @@ FoamFile
class dictionary;
object topoSetDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
actions
(
{
{
name sphere0;
type cellSet;
action new;
@ -26,9 +26,9 @@ actions
centre (0.0 0.0 2.0);
radius 0.3;
}
}
}
{
{
name sphere1;
type cellSet;
action new;
@ -38,9 +38,9 @@ actions
centre (0.0 0.0 2.0);
radius 0.45;
}
}
}
{
{
name sphere2;
type cellSet;
action new;
@ -50,9 +50,9 @@ actions
centre (0.0 0.0 2.0);
radius 0.55;
}
}
}
{
{
name sphere3;
type cellSet;
action new;
@ -62,7 +62,7 @@ actions
centre (0.0 0.0 2.0);
radius 0.65;
}
}
}
{
name sphere0;
@ -111,8 +111,6 @@ actions
set sphere3;
}
}
);
// ************************************************************************* //

70
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/DEM/DEMparameters
View File

@ -28,49 +28,49 @@ variable roll_C equal 0.5 # coefficient of rolling friction
# mean density coke = 950
# bulk density = 500
variable dens_C equal 950 # density
variable rad_C equal 0.088 # c.g. radius (cg = 4)
variable poly_C equal 0.81 # polydisp factor
variable dens_C equal 950 # density
variable rad_C equal 0.088 # c.g. radius (cg = 4)
variable poly_C equal 0.81 # polydisp factor
# DEM-parameters FINE COKE
variable yM_fC equal 5.e7 # Young's modulus
variable pR_fC equal 0.25 # Poisson ratio
variable rest_fC equal 0.5 # coefficient of restitution
variable fric_fC equal 0.9 # coefficient of friction
variable roll_fC equal 0.5 # coefficient of rolling friction
variable yM_fC equal 5.e7 # Young's modulus
variable pR_fC equal 0.25 # Poisson ratio
variable rest_fC equal 0.5 # coefficient of restitution
variable fric_fC equal 0.9 # coefficient of friction
variable roll_fC equal 0.5 # coefficient of rolling friction
# mean diameter fine coke = 0.02 m
# mean density fine coke = 950
# bulk density fine coke = 570
variable dens_fC equal 950 # density
variable rad_fC equal 0.08 # c.g. radius (cg = 8)
#variable poly_fC equal 0.92 # polydisp factor
variable poly_fC equal 1.2 # polydisp factor accounting for near-wall inaccuracy
variable dens_fC equal 950 # density
variable rad_fC equal 0.08 # c.g. radius (cg = 8)
#variable poly_fC equal 0.92 # polydisp factor
variable poly_fC equal 1.2 # polydisp factor accounting for near-wall inaccuracy
# DEM-parameters ORE (represents the mixture of sinter, pellets and actual ore)
variable yM_O equal 5.e7 # Young's modulus
variable pR_O equal 0.25 # Poisson ratio
variable rest_O equal 0.5 # coefficient of restitution
variable fric_O equal 0.8 # coefficient of friction
variable roll_O equal 0.4 # coefficient of rolling friction
variable roll_O equal 0.4 # coefficient of rolling friction
# values according to Modern Blast Furnace Ironmaking
# mean diameter ore mixture = 0.015 m
# mean density ore mixture = 3300
# bulk density = 2200
variable dens_O equal 3300 # density
variable rad_O equal 0.0825 # c.g. radius (cg = 11)
variable poly_O equal 1.03 # polydisp factor
variable dens_O equal 3300 # density
variable rad_O equal 0.0825 # c.g. radius (cg = 11)
variable poly_O equal 1.03 # polydisp factor
# DEM-parameters FINE ORE (represents the mixture of sinter, pellets and actual ore)
variable yM_fO equal 5.e7 # Young's modulus
variable pR_fO equal 0.25 # Poisson ratio
variable rest_fO equal 0.5 # coefficient of restitution
variable fric_fO equal 0.8 # coefficient of friction
variable yM_fO equal 5.e7 # Young's modulus
variable pR_fO equal 0.25 # Poisson ratio
variable rest_fO equal 0.5 # coefficient of restitution
variable fric_fO equal 0.8 # coefficient of friction
variable roll_fO equal 0.4 # coefficient of rolling friction
# values according to G. Holzinger's literature study
@ -78,10 +78,10 @@ variable roll_fO equal 0.4 # coefficient of rolling friction
# mean density fine ore mixture = 3300
# bulk density = 2200
variable dens_fO equal 3300 # density
variable rad_fO equal 0.072 # c.g. radius (cg = 24)
#variable poly_fO equal 1.03 # polydisp factor
variable poly_fO equal 1.2 # polydisp factor accounting for near-wall inaccuracy
variable dens_fO equal 3300 # density
variable rad_fO equal 0.072 # c.g. radius (cg = 24)
#variable poly_fO equal 1.03 # polydisp factor
variable poly_fO equal 1.2 # polydisp factor accounting for near-wall inaccuracy
@ -89,19 +89,19 @@ variable poly_fO equal 1.2 # polydisp factor accounting for near-wall inac
variable yM_W equal 5.e6 # Young's modulus
variable pR_W equal 0.25 # Poisson ratio
variable rest_W equal 0.5 # coefficient of restitution
variable fric_W equal 0.85 # coefficient of friction
variable fric_W equal 0.85 # coefficient of friction
variable roll_W equal 0.45 # coefficient of rolling friction
# DEM-parameters mixed
variable rest_CO equal 0.2 # coefficient of restitution
variable rest_WO equal 0.2 # coefficient of restitution
variable rest_CW equal 0.5 # coefficient of restitution
variable fric_CO equal 0.8 # coefficient of friction
variable fric_WO equal 0.2 # coefficient of friction
variable fric_CW equal 0.5 # coefficient of friction
variable roll_CO equal 0.02 # for particles and walls with friction
variable roll_WO equal 0.02 # for particles and walls with friction
variable roll_CW equal 0.02 # for particles and walls with friction
variable rest_CO equal 0.2 # coefficient of restitution
variable rest_WO equal 0.2 # coefficient of restitution
variable rest_CW equal 0.5 # coefficient of restitution
variable fric_CO equal 0.8 # coefficient of friction
variable fric_WO equal 0.2 # coefficient of friction
variable fric_CW equal 0.5 # coefficient of friction
variable roll_CO equal 0.02 # for particles and walls with friction
variable roll_WO equal 0.02 # for particles and walls with friction
variable roll_CW equal 0.02 # for particles and walls with friction
variable cV equal 0.1 # characteristic velocity
variable cV equal 0.1 # characteristic velocity

8
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/DEM/in.chemistry
View File

@ -61,17 +61,17 @@ fix Ea_chem1 fineore property/global Ea_chem1 vector 69488 73674 113859
# values for k0 and T0 according to Shen, Fuel 90 (2011) 728–738;
# values have been scaled with 1/(nA * mMol) [nA = molar concentration at standard conditions], b.c. Shen uses x P / PA while LIGGGHTS implementation has x P / (RT) * mMol = x P / PA * nA * mMol
fix chem2 coke chem/shrink/core/single speciesA O2 molMassA 0.03199 nuA 1 speciesC CO molMassC 0.02801 nuC 2 molMassB 0.01201 nuB 2 molMassD 0.001201 nuD 2 rmin 0.0075 nevery ${Nchem} screen no k 3.1e6 T 10855 Tmin 100 cg 4 layerDiffusion no heatToFluid yes shrink no
fix chem3 coke chem/shrink/core/single speciesA CO2 molMassA 0.04401 nuA 1 speciesC CO molMassC 0.02801 nuC 2 molMassB 0.01201 nuB 1 molMassD 0.001201 nuD 1 rmin 0.0075 nevery ${Nchem} screen no k 9.4e9 T 29018 Tmin 100 cg 4 layerDiffusion no heatToFluid yes shrink no
fix chem3 coke chem/shrink/core/single speciesA CO2 molMassA 0.04401 nuA 1 speciesC CO molMassC 0.02801 nuC 2 molMassB 0.01201 nuB 1 molMassD 0.001201 nuD 1 rmin 0.0075 nevery ${Nchem} screen no k 9.4e9 T 29018 Tmin 100 cg 4 layerDiffusion no heatToFluid yes shrink no
fix chem4 finecoke chem/shrink/core/single speciesA O2 molMassA 0.03199 nuA 1 speciesC CO molMassC 0.02801 nuC 2 molMassB 0.01201 nuB 2 molMassD 0.001201 nuD 2 rmin 0.0075 nevery ${Nchem} screen no k 3.1e6 T 10855 Tmin 100 cg 8 layerDiffusion no heatToFluid yes shrink no
fix chem5 finecoke chem/shrink/core/single speciesA CO2 molMassA 0.04401 nuA 1 speciesC CO molMassC 0.02801 nuC 2 molMassB 0.01201 nuB 1 molMassD 0.001201 nuD 1 rmin 0.0075 nevery ${Nchem} screen no k 9.4e9 T 29018 Tmin 100 cg 8 layerDiffusion no heatToFluid yes shrink no
fix chem5 finecoke chem/shrink/core/single speciesA CO2 molMassA 0.04401 nuA 1 speciesC CO molMassC 0.02801 nuC 2 molMassB 0.01201 nuB 1 molMassD 0.001201 nuD 1 rmin 0.0075 nevery ${Nchem} screen no k 9.4e9 T 29018 Tmin 100 cg 8 layerDiffusion no heatToFluid yes shrink no
######
# Zhuo and Shen, Applied Energy 261 (2020) 114456
# with the additional assumption that all O2 is transformed into CO upon reaction with coke
######
#fix chem2 coke chem/shrink/core/single speciesA O2 molMassA 0.03199 nuA 1 speciesC CO molMassC 0.02801 nuC 2 molMassB 0.01201 nuB 2 molMassD 0.001201 nuD 2 rmin 0.0075 nevery ${Nchem} screen no nPreFactor 1 k 0.685 T 9000 Tmin 100 cg 4 layerDiffusion no heatToFluid yes shrink no
#fix chem3 coke chem/shrink/core/single speciesA CO2 molMassA 0.04401 nuA 1 speciesC CO molMassC 0.02801 nuC 2 molMassB 0.01201 nuB 1 molMassD 0.001201 nuD 1 rmin 0.0075 nevery ${Nchem} screen no nPreFactor 1 k 589 T 26800 Tmin 100 cg 4 layerDiffusion no heatToFluid yes shrink no
#fix chem3 coke chem/shrink/core/single speciesA CO2 molMassA 0.04401 nuA 1 speciesC CO molMassC 0.02801 nuC 2 molMassB 0.01201 nuB 1 molMassD 0.001201 nuD 1 rmin 0.0075 nevery ${Nchem} screen no nPreFactor 1 k 589 T 26800 Tmin 100 cg 4 layerDiffusion no heatToFluid yes shrink no
#fix chem4 finecoke chem/shrink/core/single speciesA O2 molMassA 0.03199 nuA 1 speciesC CO molMassC 0.02801 nuC 2 molMassB 0.01201 nuB 2 molMassD 0.001201 nuD 2 rmin 0.0075 nevery ${Nchem} screen no nPreFactor 1 k 0.685 T 9000 Tmin 100 cg 8 layerDiffusion no heatToFluid yes shrink no
#fix chem5 finecoke chem/shrink/core/single speciesA CO2 molMassA 0.04401 nuA 1 speciesC CO molMassC 0.02801 nuC 2 molMassB 0.01201 nuB 1 molMassD 0.001201 nuD 1 rmin 0.0075 nevery ${Nchem} screen no nPreFactor 1 k 589 T 26800 Tmin 100 cg 8 layerDiffusion no heatToFluid yes shrink no
#fix chem5 finecoke chem/shrink/core/single speciesA CO2 molMassA 0.04401 nuA 1 speciesC CO molMassC 0.02801 nuC 2 molMassB 0.01201 nuB 1 molMassD 0.001201 nuD 1 rmin 0.0075 nevery ${Nchem} screen no nPreFactor 1 k 589 T 26800 Tmin 100 cg 8 layerDiffusion no heatToFluid yes shrink no

22
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/DEM/in.coupling_integration
View File

@ -9,18 +9,18 @@ variable Tin equal 293
#fix execSetV all execute ${Ncouple} "set group fullyDeformed property/atom vrec 0.0 0.0 ${vtrickle}" start_of_step
#####
fix cfd all couple/cfd couple_every ${Ncouple} mpi
fix cfd1 all couple/cfd/force transfer_force no transfer_torque no transfer_density yes transfer_type yes transfer_property name effvolfactor type scalar-atom
fix cfd2 all couple/cfd/recurrence transfer_force no transfer_fluctuations no
fix cfd3 all couple/cfd/convection T0 ${Tin} transfer_conduction yes max_change 1.0
fix cfd4 all couple/cfd/parttempfield T0 ${Tin}
fix cfd all couple/cfd couple_every ${Ncouple} mpi
fix cfd1 all couple/cfd/force transfer_force no transfer_torque no transfer_density yes transfer_type yes transfer_property name effvolfactor type scalar-atom
fix cfd2 all couple/cfd/recurrence transfer_force no transfer_fluctuations no
fix cfd3 all couple/cfd/convection T0 ${Tin} transfer_conduction yes max_change 1.0
fix cfd4 all couple/cfd/parttempfield T0 ${Tin}
####
# no trickling
fix cfd5 allore couple/cfd/deform rmin 0.01 monitor_heat yes use_latent_heat yes latent_heat 247000 verbose yes
fix cfd5 allore couple/cfd/deform rmin 0.01 monitor_heat yes use_latent_heat yes latent_heat 247000 verbose yes
# trickling
# fix cfd5 allore couple/cfd/deform rmin 0.01 monitor_heat yes use_latent_heat yes latent_heat 247000 delete_fully_deformed_particles no verbose yes
#fix cfd5 allore couple/cfd/deform rmin 0.01 monitor_heat yes use_latent_heat yes latent_heat 247000 delete_fully_deformed_particles no verbose yes
####
fix cfd6 all couple/cfd/chemistry n_species 4 species_names O2 CO CO2 N2 n_diff 3 diffusant_names CO O2 CO2 transfer_reactant yes init_layer_radii yes
fix cfd6 all couple/cfd/chemistry n_species 4 species_names O2 CO CO2 N2 n_diff 3 diffusant_names CO O2 CO2 transfer_reactant yes init_layer_radii yes
variable cokeBelowCZIndicator atom (f_partDeformations>0.99)*(type==1)
variable finecokeBelowCZIndicator atom (f_partDeformations>0.99)*(type==2)
@ -31,8 +31,8 @@ group finecokeBelowCZ dynamic all var finecokeBelowCZIndicator every ${Ncouple}
group aboveCZ dynamic all var aboveCZIndicator every ${Ncouple}
fix integr all nve/sphere/limit vlimit 0.05 omegalimit 1000.0
fix limV RW limit/vel 0.001
fix integr all nve/sphere/limit vlimit 0.05 omegalimit 1000.0
fix limV RW limit/vel 0.001
run 0
run 0
velocity all set 0.0 0.0 0.0

8
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/DEM/in.final
View File

@ -1,7 +1,7 @@
print 'latestTimeStep = ${currStep}' screen yes
print "variable latestTimeStep equal ${currStep}" file ../DEM/latestTimeStep
write_dump coke custom ${dumppath}/dump.relRadii1 id type x y z f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4]
write_dump finecoke custom ${dumppath}/dump.relRadii2 id type x y z f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4]
write_dump ore custom ${dumppath}/dump.relRadii3 id type x y z f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4]
write_dump fineore custom ${dumppath}/dump.relRadii4 id type x y z f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4]
write_dump coke custom ${dumppath}/dump.relRadii1 id type x y z f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4]
write_dump finecoke custom ${dumppath}/dump.relRadii2 id type x y z f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4]
write_dump ore custom ${dumppath}/dump.relRadii3 id type x y z f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4]
write_dump fineore custom ${dumppath}/dump.relRadii4 id type x y z f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4]

26
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/DEM/in.liggghts_dataDrivenCFD
View File

@ -3,23 +3,23 @@
include ../DEM/DEMparameters
###############################################
atom_style granular
atom_modify map array
communicate single vel yes
atom_style granular
atom_modify map array
communicate single vel yes
boundary f f f
newton off
units si
boundary f f f
newton off
units si
processors 2 2 8
#processors 1 1 8
#processors 1 1 8
# either (A) start new simulation and read dump or (B) restart simulation with restart file
# option A
region reg block -7.5 7.5 -7.5 7.5 -1.15 30 units box
create_box 5 reg
region reg block -7.5 7.5 -7.5 7.5 -1.15 30 units box
create_box 5 reg
fix fp all property/atom/polydispparcel effvolfactor ${poly_C} ${poly_fC} ${poly_O} ${poly_fO} 1.0
read_dump ../DEM/initFullState 0 x y z radius mass density bruteadd yes retaintimestep yes exactstep no
read_dump ../DEM/initFullState 0 x y z radius mass density bruteadd yes retaintimestep yes exactstep no
group coke type 1
group finecoke type 2
group ore type 3
@ -34,10 +34,10 @@ include ../DEM/latestTimeStep
reset_timestep ${latestTimeStep}
neighbor ${skin} bin
neigh_modify delay 0
neighbor ${skin} bin
neigh_modify delay 0
group fullyDeformed empty
neigh_modify exclude group fullyDeformed all
neigh_modify exclude group fullyDeformed all
timestep ${dt}

42
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/DEM/in.material_settings
View File

@ -2,40 +2,40 @@
# material properties required for new pair styles
##################################################
#fix m1 all property/global youngsModulus peratomtype ${yM_C} ${yM_C} ${yM_O} ${yM_O} ${yM_W}
#fix m2 all property/global poissonsRatio peratomtype ${pR_C} ${pR_C} ${pR_O} ${pR_O} ${pR_W}
#fix m3 all property/global coefficientRestitution peratomtypepair 5 ${rest_C} ${rest_C} ${rest_CO} ${rest_CO} ${rest_CW} ${rest_C} ${rest_C} ${rest_CO} ${rest_CO} ${rest_CW} ${rest_CO} ${rest_CO} ${rest_O} ${rest_O} ${rest_WO} ${rest_CO} ${rest_CO} ${rest_O} ${rest_O} ${rest_WO} ${rest_CW} ${rest_CW} ${rest_WO} ${rest_WO} ${rest_W}
#fix m4 all property/global coefficientFriction peratomtypepair 5 ${fric_C} ${fric_C} ${fric_CO} ${fric_CO} ${fric_CW} ${fric_C} ${fric_C} ${fric_CO} ${fric_CO} ${fric_CW} ${fric_CO} ${fric_CO} ${fric_O} ${fric_O} ${fric_WO} ${fric_CO} ${fric_CO} ${fric_O} ${fric_O} ${fric_WO} ${fric_CW} ${fric_CW} ${fric_WO} ${fric_WO} ${fric_W}
#fix m5 all property/global coefficientRollingFriction peratomtypepair 5 ${roll_C} ${roll_C} ${roll_CO} ${roll_CO} ${roll_CW} ${roll_C} ${roll_C} ${roll_CO} ${roll_CO} ${roll_CW} ${roll_CO} ${roll_CO} ${roll_O} ${roll_O} ${roll_WO} ${roll_CO} ${roll_CO} ${roll_O} ${roll_O} ${roll_WO} ${roll_CW} ${roll_CW} ${roll_WO} ${roll_WO} ${roll_W}
#fix m6 all property/global characteristicVelocity scalar ${cV}
#fix m1 all property/global youngsModulus peratomtype ${yM_C} ${yM_C} ${yM_O} ${yM_O} ${yM_W}
#fix m2 all property/global poissonsRatio peratomtype ${pR_C} ${pR_C} ${pR_O} ${pR_O} ${pR_W}
#fix m3 all property/global coefficientRestitution peratomtypepair 5 ${rest_C} ${rest_C} ${rest_CO} ${rest_CO} ${rest_CW} ${rest_C} ${rest_C} ${rest_CO} ${rest_CO} ${rest_CW} ${rest_CO} ${rest_CO} ${rest_O} ${rest_O} ${rest_WO} ${rest_CO} ${rest_CO} ${rest_O} ${rest_O} ${rest_WO} ${rest_CW} ${rest_CW} ${rest_WO} ${rest_WO} ${rest_W}
#fix m4 all property/global coefficientFriction peratomtypepair 5 ${fric_C} ${fric_C} ${fric_CO} ${fric_CO} ${fric_CW} ${fric_C} ${fric_C} ${fric_CO} ${fric_CO} ${fric_CW} ${fric_CO} ${fric_CO} ${fric_O} ${fric_O} ${fric_WO} ${fric_CO} ${fric_CO} ${fric_O} ${fric_O} ${fric_WO} ${fric_CW} ${fric_CW} ${fric_WO} ${fric_WO} ${fric_W}
#fix m5 all property/global coefficientRollingFriction peratomtypepair 5 ${roll_C} ${roll_C} ${roll_CO} ${roll_CO} ${roll_CW} ${roll_C} ${roll_C} ${roll_CO} ${roll_CO} ${roll_CW} ${roll_CO} ${roll_CO} ${roll_O} ${roll_O} ${roll_WO} ${roll_CO} ${roll_CO} ${roll_O} ${roll_O} ${roll_WO} ${roll_CW} ${roll_CW} ${roll_WO} ${roll_WO} ${roll_W}
#fix m6 all property/global characteristicVelocity scalar ${cV}
fix m1 all property/global kn peratomtypepair 5 ${kn} ${kn} ${kn} ${kn} ${knW} ${kn} ${kn} ${kn} ${kn} ${knW} ${kn} ${kn} ${kn} ${kn} ${knW} ${kn} ${kn} ${kn} ${kn} ${knW} ${knW} ${knW} ${knW} ${knW} ${knW}
fix m2 all property/global kt peratomtypepair 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
fix m3 all property/global gamman peratomtypepair 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
fix m4 all property/global gammat peratomtypepair 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
fix m5 all property/global coefficientFriction peratomtypepair 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
fix m6 all property/global coefficientRollingFriction peratomtypepair 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
fix m1 all property/global kn peratomtypepair 5 ${kn} ${kn} ${kn} ${kn} ${knW} ${kn} ${kn} ${kn} ${kn} ${knW} ${kn} ${kn} ${kn} ${kn} ${knW} ${kn} ${kn} ${kn} ${kn} ${knW} ${knW} ${knW} ${knW} ${knW} ${knW}
fix m2 all property/global kt peratomtypepair 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
fix m3 all property/global gamman peratomtypepair 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
fix m4 all property/global gammat peratomtypepair 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
fix m5 all property/global coefficientFriction peratomtypepair 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
fix m6 all property/global coefficientRollingFriction peratomtypepair 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
fix ftco all property/global thermalConductivity peratomtype 1.7 1.7 80 80 1.0
fix ftco all property/global thermalConductivity peratomtype 1.7 1.7 80 80 1.0
# constant Cp
#fix ftca all property/global thermalCapacity peratomtype 850 850 600 600 1.0
#fix ftca all property/global thermalCapacity peratomtype 850 850 600 600 1.0
# variable Cp
fix ftca all property/atom thermalCapacity scalar yes no no 600
fix ftca all property/atom thermalCapacity scalar yes no no 600
# do we need this? should be constructed by fix_cfd_convection; heat/gran results in heat/gran/cond with possibly undesired side effects
#fix heattransfer all heat/gran initial_temperature 293
#fix heattransfer all heat/gran initial_temperature 293
# pair style
#pair_style gran model hertz tangential history rolling_friction cdt
pair_style gran model hooke/stiffness
pair_coeff * *
fix bottom all mesh/surface file ../../CFDDEM/DEM/stl_files/surface_bottom.stl type 5
fix wall all mesh/surface file ../../CFDDEM/DEM/stl_files/surface_walls.stl type 5
fix bottom all mesh/surface file ../../CFDDEM/DEM/stl_files/surface_bottom.stl type 5
fix wall all mesh/surface file ../../CFDDEM/DEM/stl_files/surface_walls.stl type 5
#fix granwalls all wall/gran model hertz tangential history mesh n_meshes 2 meshes bottom wall
fix granwalls all wall/gran model hooke/stiffness mesh n_meshes 2 meshes bottom wall
#fix granwalls all wall/gran model hertz tangential history mesh n_meshes 2 meshes bottom wall
fix granwalls all wall/gran model hooke/stiffness mesh n_meshes 2 meshes bottom wall
#################################
# particle types and distribution
@ -49,6 +49,6 @@ group fineore type 4
####
# no need to set polydispersity if particle data is read from dump file
####
#fix fp all property/atom/polydispparcel effvolfactor ${poly_C} ${poly_fC} ${poly_O} ${poly_fO} 1.0
#fix fp all property/atom/polydispparcel effvolfactor ${poly_C} ${poly_fC} ${poly_O} ${poly_fO} 1.0
run 0

22
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/DEM/in.output
View File

@ -2,11 +2,11 @@
# output
########
#compute 1 all erotate/sphere
#thermo_style custom step atoms ke c_1 vol c_csumMT
#thermo 100
#thermo_modify lost ignore norm no
#compute_modify thermo_temp dynamic yes
#compute 1 all erotate/sphere
#thermo_style custom step atoms ke c_1 vol c_csumMT
#thermo 100
#thermo_modify lost ignore norm no
#compute_modify thermo_temp dynamic yes
# average particle temperature
@ -72,13 +72,13 @@ fix printRelRadiiCoke all print ${Nmonitor} "${time} ${rrC1} ${rrC2}" append ${d
# dump data
###########
dump dmp1 all custom/vtk ${Ndump} ${dumppath}/post/dump*.dataDrivenCFD.vtk id type radius x y z vx vy vz fx fy fz mass f_Temp[0] f_convectiveHeatFlux[0] f_fracRed[1] f_fracRed[2] f_fracRed[3] f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4] f_Aterm_chem0[1] f_Aterm_chem0[2] f_Aterm_chem0[3] f_Bterm_chem0[1] f_Bterm_chem0[2] f_Bterm_chem0[3] f_Massterm_chem0[0] f_Aterm_chem1[1] f_Aterm_chem1[2] f_Aterm_chem1[3] f_Bterm_chem1[1] f_Bterm_chem1[2] f_Bterm_chem1[3] f_Massterm_chem1[0] f_Aterm_chem2[0] f_Bterm_chem2[0] f_Massterm_chem2[0] f_Aterm_chem3[0] f_Bterm_chem3[0] f_Massterm_chem3[0] f_Aterm_chem4[0] f_Bterm_chem4[0] f_Massterm_chem4[0] f_Aterm_chem5[0] f_Bterm_chem5[0] f_Massterm_chem5[0] f_heatFlux[0] f_conductiveHeatFlux[0] v_reductionDegree f_reactionHeat[0]
dump dmp1 all custom/vtk ${Ndump} ${dumppath}/post/dump*.dataDrivenCFD.vtk id type radius x y z vx vy vz fx fy fz mass f_Temp[0] f_convectiveHeatFlux[0] f_fracRed[1] f_fracRed[2] f_fracRed[3] f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4] f_Aterm_chem0[1] f_Aterm_chem0[2] f_Aterm_chem0[3] f_Bterm_chem0[1] f_Bterm_chem0[2] f_Bterm_chem0[3] f_Massterm_chem0[0] f_Aterm_chem1[1] f_Aterm_chem1[2] f_Aterm_chem1[3] f_Bterm_chem1[1] f_Bterm_chem1[2] f_Bterm_chem1[3] f_Massterm_chem1[0] f_Aterm_chem2[0] f_Bterm_chem2[0] f_Massterm_chem2[0] f_Aterm_chem3[0] f_Bterm_chem3[0] f_Massterm_chem3[0] f_Aterm_chem4[0] f_Bterm_chem4[0] f_Massterm_chem4[0] f_Aterm_chem5[0] f_Bterm_chem5[0] f_Massterm_chem5[0] f_heatFlux[0] f_conductiveHeatFlux[0] v_reductionDegree f_reactionHeat[0]
dump dmp2 coke custom ${Nrestart} ${dumppath}/restart/dump*.relRadii1 id type x y z f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4]
dump dmp3 finecoke custom ${Nrestart} ${dumppath}/restart/dump*.relRadii2 id type x y z f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4]
dump dmp4 ore custom ${Nrestart} ${dumppath}/restart/dump*.relRadii3 id type x y z f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4]
dump dmp5 fineore custom ${Nrestart} ${dumppath}/restart/dump*.relRadii4 id type x y z f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4]
dump dmp6 all custom ${Nrestart} ${dumppath}/restart/initFullState* id type radius mass density x y z
dump dmp2 coke custom ${Nrestart} ${dumppath}/restart/dump*.relRadii1 id type x y z f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4]
dump dmp3 finecoke custom ${Nrestart} ${dumppath}/restart/dump*.relRadii2 id type x y z f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4]
dump dmp4 ore custom ${Nrestart} ${dumppath}/restart/dump*.relRadii3 id type x y z f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4]
dump dmp5 fineore custom ${Nrestart} ${dumppath}/restart/dump*.relRadii4 id type x y z f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4]
dump dmp6 all custom ${Nrestart} ${dumppath}/restart/initFullState* id type radius mass density x y z
compute 1 all erotate/sphere
thermo_style custom step atoms ke c_1 vol c_csumMT c_cm1 c_cm2 c_cm3 c_cm4 c_cTi c_cTgi c_cCOi c_cCO2i

2
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/DEM/in.testing
View File

@ -1 +1 @@
dump testdmp all custom/vtk ${Ndump} ${dumppath}/test/dump*.dataDrivenCFD.vtk x y z mass f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4] f_Aterm_chem1[1] f_Aterm_chem1[2] f_Aterm_chem1[3] f_Bterm_chem1[1] f_Bterm_chem1[2] f_Bterm_chem1[3]
dump testdmp all custom/vtk ${Ndump} ${dumppath}/test/dump*.dataDrivenCFD.vtk x y z mass f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4] f_Aterm_chem1[1] f_Aterm_chem1[2] f_Aterm_chem1[3] f_Bterm_chem1[1] f_Bterm_chem1[2] f_Bterm_chem1[3]

0
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/DEM/monitor/.gitkeep → tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/DEM/monitor/.gitignore vendored
View File

0
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/DEM/restart/.gitkeep → tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/dataDrivenCFD/DEM/restart/.gitignore vendored
View File

13
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/initLayerStructure/CFD/constant/RASProperties
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -15,11 +15,10 @@ FoamFile
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
RASModel laminar;//kEpsilon;
RASModel laminar;
turbulence off;
printCoeffs off;
// ************************************************************************* //

29
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/initLayerStructure/CFD/constant/couplingProperties
View File

@ -1,26 +1,17 @@
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.4 |
| \\ / A nd | Web: http://www.openfoam.org |
| \\/ M anipulation | |
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
root "";
case "";
instance "";
local "";
class dictionary;
object couplingProperties;
version 2.0;
format ascii;
class dictionary;
object couplingProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
#include "../../../CFDDEM/CFD/constant/couplingProperties"

11
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/initLayerStructure/CFD/constant/g
View File

@ -1,9 +1,9 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.6 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
@ -18,5 +18,4 @@ FoamFile
dimensions [0 1 -2 0 0 0 0];
value ( 0 0 -9.81 );
// ************************************************************************* //

41
tutorials/rcfdemSolverRhoSteadyPimpleChem/3D_hot_reacting/initLayerStructure/CFD/constant/liggghtsCommands
View File

@ -1,26 +1,17 @@
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.4 |
| \\ / A nd | Web: http://www.openfoam.org |
| \\/ M anipulation | |
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 6
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
root "";
case "";
instance "";
local "";
class dictionary;
object liggghtsCommands;
version 2.0;
format ascii;
class dictionary;
object liggghtsCommands;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
liggghtsCommandModels
@ -28,15 +19,5 @@ liggghtsCommandModels
runLiggghts
writeLiggghts
);
// ************************************************************************* //
/*runLiggghtsProps
{
preNo false;
}
writeLiggghtsProps
{
writeLast off;
overwrite off;
}*/

Some files were not shown because too many files have changed in this diff Show More