mirror of
https://github.com/ParticulateFlow/CFDEMcoupling-PFM.git
synced 2025-12-08 06:37:44 +00:00
Missing files for hot BF tutorial.
This commit is contained in:
tlichtenegger
@ -0,0 +1,3 @@
|
|||||||
|
rm -r postProcessing
|
||||||
|
rm -r proc*
|
||||||
|
rm -r dynamicCode
|
||||||
@ -0,0 +1,63 @@
|
|||||||
|
/*--------------------------------*- C++ -*----------------------------------*\
|
||||||
|
| ========= | |
|
||||||
|
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||||
|
| \\ / O peration | Version: 1.6 |
|
||||||
|
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||||
|
| \\/ M anipulation | |
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
FoamFile
|
||||||
|
{
|
||||||
|
version 2.0;
|
||||||
|
format ascii;
|
||||||
|
class volScalarField;
|
||||||
|
location "0";
|
||||||
|
object voidfraction;
|
||||||
|
}
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
dimensions [0 0 0 0 0 0 0];
|
||||||
|
|
||||||
|
internalField uniform 1;
|
||||||
|
|
||||||
|
boundaryField
|
||||||
|
{
|
||||||
|
bottom
|
||||||
|
{
|
||||||
|
type zeroGradient;
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
top
|
||||||
|
{
|
||||||
|
type zeroGradient;
|
||||||
|
|
||||||
|
}
|
||||||
|
|
||||||
|
wall
|
||||||
|
{
|
||||||
|
type fixedValue;
|
||||||
|
value uniform 10.0;
|
||||||
|
}
|
||||||
|
|
||||||
|
inner
|
||||||
|
{
|
||||||
|
type zeroGradient;
|
||||||
|
}
|
||||||
|
|
||||||
|
defaultFaces
|
||||||
|
{
|
||||||
|
type empty;
|
||||||
|
}
|
||||||
|
|
||||||
|
inletLeft
|
||||||
|
{
|
||||||
|
type zeroGradient;
|
||||||
|
}
|
||||||
|
inletRight
|
||||||
|
{
|
||||||
|
type zeroGradient;
|
||||||
|
}
|
||||||
|
|
||||||
|
}
|
||||||
|
|
||||||
|
// ************************************************************************* //
|
||||||
@ -0,0 +1,195 @@
|
|||||||
|
// Parametrized test case for a BF geometry
|
||||||
|
|
||||||
|
|
||||||
|
//Run using:
|
||||||
|
//m4 -P blockMeshDict.m4 > blockMeshDict
|
||||||
|
|
||||||
|
//m4 definitions:
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
//Mathematical constants:
|
||||||
|
|
||||||
|
|
||||||
|
//Geometry
|
||||||
|
|
||||||
|
// width of wedge
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
// height levels
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
// xlevels
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
//Grid points (integers!):
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
//m4_define(rGrading, 0.5)
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
| ========= | |
|
||||||
|
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||||
|
| \\ / O peration | Version: 1.4.1 |
|
||||||
|
| \\ / A nd | Web: http://www.openfoam.org |
|
||||||
|
| \\/ M anipulation | |
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
FoamFile
|
||||||
|
{
|
||||||
|
version 2.0;
|
||||||
|
format ascii;
|
||||||
|
|
||||||
|
root "";
|
||||||
|
case "";
|
||||||
|
instance "";
|
||||||
|
local "";
|
||||||
|
|
||||||
|
class dictionary;
|
||||||
|
object blockMeshDict;
|
||||||
|
}
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
convertToMeters 1;
|
||||||
|
|
||||||
|
vertices
|
||||||
|
(
|
||||||
|
(-6 0.6 -2) // Vertex V0 = 0
|
||||||
|
(0 0.6 -2) // Vertex V1 = 1
|
||||||
|
(0 0.6 0.75) // Vertex V2 = 2
|
||||||
|
(-6 0.6 0.75) // Vertex V3 = 3
|
||||||
|
|
||||||
|
(-7.4 0.6 3.9) // Vertex V4 = 4
|
||||||
|
(0 0.6 3.9) // Vertex V5 = 5
|
||||||
|
(0 0.6 6.4) // Vertex V6 = 6
|
||||||
|
(-7.4 0.6 6.4) // Vertex V7 = 7
|
||||||
|
|
||||||
|
(-4.7 0.6 24) // Vertex V8 = 8
|
||||||
|
(0 0.6 24) // Vertex V9 = 9
|
||||||
|
(0 0.6 26) // Vertex V10 = 10
|
||||||
|
(-4.7 0.6 26) // Vertex V11 = 11
|
||||||
|
|
||||||
|
// neg. y values
|
||||||
|
(-6 -0.6 -2) // Vertex V12 = 12
|
||||||
|
(0 -0.6 -2) // Vertex V13 = 13
|
||||||
|
(0 -0.6 0.75) // Vertex V14 = 14
|
||||||
|
(-6 -0.6 0.75) // Vertex V15 = 15
|
||||||
|
|
||||||
|
(-7.4 -0.6 3.9) // Vertex V16 = 16
|
||||||
|
(0 -0.6 3.9) // Vertex V17 = 17
|
||||||
|
(0 -0.6 6.4) // Vertex V18 = 18
|
||||||
|
(-7.4 -0.6 6.4) // Vertex V19 = 19
|
||||||
|
|
||||||
|
(-4.7 -0.6 24) // Vertex V20 = 20
|
||||||
|
(0 -0.6 24) // Vertex V21 = 21
|
||||||
|
(0 -0.6 26) // Vertex V22 = 22
|
||||||
|
(-4.7 -0.6 26) // Vertex V23 = 23
|
||||||
|
|
||||||
|
|
||||||
|
);
|
||||||
|
|
||||||
|
// Defining blocks:
|
||||||
|
blocks
|
||||||
|
(
|
||||||
|
hex ( 0 1 2 3 12 13 14 15 ) AB
|
||||||
|
(32 25 11)
|
||||||
|
simpleGrading (1.0 1 1)
|
||||||
|
|
||||||
|
hex ( 3 2 5 4 15 14 17 16) BC
|
||||||
|
(32 15 11)
|
||||||
|
simpleGrading (1.0 1 1)
|
||||||
|
|
||||||
|
hex ( 4 5 6 7 16 17 18 19) CD
|
||||||
|
(32 12 11)
|
||||||
|
simpleGrading (1.0 1 1)
|
||||||
|
|
||||||
|
hex ( 7 6 9 8 19 18 21 20 ) EF
|
||||||
|
(32 90 11)
|
||||||
|
simpleGrading (1.0 1 1)
|
||||||
|
|
||||||
|
hex ( 8 9 10 11 20 21 22 23 ) GH
|
||||||
|
(32 10 11)
|
||||||
|
simpleGrading (1.0 1 1)
|
||||||
|
);
|
||||||
|
|
||||||
|
// Defining patches:
|
||||||
|
boundary
|
||||||
|
(
|
||||||
|
bottom
|
||||||
|
{
|
||||||
|
type patch;
|
||||||
|
faces
|
||||||
|
(
|
||||||
|
(0 1 13 12)
|
||||||
|
);
|
||||||
|
}
|
||||||
|
wall
|
||||||
|
{
|
||||||
|
type wall;
|
||||||
|
faces
|
||||||
|
(
|
||||||
|
(0 3 15 12)
|
||||||
|
(3 4 16 15)
|
||||||
|
(4 7 19 16)
|
||||||
|
(7 8 20 19)
|
||||||
|
(8 11 23 20)
|
||||||
|
);
|
||||||
|
}
|
||||||
|
top
|
||||||
|
{
|
||||||
|
type patch;
|
||||||
|
faces
|
||||||
|
(
|
||||||
|
(11 10 22 23)
|
||||||
|
);
|
||||||
|
}
|
||||||
|
inner
|
||||||
|
{
|
||||||
|
type patch;
|
||||||
|
faces
|
||||||
|
(
|
||||||
|
(0 1 2 3)
|
||||||
|
(3 2 5 4)
|
||||||
|
(4 5 6 7)
|
||||||
|
(7 6 9 8)
|
||||||
|
(8 9 10 11)
|
||||||
|
|
||||||
|
(12 13 14 15)
|
||||||
|
(15 14 17 16)
|
||||||
|
(16 17 18 19)
|
||||||
|
(19 18 21 20)
|
||||||
|
(20 21 22 23)
|
||||||
|
);
|
||||||
|
}
|
||||||
|
);
|
||||||
|
|
||||||
|
mergePatchPairs
|
||||||
|
(
|
||||||
|
);
|
||||||
|
|
||||||
|
// ************************************************************************* //
|
||||||
@ -0,0 +1,195 @@
|
|||||||
|
// Parametrized test case for a BF geometry
|
||||||
|
|
||||||
|
|
||||||
|
//Run using:
|
||||||
|
//m4 -P blockMeshDict.m4 > blockMeshDict
|
||||||
|
|
||||||
|
//m4 definitions:
|
||||||
|
m4_changecom(//)m4_changequote([,])
|
||||||
|
m4_define(calc, [m4_esyscmd(perl -e 'use Math::Trig; printf ($1)')])
|
||||||
|
m4_define(VCOUNT, 0)
|
||||||
|
m4_define(vlabel, [[// ]Vertex $1 = VCOUNT m4_define($1, VCOUNT)m4_define([VCOUNT], m4_incr(VCOUNT))])
|
||||||
|
|
||||||
|
//Mathematical constants:
|
||||||
|
m4_define(pi, 3.1415926536)
|
||||||
|
|
||||||
|
//Geometry
|
||||||
|
|
||||||
|
// width of wedge
|
||||||
|
m4_define(y0, 0.6)
|
||||||
|
m4_define(y1, -0.6)
|
||||||
|
|
||||||
|
// height levels
|
||||||
|
m4_define(z0, -2)
|
||||||
|
m4_define(z1, 0.75)
|
||||||
|
m4_define(z2, 3.9)
|
||||||
|
m4_define(z3, 6.4)
|
||||||
|
m4_define(z4, 24)
|
||||||
|
m4_define(z5, 26)
|
||||||
|
|
||||||
|
// xlevels
|
||||||
|
m4_define(x0, -6)
|
||||||
|
m4_define(x1, -6)
|
||||||
|
m4_define(x2, -7.4)
|
||||||
|
m4_define(x3, -7.4)
|
||||||
|
m4_define(x4, -4.7)
|
||||||
|
m4_define(x5, -4.7)
|
||||||
|
|
||||||
|
|
||||||
|
//Grid points (integers!):
|
||||||
|
|
||||||
|
m4_define(xNumberOfCells, 32)
|
||||||
|
m4_define(yNumberOfCells, 11)
|
||||||
|
m4_define(zNumberOfCells0, 25)
|
||||||
|
m4_define(zNumberOfCells1, 15)
|
||||||
|
m4_define(zNumberOfCells2, 12)
|
||||||
|
m4_define(zNumberOfCells3, 90)
|
||||||
|
m4_define(zNumberOfCells4, 10)
|
||||||
|
m4_define(rGrading, 1.0)
|
||||||
|
//m4_define(rGrading, 0.5)
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
| ========= | |
|
||||||
|
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||||
|
| \\ / O peration | Version: 1.4.1 |
|
||||||
|
| \\ / A nd | Web: http://www.openfoam.org |
|
||||||
|
| \\/ M anipulation | |
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
FoamFile
|
||||||
|
{
|
||||||
|
version 2.0;
|
||||||
|
format ascii;
|
||||||
|
|
||||||
|
root "";
|
||||||
|
case "";
|
||||||
|
instance "";
|
||||||
|
local "";
|
||||||
|
|
||||||
|
class dictionary;
|
||||||
|
object blockMeshDict;
|
||||||
|
}
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
convertToMeters 1;
|
||||||
|
|
||||||
|
vertices
|
||||||
|
(
|
||||||
|
(x0 y0 z0) vlabel(V0)
|
||||||
|
(0 y0 z0) vlabel(V1)
|
||||||
|
(0 y0 z1) vlabel(V2)
|
||||||
|
(x1 y0 z1) vlabel(V3)
|
||||||
|
|
||||||
|
(x2 y0 z2) vlabel(V4)
|
||||||
|
(0 y0 z2) vlabel(V5)
|
||||||
|
(0 y0 z3) vlabel(V6)
|
||||||
|
(x3 y0 z3) vlabel(V7)
|
||||||
|
|
||||||
|
(x4 y0 z4) vlabel(V8)
|
||||||
|
(0 y0 z4) vlabel(V9)
|
||||||
|
(0 y0 z5) vlabel(V10)
|
||||||
|
(x5 y0 z5) vlabel(V11)
|
||||||
|
|
||||||
|
// neg. y values
|
||||||
|
(x0 y1 z0) vlabel(V12)
|
||||||
|
(0 y1 z0) vlabel(V13)
|
||||||
|
(0 y1 z1) vlabel(V14)
|
||||||
|
(x1 y1 z1) vlabel(V15)
|
||||||
|
|
||||||
|
(x2 y1 z2) vlabel(V16)
|
||||||
|
(0 y1 z2) vlabel(V17)
|
||||||
|
(0 y1 z3) vlabel(V18)
|
||||||
|
(x3 y1 z3) vlabel(V19)
|
||||||
|
|
||||||
|
(x4 y1 z4) vlabel(V20)
|
||||||
|
(0 y1 z4) vlabel(V21)
|
||||||
|
(0 y1 z5) vlabel(V22)
|
||||||
|
(x5 y1 z5) vlabel(V23)
|
||||||
|
|
||||||
|
|
||||||
|
);
|
||||||
|
|
||||||
|
// Defining blocks:
|
||||||
|
blocks
|
||||||
|
(
|
||||||
|
hex ( V0 V1 V2 V3 V12 V13 V14 V15 ) AB
|
||||||
|
(xNumberOfCells zNumberOfCells0 yNumberOfCells)
|
||||||
|
simpleGrading (rGrading 1 1)
|
||||||
|
|
||||||
|
hex ( V3 V2 V5 V4 V15 V14 V17 V16) BC
|
||||||
|
(xNumberOfCells zNumberOfCells1 yNumberOfCells)
|
||||||
|
simpleGrading (rGrading 1 1)
|
||||||
|
|
||||||
|
hex ( V4 V5 V6 V7 V16 V17 V18 V19) CD
|
||||||
|
(xNumberOfCells zNumberOfCells2 yNumberOfCells)
|
||||||
|
simpleGrading (rGrading 1 1)
|
||||||
|
|
||||||
|
hex ( V7 V6 V9 V8 V19 V18 V21 V20 ) EF
|
||||||
|
(xNumberOfCells zNumberOfCells3 yNumberOfCells)
|
||||||
|
simpleGrading (rGrading 1 1)
|
||||||
|
|
||||||
|
hex ( V8 V9 V10 V11 V20 V21 V22 V23 ) GH
|
||||||
|
(xNumberOfCells zNumberOfCells4 yNumberOfCells)
|
||||||
|
simpleGrading (rGrading 1 1)
|
||||||
|
);
|
||||||
|
|
||||||
|
// Defining patches:
|
||||||
|
boundary
|
||||||
|
(
|
||||||
|
bottom
|
||||||
|
{
|
||||||
|
type patch;
|
||||||
|
faces
|
||||||
|
(
|
||||||
|
(V0 V1 V13 V12)
|
||||||
|
);
|
||||||
|
}
|
||||||
|
wall
|
||||||
|
{
|
||||||
|
type wall;
|
||||||
|
faces
|
||||||
|
(
|
||||||
|
(V0 V3 V15 V12)
|
||||||
|
(V3 V4 V16 V15)
|
||||||
|
(V4 V7 V19 V16)
|
||||||
|
(V7 V8 V20 V19)
|
||||||
|
(V8 V11 V23 V20)
|
||||||
|
);
|
||||||
|
}
|
||||||
|
top
|
||||||
|
{
|
||||||
|
type patch;
|
||||||
|
faces
|
||||||
|
(
|
||||||
|
(V11 V10 V22 V23)
|
||||||
|
);
|
||||||
|
}
|
||||||
|
inner
|
||||||
|
{
|
||||||
|
type patch;
|
||||||
|
faces
|
||||||
|
(
|
||||||
|
(V0 V1 V2 V3)
|
||||||
|
(V3 V2 V5 V4)
|
||||||
|
(V4 V5 V6 V7)
|
||||||
|
(V7 V6 V9 V8)
|
||||||
|
(V8 V9 V10 V11)
|
||||||
|
|
||||||
|
(V12 V13 V14 V15)
|
||||||
|
(V15 V14 V17 V16)
|
||||||
|
(V16 V17 V18 V19)
|
||||||
|
(V19 V18 V21 V20)
|
||||||
|
(V20 V21 V22 V23)
|
||||||
|
);
|
||||||
|
}
|
||||||
|
);
|
||||||
|
|
||||||
|
mergePatchPairs
|
||||||
|
(
|
||||||
|
);
|
||||||
|
|
||||||
|
// ************************************************************************* //
|
||||||
@ -0,0 +1,8 @@
|
|||||||
|
#########
|
||||||
|
# removal
|
||||||
|
#########
|
||||||
|
|
||||||
|
# raceways
|
||||||
|
variable rCMass equal 200
|
||||||
|
fix delRW1 all remove nevery 1000 massrate ${rCMass} style delete seed 10007143 region RW1 verbose no restart_read no
|
||||||
|
fix delRW2 all remove nevery 1000 massrate ${rCMass} style delete seed 10007134 region RW2 verbose no restart_read no
|
||||||
@ -0,0 +1,22 @@
|
|||||||
|
##############
|
||||||
|
# cfd coupling
|
||||||
|
##############
|
||||||
|
|
||||||
|
variable Ncouple equal 10
|
||||||
|
variable Tin equal 293
|
||||||
|
|
||||||
|
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
|
||||||
|
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 integr all nve/sphere
|
||||||
|
fix limV all limit/vel 10.0
|
||||||
|
fix limV RW limit/vel 0.01
|
||||||
@ -0,0 +1,38 @@
|
|||||||
|
#################
|
||||||
|
# insertion setup
|
||||||
|
#################
|
||||||
|
|
||||||
|
variable outerRing equal 4.61
|
||||||
|
variable innerRing equal 4.1
|
||||||
|
region insRegTop1a block -${outerRing} -${innerRing} -0.59 0.59 26.125 26.95 units box
|
||||||
|
region insRegTop1b block -${innerRing} ${innerRing} -0.59 0.59 26.125 26.95 units box
|
||||||
|
region insRegTop1c block ${innerRing} ${outerRing} -0.59 0.59 26.125 26.95 units box
|
||||||
|
region insRegTop2 block -${outerRing} ${outerRing} -0.59 0.59 26.95 27.57 units box
|
||||||
|
|
||||||
|
region regTop union 4 insRegTop1a insRegTop1b insRegTop1c insRegTop2
|
||||||
|
variable nTop equal count(all,regTop)
|
||||||
|
|
||||||
|
group secHighestCoke empty
|
||||||
|
group secHighestOre empty
|
||||||
|
group secHighestFineMix empty
|
||||||
|
group highestFineMix empty
|
||||||
|
group highestCoke empty
|
||||||
|
group highestOre empty
|
||||||
|
|
||||||
|
# ensure reasonable initialization
|
||||||
|
group secHighestCoke union all secHighestCoke
|
||||||
|
group highestCoke union all highestCoke
|
||||||
|
group highestFineMix union all highestFineMix
|
||||||
|
variable zMaxSecHighCoke equal bound(secHighestCoke,zmax)
|
||||||
|
|
||||||
|
variable Nins equal 1000
|
||||||
|
variable counter loop 1000
|
||||||
|
variable insertedCoke equal mass(highestCoke)
|
||||||
|
variable insertedOre equal mass(highestOre)
|
||||||
|
variable insertedFineMix equal mass(highestFineMix)
|
||||||
|
fix eeIns all execute ${Nins} "../../CFDDEM/DEM/in.insert" file conditional nTop 0 0.1
|
||||||
|
|
||||||
|
variable time equal step*dt
|
||||||
|
variable np equal count(all)
|
||||||
|
fix printN all print 4000 "${time} ${np}" file ../DEM/monitor/Npart_CFDDEMflow.txt title "#time n_p"
|
||||||
|
fix printM all print 4000 "${time} ${insertedCoke} ${insertedOre} ${insertedFineMix} " file ../DEM/monitor/Mpart_CFDDEMflow.txt title "#time m_C m_O m_FM"
|
||||||
@ -0,0 +1,42 @@
|
|||||||
|
# material properties
|
||||||
|
include ../DEM/DEMparameters
|
||||||
|
|
||||||
|
|
||||||
|
###############################################
|
||||||
|
atom_style granular
|
||||||
|
atom_modify map array
|
||||||
|
communicate single vel yes
|
||||||
|
|
||||||
|
boundary f p f
|
||||||
|
newton off
|
||||||
|
units si
|
||||||
|
processors 1 1 8
|
||||||
|
|
||||||
|
read_restart ../DEM/liggghts.restart
|
||||||
|
reset_timestep 0
|
||||||
|
|
||||||
|
neighbor ${skin} bin
|
||||||
|
neigh_modify delay 0
|
||||||
|
|
||||||
|
timestep ${dt}
|
||||||
|
|
||||||
|
include ../DEM/in.material_settings
|
||||||
|
include ../DEM/in.regions
|
||||||
|
include ../DEM/in.RW
|
||||||
|
include ../DEM/in.coupling_integration
|
||||||
|
include ../DEM/in.output
|
||||||
|
|
||||||
|
run 1
|
||||||
|
|
||||||
|
include ../DEM/in.set_polydisp
|
||||||
|
include ../DEM/in.insertion_setup
|
||||||
|
|
||||||
|
# conditionally stop simulation
|
||||||
|
variable tmin equal 20.0
|
||||||
|
variable zmaxCFD equal 26.0
|
||||||
|
variable zminCFD equal 25.8
|
||||||
|
variable abort equal (v_time>=${tmin})&&(v_zMaxSecHighCoke<=${zmaxCFD})&&(v_zMaxSecHighCoke>=${zminCFD})
|
||||||
|
fix printAbort all print 4000 "${time} ${zMaxSecHighCoke} ${abort}" file ../DEM/monitor/Nabort.txt title "#time z_max_2 abort"
|
||||||
|
fix eeAbort all execute 1000 "../DEM/in.abort" file conditional abort 1 0.1
|
||||||
|
|
||||||
|
restart 200000 ../DEM/restart.A ../DEM/restart.B
|
||||||
@ -0,0 +1,48 @@
|
|||||||
|
##################################################
|
||||||
|
# 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 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
|
||||||
|
|
||||||
|
# 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 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_wall.stl type 5
|
||||||
|
|
||||||
|
fix granwalls all wall/gran model hertz tangential history mesh n_meshes 2 meshes bottom wall
|
||||||
|
|
||||||
|
|
||||||
|
#################################
|
||||||
|
# particle types and distribution
|
||||||
|
#################################
|
||||||
|
|
||||||
|
group coke type 1
|
||||||
|
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 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 fp all property/atom/polydispparcel effvolfactor ${poly_C} ${poly_fC} ${poly_O} ${poly_fO} 1.0
|
||||||
@ -0,0 +1,12 @@
|
|||||||
|
########
|
||||||
|
# output
|
||||||
|
########
|
||||||
|
|
||||||
|
#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
|
||||||
|
|
||||||
|
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]
|
||||||
@ -0,0 +1,30 @@
|
|||||||
|
#########
|
||||||
|
# regions
|
||||||
|
#########
|
||||||
|
|
||||||
|
# raceways
|
||||||
|
variable RWradFreeze equal 0.35
|
||||||
|
region RW1Freeze cylinder x 0 0 ${RWradFreeze} -6.0 -4.0
|
||||||
|
region RW2Freeze cylinder x 0 0 ${RWradFreeze} 4.0 6.0
|
||||||
|
region RW union 2 RW1Freeze RW2Freeze
|
||||||
|
group RW dynamic all region RW every 100
|
||||||
|
|
||||||
|
variable RWrad equal 0.4
|
||||||
|
region RW1 cylinder x 0 0 ${RWrad} -6.0 -4.0
|
||||||
|
region RW2 cylinder x 0 0 ${RWrad} 4.0 6.0
|
||||||
|
region RWs union 2 RW1 RW2
|
||||||
|
|
||||||
|
# deadman
|
||||||
|
variable DMheight equal 3
|
||||||
|
variable DMslope equal -0.8
|
||||||
|
region DMwall1 plane 0 0 ${DMheight} 1 0 ${DMslope}
|
||||||
|
region DMwall2 plane 0 0 ${DMheight} -1 0 ${DMslope}
|
||||||
|
region DM intersect 2 DMwall1 DMwall2
|
||||||
|
group DMcoke region DM
|
||||||
|
|
||||||
|
# hearth
|
||||||
|
variable DMheight1 equal ${DMheight}+1
|
||||||
|
region hearth block -7 7 -0.59 0.59 -2 ${DMheight1}
|
||||||
|
|
||||||
|
delete_atoms region RW1
|
||||||
|
delete_atoms region RW2
|
||||||
@ -0,0 +1,10 @@
|
|||||||
|
################
|
||||||
|
# polydispersity
|
||||||
|
################
|
||||||
|
|
||||||
|
set group coke property/atom/polydispparcel effvolfactor ${poly_C}
|
||||||
|
set group finecoke property/atom/polydispparcel effvolfactor ${poly_fC}
|
||||||
|
set group ore property/atom/polydispparcel effvolfactor ${poly_O}
|
||||||
|
set group fineore property/atom/polydispparcel effvolfactor ${poly_fO}
|
||||||
|
|
||||||
|
fix eeDM all execute ${NsetDM} "../../CFDDEM/DEM/in.DM" file
|
||||||
@ -0,0 +1 @@
|
|||||||
|
variable zInsLower equal 24.8821887469609
|
||||||
14
tutorials/rcfdemSolverRhoSteadyPimple/hotBlastFurnaceSlot/build_meshes.sh
Executable file
14
tutorials/rcfdemSolverRhoSteadyPimple/hotBlastFurnaceSlot/build_meshes.sh
Executable file
@ -0,0 +1,14 @@
|
|||||||
|
# initLayerStructure
|
||||||
|
cd initLayerStructure/CFD
|
||||||
|
bash mesh_slot.sh > mesh.log 2>&1
|
||||||
|
cd ../..
|
||||||
|
|
||||||
|
# CFDDEM
|
||||||
|
cd CFDDEM/CFD
|
||||||
|
bash mesh_slot.sh > mesh.log 2>&1
|
||||||
|
cd ../..
|
||||||
|
|
||||||
|
# dataDrivenCFD
|
||||||
|
cd dataDrivenCFD/CFD
|
||||||
|
bash mesh_slot.sh > mesh.log 2>&1
|
||||||
|
cd ../..
|
||||||
@ -0,0 +1,18 @@
|
|||||||
|
cd initLayerStructure/CFD
|
||||||
|
rm -r proc*
|
||||||
|
rm -r [0-9]*
|
||||||
|
rm -r postProcessing
|
||||||
|
rm -r dynamicCode
|
||||||
|
rm -r clockData
|
||||||
|
rm *.log
|
||||||
|
rm log.*
|
||||||
|
|
||||||
|
cd ../DEM
|
||||||
|
rm -r post*
|
||||||
|
rm -r monitor
|
||||||
|
rm init*
|
||||||
|
mkdir post
|
||||||
|
mkdir monitor
|
||||||
|
rm liggghts.restart*
|
||||||
|
rm restart.*
|
||||||
|
rm log.*
|
||||||
14
tutorials/rcfdemSolverRhoSteadyPimple/hotBlastFurnaceSlot/clean_meshes.sh
Executable file
14
tutorials/rcfdemSolverRhoSteadyPimple/hotBlastFurnaceSlot/clean_meshes.sh
Executable file
@ -0,0 +1,14 @@
|
|||||||
|
# initLayerStructure
|
||||||
|
cd initLayerStructure/CFD/constant
|
||||||
|
rm -r polyMesh
|
||||||
|
cd ../../..
|
||||||
|
|
||||||
|
# CFDDEM
|
||||||
|
cd CFDDEM/CFD/constant
|
||||||
|
rm -r polyMesh
|
||||||
|
cd ../../..
|
||||||
|
|
||||||
|
# dataDrivenCFD
|
||||||
|
cd dataDrivenCFD/CFD/constant
|
||||||
|
rm -r polyMesh
|
||||||
|
cd ../../..
|
||||||
@ -0,0 +1,3 @@
|
|||||||
|
rm -r postProcessing
|
||||||
|
rm -r proc*
|
||||||
|
rm -r dynamicCode
|
||||||
@ -0,0 +1,63 @@
|
|||||||
|
/*--------------------------------*- C++ -*----------------------------------*\
|
||||||
|
| ========= | |
|
||||||
|
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||||
|
| \\ / O peration | Version: 1.6 |
|
||||||
|
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||||
|
| \\/ M anipulation | |
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
FoamFile
|
||||||
|
{
|
||||||
|
version 2.0;
|
||||||
|
format ascii;
|
||||||
|
class volScalarField;
|
||||||
|
location "0";
|
||||||
|
object voidfraction;
|
||||||
|
}
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
dimensions [0 0 0 0 0 0 0];
|
||||||
|
|
||||||
|
internalField uniform 1;
|
||||||
|
|
||||||
|
boundaryField
|
||||||
|
{
|
||||||
|
bottom
|
||||||
|
{
|
||||||
|
type zeroGradient;
|
||||||
|
}
|
||||||
|
|
||||||
|
|
||||||
|
top
|
||||||
|
{
|
||||||
|
type zeroGradient;
|
||||||
|
|
||||||
|
}
|
||||||
|
|
||||||
|
wall
|
||||||
|
{
|
||||||
|
type fixedValue;
|
||||||
|
value uniform 10.0;
|
||||||
|
}
|
||||||
|
|
||||||
|
inner
|
||||||
|
{
|
||||||
|
type zeroGradient;
|
||||||
|
}
|
||||||
|
|
||||||
|
defaultFaces
|
||||||
|
{
|
||||||
|
type empty;
|
||||||
|
}
|
||||||
|
|
||||||
|
inletLeft
|
||||||
|
{
|
||||||
|
type zeroGradient;
|
||||||
|
}
|
||||||
|
inletRight
|
||||||
|
{
|
||||||
|
type zeroGradient;
|
||||||
|
}
|
||||||
|
|
||||||
|
}
|
||||||
|
|
||||||
|
// ************************************************************************* //
|
||||||
@ -0,0 +1,195 @@
|
|||||||
|
// Parametrized test case for a BF geometry
|
||||||
|
|
||||||
|
|
||||||
|
//Run using:
|
||||||
|
//m4 -P blockMeshDict.m4 > blockMeshDict
|
||||||
|
|
||||||
|
//m4 definitions:
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
//Mathematical constants:
|
||||||
|
|
||||||
|
|
||||||
|
//Geometry
|
||||||
|
|
||||||
|
// width of wedge
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
// height levels
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
// xlevels
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
//Grid points (integers!):
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
//m4_define(rGrading, 0.5)
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
| ========= | |
|
||||||
|
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||||
|
| \\ / O peration | Version: 1.4.1 |
|
||||||
|
| \\ / A nd | Web: http://www.openfoam.org |
|
||||||
|
| \\/ M anipulation | |
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
FoamFile
|
||||||
|
{
|
||||||
|
version 2.0;
|
||||||
|
format ascii;
|
||||||
|
|
||||||
|
root "";
|
||||||
|
case "";
|
||||||
|
instance "";
|
||||||
|
local "";
|
||||||
|
|
||||||
|
class dictionary;
|
||||||
|
object blockMeshDict;
|
||||||
|
}
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
convertToMeters 1;
|
||||||
|
|
||||||
|
vertices
|
||||||
|
(
|
||||||
|
(-6 0.6 -2) // Vertex V0 = 0
|
||||||
|
(0 0.6 -2) // Vertex V1 = 1
|
||||||
|
(0 0.6 0.75) // Vertex V2 = 2
|
||||||
|
(-6 0.6 0.75) // Vertex V3 = 3
|
||||||
|
|
||||||
|
(-7.4 0.6 3.9) // Vertex V4 = 4
|
||||||
|
(0 0.6 3.9) // Vertex V5 = 5
|
||||||
|
(0 0.6 6.4) // Vertex V6 = 6
|
||||||
|
(-7.4 0.6 6.4) // Vertex V7 = 7
|
||||||
|
|
||||||
|
(-4.7 0.6 24) // Vertex V8 = 8
|
||||||
|
(0 0.6 24) // Vertex V9 = 9
|
||||||
|
(0 0.6 26) // Vertex V10 = 10
|
||||||
|
(-4.7 0.6 26) // Vertex V11 = 11
|
||||||
|
|
||||||
|
// neg. y values
|
||||||
|
(-6 -0.6 -2) // Vertex V12 = 12
|
||||||
|
(0 -0.6 -2) // Vertex V13 = 13
|
||||||
|
(0 -0.6 0.75) // Vertex V14 = 14
|
||||||
|
(-6 -0.6 0.75) // Vertex V15 = 15
|
||||||
|
|
||||||
|
(-7.4 -0.6 3.9) // Vertex V16 = 16
|
||||||
|
(0 -0.6 3.9) // Vertex V17 = 17
|
||||||
|
(0 -0.6 6.4) // Vertex V18 = 18
|
||||||
|
(-7.4 -0.6 6.4) // Vertex V19 = 19
|
||||||
|
|
||||||
|
(-4.7 -0.6 24) // Vertex V20 = 20
|
||||||
|
(0 -0.6 24) // Vertex V21 = 21
|
||||||
|
(0 -0.6 26) // Vertex V22 = 22
|
||||||
|
(-4.7 -0.6 26) // Vertex V23 = 23
|
||||||
|
|
||||||
|
|
||||||
|
);
|
||||||
|
|
||||||
|
// Defining blocks:
|
||||||
|
blocks
|
||||||
|
(
|
||||||
|
hex ( 0 1 2 3 12 13 14 15 ) AB
|
||||||
|
(32 25 11)
|
||||||
|
simpleGrading (1.0 1 1)
|
||||||
|
|
||||||
|
hex ( 3 2 5 4 15 14 17 16) BC
|
||||||
|
(32 15 11)
|
||||||
|
simpleGrading (1.0 1 1)
|
||||||
|
|
||||||
|
hex ( 4 5 6 7 16 17 18 19) CD
|
||||||
|
(32 12 11)
|
||||||
|
simpleGrading (1.0 1 1)
|
||||||
|
|
||||||
|
hex ( 7 6 9 8 19 18 21 20 ) EF
|
||||||
|
(32 90 11)
|
||||||
|
simpleGrading (1.0 1 1)
|
||||||
|
|
||||||
|
hex ( 8 9 10 11 20 21 22 23 ) GH
|
||||||
|
(32 10 11)
|
||||||
|
simpleGrading (1.0 1 1)
|
||||||
|
);
|
||||||
|
|
||||||
|
// Defining patches:
|
||||||
|
boundary
|
||||||
|
(
|
||||||
|
bottom
|
||||||
|
{
|
||||||
|
type patch;
|
||||||
|
faces
|
||||||
|
(
|
||||||
|
(0 1 13 12)
|
||||||
|
);
|
||||||
|
}
|
||||||
|
wall
|
||||||
|
{
|
||||||
|
type wall;
|
||||||
|
faces
|
||||||
|
(
|
||||||
|
(0 3 15 12)
|
||||||
|
(3 4 16 15)
|
||||||
|
(4 7 19 16)
|
||||||
|
(7 8 20 19)
|
||||||
|
(8 11 23 20)
|
||||||
|
);
|
||||||
|
}
|
||||||
|
top
|
||||||
|
{
|
||||||
|
type patch;
|
||||||
|
faces
|
||||||
|
(
|
||||||
|
(11 10 22 23)
|
||||||
|
);
|
||||||
|
}
|
||||||
|
inner
|
||||||
|
{
|
||||||
|
type patch;
|
||||||
|
faces
|
||||||
|
(
|
||||||
|
(0 1 2 3)
|
||||||
|
(3 2 5 4)
|
||||||
|
(4 5 6 7)
|
||||||
|
(7 6 9 8)
|
||||||
|
(8 9 10 11)
|
||||||
|
|
||||||
|
(12 13 14 15)
|
||||||
|
(15 14 17 16)
|
||||||
|
(16 17 18 19)
|
||||||
|
(19 18 21 20)
|
||||||
|
(20 21 22 23)
|
||||||
|
);
|
||||||
|
}
|
||||||
|
);
|
||||||
|
|
||||||
|
mergePatchPairs
|
||||||
|
(
|
||||||
|
);
|
||||||
|
|
||||||
|
// ************************************************************************* //
|
||||||
@ -0,0 +1,27 @@
|
|||||||
|
##############
|
||||||
|
# cfd coupling
|
||||||
|
##############
|
||||||
|
|
||||||
|
variable Ncouple equal 1
|
||||||
|
variable Tin equal 293
|
||||||
|
|
||||||
|
region closeToWallLowerLeft block -6.0 -5.9 INF INF INF 0.75
|
||||||
|
region closeToWallLowerRight block 5.9 6.0 INF INF INF 0.75
|
||||||
|
region closeToWallUpperLeft prism -6.0 -5.9 INF INF 0.75 3.9 0.0 -1.4 0.0
|
||||||
|
region closeToWallUpperRight prism 5.9 6.0 INF INF 0.75 3.9 0.0 1.4 0.0
|
||||||
|
region closeToWallBottom block -6.0 6.0 INF INF -2.0 -1.9
|
||||||
|
region closeToWall union 5 closeToWallLowerLeft closeToWallLowerRight closeToWallUpperLeft closeToWallUpperRight closeToWallBottom
|
||||||
|
|
||||||
|
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 yes limit_fluctuations yes max_vfluc 0.05 relative_limit no region closeToWall
|
||||||
|
fix cfd3 all couple/cfd/convection T0 ${Tin} transfer_conduction yes
|
||||||
|
fix cfd4 all couple/cfd/parttempfield T0 ${Tin}
|
||||||
|
fix cfd5 all couple/cfd/deform rmin 0.01
|
||||||
|
|
||||||
|
fix integr all nve/sphere/limit vlimit 0.05 omegalimit 1000.0
|
||||||
|
fix limV RW limit/vel 0.01
|
||||||
|
|
||||||
|
|
||||||
|
run 0
|
||||||
|
velocity all set 0.0 0.0 0.0
|
||||||
@ -0,0 +1,9 @@
|
|||||||
|
#################
|
||||||
|
# insertion setup
|
||||||
|
#################
|
||||||
|
|
||||||
|
include ../../CFDDEM/DEM/zInsLower
|
||||||
|
region regTop block -4.61 4.61 -0.59 0.59 ${zInsLower} 26.1 units box
|
||||||
|
variable nTot equal count(all,regTop)
|
||||||
|
variable Nins equal 100
|
||||||
|
fix eeIns all execute ${Nins} "../DEM/in.insert" file conditional nTot 0 0.1
|
||||||
@ -0,0 +1,35 @@
|
|||||||
|
##################################################
|
||||||
|
# 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 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
|
||||||
|
|
||||||
|
# pair style
|
||||||
|
pair_style gran model hertz tangential history rolling_friction cdt
|
||||||
|
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_wall.stl type 5
|
||||||
|
|
||||||
|
fix granwalls all wall/gran model hertz tangential history mesh n_meshes 2 meshes bottom wall
|
||||||
|
|
||||||
|
|
||||||
|
#################################
|
||||||
|
# particle types and distribution
|
||||||
|
#################################
|
||||||
|
|
||||||
|
group coke type 1
|
||||||
|
group finecoke type 2
|
||||||
|
group ore type 3
|
||||||
|
group fineore type 4
|
||||||
|
|
||||||
|
fix fp all property/atom/polydispparcel effvolfactor ${poly_C} ${poly_fC} ${poly_O} ${poly_fO} 1.0
|
||||||
@ -0,0 +1,22 @@
|
|||||||
|
########
|
||||||
|
# output
|
||||||
|
########
|
||||||
|
|
||||||
|
compute 1 all erotate/sphere
|
||||||
|
thermo_style custom step atoms ke c_1 vol
|
||||||
|
thermo 100
|
||||||
|
thermo_modify lost ignore norm no
|
||||||
|
compute_modify thermo_temp dynamic yes
|
||||||
|
|
||||||
|
dump dmp1 all custom/vtk ${Ndump} ../DEM/post/dump*.dataDrivenCFD.vtk id type radius x y z vx vy vz fx fy fz mass f_Temp[0] f_convectiveHeatFlux[0] f_vfluc[1] f_vfluc[2] f_vfluc[3]
|
||||||
|
|
||||||
|
# calculate average particle temperature
|
||||||
|
variable time equal step*dt
|
||||||
|
variable currStep equal step
|
||||||
|
region total block INF INF INF INF INF INF units box
|
||||||
|
variable np equal count(all,total)
|
||||||
|
compute sumtemp all reduce/region total sum f_Temp[0]
|
||||||
|
variable sumT1 equal c_sumtemp
|
||||||
|
variable Tave equal v_sumT1/v_np
|
||||||
|
|
||||||
|
fix printheat all print 120 "${time} ${Tave}" append ../DEM/monitor/temp_ave.txt title "#time T_ave"
|
||||||
File diff suppressed because it is too large
Load Diff
@ -0,0 +1,222 @@
|
|||||||
|
#time N_coke_start N_coke_curr N_finecoke_start N_finecoke_curr
|
||||||
|
60 11735 11621 217 210
|
||||||
|
120 11735 11573 217 203
|
||||||
|
180 11735 11600 217 196
|
||||||
|
240 11735 11613 217 197
|
||||||
|
300 11735 11635 217 211
|
||||||
|
360 11735 11667 217 206
|
||||||
|
420 11735 11689 217 215
|
||||||
|
480 11735 11725 217 215
|
||||||
|
540 11735 11734 217 207
|
||||||
|
600 11735 11733 217 206
|
||||||
|
660 11735 11735 217 216
|
||||||
|
720 11735 11735 217 213
|
||||||
|
780 11735 11735 217 211
|
||||||
|
840 11735 11734 217 205
|
||||||
|
900 11735 11716 217 216
|
||||||
|
960 11735 11735 217 215
|
||||||
|
1020 11735 11731 217 217
|
||||||
|
1080 11735 11722 217 211
|
||||||
|
1140 11735 11735 217 215
|
||||||
|
1200 11735 11735 217 217
|
||||||
|
1260 11735 11734 217 216
|
||||||
|
1320 11735 11730 217 214
|
||||||
|
1380 11735 11735 217 212
|
||||||
|
1440 11735 11735 217 216
|
||||||
|
1500 11735 11735 217 209
|
||||||
|
1560 11735 11735 217 216
|
||||||
|
1620 11735 11734 217 216
|
||||||
|
1680 11735 11735 217 213
|
||||||
|
1740 11735 11729 217 204
|
||||||
|
1800 11735 11733 217 215
|
||||||
|
1860 11735 11733 217 212
|
||||||
|
1920 11735 11734 217 217
|
||||||
|
1980 11735 11735 217 213
|
||||||
|
2040 11735 11735 217 214
|
||||||
|
2100 11735 11735 217 217
|
||||||
|
2160 11735 11728 217 216
|
||||||
|
2220 11735 11734 217 217
|
||||||
|
2280 11735 11734 217 216
|
||||||
|
2340 11735 11735 217 210
|
||||||
|
2400 11735 11734 217 217
|
||||||
|
2460 11735 11735 217 216
|
||||||
|
2520 11735 11734 217 216
|
||||||
|
2580 11735 11726 217 214
|
||||||
|
2640 11735 11735 217 213
|
||||||
|
2700 11735 11720 217 216
|
||||||
|
2760 11735 11735 217 216
|
||||||
|
2820 11735 11723 217 214
|
||||||
|
2880 11735 11735 217 216
|
||||||
|
2940 11735 11727 217 212
|
||||||
|
3000 11735 11730 217 211
|
||||||
|
3060 11735 11734 217 216
|
||||||
|
3120 11735 11735 217 216
|
||||||
|
3180 11735 11734 217 213
|
||||||
|
3240 11735 11735 217 216
|
||||||
|
3300 11735 11735 217 211
|
||||||
|
3360 11735 11734 217 216
|
||||||
|
3420 11735 11735 217 216
|
||||||
|
3480 11735 11731 217 215
|
||||||
|
3540 11735 11735 217 213
|
||||||
|
3600 11735 11735 217 216
|
||||||
|
3660 11735 11735 217 216
|
||||||
|
3720 11735 11729 217 216
|
||||||
|
3780 11735 11735 217 215
|
||||||
|
3840 11735 11734 217 216
|
||||||
|
3900 11735 11735 217 216
|
||||||
|
3960 11735 11735 217 217
|
||||||
|
4020 11735 11735 217 216
|
||||||
|
4080 11735 11730 217 211
|
||||||
|
4140 11735 11732 217 212
|
||||||
|
4200 11735 11729 217 216
|
||||||
|
4260 11735 11734 217 215
|
||||||
|
4320 11735 11735 217 212
|
||||||
|
4380 11735 11728 217 213
|
||||||
|
4440 11735 11735 217 216
|
||||||
|
4500 11735 11733 217 217
|
||||||
|
4560 11735 11735 217 213
|
||||||
|
4620 11735 11722 217 216
|
||||||
|
4680 11735 11735 217 216
|
||||||
|
4740 11735 11735 217 216
|
||||||
|
4800 11735 11735 217 212
|
||||||
|
4860 11735 11735 217 214
|
||||||
|
4920 11735 11728 217 214
|
||||||
|
4980 11735 11735 217 212
|
||||||
|
5040 11735 11734 217 216
|
||||||
|
5100 11735 11734 217 216
|
||||||
|
5160 11735 11734 217 216
|
||||||
|
5220 11735 11726 217 217
|
||||||
|
5280 11735 11728 217 216
|
||||||
|
5340 11735 11735 217 216
|
||||||
|
5400 11735 11723 217 212
|
||||||
|
5460 11735 11734 217 214
|
||||||
|
5520 11735 11732 217 214
|
||||||
|
5580 11735 11735 217 216
|
||||||
|
5640 11735 11735 217 216
|
||||||
|
5700 11735 11735 217 216
|
||||||
|
5760 11735 11735 217 216
|
||||||
|
5820 11735 11735 217 216
|
||||||
|
5880 11735 11735 217 211
|
||||||
|
5940 11735 11735 217 212
|
||||||
|
6000 11735 11735 217 214
|
||||||
|
6060 11735 11734 217 213
|
||||||
|
6120 11735 11735 217 215
|
||||||
|
6180 11735 11734 217 211
|
||||||
|
6240 11735 11735 217 213
|
||||||
|
6300 11735 11735 217 213
|
||||||
|
6360 11735 11731 217 216
|
||||||
|
6420 11735 11734 217 213
|
||||||
|
6480 11735 11733 217 216
|
||||||
|
6540 11735 11735 217 214
|
||||||
|
6600 11735 11731 217 212
|
||||||
|
6660 11735 11735 217 216
|
||||||
|
6720 11735 11734 217 201
|
||||||
|
6780 11735 11735 217 215
|
||||||
|
6840 11735 11735 217 213
|
||||||
|
6900 11735 11724 217 208
|
||||||
|
6960 11735 11735 217 213
|
||||||
|
7020 11735 11735 217 214
|
||||||
|
7080 11735 11735 217 216
|
||||||
|
7140 11735 11735 217 213
|
||||||
|
7200 11735 11735 217 214
|
||||||
|
7260 11735 11735 217 216
|
||||||
|
7320 11735 11735 217 217
|
||||||
|
7380 11735 11733 217 211
|
||||||
|
7440 11735 11735 217 216
|
||||||
|
7500 11735 11735 217 217
|
||||||
|
7560 11735 11735 217 213
|
||||||
|
7620 11735 11735 217 217
|
||||||
|
7680 11735 11735 217 216
|
||||||
|
7740 11735 11735 217 217
|
||||||
|
7800 11735 11729 217 210
|
||||||
|
7860 11735 11734 217 211
|
||||||
|
7920 11735 11735 217 216
|
||||||
|
7980 11735 11727 217 214
|
||||||
|
8040 11735 11734 217 216
|
||||||
|
8100 11735 11729 217 214
|
||||||
|
8160 11735 11735 217 211
|
||||||
|
8220 11735 11729 217 216
|
||||||
|
8280 11735 11735 217 216
|
||||||
|
8340 11735 11734 217 217
|
||||||
|
8400 11735 11723 217 216
|
||||||
|
8460 11735 11735 217 216
|
||||||
|
8520 11735 11735 217 217
|
||||||
|
8580 11735 11721 217 217
|
||||||
|
8640 11735 11734 217 216
|
||||||
|
8700 11735 11735 217 213
|
||||||
|
8760 11735 11734 217 216
|
||||||
|
8820 11735 11729 217 217
|
||||||
|
8880 11735 11733 217 217
|
||||||
|
8940 11735 11734 217 213
|
||||||
|
9000 11735 11734 217 215
|
||||||
|
9060 11735 11735 217 216
|
||||||
|
9120 11735 11735 217 216
|
||||||
|
9180 11735 11731 217 210
|
||||||
|
9240 11735 11735 217 216
|
||||||
|
9300 11735 11730 217 216
|
||||||
|
9360 11735 11734 217 213
|
||||||
|
9420 11735 11735 217 206
|
||||||
|
9480 11735 11733 217 208
|
||||||
|
9540 11735 11727 217 213
|
||||||
|
9600 11735 11732 217 217
|
||||||
|
9660 11735 11735 217 214
|
||||||
|
9720 11735 11735 217 217
|
||||||
|
9780 11735 11730 217 213
|
||||||
|
9840 11735 11734 217 207
|
||||||
|
9900 11735 11732 217 214
|
||||||
|
9960 11735 11716 217 216
|
||||||
|
10020 11735 11735 217 216
|
||||||
|
10080 11735 11735 217 207
|
||||||
|
10140 11735 11735 217 216
|
||||||
|
10200 11735 11733 217 215
|
||||||
|
10260 11735 11735 217 215
|
||||||
|
10320 11735 11726 217 206
|
||||||
|
10380 11735 11735 217 211
|
||||||
|
10440 11735 11731 217 217
|
||||||
|
10500 11735 11734 217 216
|
||||||
|
10560 11735 11727 217 215
|
||||||
|
10620 11735 11730 217 214
|
||||||
|
10680 11735 11735 217 207
|
||||||
|
10740 11735 11734 217 207
|
||||||
|
10800 11735 11733 217 215
|
||||||
|
10860 11735 11732 217 215
|
||||||
|
10920 11735 11729 217 209
|
||||||
|
10980 11735 11733 217 213
|
||||||
|
11040 11735 11726 217 217
|
||||||
|
11100 11735 11734 217 215
|
||||||
|
11160 11735 11730 217 212
|
||||||
|
11220 11735 11735 217 212
|
||||||
|
11280 11735 11735 217 215
|
||||||
|
11340 11735 11731 217 216
|
||||||
|
11400 11735 11729 217 209
|
||||||
|
11460 11735 11733 217 210
|
||||||
|
11520 11735 11735 217 215
|
||||||
|
11580 11735 11735 217 212
|
||||||
|
11640 11735 11734 217 216
|
||||||
|
11700 11735 11729 217 216
|
||||||
|
11760 11735 11728 217 213
|
||||||
|
11820 11735 11735 217 215
|
||||||
|
11880 11735 11735 217 214
|
||||||
|
11940 11735 11721 217 217
|
||||||
|
12000 11735 11732 217 211
|
||||||
|
12060 11735 11725 217 214
|
||||||
|
12120 11735 11735 217 215
|
||||||
|
12180 11735 11734 217 208
|
||||||
|
12240 11735 11726 217 213
|
||||||
|
12300 11735 11728 217 213
|
||||||
|
12360 11735 11729 217 213
|
||||||
|
12420 11735 11735 217 215
|
||||||
|
12480 11735 11735 217 212
|
||||||
|
12540 11735 11728 217 212
|
||||||
|
12600 11735 11729 217 207
|
||||||
|
12660 11735 11732 217 216
|
||||||
|
12720 11735 11729 217 212
|
||||||
|
12780 11735 11727 217 214
|
||||||
|
12840 11735 11721 217 213
|
||||||
|
12900 11735 11731 217 214
|
||||||
|
12960 11735 11734 217 216
|
||||||
|
13020 11735 11735 217 210
|
||||||
|
13080 11735 11735 217 215
|
||||||
|
13140 11735 11730 217 211
|
||||||
|
13200 11735 11730 217 215
|
||||||
|
13260 11735 11723 217 216
|
||||||
@ -0,0 +1,222 @@
|
|||||||
|
#time n_p n_upper
|
||||||
|
60 118628 2689
|
||||||
|
120 118427 2229
|
||||||
|
180 118167 1766
|
||||||
|
240 117895 1310
|
||||||
|
300 117627 928
|
||||||
|
360 117376 640
|
||||||
|
420 117118 367
|
||||||
|
480 116871 96
|
||||||
|
540 116587 1
|
||||||
|
600 119824 2748
|
||||||
|
660 119535 2302
|
||||||
|
720 119211 1837
|
||||||
|
780 118859 1382
|
||||||
|
840 118523 973
|
||||||
|
900 118143 685
|
||||||
|
960 117802 414
|
||||||
|
1020 117475 141
|
||||||
|
1080 117129 1
|
||||||
|
1140 120344 2821
|
||||||
|
1200 120011 2368
|
||||||
|
1260 119622 1913
|
||||||
|
1320 119274 1445
|
||||||
|
1380 118881 1026
|
||||||
|
1440 118521 739
|
||||||
|
1500 118121 468
|
||||||
|
1560 117766 189
|
||||||
|
1620 117417 5
|
||||||
|
1680 120627 2891
|
||||||
|
1740 120279 2449
|
||||||
|
1800 119941 1992
|
||||||
|
1860 119590 1527
|
||||||
|
1920 119259 1087
|
||||||
|
1980 118879 794
|
||||||
|
2040 118512 517
|
||||||
|
2100 118152 241
|
||||||
|
2160 117788 15
|
||||||
|
2220 120945 2970
|
||||||
|
2280 120581 2535
|
||||||
|
2340 120212 2072
|
||||||
|
2400 119877 1606
|
||||||
|
2460 119539 1168
|
||||||
|
2520 119220 843
|
||||||
|
2580 118860 567
|
||||||
|
2640 118495 270
|
||||||
|
2700 118139 34
|
||||||
|
2760 121296 3050
|
||||||
|
2820 120931 2609
|
||||||
|
2880 120596 2138
|
||||||
|
2940 120218 1685
|
||||||
|
3000 119865 1238
|
||||||
|
3060 119529 890
|
||||||
|
3120 119171 604
|
||||||
|
3180 118776 322
|
||||||
|
3240 118429 58
|
||||||
|
3300 121597 3112
|
||||||
|
3360 121252 2688
|
||||||
|
3420 120867 2227
|
||||||
|
3480 120483 1766
|
||||||
|
3540 120119 1308
|
||||||
|
3600 119766 930
|
||||||
|
3660 119381 640
|
||||||
|
3720 119068 367
|
||||||
|
3780 118712 97
|
||||||
|
3840 118373 1
|
||||||
|
3900 121507 2748
|
||||||
|
3960 121142 2302
|
||||||
|
4020 120737 1837
|
||||||
|
4080 120390 1382
|
||||||
|
4140 120015 973
|
||||||
|
4200 119689 685
|
||||||
|
4260 119312 414
|
||||||
|
4320 118940 141
|
||||||
|
4380 118627 1
|
||||||
|
4440 121791 2821
|
||||||
|
4500 121442 2368
|
||||||
|
4560 121080 1913
|
||||||
|
4620 120689 1445
|
||||||
|
4680 120322 1026
|
||||||
|
4740 119966 739
|
||||||
|
4800 119625 468
|
||||||
|
4860 119256 189
|
||||||
|
4920 118890 5
|
||||||
|
4980 122036 2891
|
||||||
|
5040 121713 2449
|
||||||
|
5100 121338 1992
|
||||||
|
5160 120943 1527
|
||||||
|
5220 120567 1087
|
||||||
|
5280 120183 794
|
||||||
|
5340 119860 517
|
||||||
|
5400 119453 241
|
||||||
|
5460 119117 15
|
||||||
|
5520 122306 2970
|
||||||
|
5580 121942 2535
|
||||||
|
5640 121576 2072
|
||||||
|
5700 121209 1606
|
||||||
|
5760 120831 1168
|
||||||
|
5820 120438 843
|
||||||
|
5880 120083 567
|
||||||
|
5940 119730 270
|
||||||
|
6000 119347 34
|
||||||
|
6060 122531 3050
|
||||||
|
6120 122164 2609
|
||||||
|
6180 121835 2138
|
||||||
|
6240 121458 1685
|
||||||
|
6300 121071 1238
|
||||||
|
6360 120737 890
|
||||||
|
6420 120349 604
|
||||||
|
6480 119954 322
|
||||||
|
6540 119581 58
|
||||||
|
6600 122764 3112
|
||||||
|
6660 122392 2688
|
||||||
|
6720 122032 2227
|
||||||
|
6780 121701 1766
|
||||||
|
6840 121359 1308
|
||||||
|
6900 120962 930
|
||||||
|
6960 120614 640
|
||||||
|
7020 120209 367
|
||||||
|
7080 119845 97
|
||||||
|
7140 119473 1
|
||||||
|
7200 122651 2748
|
||||||
|
7260 122297 2302
|
||||||
|
7320 121959 1837
|
||||||
|
7380 121577 1382
|
||||||
|
7440 121223 973
|
||||||
|
7500 120850 685
|
||||||
|
7560 120448 414
|
||||||
|
7620 120067 141
|
||||||
|
7680 119696 1
|
||||||
|
7740 122872 2821
|
||||||
|
7800 122558 2368
|
||||||
|
7860 122212 1913
|
||||||
|
7920 121830 1445
|
||||||
|
7980 121467 1026
|
||||||
|
8040 121099 739
|
||||||
|
8100 120668 468
|
||||||
|
8160 120303 189
|
||||||
|
8220 119947 5
|
||||||
|
8280 123097 2891
|
||||||
|
8340 122728 2449
|
||||||
|
8400 122387 1992
|
||||||
|
8460 122063 1527
|
||||||
|
8520 121648 1087
|
||||||
|
8580 121302 794
|
||||||
|
8640 120909 517
|
||||||
|
8700 120513 241
|
||||||
|
8760 120160 15
|
||||||
|
8820 123305 2970
|
||||||
|
8880 122937 2535
|
||||||
|
8940 122590 2072
|
||||||
|
9000 122219 1606
|
||||||
|
9060 121852 1168
|
||||||
|
9120 121450 843
|
||||||
|
9180 121078 567
|
||||||
|
9240 120684 270
|
||||||
|
9300 120327 34
|
||||||
|
9360 123494 3050
|
||||||
|
9420 123137 2609
|
||||||
|
9480 122755 2138
|
||||||
|
9540 122378 1685
|
||||||
|
9600 122008 1238
|
||||||
|
9660 121645 890
|
||||||
|
9720 121258 604
|
||||||
|
9780 120849 322
|
||||||
|
9840 120471 58
|
||||||
|
9900 123713 3112
|
||||||
|
9960 123333 2688
|
||||||
|
10020 122941 2227
|
||||||
|
10080 122581 1766
|
||||||
|
10140 122204 1308
|
||||||
|
10200 121823 930
|
||||||
|
10260 121461 640
|
||||||
|
10320 121060 367
|
||||||
|
10380 120677 97
|
||||||
|
10440 120275 1
|
||||||
|
10500 123452 2748
|
||||||
|
10560 123089 2302
|
||||||
|
10620 122742 1837
|
||||||
|
10680 122355 1382
|
||||||
|
10740 121999 973
|
||||||
|
10800 121612 685
|
||||||
|
10860 121233 414
|
||||||
|
10920 120848 141
|
||||||
|
10980 120474 1
|
||||||
|
11040 123643 2821
|
||||||
|
11100 123227 2368
|
||||||
|
11160 122904 1913
|
||||||
|
11220 122493 1445
|
||||||
|
11280 122085 1026
|
||||||
|
11340 121742 739
|
||||||
|
11400 121335 468
|
||||||
|
11460 120940 189
|
||||||
|
11520 120581 5
|
||||||
|
11580 123779 2891
|
||||||
|
11640 123368 2449
|
||||||
|
11700 122970 1992
|
||||||
|
11760 122631 1527
|
||||||
|
11820 122275 1087
|
||||||
|
11880 121907 794
|
||||||
|
11940 121506 517
|
||||||
|
12000 121086 241
|
||||||
|
12060 120740 15
|
||||||
|
12120 123913 2970
|
||||||
|
12180 123516 2535
|
||||||
|
12240 123120 2072
|
||||||
|
12300 122775 1606
|
||||||
|
12360 122377 1168
|
||||||
|
12420 122040 843
|
||||||
|
12480 121687 567
|
||||||
|
12540 121261 270
|
||||||
|
12600 120851 34
|
||||||
|
12660 124054 3050
|
||||||
|
12720 123684 2609
|
||||||
|
12780 123274 2138
|
||||||
|
12840 122871 1685
|
||||||
|
12900 122504 1238
|
||||||
|
12960 122166 890
|
||||||
|
13020 121742 604
|
||||||
|
13080 121348 322
|
||||||
|
13140 120944 58
|
||||||
|
13200 124140 3112
|
||||||
|
13260 123807 2688
|
||||||
@ -0,0 +1,222 @@
|
|||||||
|
#time N_all N_coke N_ore
|
||||||
|
60 66241 18995 38689
|
||||||
|
120 65893 18864 38505
|
||||||
|
180 65545 18754 38289
|
||||||
|
240 65174 18670 38029
|
||||||
|
300 64796 18585 37774
|
||||||
|
360 64416 18520 37497
|
||||||
|
420 64022 18453 37206
|
||||||
|
480 63654 18386 36950
|
||||||
|
540 63283 18267 36746
|
||||||
|
600 66481 19126 38693
|
||||||
|
660 66147 18976 38553
|
||||||
|
720 65819 18839 38393
|
||||||
|
780 65444 18726 38160
|
||||||
|
840 65091 18642 37918
|
||||||
|
900 64718 18552 37658
|
||||||
|
960 64354 18474 37398
|
||||||
|
1020 63966 18401 37116
|
||||||
|
1080 63573 18331 36840
|
||||||
|
1140 66741 19241 38705
|
||||||
|
1200 66395 19104 38540
|
||||||
|
1260 66042 18946 38387
|
||||||
|
1320 65712 18796 38242
|
||||||
|
1380 65357 18657 38059
|
||||||
|
1440 64997 18559 37836
|
||||||
|
1500 64625 18483 37572
|
||||||
|
1560 64212 18413 37268
|
||||||
|
1620 63816 18364 36964
|
||||||
|
1680 66984 19304 38802
|
||||||
|
1740 66598 19198 38582
|
||||||
|
1800 66253 19071 38412
|
||||||
|
1860 65914 18904 38293
|
||||||
|
1920 65566 18725 38183
|
||||||
|
1980 65232 18597 38021
|
||||||
|
2040 64858 18490 37790
|
||||||
|
2100 64474 18418 37519
|
||||||
|
2160 64077 18366 37218
|
||||||
|
2220 67213 19304 39018
|
||||||
|
2280 66777 19260 38698
|
||||||
|
2340 66409 19183 38477
|
||||||
|
2400 66029 19046 38284
|
||||||
|
2460 65726 18885 38190
|
||||||
|
2520 65393 18677 38113
|
||||||
|
2580 65057 18523 37968
|
||||||
|
2640 64702 18424 37748
|
||||||
|
2700 64319 18341 37489
|
||||||
|
2760 67449 19286 39278
|
||||||
|
2820 67047 19251 38965
|
||||||
|
2880 66613 19213 38632
|
||||||
|
2940 66230 19138 38383
|
||||||
|
3000 65867 18993 38221
|
||||||
|
3060 65528 18835 38081
|
||||||
|
3120 65218 18635 38005
|
||||||
|
3180 64889 18484 37875
|
||||||
|
3240 64518 18366 37667
|
||||||
|
3300 67661 19278 39512
|
||||||
|
3360 67262 19222 39232
|
||||||
|
3420 66845 19183 38921
|
||||||
|
3480 66438 19148 38603
|
||||||
|
3540 66043 19074 38341
|
||||||
|
3600 65697 18930 38184
|
||||||
|
3660 65342 18747 38069
|
||||||
|
3720 65011 18564 37974
|
||||||
|
3780 64660 18415 37822
|
||||||
|
3840 64294 18326 37597
|
||||||
|
3900 67435 19251 39433
|
||||||
|
3960 67022 19193 39134
|
||||||
|
4020 66617 19152 38822
|
||||||
|
4080 66238 19109 38524
|
||||||
|
4140 65861 19012 38277
|
||||||
|
4200 65519 18837 38148
|
||||||
|
4260 65222 18659 38065
|
||||||
|
4320 64891 18457 37986
|
||||||
|
4380 64544 18331 37812
|
||||||
|
4440 67696 19256 39641
|
||||||
|
4500 67293 19214 39322
|
||||||
|
4560 66888 19176 39001
|
||||||
|
4620 66468 19139 38659
|
||||||
|
4680 66070 19108 38333
|
||||||
|
4740 65727 18948 38184
|
||||||
|
4800 65424 18755 38111
|
||||||
|
4860 65120 18550 38049
|
||||||
|
4920 64780 18364 37934
|
||||||
|
4980 67934 19303 39754
|
||||||
|
5040 67513 19258 39435
|
||||||
|
5100 67118 19224 39130
|
||||||
|
5160 66668 19185 38779
|
||||||
|
5220 66241 19131 38461
|
||||||
|
5280 65898 19047 38254
|
||||||
|
5340 65571 18840 38164
|
||||||
|
5400 65252 18639 38087
|
||||||
|
5460 64927 18452 37987
|
||||||
|
5520 68121 19340 39881
|
||||||
|
5580 67712 19277 39586
|
||||||
|
5640 67295 19228 39259
|
||||||
|
5700 66877 19192 38930
|
||||||
|
5760 66485 19160 38613
|
||||||
|
5820 66098 19070 38357
|
||||||
|
5880 65748 18916 38199
|
||||||
|
5940 65436 18727 38107
|
||||||
|
6000 65108 18539 38006
|
||||||
|
6060 68292 19411 39902
|
||||||
|
6120 67899 19336 39630
|
||||||
|
6180 67503 19271 39353
|
||||||
|
6240 67093 19223 39044
|
||||||
|
6300 66662 19184 38707
|
||||||
|
6360 66303 19106 38480
|
||||||
|
6420 65950 18957 38319
|
||||||
|
6480 65601 18777 38189
|
||||||
|
6540 65269 18594 38079
|
||||||
|
6600 68464 19474 39985
|
||||||
|
6660 68076 19382 39727
|
||||||
|
6720 67692 19317 39453
|
||||||
|
6780 67263 19260 39143
|
||||||
|
6840 66860 19214 38830
|
||||||
|
6900 66466 19161 38537
|
||||||
|
6960 66094 19004 38383
|
||||||
|
7020 65776 18834 38274
|
||||||
|
7080 65454 18649 38179
|
||||||
|
7140 65127 18504 38028
|
||||||
|
7200 68297 19419 39873
|
||||||
|
7260 67874 19362 39555
|
||||||
|
7320 67481 19305 39249
|
||||||
|
7380 67066 19263 38935
|
||||||
|
7440 66660 19224 38613
|
||||||
|
7500 66283 19108 38397
|
||||||
|
7560 65958 18948 38281
|
||||||
|
7620 65623 18750 38192
|
||||||
|
7680 65306 18557 38102
|
||||||
|
7740 68477 19448 39977
|
||||||
|
7800 68096 19374 39700
|
||||||
|
7860 67683 19324 39378
|
||||||
|
7920 67262 19274 39047
|
||||||
|
7980 66841 19236 38709
|
||||||
|
8040 66485 19156 38478
|
||||||
|
8100 66122 18988 38342
|
||||||
|
8160 65784 18797 38243
|
||||||
|
8220 65453 18600 38149
|
||||||
|
8280 68634 19474 40053
|
||||||
|
8340 68246 19395 39780
|
||||||
|
8400 67855 19340 39486
|
||||||
|
8460 67439 19289 39160
|
||||||
|
8520 67016 19244 38826
|
||||||
|
8580 66613 19186 38542
|
||||||
|
8640 66260 19038 38388
|
||||||
|
8700 65908 18864 38264
|
||||||
|
8760 65574 18664 38168
|
||||||
|
8820 68788 19518 40127
|
||||||
|
8880 68392 19427 39874
|
||||||
|
8940 67983 19359 39585
|
||||||
|
9000 67580 19309 39276
|
||||||
|
9060 67167 19263 38959
|
||||||
|
9120 66735 19213 38640
|
||||||
|
9180 66346 19088 38436
|
||||||
|
9240 66000 18927 38307
|
||||||
|
9300 65673 18732 38219
|
||||||
|
9360 68879 19566 40198
|
||||||
|
9420 68498 19452 39966
|
||||||
|
9480 68107 19372 39696
|
||||||
|
9540 67709 19322 39379
|
||||||
|
9600 67300 19276 39067
|
||||||
|
9660 66856 19227 38726
|
||||||
|
9720 66477 19119 38504
|
||||||
|
9780 66134 18960 38363
|
||||||
|
9840 65804 18782 38260
|
||||||
|
9900 69005 19592 40251
|
||||||
|
9960 68632 19476 40029
|
||||||
|
10020 68241 19396 39765
|
||||||
|
10080 67832 19341 39458
|
||||||
|
10140 67399 19289 39126
|
||||||
|
10200 66967 19245 38786
|
||||||
|
10260 66569 19141 38551
|
||||||
|
10320 66226 18998 38394
|
||||||
|
10380 65866 18804 38281
|
||||||
|
10440 65529 18618 38172
|
||||||
|
10500 68697 19491 40078
|
||||||
|
10560 68308 19412 39807
|
||||||
|
10620 67909 19347 39510
|
||||||
|
10680 67472 19293 39179
|
||||||
|
10740 67050 19253 38839
|
||||||
|
10800 66660 19169 38581
|
||||||
|
10860 66294 19016 38421
|
||||||
|
10920 65945 18839 38301
|
||||||
|
10980 65587 18635 38194
|
||||||
|
11040 68773 19507 40109
|
||||||
|
11100 68387 19417 39855
|
||||||
|
11160 68001 19358 39570
|
||||||
|
11220 67563 19302 39237
|
||||||
|
11280 67128 19255 38886
|
||||||
|
11340 66734 19188 38606
|
||||||
|
11400 66368 19039 38442
|
||||||
|
11460 66015 18859 38320
|
||||||
|
11520 65674 18666 38218
|
||||||
|
11580 68874 19529 40160
|
||||||
|
11640 68475 19431 39902
|
||||||
|
11700 68075 19361 39614
|
||||||
|
11760 67635 19305 39280
|
||||||
|
11820 67215 19259 38954
|
||||||
|
11880 66804 19201 38648
|
||||||
|
11940 66436 19065 38467
|
||||||
|
12000 66078 18894 38336
|
||||||
|
12060 65745 18702 38241
|
||||||
|
12120 68930 19546 40190
|
||||||
|
12180 68558 19455 39951
|
||||||
|
12240 68135 19372 39655
|
||||||
|
12300 67716 19319 39332
|
||||||
|
12360 67294 19266 39015
|
||||||
|
12420 66864 19216 38685
|
||||||
|
12480 66500 19093 38491
|
||||||
|
12540 66129 18927 38348
|
||||||
|
12600 65796 18736 38249
|
||||||
|
12660 68988 19564 40225
|
||||||
|
12720 68628 19470 39998
|
||||||
|
12780 68206 19377 39717
|
||||||
|
12840 67809 19328 39404
|
||||||
|
12900 67370 19272 39076
|
||||||
|
12960 66931 19227 38728
|
||||||
|
13020 66564 19121 38519
|
||||||
|
13080 66190 18965 38363
|
||||||
|
13140 65854 18770 38262
|
||||||
|
13200 69055 19588 40257
|
||||||
|
13260 68678 19478 40035
|
||||||
@ -0,0 +1,444 @@
|
|||||||
|
#time T_ave
|
||||||
|
30 889.305320789848
|
||||||
|
60 891.838353912913
|
||||||
|
90 894.473151867379
|
||||||
|
120 897.089123527059
|
||||||
|
150 899.748967166288
|
||||||
|
180 902.447780686994
|
||||||
|
210 905.086226925
|
||||||
|
240 907.775696114161
|
||||||
|
270 910.428483539357
|
||||||
|
300 913.064726630034
|
||||||
|
330 915.70513515713
|
||||||
|
360 918.246217774843
|
||||||
|
390 920.826079072607
|
||||||
|
420 923.416029814785
|
||||||
|
450 925.906925827297
|
||||||
|
480 928.371094729009
|
||||||
|
510 930.703583233878
|
||||||
|
540 932.653651992263
|
||||||
|
570 915.800388688515
|
||||||
|
600 917.675924657139
|
||||||
|
630 919.791527407083
|
||||||
|
660 921.678570246447
|
||||||
|
690 923.084918131532
|
||||||
|
720 924.987603284034
|
||||||
|
750 926.84042477452
|
||||||
|
780 928.61382875338
|
||||||
|
810 930.371335943727
|
||||||
|
840 932.323263238696
|
||||||
|
870 933.814357596226
|
||||||
|
900 935.512218921832
|
||||||
|
930 937.246752440795
|
||||||
|
960 938.988279051943
|
||||||
|
990 940.763227331367
|
||||||
|
1020 942.486625216752
|
||||||
|
1050 944.127183301109
|
||||||
|
1080 945.710111502318
|
||||||
|
1110 928.495994007816
|
||||||
|
1140 929.985412293268
|
||||||
|
1170 931.458449231017
|
||||||
|
1200 933.136981962803
|
||||||
|
1230 934.61974105503
|
||||||
|
1260 935.944069473608
|
||||||
|
1290 937.710673878627
|
||||||
|
1320 939.332976047093
|
||||||
|
1350 940.526419130537
|
||||||
|
1380 942.082660468898
|
||||||
|
1410 943.72009785569
|
||||||
|
1440 945.269307057576
|
||||||
|
1470 946.862154250492
|
||||||
|
1500 948.222207011971
|
||||||
|
1530 949.973399258398
|
||||||
|
1560 951.832681254123
|
||||||
|
1590 953.424803138778
|
||||||
|
1620 955.109280949352
|
||||||
|
1650 956.68293229609
|
||||||
|
1680 938.946499102166
|
||||||
|
1710 940.476117888105
|
||||||
|
1740 942.131112468968
|
||||||
|
1770 943.675927183837
|
||||||
|
1800 945.248592104243
|
||||||
|
1830 946.852045762838
|
||||||
|
1860 948.330640410628
|
||||||
|
1890 950.070097518372
|
||||||
|
1920 951.474877963525
|
||||||
|
1950 952.731295362448
|
||||||
|
1980 954.246147363487
|
||||||
|
2010 956.011522560744
|
||||||
|
2040 957.285214950857
|
||||||
|
2070 958.840724513129
|
||||||
|
2100 960.51545796794
|
||||||
|
2130 962.020827342256
|
||||||
|
2160 963.67778860563
|
||||||
|
2190 964.95148877482
|
||||||
|
2220 947.070653773666
|
||||||
|
2250 948.436288284839
|
||||||
|
2280 950.123992866554
|
||||||
|
2310 951.832825440757
|
||||||
|
2340 953.1068818729
|
||||||
|
2370 954.935454441704
|
||||||
|
2400 956.259231664907
|
||||||
|
2430 957.840841662049
|
||||||
|
2460 959.402621480203
|
||||||
|
2490 960.975352373107
|
||||||
|
2520 962.719581944607
|
||||||
|
2550 964.08318310568
|
||||||
|
2580 965.555026187372
|
||||||
|
2610 967.208209242457
|
||||||
|
2640 968.584872498595
|
||||||
|
2670 969.985278934755
|
||||||
|
2700 971.629008639746
|
||||||
|
2730 973.099320046426
|
||||||
|
2760 954.7299972099
|
||||||
|
2790 956.192552422768
|
||||||
|
2820 957.718616923188
|
||||||
|
2850 959.233026935006
|
||||||
|
2880 960.74979623456
|
||||||
|
2910 962.110693094275
|
||||||
|
2940 963.674895771964
|
||||||
|
2970 965.338105815732
|
||||||
|
3000 966.714203148459
|
||||||
|
3030 968.109890257642
|
||||||
|
3060 969.881506342268
|
||||||
|
3090 971.223253072856
|
||||||
|
3120 972.80191904726
|
||||||
|
3150 974.126627585077
|
||||||
|
3180 975.574697389321
|
||||||
|
3210 977.17870686291
|
||||||
|
3240 978.656912111966
|
||||||
|
3270 980.144262662635
|
||||||
|
3300 981.675117327957
|
||||||
|
3330 963.028743499226
|
||||||
|
3360 964.759278343104
|
||||||
|
3390 966.077822429413
|
||||||
|
3420 967.649293432252
|
||||||
|
3450 969.039095821758
|
||||||
|
3480 970.524736746371
|
||||||
|
3510 971.900017422494
|
||||||
|
3540 973.474071139187
|
||||||
|
3570 974.857143427599
|
||||||
|
3600 976.540405753628
|
||||||
|
3630 977.904472007536
|
||||||
|
3660 979.394936389015
|
||||||
|
3690 981.079102540988
|
||||||
|
3720 982.677331315408
|
||||||
|
3750 984.185327156242
|
||||||
|
3780 985.534593611564
|
||||||
|
3810 987.209383737842
|
||||||
|
3840 988.73363831873
|
||||||
|
3870 970.16686938216
|
||||||
|
3900 971.314573929955
|
||||||
|
3930 972.724302183748
|
||||||
|
3960 974.212922332668
|
||||||
|
3990 975.544515834565
|
||||||
|
4020 976.859118017112
|
||||||
|
4050 978.360648119758
|
||||||
|
4080 979.957339702583
|
||||||
|
4110 981.385356912583
|
||||||
|
4140 982.84002411716
|
||||||
|
4170 984.375610721969
|
||||||
|
4200 986.091807804044
|
||||||
|
4230 987.514516329627
|
||||||
|
4260 989.003684636441
|
||||||
|
4290 990.335487147411
|
||||||
|
4320 991.849159950387
|
||||||
|
4350 993.594633821583
|
||||||
|
4380 995.153380063535
|
||||||
|
4410 976.320610171536
|
||||||
|
4440 977.709556781979
|
||||||
|
4470 979.251299751233
|
||||||
|
4500 980.726650476227
|
||||||
|
4530 981.969825929014
|
||||||
|
4560 983.462323959349
|
||||||
|
4590 984.904261935389
|
||||||
|
4620 986.228387087647
|
||||||
|
4650 987.955887094135
|
||||||
|
4680 989.356593088031
|
||||||
|
4710 990.872536840402
|
||||||
|
4740 992.274082673832
|
||||||
|
4770 993.848167465842
|
||||||
|
4800 995.417122037751
|
||||||
|
4830 996.893310398074
|
||||||
|
4860 998.396402431458
|
||||||
|
4890 1000.04267149762
|
||||||
|
4920 1001.46768818773
|
||||||
|
4950 1002.93903178786
|
||||||
|
4980 983.797165239208
|
||||||
|
5010 985.59973906365
|
||||||
|
5040 987.003435276248
|
||||||
|
5070 988.441344472714
|
||||||
|
5100 989.900147336885
|
||||||
|
5130 991.132741576816
|
||||||
|
5160 992.609806124935
|
||||||
|
5190 994.136217815383
|
||||||
|
5220 995.62014546103
|
||||||
|
5250 997.009287841284
|
||||||
|
5280 998.50786957187
|
||||||
|
5310 1000.06303691281
|
||||||
|
5340 1001.85500930407
|
||||||
|
5370 1003.0761488111
|
||||||
|
5400 1004.47685551819
|
||||||
|
5430 1006.17138958271
|
||||||
|
5460 1007.8328875517
|
||||||
|
5490 1009.30979973903
|
||||||
|
5520 990.397132182681
|
||||||
|
5550 991.894930226518
|
||||||
|
5580 993.245526078119
|
||||||
|
5610 994.847021120056
|
||||||
|
5640 996.220862535871
|
||||||
|
5670 997.65653623132
|
||||||
|
5700 999.089783271254
|
||||||
|
5730 1000.34078224939
|
||||||
|
5760 1001.83028373188
|
||||||
|
5790 1003.39662822703
|
||||||
|
5820 1004.58212810037
|
||||||
|
5850 1006.32310135881
|
||||||
|
5880 1007.61828467997
|
||||||
|
5910 1009.10035206473
|
||||||
|
5940 1010.75158380796
|
||||||
|
5970 1012.29676670911
|
||||||
|
6000 1013.59187200513
|
||||||
|
6030 1015.20475047069
|
||||||
|
6060 996.085160695881
|
||||||
|
6090 997.474914404316
|
||||||
|
6120 999.042897193928
|
||||||
|
6150 1000.53336540039
|
||||||
|
6180 1002.20081168928
|
||||||
|
6210 1003.54558631296
|
||||||
|
6240 1005.05643643327
|
||||||
|
6270 1006.73904168663
|
||||||
|
6300 1007.94744986641
|
||||||
|
6330 1009.61598292334
|
||||||
|
6360 1011.02310183088
|
||||||
|
6390 1012.41493213131
|
||||||
|
6420 1013.86891257828
|
||||||
|
6450 1015.25102584135
|
||||||
|
6480 1016.73450743853
|
||||||
|
6510 1018.09536004893
|
||||||
|
6540 1019.67791726618
|
||||||
|
6570 1021.33512408337
|
||||||
|
6600 1022.91166876165
|
||||||
|
6630 1003.39673031239
|
||||||
|
6660 1004.98708167317
|
||||||
|
6690 1006.53784796054
|
||||||
|
6720 1008.02610158808
|
||||||
|
6750 1009.84430188863
|
||||||
|
6780 1011.29559846793
|
||||||
|
6810 1012.79730107325
|
||||||
|
6840 1014.35426080972
|
||||||
|
6870 1015.66620872405
|
||||||
|
6900 1017.05195913408
|
||||||
|
6930 1018.60931547925
|
||||||
|
6960 1020.16511531141
|
||||||
|
6990 1021.41352790985
|
||||||
|
7020 1022.88787311125
|
||||||
|
7050 1024.40961151166
|
||||||
|
7080 1026.00712476893
|
||||||
|
7110 1027.26060340406
|
||||||
|
7140 1028.92972135491
|
||||||
|
7170 1009.5450542415
|
||||||
|
7200 1011.05140572631
|
||||||
|
7230 1012.54055043374
|
||||||
|
7260 1013.93085600431
|
||||||
|
7290 1015.45386743186
|
||||||
|
7320 1017.08640321914
|
||||||
|
7350 1018.46575834128
|
||||||
|
7380 1019.93738096989
|
||||||
|
7410 1021.32278298765
|
||||||
|
7440 1022.84632330158
|
||||||
|
7470 1024.35678311012
|
||||||
|
7500 1025.74070015263
|
||||||
|
7530 1026.96530125747
|
||||||
|
7560 1028.48471572169
|
||||||
|
7590 1029.8477936561
|
||||||
|
7620 1031.36421327954
|
||||||
|
7650 1032.85903008868
|
||||||
|
7680 1034.29011501842
|
||||||
|
7710 1014.75694928292
|
||||||
|
7740 1016.27779758937
|
||||||
|
7770 1017.85391001844
|
||||||
|
7800 1019.54063748073
|
||||||
|
7830 1020.78831457094
|
||||||
|
7860 1022.48131995832
|
||||||
|
7890 1023.92550586627
|
||||||
|
7920 1025.268729479
|
||||||
|
7950 1026.69975504613
|
||||||
|
7980 1028.22272395673
|
||||||
|
8010 1029.66144304454
|
||||||
|
8040 1031.08238748538
|
||||||
|
8070 1032.32437634048
|
||||||
|
8100 1033.77036470161
|
||||||
|
8130 1035.28122474517
|
||||||
|
8160 1036.68224458502
|
||||||
|
8190 1038.24837831718
|
||||||
|
8220 1039.79306823109
|
||||||
|
8250 1041.23023407481
|
||||||
|
8280 1021.53416661941
|
||||||
|
8310 1023.00095966011
|
||||||
|
8340 1024.46433972313
|
||||||
|
8370 1026.07315442182
|
||||||
|
8400 1027.43919312426
|
||||||
|
8430 1028.93280691673
|
||||||
|
8460 1030.54237043718
|
||||||
|
8490 1031.88100174291
|
||||||
|
8520 1033.14967976421
|
||||||
|
8550 1034.68581022223
|
||||||
|
8580 1036.12198070573
|
||||||
|
8610 1037.60253633752
|
||||||
|
8640 1038.86602276056
|
||||||
|
8670 1040.20115855949
|
||||||
|
8700 1041.57164668452
|
||||||
|
8730 1043.08637966511
|
||||||
|
8760 1044.61599624875
|
||||||
|
8790 1046.05145732009
|
||||||
|
8820 1026.0683130667
|
||||||
|
8850 1027.57923788943
|
||||||
|
8880 1028.97317858794
|
||||||
|
8910 1030.47969902229
|
||||||
|
8940 1031.95920519325
|
||||||
|
8970 1033.50411808146
|
||||||
|
9000 1034.78928945389
|
||||||
|
9030 1036.26468681229
|
||||||
|
9060 1037.56287447789
|
||||||
|
9090 1039.01450899928
|
||||||
|
9120 1040.2381260526
|
||||||
|
9150 1041.83452248983
|
||||||
|
9180 1043.17653449451
|
||||||
|
9210 1044.69199123339
|
||||||
|
9240 1045.92154406004
|
||||||
|
9270 1047.36927622277
|
||||||
|
9300 1048.98256808961
|
||||||
|
9330 1050.35530109379
|
||||||
|
9360 1030.47542150221
|
||||||
|
9390 1031.88550608859
|
||||||
|
9420 1033.36385180094
|
||||||
|
9450 1034.60992047253
|
||||||
|
9480 1036.18737272441
|
||||||
|
9510 1037.66779762235
|
||||||
|
9540 1039.03602447527
|
||||||
|
9570 1040.4385033859
|
||||||
|
9600 1041.8711535022
|
||||||
|
9630 1043.31781710639
|
||||||
|
9660 1044.65960384285
|
||||||
|
9690 1046.16179646313
|
||||||
|
9720 1047.4841614725
|
||||||
|
9750 1048.89487348746
|
||||||
|
9780 1050.12991968562
|
||||||
|
9810 1051.66699293372
|
||||||
|
9840 1053.06362552673
|
||||||
|
9870 1054.76833174809
|
||||||
|
9900 1056.41795037063
|
||||||
|
9930 1036.25760557684
|
||||||
|
9960 1037.76342702634
|
||||||
|
9990 1039.03742527648
|
||||||
|
10020 1040.44650240585
|
||||||
|
10050 1041.96356803108
|
||||||
|
10080 1043.30732294072
|
||||||
|
10110 1044.80425873652
|
||||||
|
10140 1046.04964747654
|
||||||
|
10170 1047.58997341272
|
||||||
|
10200 1048.84928079743
|
||||||
|
10230 1050.3468646388
|
||||||
|
10260 1051.64989199182
|
||||||
|
10290 1053.01467504144
|
||||||
|
10320 1054.46777894015
|
||||||
|
10350 1055.89600321502
|
||||||
|
10380 1057.310322981
|
||||||
|
10410 1058.68981217675
|
||||||
|
10440 1060.08856242692
|
||||||
|
10470 1040.03868089938
|
||||||
|
10500 1041.52211984146
|
||||||
|
10530 1043.00004611211
|
||||||
|
10560 1044.38894913428
|
||||||
|
10590 1045.93310437418
|
||||||
|
10620 1047.39415241402
|
||||||
|
10650 1048.81324492874
|
||||||
|
10680 1050.16004776801
|
||||||
|
10710 1051.53401459213
|
||||||
|
10740 1053.02510513082
|
||||||
|
10770 1054.39162846035
|
||||||
|
10800 1055.81329442367
|
||||||
|
10830 1057.08841251989
|
||||||
|
10860 1058.63328443133
|
||||||
|
10890 1059.97660706592
|
||||||
|
10920 1061.50711154645
|
||||||
|
10950 1062.92454554407
|
||||||
|
10980 1064.35709029892
|
||||||
|
11010 1044.27981299294
|
||||||
|
11040 1045.54667927085
|
||||||
|
11070 1046.90237410343
|
||||||
|
11100 1048.2011083341
|
||||||
|
11130 1049.57902907841
|
||||||
|
11160 1051.19441205044
|
||||||
|
11190 1052.57993087606
|
||||||
|
11220 1053.83593396061
|
||||||
|
11250 1055.03608629448
|
||||||
|
11280 1056.43478103937
|
||||||
|
11310 1057.9253102247
|
||||||
|
11340 1059.41549100046
|
||||||
|
11370 1060.72299498613
|
||||||
|
11400 1062.15945572724
|
||||||
|
11430 1063.37041105969
|
||||||
|
11460 1064.85057778833
|
||||||
|
11490 1066.16741479758
|
||||||
|
11520 1067.76412763538
|
||||||
|
11550 1069.17480670993
|
||||||
|
11580 1048.99962275636
|
||||||
|
11610 1050.34395981322
|
||||||
|
11640 1051.6042498242
|
||||||
|
11670 1053.09835763596
|
||||||
|
11700 1054.39474750609
|
||||||
|
11730 1055.91483736821
|
||||||
|
11760 1057.32928763995
|
||||||
|
11790 1058.73473963456
|
||||||
|
11820 1060.32307527079
|
||||||
|
11850 1061.70718760394
|
||||||
|
11880 1063.13233513297
|
||||||
|
11910 1064.62365858764
|
||||||
|
11940 1065.8550375813
|
||||||
|
11970 1067.24230717597
|
||||||
|
12000 1068.55829727862
|
||||||
|
12030 1069.9836891129
|
||||||
|
12060 1071.51002650943
|
||||||
|
12090 1073.05900611798
|
||||||
|
12120 1052.51497046383
|
||||||
|
12150 1053.90538435209
|
||||||
|
12180 1055.21838344695
|
||||||
|
12210 1056.61278678858
|
||||||
|
12240 1057.96466119731
|
||||||
|
12270 1059.49573774374
|
||||||
|
12300 1060.96033652581
|
||||||
|
12330 1062.23276503729
|
||||||
|
12360 1063.61218557585
|
||||||
|
12390 1065.13187288275
|
||||||
|
12420 1066.61859879603
|
||||||
|
12450 1068.16855879352
|
||||||
|
12480 1069.56933995846
|
||||||
|
12510 1070.88990214989
|
||||||
|
12540 1072.11090055123
|
||||||
|
12570 1073.47849755105
|
||||||
|
12600 1074.91258307699
|
||||||
|
12630 1076.29281673476
|
||||||
|
12660 1056.05613544934
|
||||||
|
12690 1057.37785141019
|
||||||
|
12720 1058.92136696333
|
||||||
|
12750 1060.32847485315
|
||||||
|
12780 1061.48967401711
|
||||||
|
12810 1062.97791020239
|
||||||
|
12840 1064.16522967727
|
||||||
|
12870 1065.70478471559
|
||||||
|
12900 1067.0602525252
|
||||||
|
12930 1068.38166906892
|
||||||
|
12960 1070.04438970224
|
||||||
|
12990 1071.26842203683
|
||||||
|
13020 1072.58239836781
|
||||||
|
13050 1074.11004382379
|
||||||
|
13080 1075.33329850115
|
||||||
|
13110 1076.69894148889
|
||||||
|
13140 1078.16114285743
|
||||||
|
13170 1079.62413417946
|
||||||
|
13200 1081.15450725428
|
||||||
|
13230 1060.57258831942
|
||||||
|
13260 1062.20530825209
|
||||||
|
13290 1063.62060397913
|
||||||
@ -0,0 +1,3 @@
|
|||||||
|
rm -r postProcessing
|
||||||
|
rm -r proc*
|
||||||
|
rm -r dynamicCode
|
||||||
@ -0,0 +1,195 @@
|
|||||||
|
// Parametrized test case for a BF geometry
|
||||||
|
|
||||||
|
|
||||||
|
//Run using:
|
||||||
|
//m4 -P blockMeshDict.m4 > blockMeshDict
|
||||||
|
|
||||||
|
//m4 definitions:
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
//Mathematical constants:
|
||||||
|
|
||||||
|
|
||||||
|
//Geometry
|
||||||
|
|
||||||
|
// width of wedge
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
// height levels
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
// xlevels
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
//Grid points (integers!):
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
//m4_define(rGrading, 0.5)
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
|
||||||
|
/*---------------------------------------------------------------------------*\
|
||||||
|
| ========= | |
|
||||||
|
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||||
|
| \\ / O peration | Version: 1.4.1 |
|
||||||
|
| \\ / A nd | Web: http://www.openfoam.org |
|
||||||
|
| \\/ M anipulation | |
|
||||||
|
\*---------------------------------------------------------------------------*/
|
||||||
|
|
||||||
|
FoamFile
|
||||||
|
{
|
||||||
|
version 2.0;
|
||||||
|
format ascii;
|
||||||
|
|
||||||
|
root "";
|
||||||
|
case "";
|
||||||
|
instance "";
|
||||||
|
local "";
|
||||||
|
|
||||||
|
class dictionary;
|
||||||
|
object blockMeshDict;
|
||||||
|
}
|
||||||
|
|
||||||
|
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||||
|
|
||||||
|
convertToMeters 1;
|
||||||
|
|
||||||
|
vertices
|
||||||
|
(
|
||||||
|
(-6 0.6 -2) // Vertex V0 = 0
|
||||||
|
(0 0.6 -2) // Vertex V1 = 1
|
||||||
|
(0 0.6 0.75) // Vertex V2 = 2
|
||||||
|
(-6 0.6 0.75) // Vertex V3 = 3
|
||||||
|
|
||||||
|
(-7.4 0.6 3.9) // Vertex V4 = 4
|
||||||
|
(0 0.6 3.9) // Vertex V5 = 5
|
||||||
|
(0 0.6 6.4) // Vertex V6 = 6
|
||||||
|
(-7.4 0.6 6.4) // Vertex V7 = 7
|
||||||
|
|
||||||
|
(-4.7 0.6 24) // Vertex V8 = 8
|
||||||
|
(0 0.6 24) // Vertex V9 = 9
|
||||||
|
(0 0.6 26) // Vertex V10 = 10
|
||||||
|
(-4.7 0.6 26) // Vertex V11 = 11
|
||||||
|
|
||||||
|
// neg. y values
|
||||||
|
(-6 -0.6 -2) // Vertex V12 = 12
|
||||||
|
(0 -0.6 -2) // Vertex V13 = 13
|
||||||
|
(0 -0.6 0.75) // Vertex V14 = 14
|
||||||
|
(-6 -0.6 0.75) // Vertex V15 = 15
|
||||||
|
|
||||||
|
(-7.4 -0.6 3.9) // Vertex V16 = 16
|
||||||
|
(0 -0.6 3.9) // Vertex V17 = 17
|
||||||
|
(0 -0.6 6.4) // Vertex V18 = 18
|
||||||
|
(-7.4 -0.6 6.4) // Vertex V19 = 19
|
||||||
|
|
||||||
|
(-4.7 -0.6 24) // Vertex V20 = 20
|
||||||
|
(0 -0.6 24) // Vertex V21 = 21
|
||||||
|
(0 -0.6 26) // Vertex V22 = 22
|
||||||
|
(-4.7 -0.6 26) // Vertex V23 = 23
|
||||||
|
|
||||||
|
|
||||||
|
);
|
||||||
|
|
||||||
|
// Defining blocks:
|
||||||
|
blocks
|
||||||
|
(
|
||||||
|
hex ( 0 1 2 3 12 13 14 15 ) AB
|
||||||
|
(32 25 11)
|
||||||
|
simpleGrading (1.0 1 1)
|
||||||
|
|
||||||
|
hex ( 3 2 5 4 15 14 17 16) BC
|
||||||
|
(32 15 11)
|
||||||
|
simpleGrading (1.0 1 1)
|
||||||
|
|
||||||
|
hex ( 4 5 6 7 16 17 18 19) CD
|
||||||
|
(32 12 11)
|
||||||
|
simpleGrading (1.0 1 1)
|
||||||
|
|
||||||
|
hex ( 7 6 9 8 19 18 21 20 ) EF
|
||||||
|
(32 90 11)
|
||||||
|
simpleGrading (1.0 1 1)
|
||||||
|
|
||||||
|
hex ( 8 9 10 11 20 21 22 23 ) GH
|
||||||
|
(32 10 11)
|
||||||
|
simpleGrading (1.0 1 1)
|
||||||
|
);
|
||||||
|
|
||||||
|
// Defining patches:
|
||||||
|
boundary
|
||||||
|
(
|
||||||
|
bottom
|
||||||
|
{
|
||||||
|
type patch;
|
||||||
|
faces
|
||||||
|
(
|
||||||
|
(0 1 13 12)
|
||||||
|
);
|
||||||
|
}
|
||||||
|
wall
|
||||||
|
{
|
||||||
|
type wall;
|
||||||
|
faces
|
||||||
|
(
|
||||||
|
(0 3 15 12)
|
||||||
|
(3 4 16 15)
|
||||||
|
(4 7 19 16)
|
||||||
|
(7 8 20 19)
|
||||||
|
(8 11 23 20)
|
||||||
|
);
|
||||||
|
}
|
||||||
|
top
|
||||||
|
{
|
||||||
|
type patch;
|
||||||
|
faces
|
||||||
|
(
|
||||||
|
(11 10 22 23)
|
||||||
|
);
|
||||||
|
}
|
||||||
|
inner
|
||||||
|
{
|
||||||
|
type patch;
|
||||||
|
faces
|
||||||
|
(
|
||||||
|
(0 1 2 3)
|
||||||
|
(3 2 5 4)
|
||||||
|
(4 5 6 7)
|
||||||
|
(7 6 9 8)
|
||||||
|
(8 9 10 11)
|
||||||
|
|
||||||
|
(12 13 14 15)
|
||||||
|
(15 14 17 16)
|
||||||
|
(16 17 18 19)
|
||||||
|
(19 18 21 20)
|
||||||
|
(20 21 22 23)
|
||||||
|
);
|
||||||
|
}
|
||||||
|
);
|
||||||
|
|
||||||
|
mergePatchPairs
|
||||||
|
(
|
||||||
|
);
|
||||||
|
|
||||||
|
// ************************************************************************* //
|
||||||
@ -0,0 +1,2 @@
|
|||||||
|
# no end-of-run tasks
|
||||||
|
print 'no end-of-run tasks' screen yes
|
||||||
@ -0,0 +1,34 @@
|
|||||||
|
# material properties
|
||||||
|
include ../../CFDDEM/DEM/DEMparameters
|
||||||
|
|
||||||
|
|
||||||
|
###############################################
|
||||||
|
atom_style granular
|
||||||
|
atom_modify map array
|
||||||
|
communicate single vel yes
|
||||||
|
|
||||||
|
boundary f p f
|
||||||
|
newton off
|
||||||
|
units si
|
||||||
|
processors 1 1 8
|
||||||
|
|
||||||
|
read_restart ../DEM/liggghts.restart_fill
|
||||||
|
reset_timestep 0
|
||||||
|
|
||||||
|
neighbor ${skin} bin
|
||||||
|
neigh_modify delay 0
|
||||||
|
|
||||||
|
timestep ${dt}
|
||||||
|
|
||||||
|
include ../../CFDDEM/DEM/in.material_settings
|
||||||
|
include ../../CFDDEM/DEM/in.regions
|
||||||
|
include ../../CFDDEM/DEM/in.RW
|
||||||
|
include ../../CFDDEM/DEM/in.coupling_integration
|
||||||
|
include ../../CFDDEM/DEM/in.output
|
||||||
|
|
||||||
|
run 1
|
||||||
|
|
||||||
|
include ../../CFDDEM/DEM/in.set_polydisp
|
||||||
|
include ../../CFDDEM/DEM/in.insertion_setup
|
||||||
|
|
||||||
|
restart 200000 ../DEM/restart.A ../DEM/restart.B
|
||||||
Reference in New Issue
Block a user