mirror of
https://github.com/ParticulateFlow/CFDEMcoupling-PFM.git
synced 2025-12-08 06:37:44 +00:00
Reduction test case.
This commit is contained in:
Thomas Lichtenegger
@ -2,8 +2,8 @@
|
||||
|
||||
#===================================================================#
|
||||
# CFD-DEM run script for FinesColumn testcase
|
||||
# init FinesColumn
|
||||
# Thomas Lichtenegger - January 2017
|
||||
# run FinesColumn
|
||||
# Thomas Lichtenegger - June 2023
|
||||
#===================================================================#
|
||||
|
||||
#- source CFDEM env vars
|
||||
|
||||
@ -3,7 +3,7 @@
|
||||
#===================================================================#
|
||||
# DEM run script for FinesColumn testcase
|
||||
# init FinesColumn
|
||||
# Thomas Lichtenegger - January 2017
|
||||
# Thomas Lichtenegger - June 2023
|
||||
#===================================================================#
|
||||
|
||||
#- source CFDEM env vars
|
||||
|
||||
@ -2,8 +2,8 @@
|
||||
|
||||
#===================================================================#
|
||||
# CFD-DEM run script for FinesColumn testcase
|
||||
# init FinesColumn
|
||||
# Thomas Lichtenegger - January 2017
|
||||
# run FinesColumn
|
||||
# Thomas Lichtenegger - June 2023
|
||||
#===================================================================#
|
||||
|
||||
#- source CFDEM env vars
|
||||
|
||||
@ -3,7 +3,7 @@
|
||||
#===================================================================#
|
||||
# DEM run script for FinesColumn testcase
|
||||
# init FinesColumn
|
||||
# Thomas Lichtenegger - January 2017
|
||||
# Thomas Lichtenegger - June 2023
|
||||
#===================================================================#
|
||||
|
||||
#- source CFDEM env vars
|
||||
|
||||
@ -0,0 +1,43 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: 2.3.0 |
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
location "constant";
|
||||
object chemistryProperties;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
//chemistry off;
|
||||
|
||||
chemistryType
|
||||
{
|
||||
solver ode;
|
||||
}
|
||||
|
||||
chemistry on;
|
||||
|
||||
initialChemicalTimeStep 1e-07;
|
||||
|
||||
EulerImplicitCoeffs
|
||||
{
|
||||
cTauChem 0.05;
|
||||
equilibriumRateLimiter off;
|
||||
}
|
||||
|
||||
odeCoeffs
|
||||
{
|
||||
solver seulex;
|
||||
absTol 1e-8;
|
||||
relTol 1e-1;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,37 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: 2.3.0 |
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
location "constant";
|
||||
object combustionProperties;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
// OF5 -> OF6: noCombustion<> -> none (template parameters are no longer required)
|
||||
//combustionModel noCombustion<rhoThermoCombustion>; //PaSR<rhoChemistryCombustion>; //PaSR<psiChemistryCombustion>;
|
||||
combustionModel none; // OF6
|
||||
|
||||
active false; //true;
|
||||
|
||||
//noCombustionCoeffs
|
||||
noneCoeffs
|
||||
{
|
||||
}
|
||||
|
||||
|
||||
PaSRCoeffs
|
||||
{
|
||||
Cmix 1.0;
|
||||
turbulentReaction off;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,242 @@
|
||||
/*---------------------------------------------------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: 1.4 |
|
||||
| \\ / A nd | Web: http://www.openfoam.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
|
||||
root "";
|
||||
case "";
|
||||
instance "";
|
||||
local "";
|
||||
|
||||
class dictionary;
|
||||
object couplingProperties;
|
||||
}
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
//===========================================================================//
|
||||
// sub-models & settings
|
||||
|
||||
modelType "A";
|
||||
|
||||
couplingInterval 1;
|
||||
|
||||
voidFractionModel divided;
|
||||
|
||||
locateModel engine;//turboEngine;//
|
||||
|
||||
meshMotionModel noMeshMotion;
|
||||
|
||||
regionModel allRegion;
|
||||
|
||||
IOModel "off";
|
||||
|
||||
dataExchangeModel twoWayMPI;
|
||||
|
||||
averagingModel dense;
|
||||
|
||||
clockModel standardClock;
|
||||
|
||||
smoothingModel off;
|
||||
|
||||
probeModel off;
|
||||
|
||||
forceModels
|
||||
(
|
||||
freeStreaming
|
||||
BeetstraDrag
|
||||
gradPForce
|
||||
viscForce
|
||||
);
|
||||
|
||||
momCoupleModels
|
||||
(
|
||||
implicitCouple
|
||||
);
|
||||
|
||||
otherForceModels
|
||||
(
|
||||
);
|
||||
|
||||
energyModels
|
||||
(
|
||||
heatTransferGunn
|
||||
reactionHeat
|
||||
);
|
||||
|
||||
thermCondModel SyamlalThermCond;
|
||||
|
||||
chemistryModels
|
||||
(
|
||||
species
|
||||
diffusionCoefficients
|
||||
massTransferCoeff
|
||||
);
|
||||
|
||||
turbulenceModelType "turbulenceProperties";
|
||||
|
||||
coupleRecFluc false;
|
||||
coupleRecForce false;
|
||||
|
||||
//===========================================================================//
|
||||
// sub-model properties
|
||||
|
||||
freeStreamingProps
|
||||
{
|
||||
critVoidfraction 1.5;
|
||||
// density can be chosen arbitrarily
|
||||
particleDensity 667;
|
||||
ballisticVelLimit 0.0;
|
||||
// operational descent velocity
|
||||
scalingFactor 110.04;
|
||||
interpolation true;
|
||||
}
|
||||
|
||||
|
||||
implicitCoupleProps
|
||||
{
|
||||
velFieldName "U";
|
||||
granVelFieldName "Us";
|
||||
voidfractionFieldName "voidfraction";
|
||||
}
|
||||
|
||||
explicitCoupleProps
|
||||
{
|
||||
fLimit (1e3 1e3 1e3);
|
||||
}
|
||||
|
||||
dSauterProps
|
||||
{
|
||||
|
||||
}
|
||||
gradPForceProps
|
||||
{
|
||||
pFieldName "p";
|
||||
densityFieldName "rho";
|
||||
voidfractionFieldName "voidfraction";
|
||||
velocityFieldName "U";
|
||||
}
|
||||
|
||||
GidaspowDragProps
|
||||
{
|
||||
verbose false;
|
||||
velFieldName "U";
|
||||
granVelFieldName "Us";
|
||||
voidfractionFieldName "voidfraction";
|
||||
interpolation true;
|
||||
phi 1;
|
||||
}
|
||||
|
||||
heatTransferGunnProps
|
||||
{
|
||||
partTempName "Temp";
|
||||
partHeatFluxName "convectiveHeatFlux";
|
||||
calcPartTempField false;
|
||||
partRefTemp 293;
|
||||
interpolation false;
|
||||
}
|
||||
|
||||
viscForceProps
|
||||
{
|
||||
velocityFieldName "U";
|
||||
densityFieldName "rho";
|
||||
interpolation;
|
||||
}
|
||||
|
||||
BeetstraDragProps
|
||||
{
|
||||
velFieldName "U";
|
||||
granVelFieldName "Us";
|
||||
densityFieldName "rho";
|
||||
voidfractionFieldName "voidfraction";
|
||||
}
|
||||
|
||||
gravityProps
|
||||
{
|
||||
}
|
||||
|
||||
SyamlalThermCondProps
|
||||
{
|
||||
}
|
||||
|
||||
centreProps
|
||||
{
|
||||
alphaMin 0.10;
|
||||
}
|
||||
|
||||
engineProps
|
||||
{
|
||||
treeSearch true;
|
||||
}
|
||||
|
||||
turboEngineProps
|
||||
{
|
||||
treeSearch true;
|
||||
}
|
||||
|
||||
dividedProps
|
||||
{
|
||||
alphaMin 0.25;
|
||||
scaleUpVol 1.0;
|
||||
weight 1.0; //1.33;
|
||||
verbose;
|
||||
}
|
||||
|
||||
constDiffSmoothingProps
|
||||
{
|
||||
verbose;
|
||||
lowerLimit 0.25;
|
||||
upperLimit 1.0;
|
||||
smoothingLength 0.006;
|
||||
}
|
||||
|
||||
twoWayMPIProps
|
||||
{
|
||||
maxNumberOfParticles 50000;
|
||||
liggghtsPath "../DEM/in.liggghts_dataDrivenCFD";
|
||||
}
|
||||
|
||||
|
||||
particleProbeProps
|
||||
{
|
||||
particleIDsToSample "0";
|
||||
}
|
||||
|
||||
speciesProps
|
||||
{
|
||||
ChemistryFile "$FOAM_CASE/constant/foam.inp";
|
||||
Nevery 1;
|
||||
verbose false;
|
||||
interpolation false;
|
||||
}
|
||||
|
||||
diffusionCoefficientsProps
|
||||
{
|
||||
verbose false;
|
||||
interpolation false;
|
||||
ChemistryFile "$FOAM_CASE/constant/foam.inp";
|
||||
diffusantGasNames (
|
||||
CO CO2
|
||||
);
|
||||
}
|
||||
|
||||
massTransferCoeffProps
|
||||
{
|
||||
verbose false;
|
||||
interpolation false;
|
||||
}
|
||||
|
||||
reactionHeatProps
|
||||
{
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,165 @@
|
||||
/*--------------------------------*- 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;
|
||||
class dictionary;
|
||||
location "constant";
|
||||
object foam.dat;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
CO
|
||||
{
|
||||
specie
|
||||
{
|
||||
molWeight 28.0106;
|
||||
}
|
||||
thermodynamics
|
||||
{
|
||||
Tlow 200;
|
||||
Thigh 3500;
|
||||
Tcommon 1000;
|
||||
highCpCoeffs ( 2.71519 0.00206253 -9.98826e-07 2.30053e-10 -2.03648e-14 -14151.9 7.81869 );
|
||||
lowCpCoeffs ( 3.57953 -0.000610354 1.01681e-06 9.07006e-10 -9.04424e-13 -14344.1 3.50841 );
|
||||
}
|
||||
transport
|
||||
{
|
||||
As 1.67212e-06;
|
||||
Ts 170.672;
|
||||
}
|
||||
}
|
||||
|
||||
CO2
|
||||
{
|
||||
specie
|
||||
{
|
||||
molWeight 44.01;
|
||||
}
|
||||
thermodynamics
|
||||
{
|
||||
Tlow 0;
|
||||
Thigh 5000;
|
||||
Tcommon 1000;
|
||||
highCpCoeffs ( 3.85746 0.00441437 -2.21481e-06 5.2349e-10 -4.72084e-14 -48759.2 2.27164 );
|
||||
lowCpCoeffs ( 2.35677 0.0089846 -7.12356e-06 2.45919e-09 -1.437e-13 -48372 9.90105 );
|
||||
}
|
||||
transport
|
||||
{
|
||||
As 1.67212e-06;
|
||||
Ts 170.6;
|
||||
}
|
||||
}
|
||||
|
||||
CH4
|
||||
{
|
||||
specie
|
||||
{
|
||||
molWeight 16.043;
|
||||
}
|
||||
thermodynamics
|
||||
{
|
||||
Tlow 0;
|
||||
Thigh 5000;
|
||||
Tcommon 1000;
|
||||
highCpCoeffs ( 0.0748515 0.0133909 -5.73286e-06 1.22293e-09 -1.01815e-13 -9468.34 18.4373 );
|
||||
lowCpCoeffs ( 5.14988 -0.013671 4.91801e-05 -4.84743e-08 1.66694e-11 -10246.6 -4.6413 );
|
||||
}
|
||||
transport
|
||||
{
|
||||
As 1.67212e-06;
|
||||
Ts 170.672;
|
||||
}
|
||||
}
|
||||
|
||||
H2
|
||||
{
|
||||
specie
|
||||
{
|
||||
molWeight 2.01594;
|
||||
}
|
||||
thermodynamics
|
||||
{
|
||||
Tlow 0;
|
||||
Thigh 5000;
|
||||
Tcommon 1000;
|
||||
highCpCoeffs ( 3.33728 -4.94025e-05 4.99457e-07 -1.79566e-10 2.00255e-14 -950.159 -3.20502 );
|
||||
lowCpCoeffs ( 2.34433 0.00798052 -1.94782e-05 2.01572e-08 -7.37612e-12 -917.935 0.68301 );
|
||||
}
|
||||
transport
|
||||
{
|
||||
As 1.67212e-06;
|
||||
Ts 170.672;
|
||||
}
|
||||
}
|
||||
|
||||
H2O
|
||||
{
|
||||
specie
|
||||
{
|
||||
molWeight 18.0153;
|
||||
}
|
||||
thermodynamics
|
||||
{
|
||||
Tlow 0;
|
||||
Thigh 5000;
|
||||
Tcommon 1000;
|
||||
highCpCoeffs ( 3.03399 0.00217692 -1.64073e-07 -9.7042e-11 1.68201e-14 -30004.3 4.96677 );
|
||||
lowCpCoeffs ( 4.19864 -0.00203643 6.5204e-06 -5.48797e-09 1.77198e-12 -30293.7 -0.849032 );
|
||||
}
|
||||
transport
|
||||
{
|
||||
As 1.67212e-06;
|
||||
Ts 170.672;
|
||||
}
|
||||
}
|
||||
|
||||
O2
|
||||
{
|
||||
specie
|
||||
{
|
||||
molWeight 31.9988;
|
||||
}
|
||||
thermodynamics
|
||||
{
|
||||
Tlow 0;
|
||||
Thigh 5000;
|
||||
Tcommon 1000;
|
||||
highCpCoeffs ( 3.28254 0.00148309 -7.57967e-07 2.09471e-10 -2.16718e-14 -1088.46 5.45323 );
|
||||
lowCpCoeffs ( 3.78246 -0.00299673 9.8473e-06 -9.6813e-09 3.24373e-12 -1063.94 3.65768 );
|
||||
}
|
||||
transport
|
||||
{
|
||||
As 1.67212e-06;
|
||||
Ts 170.672;
|
||||
}
|
||||
}
|
||||
|
||||
N2
|
||||
{
|
||||
specie
|
||||
{
|
||||
molWeight 28.0134;
|
||||
}
|
||||
thermodynamics
|
||||
{
|
||||
Tlow 0;
|
||||
Thigh 5000;
|
||||
Tcommon 1000;
|
||||
highCpCoeffs ( 2.92664 0.00148798 -5.68476e-07 1.0097e-10 -6.75335e-15 -922.798 5.98053 );
|
||||
lowCpCoeffs ( 3.29868 0.00140824 -3.96322e-06 5.64152e-09 -2.44485e-12 -1020.9 3.95037 );
|
||||
}
|
||||
transport
|
||||
{
|
||||
As 1.67212e-06;
|
||||
Ts 170.672;
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,12 @@
|
||||
species
|
||||
(
|
||||
CO
|
||||
CO2
|
||||
N2
|
||||
);
|
||||
|
||||
reactions
|
||||
{
|
||||
}
|
||||
|
||||
|
||||
22
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/g
Executable file
22
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/g
Executable file
@ -0,0 +1,22 @@
|
||||
/*--------------------------------*- 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 uniformDimensionedVectorField;
|
||||
location "constant";
|
||||
object g;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dimensions [0 1 -2 0 0 0 0];
|
||||
value ( 0 0 -9.81 );
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,42 @@
|
||||
/*---------------------------------------------------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: 1.4 |
|
||||
| \\ / A nd | Web: http://www.openfoam.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
|
||||
root "";
|
||||
case "";
|
||||
instance "";
|
||||
local "";
|
||||
|
||||
class dictionary;
|
||||
object liggghtsCommands;
|
||||
}
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
liggghtsCommandModels
|
||||
(
|
||||
runLiggghts
|
||||
writeLiggghts
|
||||
);
|
||||
// ************************************************************************* //
|
||||
|
||||
/*runLiggghtsProps
|
||||
{
|
||||
preNo false;
|
||||
}
|
||||
|
||||
writeLiggghtsProps
|
||||
{
|
||||
writeLast off;
|
||||
overwrite off;
|
||||
}*/
|
||||
@ -0,0 +1,95 @@
|
||||
/*---------------------------------------------------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: 1.4 |
|
||||
| \\ / A nd | Web: http://www.openfoam.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
|
||||
root "";
|
||||
case "";
|
||||
instance "";
|
||||
local "";
|
||||
|
||||
class dictionary;
|
||||
object recProperties;
|
||||
}
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
//===========================================================================//
|
||||
// sub-models & settings
|
||||
|
||||
|
||||
recModel standardRecModel;
|
||||
|
||||
recNorm EuclideanNorm;
|
||||
|
||||
//recPath predefinedPath;
|
||||
//recPath MarkovPath;
|
||||
recPath noPath;
|
||||
|
||||
volScalarFields
|
||||
(
|
||||
voidfraction
|
||||
);
|
||||
|
||||
volVectorFields
|
||||
(
|
||||
Us
|
||||
);
|
||||
|
||||
surfaceScalarFields
|
||||
(
|
||||
// phi
|
||||
);
|
||||
|
||||
//verbose true;
|
||||
|
||||
//couplingSubStep 1;
|
||||
|
||||
initialRecSteps 20000; //10;
|
||||
|
||||
|
||||
//===========================================================================//
|
||||
// sub-model properties
|
||||
|
||||
standardRecModelProps
|
||||
{
|
||||
dataBases (dataBase);
|
||||
storeAveragedFields true;
|
||||
timeDilationFactor 1.0; //130.0;
|
||||
}
|
||||
|
||||
EuclideanNormProps
|
||||
{
|
||||
fieldType "volScalarField";
|
||||
// fieldName "partType";
|
||||
fieldName "voidfraction";
|
||||
}
|
||||
|
||||
MarkovPathProps
|
||||
{
|
||||
correlationSteps 0;
|
||||
meanIntervalSteps 1;
|
||||
minIntervalSteps 1;
|
||||
intervalWeights (1.0);
|
||||
}
|
||||
|
||||
predefinedPathProps
|
||||
{
|
||||
recPathName "predefRecurrencePath";
|
||||
}
|
||||
|
||||
noPathProps
|
||||
{
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,52 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: 2.3.0 |
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
location "constant";
|
||||
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
|
||||
{
|
||||
type heRhoThermo;
|
||||
mixture reactingMixture;
|
||||
transport sutherland;
|
||||
thermo janaf;
|
||||
energy sensibleInternalEnergy; //sensibleEnthalpy;
|
||||
equationOfState perfectGas;
|
||||
specie specie;
|
||||
}
|
||||
|
||||
|
||||
chemistryReader foamChemistryReader;
|
||||
|
||||
foamChemistryThermoFile "$FOAM_CASE/constant/foam.dat";
|
||||
|
||||
foamChemistryFile "$FOAM_CASE/constant/foam.inp";
|
||||
|
||||
inertSpecie N2;
|
||||
|
||||
inertLowerBound 0.2;
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,45 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: 2.3.0 |
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
location "constant";
|
||||
object transportProperties;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
transportModel Newtonian;
|
||||
|
||||
nu nu [ 0 2 -1 0 0 0 0 ] 1e-05;
|
||||
|
||||
kf kf [ 1 1 -3 -1 0 0 0 ] 0.026;
|
||||
|
||||
Cp Cp [ 0 2 -2 -1 0 0 0 ] 1005;
|
||||
|
||||
Cv Cv [ 0 2 -2 -1 0 0 0 ] 718;
|
||||
|
||||
CrossPowerLawCoeffs
|
||||
{
|
||||
nu0 nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;
|
||||
nuInf nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
|
||||
m m [ 0 0 1 0 0 0 0 ] 1;
|
||||
n n [ 0 0 0 0 0 0 0 ] 1;
|
||||
}
|
||||
|
||||
BirdCarreauCoeffs
|
||||
{
|
||||
nu0 nu0 [ 0 2 -1 0 0 0 0 ] 1e-06;
|
||||
nuInf nuInf [ 0 2 -1 0 0 0 0 ] 1e-06;
|
||||
k k [ 0 0 1 0 0 0 0 ] 0;
|
||||
n n [ 0 0 0 0 0 0 0 ] 1;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,21 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: 2.3.0 |
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
location "constant";
|
||||
object turbulenceProperties;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
simulationType laminar;
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,34 @@
|
||||
#!/bin/sh
|
||||
# Source run functions
|
||||
. $WM_PROJECT_DIR/bin/tools/RunFunctions
|
||||
|
||||
if [ $# -eq 0 ]
|
||||
then
|
||||
dBname="dataBase"
|
||||
else
|
||||
dBname=$1
|
||||
fi
|
||||
|
||||
#rm log.decomposePar
|
||||
#runApplication decomposePar -force
|
||||
|
||||
cpdirs="system constant"
|
||||
for f in $cpdirs
|
||||
do
|
||||
cp -r $f $dBname/$f
|
||||
done
|
||||
cd $dBname
|
||||
rm log.decomposePar
|
||||
runApplication decomposePar -force -time 0:
|
||||
|
||||
rm -rf $cpdirs
|
||||
rm -rf processor*/constant
|
||||
|
||||
for proc in processor*
|
||||
do
|
||||
echo "Transferring decomposed recurrence fields of $dBname/$proc to $proc."
|
||||
rm -r ../$proc/$dBname
|
||||
mkdir ../$proc/$dBname
|
||||
mv $proc/* ../$proc/$dBname/
|
||||
rm -rf $proc
|
||||
done
|
||||
42
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/CO
Executable file
42
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/CO
Executable file
@ -0,0 +1,42 @@
|
||||
/*--------------------------------*- 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 CO;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dimensions [0 0 0 0 0 0 0];
|
||||
|
||||
internalField uniform 0.4;
|
||||
|
||||
boundaryField
|
||||
{
|
||||
outlet
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
walls
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
inlet
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform 0.4;
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
42
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/CO2
Executable file
42
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/CO2
Executable file
@ -0,0 +1,42 @@
|
||||
/*--------------------------------*- 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 CO2;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dimensions [0 0 0 0 0 0 0];
|
||||
|
||||
internalField uniform 0.3;
|
||||
|
||||
boundaryField
|
||||
{
|
||||
outlet
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
walls
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
inlet
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform 0.0;
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
42
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/Ksl
Executable file
42
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/Ksl
Executable file
@ -0,0 +1,42 @@
|
||||
/*--------------------------------*- 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 Ksl;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dimensions [1 -3 -1 0 0 0 0];
|
||||
|
||||
internalField uniform 0;
|
||||
|
||||
boundaryField
|
||||
{
|
||||
walls
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
inlet
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
outlet
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
42
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/N2
Executable file
42
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/N2
Executable file
@ -0,0 +1,42 @@
|
||||
/*--------------------------------*- 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 N2;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dimensions [0 0 0 0 0 0 0];
|
||||
|
||||
internalField uniform 0.60;
|
||||
|
||||
boundaryField
|
||||
{
|
||||
outlet
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
walls
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
inlet
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform 0.60;
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
41
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/T
Executable file
41
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/T
Executable file
@ -0,0 +1,41 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: 1.6 |
|
||||
| \\ / A nd | Web: http://www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class volScalarField;
|
||||
object T;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dimensions [0 0 0 1 0 0 0];
|
||||
|
||||
internalField uniform 1223;
|
||||
|
||||
boundaryField
|
||||
{
|
||||
walls
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform 1223;
|
||||
}
|
||||
|
||||
inlet
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform 1223;
|
||||
}
|
||||
|
||||
outlet
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,46 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: 2.3.0 |
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class volVectorField;
|
||||
object U;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
dimensions [0 1 -1 0 0 0 0];
|
||||
|
||||
internalField uniform (0 0 0);
|
||||
|
||||
boundaryField
|
||||
{
|
||||
walls
|
||||
{
|
||||
//type slip;
|
||||
type uniformFixedValue;
|
||||
uniformValue constant (0 0 0);
|
||||
}
|
||||
|
||||
// ratio of crosssection to inlet currently 1.25, need to scale superficial velocity accordingly
|
||||
inlet
|
||||
{
|
||||
type flowRateInletVelocity;
|
||||
massFlowRate constant 0.00105;
|
||||
rhoInlet 0.5;
|
||||
}
|
||||
|
||||
|
||||
outlet
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
46
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/Us
Executable file
46
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/Us
Executable file
@ -0,0 +1,46 @@
|
||||
/*--------------------------------*- 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 volVectorField;
|
||||
location "0";
|
||||
object Us;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dimensions [0 1 -1 0 0 0 0];
|
||||
|
||||
internalField uniform (0 0 0);
|
||||
|
||||
boundaryField
|
||||
{
|
||||
walls
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
bottom
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
inlet
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
outlet
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,42 @@
|
||||
/*--------------------------------*- 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;
|
||||
object addSource;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dimensions [1 -1 -3 0 0 0 0];
|
||||
|
||||
internalField uniform 0;
|
||||
|
||||
boundaryField
|
||||
{
|
||||
walls
|
||||
{
|
||||
type zeroGradient;
|
||||
value uniform 0;
|
||||
}
|
||||
|
||||
inlet
|
||||
{
|
||||
type zeroGradient;
|
||||
value uniform 0;
|
||||
}
|
||||
outlet
|
||||
{
|
||||
type zeroGradient;
|
||||
value uniform 0;
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,44 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: 2.3.0 |
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class volScalarField;
|
||||
object p;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dimensions [1 -1 -2 0 0 0 0];
|
||||
|
||||
internalField uniform 1e5;
|
||||
|
||||
boundaryField
|
||||
{
|
||||
walls
|
||||
{
|
||||
// type zeroGradient;
|
||||
type fixedFluxPressure;
|
||||
}
|
||||
|
||||
inlet
|
||||
{
|
||||
type zeroGradient;
|
||||
// type fixedFluxExtrapolatedPressure;
|
||||
//type fixedFluxPressure;
|
||||
// value uniform 1.03e5;
|
||||
}
|
||||
|
||||
outlet
|
||||
{
|
||||
type uniformFixedValue;
|
||||
uniformValue constant 1e5;
|
||||
}
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,39 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: 1.6 |
|
||||
| \\ / A nd | Web: http://www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class volVectorField;
|
||||
object particleForces;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dimensions [2 -5 -2 0 0 0 0];
|
||||
|
||||
internalField uniform (0 0 0);
|
||||
|
||||
boundaryField
|
||||
{
|
||||
walls
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
inlet
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
outlet
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,43 @@
|
||||
/*--------------------------------*- 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
|
||||
{
|
||||
walls
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
inlet
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform 1.0;
|
||||
}
|
||||
|
||||
outlet
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform 1.0;
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,267 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | |
|
||||
| \\ / A nd | Copyright (C) 2016 Ehsan Madadi-Kandjani |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
object blockMeshDict;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
// General macros to create cylinder mesh
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
convertToMeters 1.0;
|
||||
|
||||
// Inner square side half
|
||||
|
||||
|
||||
// Inner square side curvature
|
||||
|
||||
|
||||
// cylinder radius
|
||||
|
||||
|
||||
// Height of cylinder
|
||||
|
||||
|
||||
// Base z
|
||||
|
||||
|
||||
// Outlet z
|
||||
|
||||
|
||||
// Number of cells at inner square
|
||||
|
||||
|
||||
// Number of cells between inner square and circle
|
||||
|
||||
|
||||
// Number of cells in the cylinder height
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
|
||||
|
||||
// 45 degree points angle
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
// Half of 45 degree points angle
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
// Inner square x and y position
|
||||
|
||||
// x
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
// y
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
// Circle x and y positions
|
||||
|
||||
// x
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
// y
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
// Inner square x and y position middle curvatures
|
||||
|
||||
// x
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
// y
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
// Circle x and y positions middle curvatures
|
||||
|
||||
// x
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
// y
|
||||
|
||||
|
||||
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
vertices
|
||||
(
|
||||
(0.02 -0.02 0) // Vertex s0b = 0
|
||||
(-0.02 -0.02 0) // Vertex s1b = 1
|
||||
(-0.02 0.02 0) // Vertex s2b = 2
|
||||
(0.02 0.02 0) // Vertex s3b = 3
|
||||
|
||||
(0.0353553390593274 -0.0353553390593274 0) // Vertex r0b = 4
|
||||
(-0.0353553390593274 -0.0353553390593274 0) // Vertex r1b = 5
|
||||
(-0.0353553390593274 0.0353553390593274 0) // Vertex r2b = 6
|
||||
(0.0353553390593274 0.0353553390593274 0) // Vertex r3b = 7
|
||||
|
||||
(0.02 -0.02 0.2) // Vertex s0t = 8
|
||||
(-0.02 -0.02 0.2) // Vertex s1t = 9
|
||||
(-0.02 0.02 0.2) // Vertex s2t = 10
|
||||
(0.02 0.02 0.2) // Vertex s3t = 11
|
||||
|
||||
(0.0353553390593274 -0.0353553390593274 0.2) // Vertex r0t = 12
|
||||
(-0.0353553390593274 -0.0353553390593274 0.2) // Vertex r1t = 13
|
||||
(-0.0353553390593274 0.0353553390593274 0.2) // Vertex r2t = 14
|
||||
(0.0353553390593274 0.0353553390593274 0.2) // Vertex r3t = 15
|
||||
);
|
||||
|
||||
blocks
|
||||
(
|
||||
//block0
|
||||
hex (1 0 3 2 9 8 11 10)
|
||||
square
|
||||
(6 6 20)
|
||||
simpleGrading (1 1 1)
|
||||
|
||||
//block1
|
||||
hex (0 4 7 3 8 12 15 11)
|
||||
innerCircle
|
||||
(4 6 20)
|
||||
simpleGrading (0.5 1 1)
|
||||
|
||||
//block2
|
||||
hex (3 7 6 2 11 15 14 10)
|
||||
innerCircle
|
||||
(4 6 20)
|
||||
simpleGrading (0.5 1 1)
|
||||
|
||||
//block3
|
||||
hex (2 6 5 1 10 14 13 9)
|
||||
innerCircle
|
||||
(4 6 20)
|
||||
simpleGrading (0.5 1 1)
|
||||
|
||||
//block4
|
||||
hex (1 5 4 0 9 13 12 8)
|
||||
innerCircle
|
||||
(4 6 20)
|
||||
simpleGrading (0.5 1 1)
|
||||
);
|
||||
|
||||
edges
|
||||
(
|
||||
//Circle edges
|
||||
arc 7 4 (0.05 0 0)
|
||||
arc 4 5 (3.06161699786839e-18 -0.05 0)
|
||||
arc 5 6 (-0.05 6.12323399573675e-18 0)
|
||||
arc 6 7 (3.06161699786839e-18 0.05 0)
|
||||
|
||||
//Circle edges
|
||||
arc 15 12 (0.05 0 0.2)
|
||||
arc 12 13 (3.06161699786839e-18 -0.05 0.2)
|
||||
arc 13 14 (-0.05 6.12323399573675e-18 0.2)
|
||||
arc 14 15 (3.06161699786839e-18 0.05 0.2)
|
||||
|
||||
arc 3 0 (0.025 0 0)
|
||||
arc 0 1 (0 -0.025 0)
|
||||
arc 1 2 (-0.025 0 0)
|
||||
arc 2 3 (0 0.025 0)
|
||||
|
||||
arc 11 8 (0.025 0 0.2)
|
||||
arc 8 9 (0 -0.025 0.2)
|
||||
arc 9 10 (-0.025 0 0.2)
|
||||
arc 10 11 (0 0.025 0.2)
|
||||
|
||||
);
|
||||
|
||||
patches
|
||||
(
|
||||
wall walls
|
||||
(
|
||||
(4 7 15 12)
|
||||
(5 4 12 13)
|
||||
(6 5 13 14)
|
||||
(7 6 14 15)
|
||||
)
|
||||
|
||||
patch inlet
|
||||
(
|
||||
(3 0 1 2)
|
||||
(3 7 4 0)
|
||||
(2 6 7 3)
|
||||
(1 5 6 2)
|
||||
(0 4 5 1)
|
||||
)
|
||||
|
||||
patch outlet
|
||||
(
|
||||
(11 10 9 8)
|
||||
(11 8 12 15)
|
||||
(10 11 15 14)
|
||||
(9 10 14 13)
|
||||
(8 9 13 12)
|
||||
)
|
||||
);
|
||||
|
||||
mergePatchPairs
|
||||
(
|
||||
);
|
||||
@ -0,0 +1,267 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | |
|
||||
| \\ / A nd | Copyright (C) 2016 Ehsan Madadi-Kandjani |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
`format' ascii;
|
||||
class dictionary;
|
||||
object blockMeshDict;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
// General macros to create cylinder mesh
|
||||
|
||||
changecom(//)changequote([,])
|
||||
define(calc, [esyscmd(perl -e 'use Math::Trig; print ($1)')])
|
||||
define(VCOUNT, 0)
|
||||
define(vlabel, [[// ]Vertex $1 = VCOUNT define($1, VCOUNT)define([VCOUNT], incr(VCOUNT))])
|
||||
|
||||
define(hex2D, hex ($1b $2b $3b $4b $1t $2t $3t $4t))
|
||||
define(btQuad, ($1b $2b $2t $1t))
|
||||
define(topQuad, ($1t $4t $3t $2t))
|
||||
define(bottomQuad, ($1b $2b $3b $4b))
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
convertToMeters 1.0;
|
||||
|
||||
// Inner square side half
|
||||
define(s, 0.02)
|
||||
|
||||
// Inner square side curvature
|
||||
define(sc, 0.025)
|
||||
|
||||
// cylinder radius
|
||||
define(r, 0.05)
|
||||
|
||||
// Height of cylinder
|
||||
define(z, 0.2)
|
||||
|
||||
// Base z
|
||||
define(Zb, 0)
|
||||
|
||||
// Outlet z
|
||||
define(Zt, calc(Zb + z))
|
||||
|
||||
// Number of cells at inner square
|
||||
define(Ns, 6)
|
||||
|
||||
// Number of cells between inner square and circle
|
||||
define(Ni, 4)
|
||||
|
||||
// Number of cells in the cylinder height
|
||||
define(Nz, 20)
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
define(vert, (x$1$2 y$1$2 $3))
|
||||
define(evert, (ex$1$2 ey$1$2 $3))
|
||||
|
||||
// 45 degree points angle
|
||||
define(a0, -45)
|
||||
define(a1, -135)
|
||||
define(a2, 135)
|
||||
define(a3, 45)
|
||||
|
||||
// Half of 45 degree points angle
|
||||
define(ea0, 0)
|
||||
define(ea1, -90)
|
||||
define(ea2, 180)
|
||||
define(ea3, 90)
|
||||
|
||||
define(ca0, calc(cos((pi/180)*a0)))
|
||||
define(ca1, calc(cos((pi/180)*a1)))
|
||||
define(ca2, calc(cos((pi/180)*a2)))
|
||||
define(ca3, calc(cos((pi/180)*a3)))
|
||||
|
||||
define(sa0, calc(sin((pi/180)*a0)))
|
||||
define(sa1, calc(sin((pi/180)*a1)))
|
||||
define(sa2, calc(sin((pi/180)*a2)))
|
||||
define(sa3, calc(sin((pi/180)*a3)))
|
||||
|
||||
define(cea0, calc(cos((pi/180)*ea0)))
|
||||
define(cea1, calc(cos((pi/180)*ea1)))
|
||||
define(cea2, calc(cos((pi/180)*ea2)))
|
||||
define(cea3, calc(cos((pi/180)*ea3)))
|
||||
|
||||
define(sea0, calc(sin((pi/180)*ea0)))
|
||||
define(sea1, calc(sin((pi/180)*ea1)))
|
||||
define(sea2, calc(sin((pi/180)*ea2)))
|
||||
define(sea3, calc(sin((pi/180)*ea3)))
|
||||
|
||||
// Inner square x and y position
|
||||
|
||||
// x
|
||||
define(x00, s)
|
||||
define(x01, calc(-1.0*s))
|
||||
define(x02, calc(-1.0*s))
|
||||
define(x03, s)
|
||||
|
||||
// y
|
||||
define(y00, calc(-1.0*s))
|
||||
define(y01, calc(-1.0*s))
|
||||
define(y02, s)
|
||||
define(y03, s)
|
||||
|
||||
// Circle x and y positions
|
||||
|
||||
// x
|
||||
define(x10, calc(r*ca0))
|
||||
define(x11, calc(r*ca1))
|
||||
define(x12, calc(r*ca2))
|
||||
define(x13, calc(r*ca3))
|
||||
|
||||
// y
|
||||
define(y10, calc(r*sa0))
|
||||
define(y11, calc(r*sa1))
|
||||
define(y12, calc(r*sa2))
|
||||
define(y13, calc(r*sa3))
|
||||
|
||||
// Inner square x and y position middle curvatures
|
||||
|
||||
// x
|
||||
define(ex00, sc)
|
||||
define(ex01, 0)
|
||||
define(ex02, calc(-1.0*sc))
|
||||
define(ex03, 0)
|
||||
|
||||
// y
|
||||
define(ey00, 0)
|
||||
define(ey01, calc(-1.0*sc))
|
||||
define(ey02, 0)
|
||||
define(ey03, sc)
|
||||
|
||||
// Circle x and y positions middle curvatures
|
||||
|
||||
// x
|
||||
define(ex10, calc(r*cea0))
|
||||
define(ex11, calc(r*cea1))
|
||||
define(ex12, calc(r*cea2))
|
||||
define(ex13, calc(r*cea3))
|
||||
|
||||
// y
|
||||
define(ey10, calc(r*sea0))
|
||||
define(ey11, calc(r*sea1))
|
||||
define(ey12, calc(r*sea2))
|
||||
define(ey13, calc(r*sea3))
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
vertices
|
||||
(
|
||||
vert(0, 0, Zb) vlabel(s0b)
|
||||
vert(0, 1, Zb) vlabel(s1b)
|
||||
vert(0, 2, Zb) vlabel(s2b)
|
||||
vert(0, 3, Zb) vlabel(s3b)
|
||||
|
||||
vert(1, 0, Zb) vlabel(r0b)
|
||||
vert(1, 1, Zb) vlabel(r1b)
|
||||
vert(1, 2, Zb) vlabel(r2b)
|
||||
vert(1, 3, Zb) vlabel(r3b)
|
||||
|
||||
vert(0, 0, Zt) vlabel(s0t)
|
||||
vert(0, 1, Zt) vlabel(s1t)
|
||||
vert(0, 2, Zt) vlabel(s2t)
|
||||
vert(0, 3, Zt) vlabel(s3t)
|
||||
|
||||
vert(1, 0, Zt) vlabel(r0t)
|
||||
vert(1, 1, Zt) vlabel(r1t)
|
||||
vert(1, 2, Zt) vlabel(r2t)
|
||||
vert(1, 3, Zt) vlabel(r3t)
|
||||
);
|
||||
|
||||
blocks
|
||||
(
|
||||
//block0
|
||||
hex2D(s1, s0, s3, s2)
|
||||
square
|
||||
(Ns Ns Nz)
|
||||
simpleGrading (1 1 1)
|
||||
|
||||
//block1
|
||||
hex2D(s0, r0, r3, s3)
|
||||
innerCircle
|
||||
(Ni Ns Nz)
|
||||
simpleGrading (0.5 1 1)
|
||||
|
||||
//block2
|
||||
hex2D(s3, r3, r2, s2)
|
||||
innerCircle
|
||||
(Ni Ns Nz)
|
||||
simpleGrading (0.5 1 1)
|
||||
|
||||
//block3
|
||||
hex2D(s2, r2, r1, s1)
|
||||
innerCircle
|
||||
(Ni Ns Nz)
|
||||
simpleGrading (0.5 1 1)
|
||||
|
||||
//block4
|
||||
hex2D(s1, r1, r0, s0)
|
||||
innerCircle
|
||||
(Ni Ns Nz)
|
||||
simpleGrading (0.5 1 1)
|
||||
);
|
||||
|
||||
edges
|
||||
(
|
||||
//Circle edges
|
||||
arc r3b r0b evert(1, 0, Zb)
|
||||
arc r0b r1b evert(1, 1, Zb)
|
||||
arc r1b r2b evert(1, 2, Zb)
|
||||
arc r2b r3b evert(1, 3, Zb)
|
||||
|
||||
//Circle edges
|
||||
arc r3t r0t evert(1, 0, Zt)
|
||||
arc r0t r1t evert(1, 1, Zt)
|
||||
arc r1t r2t evert(1, 2, Zt)
|
||||
arc r2t r3t evert(1, 3, Zt)
|
||||
|
||||
arc s3b s0b evert(0, 0, Zb)
|
||||
arc s0b s1b evert(0, 1, Zb)
|
||||
arc s1b s2b evert(0, 2, Zb)
|
||||
arc s2b s3b evert(0, 3, Zb)
|
||||
|
||||
arc s3t s0t evert(0, 0, Zt)
|
||||
arc s0t s1t evert(0, 1, Zt)
|
||||
arc s1t s2t evert(0, 2, Zt)
|
||||
arc s2t s3t evert(0, 3, Zt)
|
||||
|
||||
);
|
||||
|
||||
patches
|
||||
(
|
||||
wall walls
|
||||
(
|
||||
btQuad(r0, r3)
|
||||
btQuad(r1, r0)
|
||||
btQuad(r2, r1)
|
||||
btQuad(r3, r2)
|
||||
)
|
||||
|
||||
patch inlet
|
||||
(
|
||||
bottomQuad(s3, s0, s1, s2)
|
||||
bottomQuad(s3, r3, r0, s0)
|
||||
bottomQuad(s2, r2, r3, s3)
|
||||
bottomQuad(s1, r1, r2, s2)
|
||||
bottomQuad(s0, r0, r1, s1)
|
||||
)
|
||||
|
||||
patch outlet
|
||||
(
|
||||
topQuad(s3, s0, s1, s2)
|
||||
topQuad(s3, r3, r0, s0)
|
||||
topQuad(s2, r2, r3, s3)
|
||||
topQuad(s1, r1, r2, s2)
|
||||
topQuad(s0, r0, r1, s1)
|
||||
)
|
||||
);
|
||||
|
||||
mergePatchPairs
|
||||
(
|
||||
);
|
||||
@ -0,0 +1,57 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: 2.3.0 |
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
location "system";
|
||||
object controlDict;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
application potentialFoam;
|
||||
|
||||
startFrom startTime;
|
||||
|
||||
startTime 0;
|
||||
|
||||
stopAt endTime;
|
||||
|
||||
endTime 8000;
|
||||
|
||||
deltaT 0.125;
|
||||
|
||||
writeControl timeStep;
|
||||
|
||||
writeInterval 100000000;
|
||||
|
||||
purgeWrite 0;
|
||||
|
||||
writeFormat ascii;
|
||||
|
||||
writePrecision 6;
|
||||
|
||||
writeCompression off;
|
||||
|
||||
timeFormat general;
|
||||
|
||||
timePrecision 6;
|
||||
|
||||
runTimeModifiable true;
|
||||
|
||||
|
||||
functions
|
||||
{
|
||||
moleFrac
|
||||
{
|
||||
type rhoReactionThermoMoleFractions;
|
||||
}
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,31 @@
|
||||
/*-------------------------------*- C++ -*---------------------------------*\
|
||||
| ========= |
|
||||
| \\ / OpenFOAM |
|
||||
| \\ / |
|
||||
| \\ / The Open Source CFD Toolbox |
|
||||
| \\/ http://www.OpenFOAM.org |
|
||||
\*-------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
note "mesh decomposition control dictionary";
|
||||
location "system";
|
||||
object decomposeParDict;
|
||||
}
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
numberOfSubdomains 8;
|
||||
|
||||
method scotch;
|
||||
|
||||
simpleCoeffs
|
||||
{
|
||||
n (1 1 8);
|
||||
delta 0.001;
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
|
||||
@ -0,0 +1,32 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: 2.3.0 |
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
location "system";
|
||||
object fvOptions;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
// OF 6
|
||||
|
||||
source1
|
||||
{
|
||||
type limitTemperature;
|
||||
active yes;
|
||||
limitTemperatureCoeffs
|
||||
{
|
||||
active yes;
|
||||
selectionMode all;
|
||||
min 1222;
|
||||
max 1224;
|
||||
}
|
||||
}
|
||||
|
||||
@ -0,0 +1,93 @@
|
||||
/*--------------------------------*- 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 dictionary;
|
||||
location "system";
|
||||
object fvSchemes;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
ddtSchemes
|
||||
{
|
||||
default Euler;
|
||||
}
|
||||
|
||||
gradSchemes
|
||||
{
|
||||
default Gauss linear;
|
||||
grad(p) cellLimited Gauss linear 1;
|
||||
// grad(U) Gauss linear;
|
||||
grad(U) cellLimited Gauss linear 1;
|
||||
}
|
||||
|
||||
divSchemes
|
||||
{
|
||||
default Gauss linear;
|
||||
div(phi,U) Gauss upwind;
|
||||
div(phi,k) Gauss limitedLinear 1;
|
||||
div(phi,epsilon) Gauss limitedLinear 1;
|
||||
div(phi,R) Gauss limitedLinear 1;
|
||||
div(R) Gauss linear;
|
||||
div(phi,nuTilda) Gauss limitedLinear 1;
|
||||
div((viscousTerm*dev(grad(U).T()))) Gauss linear;
|
||||
div((nu*dev(grad(U).T()))) Gauss linear;
|
||||
div(phi,Yi_h) Gauss multivariateSelection
|
||||
{
|
||||
CO limitedLinear01 1;
|
||||
H2 limitedLinear01 1;
|
||||
N2 limitedLinear01 1;
|
||||
H2O limitedLinear01 1;
|
||||
CO2 limitedLinear01 1;
|
||||
O2 limitedLinear01 1;
|
||||
h limitedLinear 1;
|
||||
e limitedLinear 1;
|
||||
}
|
||||
div((nuEff*dev(grad(U).T()))) Gauss linear;
|
||||
|
||||
}
|
||||
|
||||
laplacianSchemes
|
||||
{
|
||||
/*
|
||||
default Gauss linear corrected;
|
||||
laplacian(viscousTerm,U) Gauss linear corrected;
|
||||
laplacian(nu,U) Gauss linear corrected;
|
||||
laplacian(nuEff,U) Gauss linear corrected;
|
||||
laplacian((1|A(U)),p) Gauss linear corrected;
|
||||
laplacian((voidfraction2|A(U)),p) Gauss linear corrected;
|
||||
laplacian(DkEff,k) Gauss linear corrected;
|
||||
laplacian(DepsilonEff,epsilon) Gauss linear corrected;
|
||||
laplacian(DREff,R) Gauss linear corrected;
|
||||
laplacian(DnuTildaEff,nuTilda) Gauss linear corrected;
|
||||
*/
|
||||
default Gauss linear limited 0.333;
|
||||
}
|
||||
|
||||
interpolationSchemes
|
||||
{
|
||||
default linear;
|
||||
interpolate(U) linear;
|
||||
}
|
||||
|
||||
snGradSchemes
|
||||
{
|
||||
// default corrected;
|
||||
default limited 0.333;
|
||||
}
|
||||
|
||||
fluxRequired
|
||||
{
|
||||
default no;
|
||||
p ;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,150 @@
|
||||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: 2.3.0 |
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
location "system";
|
||||
object fvSolution;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
solvers
|
||||
{
|
||||
"(p|rho|correctedField|partTemp)"
|
||||
{
|
||||
solver PCG;
|
||||
preconditioner DIC;
|
||||
tolerance 1e-7;
|
||||
relTol 1e-2;
|
||||
}
|
||||
|
||||
"(p|rho|correctedField|partTemp)Final"
|
||||
{
|
||||
$p;
|
||||
relTol 0;
|
||||
}
|
||||
|
||||
"(U|k|e|epsilon|R|nuTilda|c)"
|
||||
{
|
||||
solver smoothSolver;
|
||||
smoother GaussSeidel;
|
||||
tolerance 1e-8;
|
||||
relTol 1e-2;
|
||||
}
|
||||
|
||||
"(U|e|k|nuTilda)Final"
|
||||
{
|
||||
$U;
|
||||
relTol 0;
|
||||
}
|
||||
|
||||
"(Yi|CO|CO2|H2|H2O|O2|N2)"
|
||||
{
|
||||
solver smoothSolver;
|
||||
smoother symGaussSeidel;
|
||||
tolerance 1e-07;
|
||||
relTol 0.01;
|
||||
maxIter 100;
|
||||
}
|
||||
|
||||
"(Yi|CO|CO2|H2|H2O|O2|N2)Final"
|
||||
{
|
||||
$Yi;
|
||||
tolerance 1e-06;
|
||||
relTol 0;
|
||||
maxIter 100;
|
||||
}
|
||||
|
||||
T
|
||||
{
|
||||
solver PBiCG;
|
||||
preconditioner DILU;
|
||||
tolerance 1e-7;
|
||||
relTol 0;
|
||||
}
|
||||
|
||||
|
||||
|
||||
"(voidfraction|Ksl|UsNext|voidfractionNext|sSmoothField|vSmoothField|fSmooth)"
|
||||
{
|
||||
solver PCG;
|
||||
preconditioner DIC;
|
||||
tolerance 1e-09;
|
||||
relTol 1e-06;
|
||||
}
|
||||
}
|
||||
|
||||
PIMPLE
|
||||
{
|
||||
momentumPredictor yes;
|
||||
nOuterCorrectors 10;
|
||||
nCorrectors 1;
|
||||
nNonOrthogonalCorrectors 0;
|
||||
rhoMin rhoMin [ 1 -3 0 0 0 ] 0.1;
|
||||
rhoMax rhoMax [ 1 -3 0 0 0 ] 5.0;
|
||||
|
||||
pMin pMin [ 1 -1 -2 0 0 ] 90000;
|
||||
pMax pMax [ 1 -1 -2 0 0 ] 500000;
|
||||
|
||||
// pMin 90000;
|
||||
// pMax 500000;
|
||||
|
||||
UMax UMax [ 0 1 -1 0 0 ] 250;
|
||||
}
|
||||
|
||||
PISO
|
||||
{
|
||||
nCorrectors 4;
|
||||
nNonOrthogonalCorrectors 1;
|
||||
pRefCell 0;
|
||||
pRefValue 0;
|
||||
}
|
||||
|
||||
relaxationFactors
|
||||
{
|
||||
fields
|
||||
{
|
||||
U 0.3;
|
||||
p 0.3;
|
||||
UFinal 0.3;
|
||||
pFinal 0.3;
|
||||
CO 0.3;
|
||||
COFinal 0.3;
|
||||
CO2 0.3;
|
||||
CO2Final 0.3;
|
||||
O2 0.3;
|
||||
O2Final 0.3;
|
||||
N2 0.3;
|
||||
N2Final 0.3;
|
||||
e 0.3;
|
||||
he 0.3;
|
||||
eFinal 0.3;
|
||||
heFinal 0.3;
|
||||
}
|
||||
equations
|
||||
{
|
||||
// CO 0.3;
|
||||
// COFinal 0.3;
|
||||
// CO2 0.3;
|
||||
// CO2Final 0.3;
|
||||
// O2 0.3;
|
||||
// O2Final 0.3;
|
||||
// N2 0.3;
|
||||
// N2Final 0.3;
|
||||
// e 0.3;
|
||||
// he 0.3;
|
||||
// eFinal 0.3;
|
||||
// heFinal 0.3;
|
||||
".*" 0.3;
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
@ -0,0 +1,30 @@
|
||||
# DEM-parameters
|
||||
variable partDens equal 3300 # density
|
||||
variable partDensHeavy equal 8000 # density
|
||||
#variable partRad equal 0.001 # particle radius
|
||||
variable partRad equal 0.0055 # particle radius
|
||||
variable partVol equal 4*3.14/3*${partRad}^3
|
||||
variable yM equal 5.e6 # Young's modulus
|
||||
variable pR equal 0.25 # Poisson ratio
|
||||
variable rest equal 0.66 # coefficient of restitution
|
||||
variable fric equal 0.95 # coefficient of friction
|
||||
variable roll equal 0.15 # for particles and walls with friction
|
||||
variable noRoll equal 0 # for frictionless walls
|
||||
variable cV equal 0.1 # characteristic velocity
|
||||
|
||||
variable Tpart equal 298
|
||||
|
||||
variable dt equal 0.125
|
||||
variable skin equal 0.005
|
||||
|
||||
# particle settling loop parameters
|
||||
variable vthre equal 0.0001
|
||||
|
||||
# Geometry
|
||||
variable cylrad equal 0.0415
|
||||
|
||||
variable fillH equal 0.35
|
||||
variable fillH1 equal ${fillH}-0.025
|
||||
variable fillH2 equal ${fillH1}+2*0.0024
|
||||
variable fillV equal ${fillH}*${cylrad}^2*3.14
|
||||
variable nPart equal ceil(${fillV}*0.7/${partVol})
|
||||
@ -0,0 +1,32 @@
|
||||
# particle properties
|
||||
#####################
|
||||
|
||||
# ore
|
||||
# obtain H volume fraction from initial particle density (fully unreduced); then compute volume fractions for inner layers according to Kinaci et al. CES 227 (2020) [ATTENTION: correct conversion factors are q_HM = 0.9666, q_MW = 0.9309, q_WFe = 0.7773]; set porosities = 1 - volume fraction
|
||||
|
||||
#fix porosity_ore all property/atom porosity_ore vector yes no no 0.61 0.34 0.19 0.17 # values for ore in Kinaci et al. CES 227 (2020), corresponding to rho = 4300 for hematite
|
||||
fix porosity_ore all property/global porosity_ore vector 0.71 0.48 0.38 0.37
|
||||
fix tortuosity_ore all property/global tortuosity_ore scalar 3
|
||||
fix pore_diameter_ore all property/global pore_diameter_ore vector 3.e-07 1.75e-05 1.1e-05
|
||||
fix layerDensities_ore all property/global density_ore vector 7870. 5740. 5170. 5240.
|
||||
|
||||
|
||||
# define layer properties
|
||||
fix LayerRelRadii all property/atom relRadii vector yes no no 1.0 0.998 0.995 0.98
|
||||
|
||||
# define fix for rho_eff and fracRed - initialize as zero
|
||||
fix fracRed all property/atom fracRed vector yes no no 0. 0. 0.
|
||||
fix rhoeff all property/atom rhoeff vector yes no no 0. 0. 0. 0.
|
||||
|
||||
# define fix for mass layer - initial testing
|
||||
fix LayerMasses all property/atom massLayer vector yes no no 0. 0. 0. 0.
|
||||
|
||||
|
||||
# ore reactions
|
||||
###############
|
||||
# Activate for 3-layer unreacted core shrink model
|
||||
fix chem0 ore chem/shrink/core speciesA CO molMassA 0.02801 speciesC CO2 molMassC 0.04401 nevery 1 screen no cg 1 limit_reactant_consumption no
|
||||
|
||||
# Material properties for chemical reaction
|
||||
fix k0_chem0 ore property/global k0_chem0 vector 10 18 3700
|
||||
fix Ea_chem0 ore property/global Ea_chem0 vector 69488 73674 113859
|
||||
@ -0,0 +1,23 @@
|
||||
##############
|
||||
# cfd coupling
|
||||
##############
|
||||
variable Tin equal 1223
|
||||
|
||||
####
|
||||
# uncomment if molten iron to be pictured as trickling particles
|
||||
#variable vtrickle equal -0.05
|
||||
#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 1 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 cfd1b all couple/cfd/recurrence transfer_force no transfer_fluctuations no
|
||||
fix cfd2 all couple/cfd/convection T0 ${Tin} transfer_conduction no max_change 1.0
|
||||
fix cfd3 all couple/cfd/parttempfield T0 ${Tin}
|
||||
fix cfd4 all couple/cfd/chemistry n_species 3 species_names CO CO2 N2 n_diff 2 diffusant_names CO CO2 transfer_reactant no init_layer_radii no
|
||||
|
||||
|
||||
#fix integr all nve/sphere/limit vlimit 0.05 omegalimit 100
|
||||
|
||||
run 0
|
||||
velocity all set 0.0 0.0 0.0
|
||||
@ -0,0 +1,42 @@
|
||||
|
||||
## material properties
|
||||
include ../DEM/DEMparameters
|
||||
|
||||
###############################################
|
||||
atom_style granular
|
||||
atom_modify map array
|
||||
communicate single vel yes
|
||||
|
||||
boundary f f f
|
||||
newton off
|
||||
units si
|
||||
|
||||
processors 1 1 8
|
||||
|
||||
# option B
|
||||
read_restart ../DEM/liggghts.restart
|
||||
|
||||
group dummy type 1
|
||||
group ore type 2
|
||||
|
||||
reset_timestep 0
|
||||
|
||||
|
||||
neighbor ${skin} bin
|
||||
neigh_modify delay 0
|
||||
group fullyDeformed empty
|
||||
neigh_modify exclude group fullyDeformed all
|
||||
|
||||
timestep 0.125
|
||||
|
||||
include ../DEM/in.material_parameters
|
||||
include ../DEM/in.chemistry
|
||||
include ../DEM/in.coupling_integration
|
||||
#include ../DEM/in.output
|
||||
include ../DEM/in.monitor
|
||||
|
||||
########
|
||||
# limits
|
||||
########
|
||||
# particle temperatures
|
||||
fix limTpart all limit/property/atom Temp min 1222 max 1224
|
||||
93
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.liggghts_init
Executable file
93
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.liggghts_init
Executable file
@ -0,0 +1,93 @@
|
||||
###############################################
|
||||
# material properties
|
||||
include ../DEM/DEMparameters
|
||||
|
||||
###############################################
|
||||
|
||||
atom_style granular
|
||||
atom_modify map array
|
||||
communicate single vel yes
|
||||
|
||||
boundary f f f
|
||||
newton off
|
||||
units si
|
||||
#processors 1 1 8
|
||||
|
||||
region reg block -${cylrad} ${cylrad} -${cylrad} ${cylrad} 0.0 0.45 units box
|
||||
create_box 2 reg
|
||||
|
||||
neighbor 0.001 bin
|
||||
neigh_modify delay 0
|
||||
|
||||
#Material properties required for granular pair styles
|
||||
|
||||
fix m1 all property/global youngsModulus peratomtype ${yM} ${yM}
|
||||
fix m2 all property/global poissonsRatio peratomtype ${pR} ${pR}
|
||||
fix m3 all property/global coefficientRestitution peratomtypepair 2 ${rest} ${rest} ${rest} ${rest}
|
||||
fix m4 all property/global coefficientFriction peratomtypepair 2 ${fric} ${fric} ${fric} ${fric}
|
||||
fix m5 all property/global coefficientRollingFriction peratomtypepair 2 ${roll} ${roll} ${roll} ${roll}
|
||||
fix m6 all property/global characteristicVelocity scalar ${cV}
|
||||
|
||||
#pair style
|
||||
pair_style gran model hertz tangential history rolling_friction cdt
|
||||
pair_coeff * *
|
||||
|
||||
#timestep, gravity
|
||||
timestep 0.00001
|
||||
fix gravi all gravity 9.81 vector 0.0 0.0 -1.0
|
||||
|
||||
fix zwall1 all wall/gran model hertz tangential history primitive type 1 zplane 0.0
|
||||
#fix zwall2 all wall/gran model hertz tangential history primitive type 1 zplane 1.3
|
||||
fix cylwalls all wall/gran model hertz tangential history primitive type 1 zcylinder ${cylrad} 0. 0.
|
||||
|
||||
fix integr all nve/sphere
|
||||
|
||||
#fix ts_check all check/timestep/gran 10000 0.1 0.1
|
||||
|
||||
#screen output
|
||||
compute rke all erotate/sphere
|
||||
thermo_style custom step atoms ke c_rke vol
|
||||
thermo 5000
|
||||
thermo_modify lost ignore norm no
|
||||
compute_modify thermo_temp dynamic yes
|
||||
|
||||
fix pts1 all particletemplate/sphere 1 atom_type 1 density constant ${partDens} radius constant ${partRad}
|
||||
fix pdd1 all particledistribution/discrete 1.0 1 pts1 1.0
|
||||
fix pts2 all particletemplate/sphere 1 atom_type 2 density constant ${partDens} radius constant ${partRad}
|
||||
fix pdd2 all particledistribution/discrete 1.0 1 pts2 1.0
|
||||
|
||||
|
||||
region insReg0 cylinder z 0. 0. ${cylrad} 0.0 0.05 units box
|
||||
fix ins0 all insert/pack/dense seed 123 region insReg0 distributiontemplate pdd1 volumefraction_region 0.57
|
||||
run 1
|
||||
|
||||
region insReg1 cylinder z 0. 0. ${cylrad} 0.055 0.25 units box
|
||||
fix ins1 all insert/pack seed 1001 distributiontemplate pdd2 insert_every once overlapcheck yes mass_in_region 0.5 region insReg1 ntry_mc 10000
|
||||
run 1
|
||||
|
||||
region insReg2 cylinder z 0. 0. ${cylrad} 0.26 0.4 units box
|
||||
fix ins2 all insert/pack/dense seed 123 region insReg2 distributiontemplate pdd1 volumefraction_region 0.57
|
||||
run 1
|
||||
|
||||
|
||||
#======================================================================
|
||||
#- settling
|
||||
|
||||
variable vcmx equal vcm(all,x)
|
||||
variable vcmy equal vcm(all,y)
|
||||
variable vcmz equal vcm(all,z)
|
||||
|
||||
variable vave equal sqrt(v_vcmx*v_vcmx+v_vcmy*v_vcmy+v_vcmz*v_vcmz)
|
||||
print ${vave}
|
||||
|
||||
run 1000000
|
||||
|
||||
#======================================================================
|
||||
region delReg block -${cylrad} ${cylrad} -${cylrad} ${cylrad} 0.2 INF units box
|
||||
run 1 every 1 "delete_atoms region delReg"
|
||||
fix zwall2 all wall/gran model hertz tangential history primitive type 1 zplane 0.2
|
||||
run 5000
|
||||
run 1 every 1 "delete_atoms region delReg"
|
||||
run 5000
|
||||
write_restart liggghts.restart
|
||||
write_dump all custom/vtk ../DEM/finalState.vtk id type radius mass density x y z
|
||||
@ -0,0 +1,12 @@
|
||||
fix m1 all property/global youngsModulus peratomtype ${yM} ${yM}
|
||||
fix m2 all property/global poissonsRatio peratomtype ${pR} ${pR}
|
||||
fix m3 all property/global coefficientRestitution peratomtypepair 2 ${rest} ${rest} ${rest} ${rest}
|
||||
fix m4 all property/global coefficientFriction peratomtypepair 2 ${fric} ${fric} ${fric} ${fric}
|
||||
fix m5 all property/global coefficientRollingFriction peratomtypepair 2 ${roll} ${roll} ${roll} ${roll}
|
||||
fix m6 all property/global characteristicVelocity scalar ${cV}
|
||||
|
||||
fix ftco all property/global thermalConductivity peratomtype 1.7 80
|
||||
# constant Cp
|
||||
#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
|
||||
@ -0,0 +1,15 @@
|
||||
variable time equal time
|
||||
|
||||
variable m1 equal mass(ore)
|
||||
print 'initial mass = ${m1}'
|
||||
compute layermasses all reduce sum f_LayerMasses[1] f_LayerMasses[2] f_LayerMasses[3] f_LayerMasses[4]
|
||||
fix layermasses all ave/time 1 1 1 c_layermasses[1] c_layermasses[2] c_layermasses[3] c_layermasses[4]
|
||||
variable lmass1 equal f_layermasses[1]
|
||||
variable lmass2 equal f_layermasses[2]
|
||||
variable lmass3 equal f_layermasses[3]
|
||||
variable lmass4 equal f_layermasses[4]
|
||||
variable lmasstot equal f_layermasses[1]+f_layermasses[2]+f_layermasses[3]+f_layermasses[4]
|
||||
|
||||
fix printmasslayer all print 80 "${time} ${m1} ${lmass1} ${lmass2} ${lmass3} ${lmass4} ${lmasstot}" &
|
||||
file ../DEM/monitor/mass_layer.dat title "#time mass lmass1 lmass2 lmass3 lmass4 lmass_{tot}" screen no
|
||||
|
||||
@ -0,0 +1,93 @@
|
||||
########
|
||||
# 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
|
||||
|
||||
|
||||
# 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
|
||||
|
||||
# mass-weighted temperature
|
||||
###########################
|
||||
|
||||
variable mT atom mass*f_Temp
|
||||
compute csumMT all reduce sum v_mT
|
||||
variable sumMT equal c_csumMT
|
||||
variable totMass equal mass(all)
|
||||
variable TaveM equal v_sumMT/v_totMass
|
||||
|
||||
fix printheat all print ${Nmonitor} "${time} ${Tave} ${TaveM}" append ${dumppath}/monitor/temp_ave.txt title "#time T_ave T_ave_mass"
|
||||
|
||||
|
||||
# chemistry output
|
||||
##################
|
||||
|
||||
variable nore equal count(ore)
|
||||
variable ncoke equal count(coke)
|
||||
|
||||
## Check Layer Radii and Write them to files ###
|
||||
compute layerRadO1 ore reduce sum f_LayerRelRadii[1]
|
||||
fix redRadO1 ore ave/time 1 1 1 c_layerRadO1
|
||||
variable rrO1 equal f_redRadO1/v_nore
|
||||
|
||||
compute layerRadO2 ore reduce sum f_LayerRelRadii[2]
|
||||
fix redRadO2 ore ave/time 1 1 1 c_layerRadO2
|
||||
variable rrO2 equal f_redRadO2/v_nore
|
||||
|
||||
compute layerRadO3 ore reduce sum f_LayerRelRadii[3]
|
||||
fix redRadO3 ore ave/time 1 1 1 c_layerRadO3
|
||||
variable rrO3 equal f_redRadO3/v_nore
|
||||
|
||||
compute layerRadO4 ore reduce sum f_LayerRelRadii[4]
|
||||
fix redRadO4 ore ave/time 1 1 1 c_layerRadO4
|
||||
variable rrO4 equal f_redRadO4/v_nore
|
||||
|
||||
compute layerRadC1 coke reduce sum f_LayerRelRadii[1]
|
||||
fix redRadC1 coke ave/time 1 1 1 c_layerRadC1
|
||||
variable rrC1 equal f_redRadC1/v_ncoke
|
||||
|
||||
compute layerRadC2 coke reduce sum f_LayerRelRadii[2]
|
||||
fix redRadC2 coke ave/time 1 1 1 c_layerRadC2
|
||||
variable rrC2 equal f_redRadC2/v_ncoke
|
||||
|
||||
variable reductionDegree atom 1./9.*f_fracRed[1]+2./9.*f_fracRed[2]+6./9.*f_fracRed[3]
|
||||
|
||||
|
||||
fix printRelRadiiOre all print ${Nmonitor} "${time} ${rrO1} ${rrO2} ${rrO3} ${rrO4}" append ${dumppath}/post/relRadiiOre.txt title "#time relRadO_1 relRadO_2 relRadO_3 relRadO_4"
|
||||
fix printRelRadiiCoke all print ${Nmonitor} "${time} ${rrC1} ${rrC2}" append ${dumppath}/post/relRadiiCoke.txt title "#time relRadC_1 relRadC_2"
|
||||
|
||||
# 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 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
|
||||
thermo 100
|
||||
thermo_modify lost ignore norm no
|
||||
compute_modify thermo_temp dynamic yes
|
||||
|
||||
|
||||
#dump dmp7 ore custom ${Ndump} ${dumppath}/post/dump*.ore id type radius x y z mass density f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4] f_LayerMasses[1] f_LayerMasses[2] f_LayerMasses[3] f_LayerMasses[4] f_fp[0] f_ftca[0]
|
||||
#dump dmp8 coke custom ${Ndump} ${dumppath}/post/dump*.coke id type radius x y z mass density f_fp[0] f_ftca[0]
|
||||
|
||||
#dump dmp9 fineore custom ${Ndump} ${dumppath}/post/dump*.fineore id type radius x y z mass density f_LayerRelRadii[1] f_LayerRelRadii[2] f_LayerRelRadii[3] f_LayerRelRadii[4] f_LayerMasses[1] f_LayerMasses[2] f_LayerMasses[3] f_LayerMasses[4] f_fp[0] f_ftca[0] f_reactionHeat[0]
|
||||
Binary file not shown.
@ -0,0 +1,10 @@
|
||||
Test case for ore reduction in a small cylinder
|
||||
|
||||
1) fill a cylinder with sperical particles (parDEMrun.sh)
|
||||
2) create CFD mesh (mesh.sh)
|
||||
3) introduce gas flow to reduce ore (parCFDDEMrun.sh)
|
||||
|
||||
|
||||
Comments:
|
||||
|
||||
After particles have settled, they remain frozen to their positions.
|
||||
@ -0,0 +1,801 @@
|
||||
#time mass lmass1 lmass2 lmass3 lmass4 lmass_{tot}
|
||||
10 0.493434041863914 0.00205514989592862 0.00660589076126408 0.0424146682964906 0.442358332910231 0.493434041863914
|
||||
20 0.49245951911578 0.00205514989592862 0.00950580425245368 0.0612724294457584 0.419626135521639 0.492459519115779
|
||||
30 0.491533853747881 0.00205514989592862 0.0126014178803977 0.0780855697545337 0.398791716217021 0.491533853747881
|
||||
40 0.490649936786299 0.0020551948532625 0.0158465784505578 0.0932080011432098 0.379540162339268 0.490649936786298
|
||||
50 0.48980226557356 0.00205651826396607 0.0192057421566188 0.106898176148744 0.361641829004232 0.48980226557356
|
||||
60 0.488986100451775 0.00206271221043687 0.0226469919108037 0.1193514445094 0.344924951821135 0.488986100451775
|
||||
70 0.488197674510037 0.00207593629857599 0.0261471189292331 0.130724756372543 0.329249862909684 0.488197674510037
|
||||
80 0.487434187073883 0.00209646668971253 0.0296903546886684 0.141146545826279 0.314500819869223 0.487434187073883
|
||||
90 0.486693464306797 0.00212374561405316 0.0332654287398584 0.150722982002698 0.300581307950187 0.486693464306797
|
||||
100 0.485973426076247 0.00215857104586299 0.0368608469436771 0.159541381689852 0.287412626396855 0.485973426076247
|
||||
110 0.485272086812471 0.00220272771429409 0.0404648940412023 0.167674765311812 0.274929699745162 0.485272086812471
|
||||
120 0.484587966175177 0.0022564027638926 0.0440703459997776 0.17518692735743 0.263074290054077 0.484587966175177
|
||||
130 0.483919896378486 0.0023191383558737 0.0476722750226317 0.182133238609454 0.251795244390526 0.483919896378486
|
||||
140 0.483266855899709 0.00239048097776341 0.0512666856408547 0.188561778584629 0.241047910696461 0.483266855899709
|
||||
150 0.482627938946594 0.00246979716551087 0.054851016201419 0.194514894919727 0.230792230659937 0.482627938946594
|
||||
160 0.482002380434478 0.00255647440403806 0.0584230043346048 0.200029638095666 0.220993263600169 0.482002380434478
|
||||
170 0.481389418337021 0.00265033116508886 0.0619800798950085 0.205138601749306 0.211620405527618 0.481389418337021
|
||||
180 0.480788333898775 0.00275131939578165 0.0655198900838632 0.209870766780944 0.202646357638186 0.480788333898775
|
||||
190 0.480198511647239 0.00285932489569668 0.0690403766546251 0.214251998485013 0.194046811611904 0.480198511647239
|
||||
200 0.479619438333945 0.00297398109578522 0.0725403101796411 0.218305671485241 0.185799475573278 0.479619438333945
|
||||
210 0.479050676506718 0.00309483779406924 0.0760188054338961 0.222052897735656 0.177884135543097 0.479050676506718
|
||||
220 0.478491751332173 0.00322175272220709 0.0794747520097794 0.225512766928222 0.170282479671965 0.478491751332173
|
||||
230 0.477942206696138 0.00335473745553208 0.0829068287565418 0.228702592806263 0.162978047677801 0.477942206696138
|
||||
240 0.477401606492865 0.003493776230328 0.0863141075730743 0.231638183425298 0.155955539264164 0.477401606492865
|
||||
250 0.476869543823169 0.00363889601366458 0.0896956740540771 0.234333926011858 0.149201047743569 0.476869543823169
|
||||
260 0.476345684967191 0.0037899617860583 0.0930509463489371 0.236803032879381 0.142701743952814 0.476345684967191
|
||||
270 0.475829704693984 0.00394691925711777 0.096379276971826 0.239057678159082 0.136445830305958 0.475829704693984
|
||||
280 0.47532121341203 0.00411004085599288 0.0996796423624596 0.241109047141644 0.130422483051934 0.47532121341203
|
||||
290 0.474819850835155 0.00427959933027457 0.102951012887075 0.242967492620945 0.12462174599686 0.474819850835155
|
||||
300 0.474325375881525 0.00445542624577653 0.106193199873481 0.244642694439716 0.119034055322551 0.474325375881525
|
||||
310 0.47383759065302 0.00463723078969086 0.109406320375662 0.246143619467683 0.113650420019984 0.47383759065302
|
||||
320 0.473356238918391 0.00482505187395772 0.112589940927639 0.247478596093821 0.108462650022973 0.473356238918391
|
||||
330 0.472881039426415 0.00501906803280593 0.115743491145508 0.248655381549098 0.103463098699003 0.472881039426415
|
||||
340 0.472411808038782 0.00521911981755022 0.118866980830564 0.249681256348929 0.0986444510417402 0.472411808038783
|
||||
350 0.471948356462611 0.00542512585890607 0.121960308936916 0.250563007389633 0.0939999142771559 0.471948356462611
|
||||
360 0.471490488214803 0.00563706022735517 0.125023369804103 0.251306983962763 0.0895230742205818 0.471490488214803
|
||||
370 0.471038022703491 0.00585485525544834 0.128056203834834 0.251919123848715 0.0852078397644945 0.471038022703491
|
||||
380 0.470590819410091 0.00607834976417837 0.131058963010939 0.252405008515594 0.0810484981193789 0.470590819410091
|
||||
390 0.470148730547618 0.00630742959128246 0.134031789968648 0.252769875152351 0.0770396358353364 0.470148730547618
|
||||
400 0.469711597245998 0.00654206809128671 0.136974686672405 0.253018646583775 0.0731761958985304 0.469711597245998
|
||||
410 0.46927927684286 0.00678220753438163 0.139887721012722 0.253155950717722 0.0694533975780347 0.469279276842859
|
||||
420 0.468851642570818 0.00702773843129973 0.142771104728015 0.253186135618856 0.0658666637926472 0.468851642570818
|
||||
430 0.468428551206996 0.00727864125908593 0.145624925798717 0.253113305495621 0.0624116786535722 0.468428551206996
|
||||
440 0.468009898011481 0.0075347862065936 0.148449419994691 0.252941346227296 0.0590843455829006 0.468009898011481
|
||||
450 0.467595557425268 0.00779615364866652 0.151244648878321 0.252673937229072 0.0558808176692078 0.467595557425268
|
||||
460 0.467185405363116 0.00806273740858438 0.154010671496229 0.252314556114025 0.0527974403442773 0.467185405363116
|
||||
470 0.466779350703461 0.00833441891072187 0.156747743814719 0.251866489819119 0.0498306981589008 0.466779350703461
|
||||
480 0.466377292300248 0.00861112981130234 0.159456074605791 0.251332842750522 0.0469772451326322 0.466377292300248
|
||||
490 0.465979154314444 0.00889273031895409 0.162135972093201 0.25071655501638 0.0442338968859091 0.465979154314444
|
||||
500 0.46558484460185 0.00917916495190843 0.164787610395181 0.250020422216833 0.0415976470379266 0.465584844601849
|
||||
510 0.465194265681034 0.00947041297019442 0.167411122035056 0.249247096374911 0.0390656343008725 0.465194265681034
|
||||
520 0.464807335507448 0.00976641200777846 0.170006703394544 0.248399094675386 0.0366351254297398 0.464807335507447
|
||||
530 0.464423988030738 0.0100670570760835 0.17257461163946 0.247478805636361 0.0343035136788331 0.464423988030738
|
||||
540 0.464044171637892 0.0103721936647314 0.175115207352853 0.246488482580738 0.0320682880395706 0.464044171637892
|
||||
550 0.463667824660517 0.0106817205157546 0.177628769782874 0.245430269650257 0.0299270647116316 0.463667824660517
|
||||
560 0.463294874258445 0.0109955952793715 0.180115481358623 0.244306213673118 0.0278775839473317 0.463294874258445
|
||||
570 0.46292526378887 0.0113137278611229 0.182575595944244 0.243118256311083 0.0259176836724196 0.46292526378887
|
||||
580 0.462558925722485 0.0116360676373137 0.185009349372034 0.241868221498014 0.0240452872151237 0.462558925722485
|
||||
590 0.462195807551822 0.0119625273276453 0.187417010436372 0.240557848103057 0.0222584216847481 0.462195807551823
|
||||
600 0.461835831805752 0.0122931176765889 0.189798721943154 0.239188784334212 0.0205552078517971 0.461835831805752
|
||||
610 0.461478924428081 0.0126278527438526 0.19215460690415 0.237762602127157 0.0189338626529216 0.461478924428081
|
||||
620 0.461125027312407 0.012966701024727 0.194484845791486 0.236280793335352 0.0173926871608427 0.461125027312407
|
||||
630 0.460774103626442 0.013309551572722 0.196789766529137 0.23474474696556 0.0159300385590222 0.460774103626442
|
||||
640 0.460426098541204 0.0136563692598133 0.199069588534976 0.233155779576855 0.0145443611695601 0.460426098541204
|
||||
650 0.460080974273298 0.0140070661699124 0.20132460621603 0.231515129285039 0.013234172602316 0.460080974273298
|
||||
660 0.459738704749152 0.0143615227710436 0.203555152774049 0.22982396350446 0.0119980656995999 0.459738704749153
|
||||
670 0.459399266560965 0.0147196196602568 0.205761558440754 0.22808338046766 0.0108347079922936 0.459399266560964
|
||||
680 0.459062629169544 0.0150812644044363 0.207944139178809 0.226294395352593 0.00974283023370586 0.459062629169544
|
||||
690 0.458728759874213 0.0154463787916989 0.2101031980732 0.224457950522122 0.00872123248719193 0.458728759874213
|
||||
700 0.458397628687537 0.0158148863604327 0.212239024265503 0.22257493298302 0.00776878507858055 0.458397628687536
|
||||
710 0.458069199886172 0.0161867401711689 0.214351861743143 0.22064617408971 0.00688442388214942 0.458069199886172
|
||||
720 0.457743441357678 0.0165618873801316 0.216441961446955 0.218672451719277 0.00606714081131383 0.457743441357677
|
||||
730 0.457420310826453 0.0169403230634692 0.218509487651847 0.216654521186756 0.00531597892438057 0.457420310826453
|
||||
740 0.457099771379877 0.0173220177209467 0.220554662248781 0.214593102581762 0.00462998882838667 0.457099771379877
|
||||
750 0.456781795768988 0.0177069069204143 0.222577742072926 0.212488958159533 0.0040081886161141 0.456781795768988
|
||||
760 0.456466350988364 0.0180949429805701 0.224578911337735 0.210343042472018 0.00344945419804175 0.456466350988364
|
||||
770 0.456153409005604 0.0184860233022929 0.226558351805158 0.208156857062682 0.00295217683547067 0.456153409005604
|
||||
780 0.455842932417995 0.0188799618813823 0.228516070768917 0.205933774433245 0.00251312533445073 0.455842932417995
|
||||
790 0.455534871188981 0.0192764549621994 0.230451832873147 0.203679778564795 0.00212680478883927 0.455534871188981
|
||||
800 0.455229177971096 0.0196752323121334 0.232365456188619 0.201400187151041 0.00178830231930146 0.455229177971096
|
||||
810 0.454925797579346 0.0200760014344116 0.234256675881942 0.199101088326314 0.00149203193667854 0.454925797579346
|
||||
820 0.454624686568758 0.0204785741312468 0.236125437947495 0.196786254018512 0.00123442047150409 0.454624686568758
|
||||
830 0.454325805727925 0.0208828155802221 0.2379717787189 0.194458378905758 0.00101283252304556 0.454325805727925
|
||||
840 0.454029119016436 0.0212885571973765 0.239795696991469 0.192120606400478 0.000824258427111858 0.454029119016436
|
||||
850 0.453734577280887 0.0216956787724586 0.241597174720594 0.189775912368609 0.000665811419225057 0.453734577280887
|
||||
860 0.453442131318399 0.0221040191394666 0.243376120161609 0.187428020814199 0.000533971203124198 0.453442131318399
|
||||
870 0.453151721146744 0.0225134320582616 0.245132338566594 0.185081113969924 0.000424836551964929 0.453151721146744
|
||||
880 0.452863280163598 0.0229238340411649 0.246865653774602 0.182738725432085 0.000335066915745356 0.452863280163598
|
||||
890 0.452576760940244 0.0233351020733977 0.24857601732094 0.180403913860032 0.000261727685874043 0.452576760940244
|
||||
900 0.452292123544132 0.0237471303217468 0.250263460302844 0.178079050983261 0.000202481936279968 0.452292123544132
|
||||
910 0.452009330359872 0.024159819492436 0.251928030320685 0.175766306944584 0.000155173602166692 0.452009330359872
|
||||
920 0.4517283429453 0.0245730759549963 0.253569772675869 0.173467812725022 0.000117681589412949 0.451728342945301
|
||||
930 0.451449125684644 0.0249868221074997 0.255188764980451 0.171185314082932 8.82245137614295e-05 0.451449125684644
|
||||
940 0.45117164690618 0.025400993425205 0.256785123406956 0.168920155783409 6.5374290610093e-05 0.45117164690618
|
||||
950 0.450895877035965 0.025815530354557 0.258358975742025 0.166673522497199 4.78484421841628e-05 0.450895877035965
|
||||
960 0.450621788863497 0.0262303860101279 0.259910477073756 0.164446285052249 3.46407273632038e-05 0.450621788863497
|
||||
970 0.450349354494856 0.0266455111020029 0.261439768429187 0.162239409324068 2.46656395986892e-05 0.450349354494856
|
||||
980 0.450078549463427 0.0270608557511036 0.262947005692316 0.160053721275328 1.69667446792725e-05 0.450078549463427
|
||||
990 0.449809351910175 0.027476401059116 0.264432397663347 0.157889415688769 1.11374989427591e-05 0.449809351910175
|
||||
1000 0.449541744782578 0.0278921303729358 0.265896180373471 0.155746432579906 7.00145626489185e-06 0.449541744782578
|
||||
1010 0.449275709388453 0.0283080190255109 0.26733856262303 0.153624946817497 4.18092241537987e-06 0.449275709388453
|
||||
1020 0.449011225516672 0.0287240412285322 0.268759734686927 0.151525250831332 2.19876988086144e-06 0.449011225516671
|
||||
1030 0.448748277359769 0.0291401805601684 0.270159920857384 0.149447283528172 8.92414044661887e-07 0.448748277359769
|
||||
1040 0.448486850095637 0.029556425635008 0.271539355379212 0.147390860594663 2.0848675499274e-07 0.448486850095637
|
||||
1050 0.448226927744799 0.0299727640758805 0.272898256832537 0.145355894032952 1.28034289795913e-08 0.448226927744799
|
||||
1060 0.447968492526893 0.0303891783190104 0.27423681770961 0.143342496498273 0 0.447968492526893
|
||||
1070 0.447711527805661 0.030805657332688 0.275555242360397 0.141350628112575 0 0.44771152780566
|
||||
1080 0.447456020700517 0.031222190223118 0.276853754594765 0.139380075882634 0 0.447456020700517
|
||||
1090 0.447201957716433 0.0316387699994147 0.278132572174002 0.137430615543017 0 0.447201957716433
|
||||
1100 0.446949325553582 0.0320553895415832 0.279391910894222 0.135502025117776 0 0.446949325553582
|
||||
1110 0.44669811119467 0.0324720419999529 0.280631981354591 0.133594087840126 0 0.44669811119467
|
||||
1120 0.446448301840927 0.0328887203572838 0.281852992072866 0.131706589410778 0 0.446448301840927
|
||||
1130 0.446199884981126 0.0333054176346315 0.283055147496335 0.12983931985016 0 0.446199884981126
|
||||
1140 0.445952848158957 0.0337221275532908 0.284238647729277 0.127992072876389 0 0.445952848158957
|
||||
1150 0.445707178718869 0.0341388450263185 0.285403689427406 0.126164644265145 0 0.445707178718869
|
||||
1160 0.445462864434994 0.0345555644626056 0.286550466052608 0.12435683391978 0 0.445462864434994
|
||||
1170 0.445219893085252 0.0349722808424508 0.287679168075415 0.122568444167386 0 0.445219893085252
|
||||
1180 0.444978252547252 0.0353889896909243 0.288789981820261 0.120799281036066 0 0.444978252547252
|
||||
1190 0.44473793122682 0.0358056854515252 0.289883092061517 0.119049153713779 0 0.44473793122682
|
||||
1200 0.444498917684854 0.0362223627893441 0.290958680344431 0.117317874551079 0 0.444498917684854
|
||||
1210 0.444261200486622 0.0366390170424661 0.292016924692286 0.11560525875187 0 0.444261200486622
|
||||
1220 0.444024768320574 0.0370556438561171 0.293057999888579 0.113911124575878 0 0.444024768320574
|
||||
1230 0.443789609977493 0.0374722388817395 0.294082079306311 0.112235291789443 0 0.443789609977493
|
||||
1240 0.443555714658514 0.0378887975142117 0.295089332110967 0.110577585033335 0 0.443555714658514
|
||||
1250 0.443323071473129 0.0383053157334633 0.296079925691969 0.108937830047696 0 0.443323071473129
|
||||
1260 0.443091669584401 0.0387217900477453 0.297054023998551 0.107315855538105 0 0.443091669584401
|
||||
1270 0.442861498710032 0.039138215855163 0.298011789221018 0.105711493633851 0 0.442861498710032
|
||||
1280 0.442632548330566 0.0395545897401235 0.298953380640011 0.104124577950431 0 0.442632548330566
|
||||
1290 0.44240480831711 0.0399709076633717 0.299878955484169 0.102554945169569 0 0.44240480831711
|
||||
1300 0.442178268132117 0.0403871673476347 0.300788667420356 0.101002433364127 0 0.442178268132117
|
||||
1310 0.441952917855166 0.0408033651930379 0.301682668608025 0.0994668840541024 0 0.441952917855166
|
||||
1320 0.441728747598508 0.0412194980328398 0.302561108537709 0.097948141027959 0 0.441728747598508
|
||||
1330 0.441505747393386 0.0416355632984749 0.303424134713163 0.0964460493817483 0 0.441505747393386
|
||||
1340 0.441283907680808 0.0420515577999008 0.304271892316228 0.0949604575646792 0 0.441283907680808
|
||||
1350 0.441063218725744 0.0424674792199563 0.305104524408512 0.0934912150972752 0 0.441063218725744
|
||||
1360 0.440843671441774 0.0428833238767342 0.305922172336961 0.0920381752280791 0 0.440843671441774
|
||||
1370 0.440625256291648 0.0432990898112337 0.306724974655946 0.0906011918244685 0 0.440625256291648
|
||||
1380 0.440407964082792 0.0437147744579734 0.307513068397574 0.0891801212272437 0 0.440407964082791
|
||||
1390 0.440191785694851 0.0441303754953021 0.308286588366425 0.0877748218331236 0 0.44019178569485
|
||||
1400 0.439976712018596 0.0445458909815439 0.309045667261952 0.0863851537751001 0 0.439976712018596
|
||||
1410 0.439762733150696 0.0449613219390581 0.309790433228675 0.0850109779829631 0 0.439762733150696
|
||||
1420 0.439549841306868 0.0453766631431862 0.310521018011359 0.0836521601523232 0 0.439549841306868
|
||||
1430 0.439338027884298 0.0457919123530222 0.311237549027423 0.0823085665038531 0 0.439338027884298
|
||||
1440 0.439127284465263 0.0462070670464795 0.311940152647833 0.0809800647709502 0 0.439127284465263
|
||||
1450 0.438917602581491 0.0466221254036239 0.312628952311436 0.0796665248664306 0 0.438917602581491
|
||||
1460 0.438708973787445 0.047037085740404 0.313304070451614 0.0783678175954265 0 0.438708973787445
|
||||
1470 0.438501389781314 0.0474519464756628 0.313965627036148 0.0770838162695032 0 0.438501389781314
|
||||
1480 0.438294842225595 0.0478667063106975 0.314613741060234 0.0758143948546633 0 0.438294842225595
|
||||
1490 0.438089323235239 0.0482813631979058 0.315248529277766 0.074559430759567 0 0.438089323235239
|
||||
1500 0.437884824884554 0.0486959153414541 0.315870107694613 0.0733188018484869 0 0.437884824884554
|
||||
1510 0.437681339066706 0.0491103618918775 0.316478589887018 0.0720923872878106 0 0.437681339066706
|
||||
1520 0.437478857926716 0.0495247015518446 0.317074088343141 0.0708800680317303 0 0.437478857926716
|
||||
1530 0.437277373561739 0.0499389335297977 0.317656713593349 0.0696817264385921 0 0.437277373561739
|
||||
1540 0.437076878271982 0.0503530567224081 0.318226575037053 0.0684972465125212 0 0.437076878271982
|
||||
1550 0.436877364629293 0.0507670694908236 0.318783781170107 0.067326513968362 0 0.436877364629293
|
||||
1560 0.436678825016339 0.0511809710801582 0.319328438488115 0.066169415448066 0 0.436678825016339
|
||||
1570 0.436481251990724 0.0515947605645857 0.319860652009781 0.0650258394163574 0 0.436481251990724
|
||||
1580 0.436284638262365 0.0520084368018543 0.320380525815688 0.0638956756448236 0 0.436284638262366
|
||||
1590 0.436088976527405 0.0524219989996761 0.32088816237117 0.0627788151565587 0 0.436088976527405
|
||||
1600 0.435894259566785 0.0528354463335026 0.321383663165252 0.0616751500680304 0 0.435894259566785
|
||||
1610 0.435700480404787 0.0532487775788413 0.321867128466958 0.0605845743589876 0 0.435700480404787
|
||||
1620 0.435507631927548 0.0536619921930414 0.322338656897719 0.0595069828367873 0 0.435507631927548
|
||||
1630 0.435315707288905 0.0540750890725945 0.322798346358803 0.0584422718575077 0 0.435315707288905
|
||||
1640 0.435124699525443 0.0544880677104212 0.323246293261923 0.0573903385530991 0 0.435124699525443
|
||||
1650 0.434934601875669 0.0549009272854378 0.323682592914532 0.0563510816756992 0 0.434934601875669
|
||||
1660 0.434745407744692 0.0553136667064795 0.324107339769578 0.0553244012686353 0 0.434745407744693
|
||||
1670 0.434557110017809 0.0557262868402225 0.324520625207785 0.054310197969802 0 0.434557110017809
|
||||
1680 0.434369702403047 0.056138785997629 0.324922542678042 0.0533083737273764 0 0.434369702403047
|
||||
1690 0.434183176371174 0.0565511701809148 0.32531317618982 0.0523188300004393 0 0.434183176371174
|
||||
1700 0.433997527391334 0.0569634325555453 0.325692621945697 0.0513414728900919 0 0.433997527391334
|
||||
1710 0.433812749451541 0.0573755714112144 0.32606096973242 0.0503762083079063 0 0.43381274945154
|
||||
1720 0.433628836364504 0.0577875858477775 0.326418307529282 0.0494229429874447 0 0.433628836364504
|
||||
1730 0.433445781923718 0.0581994751475407 0.326764722635658 0.048481584140519 0 0.433445781923718
|
||||
1740 0.433263579886109 0.0586112390431259 0.327100300526251 0.0475520403167324 0 0.433263579886109
|
||||
1750 0.433082224084558 0.0590228771889987 0.327425126054812 0.0466342208407476 0 0.433082224084559
|
||||
1760 0.432901708432193 0.0594343891539749 0.327739283434843 0.0457280358433745 0 0.432901708432193
|
||||
1770 0.432722027119511 0.0598457739462246 0.328042856028359 0.0448333971449276 0 0.432722027119511
|
||||
1780 0.432543174409041 0.0602570305393375 0.328335926405016 0.0439502174646884 0 0.432543174409041
|
||||
1790 0.432365144606493 0.0606681579351838 0.328618576407274 0.0430784102640361 0 0.432365144606494
|
||||
1800 0.432187932101585 0.0610791550123792 0.328890887380175 0.0422178897090301 0 0.432187932101585
|
||||
1810 0.432011531209591 0.0614900211651439 0.32915293901837 0.0413685710260772 0 0.432011531209591
|
||||
1820 0.431835935594994 0.0619007580086444 0.329404808354879 0.0405303692314713 0 0.431835935594994
|
||||
1830 0.431661139868729 0.0623113643762961 0.329646573987351 0.0397032015050819 0 0.431661139868729
|
||||
1840 0.431487138399997 0.0627218400540741 0.329878312868175 0.0388869854777481 0 0.431487138399997
|
||||
1850 0.431313925653634 0.0631321846454273 0.330100101600704 0.0380816394075022 0 0.431313925653634
|
||||
1860 0.431141496344355 0.0635423970924848 0.330312016825194 0.0372870824266764 0 0.431141496344355
|
||||
1870 0.430969843028812 0.063952483823277 0.330514125732306 0.0365032334732286 0 0.430969843028812
|
||||
1880 0.430798961837283 0.0643624393859325 0.330706508503889 0.0357300139474611 0 0.430798961837283
|
||||
1890 0.430628847958432 0.0647722618067328 0.330889240108421 0.0349673460432787 0 0.430628847958432
|
||||
1900 0.430459495915328 0.0651819514145952 0.331062392631777 0.034215151868956 0 0.430459495915328
|
||||
1910 0.430290900853697 0.0655915065984029 0.331226039755012 0.033473354500282 0 0.430290900853697
|
||||
1920 0.430123057763199 0.0660009265257128 0.331380253254683 0.0327418779828034 0 0.430123057763199
|
||||
1930 0.429955961533266 0.0664102108269541 0.331525103875672 0.0320206468306393 0 0.429955961533266
|
||||
1940 0.429789606938626 0.0668193595133978 0.331660661948321 0.0313095854769068 0 0.429789606938626
|
||||
1950 0.429623988885572 0.0672283725117359 0.331786996465892 0.0306086199079445 0 0.429623988885572
|
||||
1960 0.429459103198818 0.0676372467833704 0.331904179308414 0.0299176771070336 0 0.429459103198817
|
||||
1970 0.429294944202043 0.0680459843449261 0.332012276563166 0.0292366832939512 0 0.429294944202043
|
||||
1980 0.429131506707577 0.0684545857985981 0.332111355013262 0.0285655658957173 0 0.429131506707577
|
||||
1990 0.428968787144363 0.0688630465265093 0.332201486525793 0.0279042540920603 0 0.428968787144363
|
||||
2000 0.428806779452763 0.0692713702492374 0.33228273358358 0.0272526756199457 0 0.428806779452763
|
||||
2010 0.42864547958226 0.0696795542481675 0.332355164739362 0.0266107605947307 0 0.42864547958226
|
||||
2020 0.428484883245173 0.0700875966770316 0.332418847251794 0.0259784393163471 0 0.428484883245173
|
||||
2030 0.428324984958427 0.0704954996769224 0.332473843964522 0.0253556413169829 0 0.428324984958427
|
||||
2040 0.428165780984084 0.0709032598415467 0.332520222758224 0.0247422983843134 0 0.428165780984084
|
||||
2050 0.428007266980916 0.0713108759054014 0.332558048624407 0.0241383424511075 0 0.428007266980916
|
||||
2060 0.427849437651967 0.0717183498262314 0.332587382843043 0.0235437049826919 0 0.427849437651967
|
||||
2070 0.427692288529164 0.0721256809215442 0.33260828911432 0.0229583184933 0 0.427692288529164
|
||||
2080 0.427535815412773 0.0725328677090307 0.332620831639374 0.0223821160643689 0 0.427535815412773
|
||||
2090 0.427380013667213 0.072939910302017 0.332625072279799 0.0218150310853972 0 0.427380013667214
|
||||
2100 0.427224879289353 0.0733468068680316 0.332621074390164 0.0212569980311571 0 0.427224879289353
|
||||
2110 0.427070407575773 0.0737535578816815 0.33260889888685 0.0207079508072421 0 0.427070407575773
|
||||
2120 0.42691659425079 0.0741601625828246 0.332588607278011 0.0201678243899541 0 0.42691659425079
|
||||
2130 0.42676343505056 0.0745666202018714 0.332560260960893 0.0196365538877951 0 0.42676343505056
|
||||
2140 0.426610925641703 0.0749729303410152 0.332523920356386 0.0191140749443021 0 0.426610925641703
|
||||
2150 0.42645906264756 0.0753790894986689 0.332479648968799 0.0186003241800921 0 0.42645906264756
|
||||
2160 0.426307842073313 0.0757850962918447 0.332427507727769 0.0180952380536992 0 0.426307842073313
|
||||
2170 0.426157259027571 0.0761909524306535 0.332367553724285 0.0175987528726324 0 0.426157259027571
|
||||
2180 0.426007309472312 0.0765966568453997 0.332299846857847 0.0171108057690658 0 0.426007309472312
|
||||
2190 0.42585798928455 0.0770022088066036 0.332224446421899 0.0166313340560473 0 0.42585798928455
|
||||
2200 0.425709294539754 0.0774076070470801 0.332141411803367 0.0161602756893071 0 0.425709294539754
|
||||
2210 0.425561221174843 0.0778128507763156 0.33205080180278 0.0156975685957471 0 0.425561221174843
|
||||
2220 0.425413765287399 0.078217938751326 0.331952675384646 0.0152431511514277 0 0.4254137652874
|
||||
2230 0.425266922922551 0.0786228699411273 0.33184709112924 0.0147969618521827 0 0.425266922922551
|
||||
2240 0.425120690244954 0.079027643010299 0.331734107573185 0.0143589396614701 0 0.425120690244954
|
||||
2250 0.42497506337792 0.0794322568238804 0.331613782778593 0.0139290237754467 0 0.42497506337792
|
||||
2260 0.424830038552614 0.0798367099232419 0.331486175020169 0.0135071536092034 0 0.424830038552614
|
||||
2270 0.424685611915437 0.0802410011502227 0.331351342168376 0.0130932685968386 0 0.424685611915437
|
||||
2280 0.424541779626556 0.0806451293724847 0.331209341775315 0.0126873084787565 0 0.424541779626556
|
||||
2290 0.424398537529335 0.0810490945417432 0.331060230082451 0.0122892129051406 0 0.424398537529335
|
||||
2300 0.42425588218223 0.0814528942864597 0.33090406565324 0.0118989222425304 0 0.42425588218223
|
||||
2310 0.424113810080523 0.081856526452301 0.330740906784438 0.0115163768437834 0 0.424113810080523
|
||||
2320 0.423972317409755 0.0822599899320771 0.330570810556958 0.01114151692072 0 0.423972317409755
|
||||
2330 0.42383140056019 0.082663282982101 0.330393834569157 0.0107742830089325 0 0.42383140056019
|
||||
2340 0.423691055935424 0.0830664038297741 0.330210036388609 0.0104146157170405 0 0.423691055935424
|
||||
2350 0.423551279980931 0.0834693505713213 0.330019473693713 0.0100624557158967 0 0.423551279980931
|
||||
2360 0.42341206894783 0.083872121955803 0.329822203421704 0.00971774357032254 0 0.423412068947829
|
||||
2370 0.423273419671638 0.0842747147988716 0.329618284586271 0.00938042028649523 0 0.423273419671638
|
||||
2380 0.423135328423159 0.0846771278000725 0.329407774114258 0.00905042650882896 0 0.423135328423159
|
||||
2390 0.422997791734857 0.0850793587699245 0.329190729982432 0.00872770298250057 0 0.422997791734857
|
||||
2400 0.422860806027855 0.0854814058716393 0.328967209852539 0.00841219030367653 0 0.422860806027855
|
||||
2410 0.422724368029307 0.0858832662796478 0.328737272351844 0.00810382939781493 0 0.422724368029307
|
||||
2420 0.422588474105974 0.0862849383062322 0.328500975107319 0.0078025606924231 0 0.422588474105974
|
||||
2430 0.422453120740654 0.0866864198714675 0.328258376201023 0.00750832466816376 0 0.422453120740654
|
||||
2440 0.42231830454242 0.0870877084507361 0.328009534302234 0.0072210617894502 0 0.422318304542421
|
||||
2450 0.422184021131636 0.0874888048245077 0.327754504391891 0.00694071191523704 0 0.422184021131636
|
||||
2460 0.42205026764135 0.0878897045990795 0.327493347698696 0.00666721534357484 0 0.42205026764135
|
||||
2470 0.421917041122617 0.0882904036235475 0.327226125306553 0.00640051219251615 0 0.421917041122617
|
||||
2480 0.42178433815371 0.0886908992841975 0.326952896758109 0.00614054211140294 0 0.42178433815371
|
||||
2490 0.421652155373411 0.0890911886778746 0.326673722270312 0.00588724442522354 0 0.42165215537341
|
||||
2500 0.421520489532435 0.0894912684669756 0.32638866278781 0.00564055827764975 0 0.421520489532435
|
||||
2510 0.42138933727417 0.0898911356329241 0.326097779059446 0.00540042258180018 0 0.42138933727417
|
||||
2520 0.421258695199144 0.0902907871537431 0.325801132428879 0.00516677561652211 0 0.421258695199144
|
||||
2530 0.421128560099796 0.0906902193245424 0.32549878516869 0.00493955560656335 0 0.421128560099796
|
||||
2540 0.420998928570326 0.0910894289800615 0.325190799445076 0.00471870014518778 0 0.420998928570326
|
||||
2550 0.420869797335073 0.0914884124176672 0.324877238436346 0.00450414648105936 0 0.420869797335073
|
||||
2560 0.420741163099989 0.0918871658584301 0.324558165961107 0.00429583128045217 0 0.420741163099989
|
||||
2570 0.420613022800991 0.0922856846108202 0.324233647436127 0.00409369075404368 0 0.420613022800991
|
||||
2580 0.420485373046611 0.0926839649755392 0.323903747582575 0.00389766048849754 0 0.420485373046611
|
||||
2590 0.420358210393066 0.0930820032656116 0.323568531762922 0.00370767536453182 0 0.420358210393066
|
||||
2600 0.420231531625953 0.0934797948516396 0.323228067099134 0.00352366967517947 0 0.420231531625953
|
||||
2610 0.420105333427577 0.0938773352817929 0.32288242119093 0.00334557695485399 0 0.420105333427577
|
||||
2620 0.41997961262306 0.0942746194368676 0.322531663148258 0.00317333003793499 0 0.41997961262306
|
||||
2630 0.419854366007485 0.0946716420896807 0.322175863041357 0.00300686087644755 0 0.419854366007485
|
||||
2640 0.419729590105994 0.0950683986826737 0.321815091131083 0.00284610029223673 0 0.419729590105994
|
||||
2650 0.419605281856922 0.0954648830657219 0.321449420307144 0.00269097848405556 0 0.419605281856922
|
||||
2660 0.419481437817991 0.0958610901095617 0.321078923333343 0.00254142437508629 0 0.419481437817991
|
||||
2670 0.419358054796153 0.0962570135850767 0.320703675275689 0.00239736593538706 0 0.419358054796153
|
||||
2680 0.419235129405077 0.096652647623024 0.320323751952091 0.00225872982996163 0 0.419235129405077
|
||||
2690 0.41911265828296 0.0970479859915189 0.319939230716029 0.00212544157541195 0 0.41911265828296
|
||||
2700 0.418990636854966 0.0974430263003636 0.319550185493289 0.00199742506131319 0 0.418990636854966
|
||||
2710 0.418869062999154 0.0978377573566711 0.319156702463406 0.00187460317907635 0 0.418869062999154
|
||||
2720 0.418747933330545 0.0982321719327782 0.318758864416092 0.00175689698167464 0 0.418747933330545
|
||||
2730 0.418627244486197 0.0986262623464631 0.318356756183247 0.00164422595648694 0 0.418627244486197
|
||||
2740 0.418506993106155 0.0990200204717015 0.317950464840184 0.00153650779426979 0 0.418506993106155
|
||||
2750 0.418387175781639 0.0994134379001719 0.31754007957259 0.0014336583088766 0 0.418387175781638
|
||||
2760 0.418267789179338 0.0998065054770269 0.317125692390877 0.00133559131143398 0 0.418267789179338
|
||||
2770 0.418148829706359 0.100199214424041 0.316707396864283 0.00124221841803558 0 0.418148829706359
|
||||
2780 0.41803029388017 0.100591555033082 0.31628528986043 0.001153448986657 0 0.418030293880169
|
||||
2790 0.417912178219791 0.100983516985691 0.315859471366581 0.00106918986751843 0 0.41791217821979
|
||||
2800 0.417794479158506 0.10137508962158 0.315430044283882 0.000989345253044353 0 0.417794479158506
|
||||
2810 0.417677193275433 0.101766261103678 0.314997115604415 0.00091381656734108 0 0.417677193275433
|
||||
2820 0.417560316783121 0.102157020119208 0.314560794560318 0.000842502103594881 0 0.417560316783121
|
||||
2830 0.417443845818292 0.102547354838851 0.31412119406315 0.000775296916290516 0 0.417443845818292
|
||||
2840 0.417327776907825 0.102937251239884 0.313678433055426 0.000712092612514803 0 0.417327776907825
|
||||
2850 0.417212106305962 0.103326695330026 0.313232633995625 0.000652776980310418 0 0.417212106305962
|
||||
2860 0.4170968302574 0.103715672178741 0.312783924298786 0.000597233779872221 0 0.4170968302574
|
||||
2870 0.41698194510007 0.104104165483933 0.312332437181691 0.000545342434446151 0 0.41698194510007
|
||||
2880 0.416867446661543 0.104492159562787 0.311878309549417 0.000496977549338103 0 0.416867446661543
|
||||
2890 0.416753330904613 0.104879637044229 0.311421685172748 0.000452008687636219 0 0.416753330904613
|
||||
2900 0.416639593643823 0.105266579748876 0.310962713978911 0.00041029991603621 0 0.416639593643823
|
||||
2910 0.416526230722781 0.105652967959807 0.310501553415289 0.00037170934768457 0 0.416526230722781
|
||||
2920 0.416413238004462 0.1060387803274 0.310038369069811 0.000336088607251058 0 0.416413238004462
|
||||
2930 0.416300611072742 0.106423994776527 0.309573334023322 0.000303282272892662 0 0.416300611072742
|
||||
2940 0.416188345581263 0.10680858715462 0.309106631233296 0.000273127193346156 0 0.416188345581263
|
||||
2950 0.416076437050413 0.107192532227104 0.3086384508217 0.000245454001609847 0 0.416076437050413
|
||||
2960 0.415964880269379 0.107575808593785 0.308168972966002 0.000220098709592496 0 0.415964880269379
|
||||
2970 0.415853671235404 0.107958393256096 0.307698370168411 0.000196907810897002 0 0.415853671235404
|
||||
2980 0.415742805516276 0.10834026616472 0.307226805042486 0.000175734309069943 0 0.415742805516276
|
||||
2990 0.415632279188415 0.108721409738738 0.306754420945037 0.000156448504640386 0 0.415632279188415
|
||||
3000 0.415522088366759 0.109101809348928 0.306281345515535 0.000138933502295226 0 0.415522088366759
|
||||
3010 0.415412229380687 0.109481450634594 0.305807702121707 0.00012307662438593 0 0.415412229380687
|
||||
3020 0.415302698419087 0.109860319800219 0.305333613129143 0.000108765489725053 0 0.415302698419087
|
||||
3030 0.415193491775501 0.110238402563937 0.304859201956855 9.58872547082936e-05 0 0.415193491775501
|
||||
3040 0.415084605585286 0.11061568600342 0.304384587107219 8.43324746464453e-05 0 0.415084605585285
|
||||
3050 0.414976036347436 0.110992156627023 0.303909885202889 7.39945175237477e-05 0 0.414976036347436
|
||||
3060 0.414867780344681 0.111367802017269 0.30343521024491 6.47680825019088e-05 0 0.414867780344682
|
||||
3070 0.414759834021692 0.111742610084615 0.302960672295048 5.65516420290873e-05 0 0.414759834021692
|
||||
3080 0.414652193814798 0.112116570146868 0.3024863743354 4.92493325302874e-05 0 0.414652193814798
|
||||
3090 0.414544856497393 0.112489672600995 0.302012409145575 4.27747508232629e-05 0 0.414544856497394
|
||||
3100 0.414437819026588 0.112861909327502 0.301538859688684 3.70500104018557e-05 0 0.414437819026588
|
||||
3110 0.414331078335843 0.113233273044106 0.30106580455323 3.20007385075534e-05 0 0.414331078335843
|
||||
3120 0.414224631587837 0.113603757659687 0.300593313430419 2.75604977315707e-05 0 0.414224631587837
|
||||
3130 0.414118475927152 0.113973357963244 0.300121451511615 2.3666452292529e-05 0 0.414118475927152
|
||||
3140 0.41401260869169 0.114342069334729 0.299650278519481 2.02608374794543e-05 0 0.41401260869169
|
||||
3150 0.413907027230923 0.114709888007416 0.299179849502336 1.72897211717241e-05 0 0.413907027230923
|
||||
3160 0.413801728965665 0.115076811131047 0.298710213846041 1.4703988577817e-05 0 0.413801728965665
|
||||
3170 0.413696711544375 0.115442836375544 0.29824141526885 1.24598999817276e-05 0 0.413696711544375
|
||||
3180 0.413591972674414 0.11580796238768 0.297773491551356 1.05187353782176e-05 0 0.413591972674414
|
||||
3190 0.413487510215722 0.116172188220154 0.29730647643584 8.8455597275486e-06 0 0.413487510215722
|
||||
3200 0.413383322003293 0.1165355136341 0.296840400329149 7.40804004285742e-06 0 0.413383322003293
|
||||
3210 0.413279406115707 0.116897938174463 0.296375291416396 6.17652484853402e-06 0 0.413279406115707
|
||||
3220 0.413175760826759 0.117259461672901 0.295911173830975 5.12532288262703e-06 0 0.413175760826759
|
||||
3230 0.413072383751012 0.117620086751079 0.295448066253704 4.23074622947696e-06 0 0.413072383751012
|
||||
3240 0.412969272946707 0.1179798148063 0.294985987612721 3.4705276871283e-06 0 0.412969272946707
|
||||
3250 0.412866426585174 0.118338647530256 0.294524953875561 2.8251793568099e-06 0 0.412866426585174
|
||||
3260 0.412763842680857 0.118696587668486 0.294064977635029 2.27737734176786e-06 0 0.412763842680857
|
||||
3270 0.412661519676269 0.119053637454174 0.293606068356779 1.8138653153911e-06 0 0.412661519676269
|
||||
3280 0.412559455932274 0.119409800233014 0.293148230903881 1.42479537889906e-06 0 0.412559455932274
|
||||
3290 0.412457649796786 0.11976507992958 0.292691467362068 1.10250513716795e-06 0 0.412457649796786
|
||||
3300 0.412356100280447 0.120119478200945 0.292235782531923 8.39547578937094e-07 0 0.412356100280447
|
||||
3310 0.412254805309597 0.120473000385631 0.291781176888916 6.28035050211622e-07 0 0.412254805309597
|
||||
3320 0.41215376332527 0.120825649984162 0.2913276533963 4.59944808339263e-07 0 0.41215376332527
|
||||
3330 0.412052972779358 0.121177430587842 0.290875214477922 3.27713594273931e-07 0 0.412052972779358
|
||||
3340 0.411952432288565 0.121528345285595 0.290423862810768 2.24192201907604e-07 0 0.411952432288565
|
||||
3350 0.411852140646789 0.121878397180341 0.289973598688058 1.44778390113921e-07 0 0.411852140646789
|
||||
3360 0.41175209643872 0.122227590477876 0.289524419545441 8.64154038898031e-08 0 0.41175209643872
|
||||
3370 0.411652298138544 0.122575929789613 0.289076322207755 4.61411759940828e-08 0 0.411652298138544
|
||||
3380 0.411552744404293 0.122923419042589 0.288629304544892 2.08168115689801e-08 0 0.411552744404293
|
||||
3390 0.41145343386228 0.123270062227585 0.288183364526026 7.10866898319633e-09 0 0.41145343386228
|
||||
3400 0.411354365221949 0.123615862992888 0.287738500762853 1.46620809244937e-09 0 0.411354365221949
|
||||
3410 0.411255537135353 0.123960825137831 0.287294711888615 1.08907034074468e-10 0 0.411255537135353
|
||||
3420 0.411156948251663 0.124304952635798 0.286851995615865 0 0 0.411156948251663
|
||||
3430 0.411058597348801 0.124648249477496 0.286410347871304 0 0 0.411058597348801
|
||||
3440 0.410960483188279 0.124990719962796 0.285969763225483 0 0 0.410960483188279
|
||||
3450 0.410862604566028 0.125332368294543 0.285530236271485 0 0 0.410862604566028
|
||||
3460 0.410764960251866 0.12567319876672 0.285091761485146 0 0 0.410764960251867
|
||||
3470 0.410667549048317 0.12601321555916 0.284654333489157 0 0 0.410667549048316
|
||||
3480 0.410570369764851 0.126352422827433 0.284217946937418 0 0 0.410570369764851
|
||||
3490 0.410473421223892 0.126690824681918 0.283782596541973 0 0 0.410473421223892
|
||||
3500 0.410376702305437 0.12702842503201 0.283348277273427 0 0 0.410376702305437
|
||||
3510 0.410280211924028 0.127365227666534 0.282914984257494 0 0 0.410280211924028
|
||||
3520 0.410183948908189 0.127701236674567 0.282482712233623 0 0 0.410183948908189
|
||||
3530 0.410087912109298 0.128036456065409 0.282051456043889 0 0 0.410087912109298
|
||||
3540 0.409992100315292 0.128370890069801 0.281621210245491 0 0 0.409992100315292
|
||||
3550 0.409896512457908 0.128704542416547 0.281191970041361 0 0 0.409896512457908
|
||||
3560 0.409801147391164 0.129037417105738 0.280763730285426 0 0 0.409801147391164
|
||||
3570 0.409706003939975 0.129369518239035 0.28033648570094 0 0 0.409706003939975
|
||||
3580 0.409611081104734 0.129700849305603 0.279910231799132 0 0 0.409611081104734
|
||||
3590 0.409516377766489 0.130031414211181 0.279484963555308 0 0 0.409516377766489
|
||||
3600 0.40942189281011 0.130361216848156 0.279060675961954 0 0 0.40942189281011
|
||||
3610 0.409327625257845 0.130690260629406 0.27863736462844 0 0 0.409327625257845
|
||||
3620 0.409233573928233 0.131018549678858 0.278215024249375 0 0 0.409233573928233
|
||||
3630 0.409139737889252 0.131346087249764 0.277793650639487 0 0 0.409139737889252
|
||||
3640 0.409046116278579 0.131672876352085 0.277373239926493 0 0 0.409046116278579
|
||||
3650 0.408952707984844 0.131998920865093 0.276953787119751 0 0 0.408952707984844
|
||||
3660 0.40885951201131 0.132324224267932 0.276535287743377 0 0 0.40885951201131
|
||||
3670 0.40876652733686 0.132648790124837 0.276117737212024 0 0 0.40876652733686
|
||||
3680 0.40867375288783 0.132972622183472 0.275701130704358 0 0 0.40867375288783
|
||||
3690 0.408581187678597 0.13329572388418 0.275285463794417 0 0 0.408581187678597
|
||||
3700 0.408488830787303 0.133618098444733 0.27487073234257 0 0 0.408488830787303
|
||||
3710 0.408396681106564 0.133939749730492 0.274456931376072 0 0 0.408396681106564
|
||||
3720 0.408304737607914 0.134260681331351 0.274044056276563 0 0 0.408304737607914
|
||||
3730 0.408212999452728 0.13458089617455 0.273632103278178 0 0 0.408212999452728
|
||||
3740 0.4081214656493 0.134900397721674 0.273221067927626 0 0 0.4081214656493
|
||||
3750 0.408030135270756 0.135219189208004 0.272810946062753 0 0 0.408030135270756
|
||||
3760 0.407939007370519 0.135537273937596 0.272401733432924 0 0 0.407939007370519
|
||||
3770 0.407848080881393 0.13585465563553 0.271993425245863 0 0 0.407848080881393
|
||||
3780 0.407757355088113 0.13617133679846 0.271586018289653 0 0 0.407757355088113
|
||||
3790 0.407666828991041 0.136487320915655 0.271179508075386 0 0 0.407666828991041
|
||||
3800 0.407576501655722 0.136802611248854 0.270773890406869 0 0 0.407576501655722
|
||||
3810 0.407486372210025 0.137117210842264 0.270369161367762 0 0 0.407486372210025
|
||||
3820 0.407396439746566 0.137431122863131 0.269965316883435 0 0 0.407396439746566
|
||||
3830 0.407306703320967 0.137744350607839 0.269562352713129 0 0 0.407306703320967
|
||||
3840 0.407217162129167 0.138056896882986 0.269160265246181 0 0 0.407217162129167
|
||||
3850 0.407127815248932 0.13836876490766 0.268759050341272 0 0 0.407127815248932
|
||||
3860 0.40703866174143 0.138679957958875 0.268358703782555 0 0 0.40703866174143
|
||||
3870 0.40694970083782 0.138990478720305 0.267959222117515 0 0 0.40694970083782
|
||||
3880 0.406860931515711 0.139300330760634 0.267560600755077 0 0 0.406860931515711
|
||||
3890 0.406772353011059 0.139609516746774 0.267162836264285 0 0 0.406772353011059
|
||||
3900 0.406683964538258 0.139918039420915 0.266765925117343 0 0 0.406683964538258
|
||||
3910 0.406595765253824 0.140225901727261 0.266369863526563 0 0 0.406595765253824
|
||||
3920 0.40650775427838 0.140533106735302 0.265974647543078 0 0 0.40650775427838
|
||||
3930 0.406419930857464 0.140839657078514 0.26558027377895 0 0 0.406419930857464
|
||||
3940 0.40633229413952 0.14114555572928 0.265186738410239 0 0 0.40633229413952
|
||||
3950 0.406244843253225 0.141450805728975 0.26479403752425 0 0 0.406244843253225
|
||||
3960 0.406157577540896 0.141755409373269 0.264402168167627 0 0 0.406157577540896
|
||||
3970 0.406070496074552 0.1420593699013 0.264011126173252 0 0 0.406070496074552
|
||||
3980 0.405983598058491 0.142362690090492 0.263620907967999 0 0 0.405983598058491
|
||||
3990 0.405896882751168 0.142665372529236 0.263231510221932 0 0 0.405896882751168
|
||||
4000 0.405810349414148 0.14296741979506 0.262842929619087 0 0 0.405810349414148
|
||||
4010 0.405723997289505 0.143268834533524 0.262455162755981 0 0 0.405723997289505
|
||||
4020 0.405637825618767 0.1435696193921 0.262068206226667 0 0 0.405637825618767
|
||||
4030 0.405551833525401 0.143869777430353 0.261682056095048 0 0 0.405551833525401
|
||||
4040 0.405466020328269 0.144169311025817 0.261296709302452 0 0 0.405466020328269
|
||||
4050 0.405380385161499 0.144468223200845 0.260912161960654 0 0 0.405380385161499
|
||||
4060 0.405294927396031 0.144766516151188 0.260528411244844 0 0 0.405294927396032
|
||||
4070 0.405209646212375 0.145064192737305 0.260145453475071 0 0 0.405209646212375
|
||||
4080 0.405124540794272 0.14536125580837 0.259763284985902 0 0 0.405124540794272
|
||||
4090 0.405039610506577 0.145657707581377 0.2593819029252 0 0 0.405039610506577
|
||||
4100 0.404954854637477 0.145953550540929 0.259001304096549 0 0 0.404954854637477
|
||||
4110 0.404870272498613 0.146248787089769 0.258621485408843 0 0 0.404870272498613
|
||||
4120 0.404785863321037 0.146543419911928 0.258242443409109 0 0 0.404785863321037
|
||||
4130 0.404701626382158 0.146837451529632 0.257864174852526 0 0 0.404701626382159
|
||||
4140 0.404617560910813 0.147130884634655 0.257486676276158 0 0 0.404617560910813
|
||||
4150 0.40453366630791 0.147423721318143 0.257109944989767 0 0 0.404533666307911
|
||||
4160 0.404449941836091 0.147715964153869 0.256733977682222 0 0 0.404449941836091
|
||||
4170 0.40436638682749 0.148007615473035 0.256358771354455 0 0 0.404366386827489
|
||||
4180 0.404283000577268 0.148298677735903 0.255984322841364 0 0 0.404283000577268
|
||||
4190 0.404199782422479 0.148589153256509 0.25561062916597 0 0 0.404199782422479
|
||||
4200 0.404116731586843 0.148879044744482 0.255237686842361 0 0 0.404116731586843
|
||||
4210 0.404033847470961 0.149168354292045 0.254865493178916 0 0 0.404033847470961
|
||||
4220 0.403951129309858 0.149457084569361 0.254494044740497 0 0 0.403951129309858
|
||||
4230 0.403868576562623 0.1497452374645 0.254123339098123 0 0 0.403868576562623
|
||||
4240 0.403786188525562 0.150032815433728 0.253753373091834 0 0 0.403786188525562
|
||||
4250 0.403703964585247 0.150319820618234 0.253384143967013 0 0 0.403703964585247
|
||||
4260 0.40362190399064 0.150606255639537 0.253015648351102 0 0 0.40362190399064
|
||||
4270 0.40354000612796 0.150892122640057 0.252647883487903 0 0 0.40354000612796
|
||||
4280 0.403458270358216 0.151177423850218 0.252280846507999 0 0 0.403458270358216
|
||||
4290 0.403376696016535 0.151462161590775 0.25191453442576 0 0 0.403376696016535
|
||||
4300 0.40329528251513 0.151746337913421 0.251548944601708 0 0 0.40329528251513
|
||||
4310 0.403214029141097 0.152029955306561 0.251184073834536 0 0 0.403214029141097
|
||||
4320 0.40313293533254 0.152313015731516 0.250819919601024 0 0 0.40313293533254
|
||||
4330 0.403052000475751 0.152595521330452 0.2504564791453 0 0 0.403052000475752
|
||||
4340 0.402971223885218 0.15287747449617 0.250093749389048 0 0 0.402971223885218
|
||||
4350 0.402890604996796 0.153158877197836 0.24973172779896 0 0 0.402890604996796
|
||||
4360 0.402810143186352 0.153439731614 0.249370411572352 0 0 0.402810143186352
|
||||
4370 0.402729837787951 0.153720040069137 0.249009797718813 0 0 0.402729837787951
|
||||
4380 0.402649688246253 0.153999804501669 0.248649883744584 0 0 0.402649688246253
|
||||
4390 0.402569693992548 0.154279026896705 0.248290667095843 0 0 0.402569693992548
|
||||
4400 0.402489854416061 0.15455770938617 0.24793214502989 0 0 0.40248985441606
|
||||
4410 0.402410168901717 0.154835854116999 0.247574314784717 0 0 0.402410168901717
|
||||
4420 0.402330636906389 0.155113462984996 0.247217173921393 0 0 0.402330636906389
|
||||
4430 0.402251257801679 0.155390538183605 0.246860719618074 0 0 0.402251257801679
|
||||
4440 0.402172030969962 0.155667081868661 0.246504949101301 0 0 0.402172030969962
|
||||
4450 0.402092955901812 0.155943095818335 0.246149860083477 0 0 0.402092955901812
|
||||
4460 0.402014031980614 0.156218582184939 0.245795449795675 0 0 0.402014031980614
|
||||
4470 0.401935258693735 0.156493542757834 0.245441715935902 0 0 0.401935258693735
|
||||
4480 0.401856635454439 0.156767979585042 0.245088655869397 0 0 0.401856635454439
|
||||
4490 0.401778161665874 0.157041894749892 0.244736266915982 0 0 0.401778161665873
|
||||
4500 0.401699836796523 0.157315290107656 0.244384546688867 0 0 0.401699836796523
|
||||
4510 0.401621660254169 0.157588167725494 0.244033492528675 0 0 0.401621660254169
|
||||
4520 0.401543631623859 0.157860529051806 0.243683102572053 0 0 0.401543631623859
|
||||
4530 0.401465750249056 0.158132376378257 0.243333373870799 0 0 0.401465750249055
|
||||
4540 0.401388015578373 0.158403711629476 0.242984303948897 0 0 0.401388015578373
|
||||
4550 0.401310427034272 0.158674536821382 0.24263589021289 0 0 0.401310427034272
|
||||
4560 0.401232984124127 0.158944853673501 0.242288130450626 0 0 0.401232984124127
|
||||
4570 0.401155686309227 0.159214664066222 0.241941022243004 0 0 0.401155686309226
|
||||
4580 0.40107853306554 0.159483969828686 0.241594563236854 0 0 0.40107853306554
|
||||
4590 0.401001523870355 0.159752772785439 0.241248751084916 0 0 0.401001523870355
|
||||
4600 0.400924658187847 0.160021074806797 0.24090358338105 0 0 0.400924658187847
|
||||
4610 0.400847935508098 0.160288877672635 0.240559057835463 0 0 0.400847935508098
|
||||
4620 0.400771355296941 0.16055618324749 0.24021517204945 0 0 0.400771355296941
|
||||
4630 0.400694917050819 0.160822993289033 0.239871923761786 0 0 0.400694917050819
|
||||
4640 0.400618620262536 0.161089309567656 0.239529310694879 0 0 0.400618620262536
|
||||
4650 0.400542464399523 0.161355133942302 0.239187330457221 0 0 0.400542464399523
|
||||
4660 0.400466448954751 0.161620468182772 0.238845980771978 0 0 0.400466448954751
|
||||
4670 0.4003905734525 0.161885313949584 0.238505259502915 0 0 0.4003905734525
|
||||
4680 0.400314837342236 0.162149673164389 0.238165164177847 0 0 0.400314837342236
|
||||
4690 0.400239240198398 0.162413547312626 0.237825692885772 0 0 0.400239240198398
|
||||
4700 0.400163781502947 0.162676938202524 0.237486843300423 0 0 0.400163781502947
|
||||
4710 0.40008846075303 0.162939847589303 0.237148613163727 0 0 0.40008846075303
|
||||
4720 0.400013277469465 0.163202277145569 0.236811000323896 0 0 0.400013277469465
|
||||
4730 0.399938231150961 0.16346422862109 0.236474002529871 0 0 0.399938231150961
|
||||
4740 0.399863321345915 0.1637257035922 0.236137617753715 0 0 0.399863321345915
|
||||
4750 0.399788547581553 0.163986703709132 0.235801843872421 0 0 0.399788547581553
|
||||
4760 0.39971390935602 0.164247230723627 0.235466678632393 0 0 0.39971390935602
|
||||
4770 0.399639406216988 0.164507286214548 0.235132120002439 0 0 0.399639406216988
|
||||
4780 0.399565037736338 0.164766871676257 0.234798166060081 0 0 0.399565037736338
|
||||
4790 0.399490803371791 0.165025989001599 0.234464814370191 0 0 0.399490803371791
|
||||
4800 0.399416702677384 0.165284639747218 0.234132062930166 0 0 0.399416702677384
|
||||
4810 0.399342735225351 0.165542825406246 0.233799909819104 0 0 0.39934273522535
|
||||
4820 0.399268900538388 0.165800547644726 0.233468352893662 0 0 0.399268900538388
|
||||
4830 0.399195198166963 0.166057808031765 0.233137390135198 0 0 0.399195198166963
|
||||
4840 0.399121627674827 0.166314608090102 0.232807019584725 0 0 0.399121627674827
|
||||
4850 0.399048188583675 0.166570949489275 0.2324772390944 0 0 0.399048188583675
|
||||
4860 0.39897488043949 0.166826833814045 0.232148046625445 0 0 0.39897488043949
|
||||
4870 0.398901702702714 0.167082262947758 0.231819439754955 0 0 0.398901702702714
|
||||
4880 0.398828654971603 0.167337238292712 0.231491416678891 0 0 0.398828654971603
|
||||
4890 0.398755736768907 0.167591761514758 0.231163975254149 0 0 0.398755736768907
|
||||
4900 0.39868294764855 0.167845834170947 0.230837113477603 0 0 0.39868294764855
|
||||
4910 0.39861028722809 0.168099457596198 0.230510829631891 0 0 0.39861028722809
|
||||
4920 0.398537755052297 0.168352633379507 0.23018512167279 0 0 0.398537755052297
|
||||
4930 0.398465350664043 0.168605363116491 0.229859987547552 0 0 0.398465350664043
|
||||
4940 0.398393073678066 0.168857648151924 0.229535425526142 0 0 0.398393073678066
|
||||
4950 0.398320923703047 0.169109489851714 0.229211433851333 0 0 0.398320923703047
|
||||
4960 0.398248900242342 0.169360889949413 0.228888010292929 0 0 0.398248900242342
|
||||
4970 0.398177002914658 0.169611849775935 0.228565153138723 0 0 0.398177002914658
|
||||
4980 0.398105231277343 0.16986237087637 0.228242860400973 0 0 0.398105231277343
|
||||
4990 0.398033584926335 0.170112454661106 0.227921130265229 0 0 0.398033584926335
|
||||
5000 0.397962063414562 0.170362102690671 0.227599960723891 0 0 0.397962063414562
|
||||
5010 0.397890666311057 0.170611316469357 0.2272793498417 0 0 0.397890666311057
|
||||
5020 0.397819393278718 0.170860097173838 0.226959296104879 0 0 0.397819393278718
|
||||
5030 0.397748243889958 0.171108446296614 0.226639797593344 0 0 0.397748243889958
|
||||
5040 0.39767721776999 0.171356365145893 0.226320852624097 0 0 0.39767721776999
|
||||
5050 0.397606314500268 0.171603855182623 0.226002459317644 0 0 0.397606314500268
|
||||
5060 0.397535533660397 0.171850917874202 0.225684615786195 0 0 0.397535533660397
|
||||
5070 0.397464874856412 0.17209755459578 0.225367320260632 0 0 0.397464874856412
|
||||
5080 0.39739433769164 0.172343766731959 0.225050570959681 0 0 0.39739433769164
|
||||
5090 0.39732392173867 0.172589555774624 0.224734365964046 0 0 0.39732392173867
|
||||
5100 0.397253626653402 0.17283492292487 0.224418703728533 0 0 0.397253626653402
|
||||
5110 0.39718345205372 0.173079869516485 0.224103582537235 0 0 0.39718345205372
|
||||
5120 0.397113397509107 0.173324397052201 0.223789000456905 0 0 0.397113397509107
|
||||
5130 0.397043462653375 0.173568506810208 0.223474955843166 0 0 0.397043462653375
|
||||
5140 0.396973647126326 0.17381220004779 0.223161447078536 0 0 0.396973647126325
|
||||
5150 0.396903950569434 0.174055478016389 0.222848472553044 0 0 0.396903950569434
|
||||
5160 0.396834372575026 0.174298342139004 0.222536030436022 0 0 0.396834372575026
|
||||
5170 0.396764912790511 0.174540793646361 0.222224119144151 0 0 0.396764912790511
|
||||
5180 0.396695570757795 0.174782834137458 0.221912736620337 0 0 0.396695570757794
|
||||
5190 0.396626346161653 0.17502446471259 0.221601881449063 0 0 0.396626346161653
|
||||
5200 0.396557238658374 0.175265686571493 0.221291552086882 0 0 0.396557238658374
|
||||
5210 0.396488247863349 0.17550650105666 0.22098174680669 0 0 0.396488247863349
|
||||
5220 0.396419373386175 0.175746909530803 0.220672463855372 0 0 0.396419373386175
|
||||
5230 0.396350614873062 0.175986913228832 0.22036370164423 0 0 0.396350614873062
|
||||
5240 0.396281971965888 0.176226513400789 0.220055458565099 0 0 0.396281971965888
|
||||
5250 0.39621344432538 0.176465711230914 0.219747733094465 0 0 0.39621344432538
|
||||
5260 0.396145031556698 0.176704508097411 0.219440523459287 0 0 0.396145031556697
|
||||
5270 0.396076733311134 0.176942905217448 0.219133828093686 0 0 0.396076733311134
|
||||
5280 0.396008549263866 0.177180903724834 0.218827645539032 0 0 0.396008549263866
|
||||
5290 0.395940479055354 0.177418504874553 0.218521974180801 0 0 0.395940479055354
|
||||
5300 0.395872522295659 0.177655710027698 0.218216812267961 0 0 0.395872522295659
|
||||
5310 0.395804678670421 0.177892520281552 0.21791215838887 0 0 0.395804678670422
|
||||
5320 0.395736947826335 0.17812893686934 0.217608010956995 0 0 0.395736947826335
|
||||
5330 0.395669329416645 0.178364961001418 0.217304368415227 0 0 0.395669329416645
|
||||
5340 0.395601823060774 0.178600594006206 0.217001229054569 0 0 0.395601823060774
|
||||
5350 0.395534428464588 0.178835836910391 0.216698591554197 0 0 0.395534428464588
|
||||
5360 0.395467145289271 0.179070690896613 0.216396454392657 0 0 0.395467145289271
|
||||
5370 0.395399973142679 0.179305157333667 0.216094815809012 0 0 0.395399973142679
|
||||
5380 0.395332911770061 0.17953923711077 0.215793674659291 0 0 0.39533291177006
|
||||
5390 0.39526596077442 0.179772931613649 0.215493029160771 0 0 0.39526596077442
|
||||
5400 0.395199119840774 0.180006241941758 0.215192877899016 0 0 0.395199119840774
|
||||
5410 0.395132388676125 0.180239169117821 0.214893219558305 0 0 0.395132388676125
|
||||
5420 0.395065766910348 0.180471714433768 0.214594052476581 0 0 0.395065766910348
|
||||
5430 0.394999254186499 0.180703879135521 0.214295375050978 0 0 0.394999254186499
|
||||
5440 0.394932850220389 0.18093566421505 0.213997186005339 0 0 0.39493285022039
|
||||
5450 0.394866554672414 0.18116707085776 0.213699483814654 0 0 0.394866554672413
|
||||
5460 0.394800367166611 0.181398100375951 0.21340226679066 0 0 0.394800367166611
|
||||
5470 0.394734287471307 0.181628753578291 0.213105533893016 0 0 0.394734287471307
|
||||
5480 0.394668315221202 0.18185903173987 0.212809283481332 0 0 0.394668315221202
|
||||
5490 0.394602450097176 0.182088935974588 0.212513514122588 0 0 0.394602450097176
|
||||
5500 0.394536691799849 0.182318467327439 0.212218224472409 0 0 0.394536691799849
|
||||
5510 0.394471039992296 0.182547626974465 0.21192341301783 0 0 0.394471039992296
|
||||
5520 0.394405494327496 0.182776416126959 0.211629078200537 0 0 0.394405494327496
|
||||
5530 0.394340054577637 0.183004835580105 0.211335218997532 0 0 0.394340054577637
|
||||
5540 0.394274720331875 0.183232886767965 0.21104183356391 0 0 0.394274720331875
|
||||
5550 0.39420949136437 0.183460570478841 0.21074892088553 0 0 0.39420949136437
|
||||
5560 0.39414436734137 0.183687887877705 0.210456479463665 0 0 0.39414436734137
|
||||
5570 0.394079347915998 0.183914840175344 0.210164507740654 0 0 0.394079347915998
|
||||
5580 0.394014432836715 0.184141428249758 0.209873004586957 0 0 0.394014432836715
|
||||
5590 0.39394962177644 0.184367653242635 0.209581968533805 0 0 0.39394962177644
|
||||
5600 0.39388491442163 0.18459351624841 0.209291398173219 0 0 0.39388491442163
|
||||
5610 0.393820310512973 0.184819018172214 0.209001292340759 0 0 0.393820310512973
|
||||
5620 0.393755809692289 0.18504416026429 0.208711649427999 0 0 0.393755809692289
|
||||
5630 0.393691411711266 0.185268943391372 0.208422468319894 0 0 0.393691411711266
|
||||
5640 0.393627116272729 0.185493368590766 0.208133747681963 0 0 0.393627116272729
|
||||
5650 0.393562922829392 0.185717437772788 0.207845485056604 0 0 0.393562922829392
|
||||
5660 0.393498831196481 0.185941151582395 0.207557679614086 0 0 0.393498831196481
|
||||
5670 0.39343484119512 0.186164510643961 0.20727033055116 0 0 0.393434841195121
|
||||
5680 0.393370952571083 0.186387515844871 0.206983436726212 0 0 0.393370952571083
|
||||
5690 0.393307164995937 0.186610168331529 0.206696996664408 0 0 0.393307164995937
|
||||
5700 0.393243478219079 0.186832468978672 0.206411009240406 0 0 0.393243478219079
|
||||
5710 0.3931798919299 0.187054418870488 0.206125473059412 0 0 0.3931798919299
|
||||
5720 0.393116405828643 0.187276019053295 0.205840386775347 0 0 0.393116405828643
|
||||
5730 0.393053019679551 0.187497270350002 0.20555574932955 0 0 0.393053019679552
|
||||
5740 0.392989733138769 0.187718173960853 0.205271559177916 0 0 0.392989733138769
|
||||
5750 0.392926546007643 0.187938730579248 0.204987815428396 0 0 0.392926546007644
|
||||
5760 0.392863457944229 0.188158941398759 0.20470451654547 0 0 0.392863457944229
|
||||
5770 0.392800468692103 0.188378807314442 0.20442166137766 0 0 0.392800468692103
|
||||
5780 0.392737577954531 0.188598329362051 0.20413924859248 0 0 0.392737577954531
|
||||
5790 0.392674785487063 0.18881750839485 0.203857277092213 0 0 0.392674785487063
|
||||
5800 0.39261209101454 0.189036345373288 0.203575745641252 0 0 0.39261209101454
|
||||
5810 0.392549494174671 0.189254841561955 0.203294652612716 0 0 0.392549494174671
|
||||
5820 0.392486994774977 0.189472997632707 0.20301399714227 0 0 0.392486994774976
|
||||
5830 0.392424592500414 0.189690814685208 0.202733777815206 0 0 0.392424592500413
|
||||
5840 0.392362287142814 0.189908293446077 0.202453993696737 0 0 0.392362287142814
|
||||
5850 0.392300078399363 0.190125434972297 0.202174643427066 0 0 0.392300078399363
|
||||
5860 0.392237966032626 0.190342240092646 0.20189572593998 0 0 0.392237966032626
|
||||
5870 0.39217594975606 0.190558709807303 0.201617239948757 0 0 0.39217594975606
|
||||
5880 0.39211402934586 0.190774844897475 0.201339184448384 0 0 0.392114029345859
|
||||
5890 0.392052204498126 0.190990646423925 0.201061558074201 0 0 0.392052204498126
|
||||
5900 0.391990474958666 0.191206115273924 0.200784359684742 0 0 0.391990474958666
|
||||
5910 0.391928840411132 0.191421252551697 0.200507587859435 0 0 0.391928840411131
|
||||
5920 0.391867300676651 0.191636058881595 0.200231241795056 0 0 0.391867300676651
|
||||
5930 0.391805855491565 0.191850535183935 0.199955320307631 0 0 0.391805855491565
|
||||
5940 0.391744504569936 0.192064682456789 0.199679822113147 0 0 0.391744504569936
|
||||
5950 0.391683247691406 0.192278501469319 0.199404746222087 0 0 0.391683247691406
|
||||
5960 0.391622084538459 0.192491993329832 0.199130091208627 0 0 0.391622084538459
|
||||
5970 0.391561014935471 0.192705158651345 0.198855856284125 0 0 0.391561014935471
|
||||
5980 0.391500038569907 0.192917998524769 0.198582040045138 0 0 0.391500038569907
|
||||
5990 0.39143915522067 0.193130513721853 0.198308641498817 0 0 0.39143915522067
|
||||
6000 0.391378364631441 0.193342705137284 0.198035659494157 0 0 0.391378364631441
|
||||
6010 0.391317666563995 0.193554573602596 0.1977630929614 0 0 0.391317666563995
|
||||
6020 0.391257060728792 0.193766120128434 0.197490940600358 0 0 0.391257060728792
|
||||
6030 0.391196546917149 0.193977345443219 0.19721920147393 0 0 0.391196546917149
|
||||
6040 0.391136124842301 0.194188250547904 0.196947874294398 0 0 0.391136124842301
|
||||
6050 0.391075794310259 0.194398836119623 0.196676958190636 0 0 0.391075794310259
|
||||
6060 0.391015555046624 0.194609103116164 0.19640645193046 0 0 0.391015555046624
|
||||
6070 0.39095540681571 0.194819052360205 0.196136354455505 0 0 0.39095540681571
|
||||
6080 0.390895349388086 0.195028684652577 0.195866664735509 0 0 0.390895349388086
|
||||
6090 0.390835382515949 0.195238000858246 0.195597381657702 0 0 0.390835382515949
|
||||
6100 0.390775505965668 0.195447001792708 0.19532850417296 0 0 0.390775505965667
|
||||
6110 0.390715719500725 0.195655688281535 0.195060031219189 0 0 0.390715719500724
|
||||
6120 0.390656022874485 0.195864061185611 0.194791961688875 0 0 0.390656022874486
|
||||
6130 0.390596415855715 0.196072121312072 0.194524294543642 0 0 0.390596415855715
|
||||
6140 0.390536898208406 0.196279869484706 0.1942570287237 0 0 0.390536898208406
|
||||
6150 0.390477469706518 0.196487306492514 0.193990163214005 0 0 0.390477469706518
|
||||
6160 0.390418130139052 0.196694433071997 0.193723697067055 0 0 0.390418130139052
|
||||
6170 0.390358879221931 0.196901250214731 0.193457629007201 0 0 0.390358879221931
|
||||
6180 0.390299716775208 0.19710775854883 0.193191958226378 0 0 0.390299716775208
|
||||
6190 0.390240642527452 0.197313959021733 0.192926683505719 0 0 0.390240642527451
|
||||
6200 0.390181656281297 0.197519852322346 0.19266180395895 0 0 0.390181656281297
|
||||
6210 0.390122757802992 0.197725439266587 0.192397318536406 0 0 0.390122757802992
|
||||
6220 0.390063946881306 0.197930720591767 0.19213322628954 0 0 0.390063946881307
|
||||
6230 0.390005223249445 0.198135697229136 0.191869526020309 0 0 0.390005223249445
|
||||
6240 0.389946586726839 0.198340369808982 0.191606216917857 0 0 0.389946586726839
|
||||
6250 0.389888037039062 0.198544739289195 0.191343297749867 0 0 0.389888037039062
|
||||
6260 0.389829574019786 0.198748806250345 0.191080767769442 0 0 0.389829574019786
|
||||
6270 0.389771197427905 0.198952571534028 0.190818625893877 0 0 0.389771197427904
|
||||
6280 0.389712907026917 0.199156035965754 0.190556871061163 0 0 0.389712907026916
|
||||
6290 0.389654702614412 0.199359200252042 0.190295502362369 0 0 0.389654702614411
|
||||
6300 0.389596583963656 0.199562065184314 0.190034518779342 0 0 0.389596583963656
|
||||
6310 0.389538550885011 0.199764631424513 0.189773919460499 0 0 0.389538550885011
|
||||
6320 0.389480603111094 0.199966899905943 0.189513703205151 0 0 0.389480603111094
|
||||
6330 0.389422740465482 0.200168871244414 0.189253869221069 0 0 0.389422740465482
|
||||
6340 0.38936496271008 0.200370546271005 0.188994416439075 0 0 0.38936496271008
|
||||
6350 0.389307269633809 0.200571925722492 0.188735343911317 0 0 0.389307269633809
|
||||
6360 0.389249661021055 0.200773010351484 0.188476650669571 0 0 0.389249661021055
|
||||
6370 0.389192136684513 0.200973800811774 0.188218335872739 0 0 0.389192136684514
|
||||
6380 0.389134696370393 0.201174297989225 0.187960398381168 0 0 0.389134696370393
|
||||
6390 0.389077339891513 0.201374502537203 0.187702837354309 0 0 0.389077339891513
|
||||
6400 0.389020067056172 0.20157441512484 0.187445651931332 0 0 0.389020067056172
|
||||
6410 0.388962877629737 0.201774036571135 0.187188841058602 0 0 0.388962877629737
|
||||
6420 0.38890577142854 0.201973367517185 0.186932403911355 0 0 0.38890577142854
|
||||
6430 0.388848748220814 0.20217240877198 0.186676339448834 0 0 0.388848748220814
|
||||
6440 0.388791807820019 0.202371160986642 0.186420646833376 0 0 0.388791807820019
|
||||
6450 0.388734950002821 0.202569624940729 0.186165325062092 0 0 0.38873495000282
|
||||
6460 0.388678174566162 0.202767801343016 0.185910373223146 0 0 0.388678174566162
|
||||
6470 0.388621481363515 0.202965690704965 0.18565579065855 0 0 0.388621481363515
|
||||
6480 0.388564870214348 0.203163293656728 0.18540157655762 0 0 0.388564870214348
|
||||
6490 0.388508340831072 0.203360611202138 0.185147729628934 0 0 0.388508340831072
|
||||
6500 0.388451893041943 0.203557643940678 0.184894249101266 0 0 0.388451893041943
|
||||
6510 0.388395526615951 0.203754392678696 0.184641133937254 0 0 0.388395526615951
|
||||
6520 0.388339241370328 0.203950858054143 0.184388383316185 0 0 0.388339241370329
|
||||
6530 0.38828303708697 0.204147040828328 0.184135996258642 0 0 0.38828303708697
|
||||
6540 0.388226913589466 0.204342941617008 0.183883971972458 0 0 0.388226913589466
|
||||
6550 0.388170870673669 0.204538561132769 0.1836323095409 0 0 0.388170870673669
|
||||
6560 0.388114908134925 0.204733900089963 0.183381008044962 0 0 0.388114908134925
|
||||
6570 0.388059025802203 0.204928959085581 0.183130066716622 0 0 0.388059025802203
|
||||
6580 0.388003223462422 0.205123738863385 0.182879484599037 0 0 0.388003223462422
|
||||
6590 0.387947500920525 0.205318240104228 0.182629260816297 0 0 0.387947500920525
|
||||
6600 0.387891857996009 0.20551246343816 0.182379394557849 0 0 0.387891857996009
|
||||
6610 0.387836294470186 0.205706409628522 0.182129884841664 0 0 0.387836294470186
|
||||
6620 0.387780810170305 0.205900079278302 0.181880730892003 0 0 0.387780810170305
|
||||
6630 0.38772540492896 0.206093472971839 0.181631931957121 0 0 0.38772540492896
|
||||
6640 0.387670078507299 0.206286591542856 0.181383486964443 0 0 0.387670078507299
|
||||
6650 0.387614830765597 0.206479435479059 0.181135395286538 0 0 0.387614830765597
|
||||
6660 0.387559661495397 0.206672005508081 0.180887655987315 0 0 0.387559661495397
|
||||
6670 0.387504570484759 0.206864302369701 0.180640268115058 0 0 0.387504570484759
|
||||
6680 0.387449557600169 0.207056326529959 0.18039323107021 0 0 0.387449557600169
|
||||
6690 0.387394622608669 0.207248078801994 0.180146543806675 0 0 0.387394622608669
|
||||
6700 0.387339765354096 0.207439559730904 0.179900205623193 0 0 0.387339765354096
|
||||
6710 0.387284985646796 0.20763076997868 0.179654215668116 0 0 0.387284985646796
|
||||
6720 0.387230283288683 0.207821710236747 0.179408573051936 0 0 0.387230283288683
|
||||
6730 0.387175658105897 0.20801238111197 0.179163276993928 0 0 0.387175658105897
|
||||
6740 0.387121109912725 0.208202783252583 0.178918326660142 0 0 0.387121109912725
|
||||
6750 0.387066638536333 0.208392917261873 0.17867372127446 0 0 0.387066638536333
|
||||
6760 0.387012243803994 0.208582783742745 0.17842946006125 0 0 0.387012243803994
|
||||
6770 0.386957925541691 0.208772383302613 0.178185542239078 0 0 0.386957925541691
|
||||
6780 0.386903683528485 0.20896171671267 0.177941966815814 0 0 0.386903683528485
|
||||
6790 0.386849517623719 0.209150784463882 0.177698733159838 0 0 0.386849517623719
|
||||
6800 0.386795427644255 0.2093395871955 0.177455840448755 0 0 0.386795427644255
|
||||
6810 0.386741413406428 0.209528125548613 0.177213287857815 0 0 0.386741413406428
|
||||
6820 0.386687474736345 0.209716400130199 0.176971074606146 0 0 0.386687474736345
|
||||
6830 0.386633611444504 0.20990441160172 0.176729199842784 0 0 0.386633611444504
|
||||
6840 0.386579823351734 0.210092160588578 0.176487662763156 0 0 0.386579823351734
|
||||
6850 0.386526110328109 0.210279647544283 0.176246462783826 0 0 0.386526110328109
|
||||
6860 0.386472472154912 0.210466873232276 0.176005598922636 0 0 0.386472472154912
|
||||
6870 0.386418908718552 0.210653838049048 0.175765070669503 0 0 0.386418908718552
|
||||
6880 0.386365419841955 0.210840542612678 0.175524877229277 0 0 0.386365419841955
|
||||
6890 0.386312005314127 0.211026987659655 0.175285017654472 0 0 0.386312005314127
|
||||
6900 0.386258664978216 0.211213173737468 0.175045491240748 0 0 0.386258664978216
|
||||
6910 0.386205398656989 0.21139910146476 0.174806297192229 0 0 0.386205398656989
|
||||
6920 0.386152206166405 0.211584771483927 0.174567434682478 0 0 0.386152206166405
|
||||
6930 0.38609908735592 0.21177018432045 0.17432890303547 0 0 0.38609908735592
|
||||
6940 0.386046042173251 0.211955340156825 0.174090702016426 0 0 0.386046042173251
|
||||
6950 0.385993070353586 0.212140239917385 0.173852830436201 0 0 0.385993070353586
|
||||
6960 0.385940171748121 0.212324884121536 0.173615287626585 0 0 0.385940171748121
|
||||
6970 0.385887346165198 0.212509273438262 0.173378072726937 0 0 0.385887346165198
|
||||
6980 0.385834593451157 0.212693408403925 0.173141185047231 0 0 0.385834593451157
|
||||
6990 0.385781913408419 0.212877289708173 0.172904623700247 0 0 0.385781913408419
|
||||
7000 0.385729305925029 0.213060917741796 0.172668388183233 0 0 0.385729305925029
|
||||
7010 0.385676770798561 0.213244293211364 0.172432477587197 0 0 0.385676770798561
|
||||
7020 0.385624307872098 0.213427416664601 0.172196891207497 0 0 0.385624307872098
|
||||
7030 0.38557191698679 0.213610288655975 0.171961628330815 0 0 0.38557191698679
|
||||
7040 0.385519597960242 0.213792909822145 0.171726688138097 0 0 0.385519597960242
|
||||
7050 0.38546735069745 0.213975280494728 0.171492070202722 0 0 0.38546735069745
|
||||
7060 0.385415174936738 0.214157401587106 0.171257773349632 0 0 0.385415174936738
|
||||
7070 0.3853630705537 0.214339273533523 0.171023797020177 0 0 0.385363070553701
|
||||
7080 0.385311037426852 0.21452089675803 0.170790140668822 0 0 0.385311037426852
|
||||
7090 0.385259075463875 0.214702271582856 0.170556803881019 0 0 0.385259075463875
|
||||
7100 0.385207184420708 0.214883398859917 0.170323785560791 0 0 0.385207184420708
|
||||
7110 0.385155364115791 0.215064279222944 0.170091084892847 0 0 0.385155364115791
|
||||
7120 0.385103614461154 0.215244912979005 0.169858701482148 0 0 0.385103614461154
|
||||
7130 0.385051935242222 0.215425300877076 0.169626634365147 0 0 0.385051935242222
|
||||
7140 0.385000326302344 0.215605443463956 0.169394882838388 0 0 0.385000326302344
|
||||
7150 0.384948787508231 0.215785341204893 0.169163446303338 0 0 0.384948787508231
|
||||
7160 0.38489731871724 0.215964994597786 0.168932324119454 0 0 0.38489731871724
|
||||
7170 0.384845919725317 0.216144404354891 0.168701515370427 0 0 0.384845919725317
|
||||
7180 0.384794590397861 0.216323570946039 0.168471019451822 0 0 0.384794590397861
|
||||
7190 0.384743330648765 0.216502494671788 0.168240835976978 0 0 0.384743330648765
|
||||
7200 0.384692140264939 0.216681176275936 0.168010963989004 0 0 0.384692140264939
|
||||
7210 0.384641019094317 0.216859616289275 0.167781402805042 0 0 0.384641019094317
|
||||
7220 0.384589966980159 0.217037815258906 0.167552151721253 0 0 0.384589966980159
|
||||
7230 0.384538983883486 0.217215773320892 0.167323210562594 0 0 0.384538983883486
|
||||
7240 0.384488069549463 0.217393491364742 0.16709457818472 0 0 0.384488069549462
|
||||
7250 0.384437223610433 0.217570970673769 0.166866252936663 0 0 0.384437223610432
|
||||
7260 0.384386446084141 0.217748211186034 0.166638234898108 0 0 0.384386446084141
|
||||
7270 0.3843357368916 0.217925213177252 0.166410523714348 0 0 0.3843357368916
|
||||
7280 0.384285095867745 0.218101977223578 0.166183118644167 0 0 0.384285095867745
|
||||
7290 0.38423452290443 0.218278503702502 0.165956019201929 0 0 0.38423452290443
|
||||
7300 0.384184017796269 0.218454793330934 0.165729224465335 0 0 0.384184017796269
|
||||
7310 0.384133580430773 0.218630846501513 0.16550273392926 0 0 0.384133580430773
|
||||
7320 0.384083210708316 0.218806663561994 0.165276547146322 0 0 0.384083210708316
|
||||
7330 0.384032908426618 0.218982245218435 0.165050663208183 0 0 0.384032908426618
|
||||
7340 0.383982673423143 0.219157592038179 0.164825081384964 0 0 0.383982673423143
|
||||
7350 0.383932505576771 0.219332704443995 0.164599801132775 0 0 0.383932505576771
|
||||
7360 0.383882404728874 0.219507582989577 0.164374821739297 0 0 0.383882404728874
|
||||
7370 0.383832370771199 0.219682228052794 0.164150142718405 0 0 0.383832370771199
|
||||
7380 0.383782403537087 0.219856640215365 0.163925763321722 0 0 0.383782403537087
|
||||
7390 0.383732502862036 0.220030820051494 0.163701682810542 0 0 0.383732502862036
|
||||
7400 0.38368266861773 0.220204768009069 0.163477900608661 0 0 0.38368266861773
|
||||
7410 0.383632900665694 0.220378484571445 0.163254416094249 0 0 0.383632900665694
|
||||
7420 0.383583198896237 0.220551970121499 0.163031228774738 0 0 0.383583198896237
|
||||
7430 0.383533563148466 0.220725225220835 0.162808337927631 0 0 0.383533563148466
|
||||
7440 0.383483993277941 0.220898250373625 0.162585742904316 0 0 0.383483993277941
|
||||
7450 0.383434489154241 0.221071046035105 0.162363443119137 0 0 0.383434489154241
|
||||
7460 0.383385050657527 0.221243612623577 0.16214143803395 0 0 0.383385050657527
|
||||
7470 0.383335677636366 0.221415950667628 0.161919726968738 0 0 0.383335677636366
|
||||
7480 0.383286369944183 0.221588060678877 0.161698309265306 0 0 0.383286369944183
|
||||
7490 0.383237127433505 0.221759943172083 0.161477184261422 0 0 0.383237127433505
|
||||
7500 0.383187949965154 0.221931598633052 0.161256351332102 0 0 0.383187949965154
|
||||
7510 0.383138837433526 0.222103027430401 0.161035810003125 0 0 0.383138837433526
|
||||
7520 0.38308978969063 0.222274230080693 0.160815559609937 0 0 0.38308978969063
|
||||
7530 0.383040806616588 0.222445207002368 0.160595599614219 0 0 0.383040806616588
|
||||
7540 0.382991888048909 0.222615958762604 0.160375929286306 0 0 0.382991888048909
|
||||
7550 0.382943033869436 0.222786485773833 0.160156548095604 0 0 0.382943033869437
|
||||
7560 0.38289424393332 0.222956788541658 0.159937455391662 0 0 0.38289424393332
|
||||
7570 0.382845518120335 0.223126867485724 0.159718650634612 0 0 0.382845518120335
|
||||
7580 0.382796856299152 0.223296723064447 0.159500133234705 0 0 0.382796856299152
|
||||
7590 0.382748258328301 0.223466355771625 0.159281902556676 0 0 0.382748258328302
|
||||
7600 0.382699724075408 0.22363576606932 0.159063958006088 0 0 0.382699724075408
|
||||
7610 0.382651253420638 0.223804954375818 0.158846299044819 0 0 0.382651253420638
|
||||
7620 0.382602846219231 0.223973921196403 0.158628925022828 0 0 0.382602846219231
|
||||
7630 0.382554502358038 0.224142666926032 0.158411835432006 0 0 0.382554502358038
|
||||
7640 0.382506221702866 0.224311192033106 0.15819502966976 0 0 0.382506221702866
|
||||
7650 0.382458004101919 0.224479497047474 0.157978507054445 0 0 0.382458004101919
|
||||
7660 0.382409849464752 0.224647582284842 0.15776226717991 0 0 0.382409849464751
|
||||
7670 0.382361757642396 0.224815448265179 0.157546309377217 0 0 0.382361757642396
|
||||
7680 0.382313728508627 0.224983095429087 0.15733063307954 0 0 0.382313728508627
|
||||
7690 0.382265761939303 0.225150524209878 0.157115237729426 0 0 0.382265761939304
|
||||
7700 0.382217857797598 0.225317735085154 0.156900122712444 0 0 0.382217857797598
|
||||
7710 0.382170015983447 0.225484728404188 0.156685287579259 0 0 0.382170015983447
|
||||
7720 0.382122236355584 0.225651504660082 0.156470731695501 0 0 0.382122236355584
|
||||
7730 0.382074518791155 0.225818064281655 0.1562564545095 0 0 0.382074518791155
|
||||
7740 0.382026863152581 0.225984407749131 0.15604245540345 0 0 0.382026863152581
|
||||
7750 0.381979269326678 0.226150535457584 0.155828733869094 0 0 0.381979269326678
|
||||
7760 0.381931737190082 0.226316447837622 0.15561528935246 0 0 0.381931737190082
|
||||
7770 0.38188426660932 0.226482145355134 0.155402121254186 0 0 0.381884266609319
|
||||
7780 0.381836857484616 0.226647628358389 0.155189229126227 0 0 0.381836857484616
|
||||
7790 0.381789509647377 0.226812897435868 0.154976612211509 0 0 0.381789509647377
|
||||
7800 0.381742223016644 0.226977952870163 0.154764270146481 0 0 0.381742223016644
|
||||
7810 0.381694997451365 0.22714279515362 0.154552202297745 0 0 0.381694997451365
|
||||
7820 0.381647832833362 0.22730742469874 0.154340408134622 0 0 0.381647832833362
|
||||
7830 0.381600729052666 0.227471841889375 0.15412888716329 0 0 0.381600729052666
|
||||
7840 0.381553685977699 0.227636047184802 0.153917638792897 0 0 0.381553685977699
|
||||
7850 0.381506703474806 0.227800041051545 0.153706662423261 0 0 0.381506703474806
|
||||
7860 0.381459781501152 0.227963823639128 0.153495957862024 0 0 0.381459781501152
|
||||
7870 0.381412919880629 0.228127395562256 0.153285524318373 0 0 0.381412919880629
|
||||
7880 0.381366118440433 0.228290757424106 0.153075361016327 0 0 0.381366118440433
|
||||
7890 0.381319377043197 0.228453909704164 0.152865467339033 0 0 0.381319377043197
|
||||
7900 0.381272695704487 0.228616852348095 0.152655843356392 0 0 0.381272695704487
|
||||
7910 0.381226074223305 0.22877958605749 0.152446488165815 0 0 0.381226074223305
|
||||
7920 0.381179512510193 0.228942111144604 0.152237401365589 0 0 0.381179512510193
|
||||
7930 0.38113301042499 0.229104428098673 0.152028582326316 0 0 0.38113301042499
|
||||
7940 0.381086567851745 0.229266537324428 0.151820030527318 0 0 0.381086567851745
|
||||
7950 0.381040184712171 0.229428439095132 0.151611745617039 0 0 0.381040184712171
|
||||
7960 0.380993860835438 0.229590134007076 0.151403726828362 0 0 0.380993860835438
|
||||
7970 0.38094759615096 0.229751622306635 0.151195973844325 0 0 0.38094759615096
|
||||
7980 0.380901390507432 0.229912904521952 0.15098848598548 0 0 0.380901390507432
|
||||
7990 0.380855243825462 0.230073980930138 0.150781262895324 0 0 0.380855243825463
|
||||
8000 0.380809155967371 0.230234852011776 0.150574303955595 0 0 0.380809155967371
|
||||
8
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/mesh.sh
Executable file
8
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/mesh.sh
Executable file
@ -0,0 +1,8 @@
|
||||
cd CFD/
|
||||
rm -r 0
|
||||
rm -r constant/polyMesh
|
||||
m4 system/blockMeshDict.m4 > system/blockMeshDict
|
||||
blockMesh
|
||||
cp -r orig.0 0
|
||||
decomposePar
|
||||
./decomposeParDB.sh
|
||||
27
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/parCFDDEMrun.sh
Executable file
27
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/parCFDDEMrun.sh
Executable file
@ -0,0 +1,27 @@
|
||||
#!/bin/bash
|
||||
|
||||
#===================================================================#
|
||||
# CFD-DEM run script for ore reduction testcase
|
||||
# Thomas Lichtenegger - January 2023
|
||||
#===================================================================#
|
||||
|
||||
#- source CFDEM env vars
|
||||
. ~/.bashrc
|
||||
|
||||
#- include functions
|
||||
source $CFDEM_PROJECT_DIR/etc/functions.sh
|
||||
|
||||
#--------------------------------------------------------------------------------#
|
||||
#- define variables
|
||||
casePath="$(dirname "$(readlink -f ${BASH_SOURCE[0]})")"
|
||||
logpath=$casePath
|
||||
headerText="run_parallel_rcfdemSolverRhoSteadyPimpleChem_reductionTest_CFDDEM"
|
||||
logfileName="log_$headerText"
|
||||
solverName="rcfdemSolverRhoSteadyPimpleChem"
|
||||
nrProcs="8"
|
||||
machineFileName="none" # yourMachinefileName | none
|
||||
debugMode="off" # on | off| strict
|
||||
#--------------------------------------------------------------------------------#
|
||||
|
||||
#- call function to run a parallel CFD-DEM case
|
||||
parCFDDEMrun $logpath $logfileName $casePath $headerText $solverName $nrProcs $machineFileName $debugMode
|
||||
29
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/parDEMrun.sh
Executable file
29
tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/parDEMrun.sh
Executable file
@ -0,0 +1,29 @@
|
||||
#!/bin/bash
|
||||
|
||||
#===================================================================#
|
||||
# DEM run script for ore reduction testcase
|
||||
# Thomas Lichtenegger - January 2023
|
||||
#===================================================================#
|
||||
|
||||
#- source CFDEM env vars
|
||||
. ~/.bashrc
|
||||
|
||||
#- include functions
|
||||
source $CFDEM_PROJECT_DIR/etc/functions.sh
|
||||
|
||||
echo "starting DEM run in parallel..."
|
||||
#--------------------------------------------------------------------------------#
|
||||
#- define variables
|
||||
casePath="$(dirname "$(readlink -f ${BASH_SOURCE[0]})")"
|
||||
logpath="$casePath"
|
||||
headerText="run_liggghts_init_DEM"
|
||||
logfileName="log_$headerText"
|
||||
solverName="in.liggghts_init"
|
||||
nrProcs=8
|
||||
machineFileName="none"
|
||||
debugMode="off"
|
||||
#--------------------------------------------------------------------------------#
|
||||
|
||||
#- call function to run DEM case
|
||||
parDEMrun $logpath $logfileName $casePath $headerText $solverName $nrProcs $machineFileName $debugMode
|
||||
|
||||
File diff suppressed because it is too large
Load Diff
Reference in New Issue
Block a user