diff --git a/tutorials/cfdemSolverRhoPimple/FinesColumnStickyBidisperse/parCFDDEMrun.sh b/tutorials/cfdemSolverRhoPimple/FinesColumnStickyBidisperse/parCFDDEMrun.sh index 6514ca69..114a8a03 100755 --- a/tutorials/cfdemSolverRhoPimple/FinesColumnStickyBidisperse/parCFDDEMrun.sh +++ b/tutorials/cfdemSolverRhoPimple/FinesColumnStickyBidisperse/parCFDDEMrun.sh @@ -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 diff --git a/tutorials/cfdemSolverRhoPimple/FinesColumnStickyBidisperse/parDEMrun.sh b/tutorials/cfdemSolverRhoPimple/FinesColumnStickyBidisperse/parDEMrun.sh index 43f0cadf..85e05ab2 100755 --- a/tutorials/cfdemSolverRhoPimple/FinesColumnStickyBidisperse/parDEMrun.sh +++ b/tutorials/cfdemSolverRhoPimple/FinesColumnStickyBidisperse/parDEMrun.sh @@ -3,7 +3,7 @@ #===================================================================# # DEM run script for FinesColumn testcase # init FinesColumn -# Thomas Lichtenegger - January 2017 +# Thomas Lichtenegger - June 2023 #===================================================================# #- source CFDEM env vars diff --git a/tutorials/cfdemSolverRhoPimple/FinesColumnStickyMonodisperse/parCFDDEMrun.sh b/tutorials/cfdemSolverRhoPimple/FinesColumnStickyMonodisperse/parCFDDEMrun.sh index b584c904..5514f4b7 100755 --- a/tutorials/cfdemSolverRhoPimple/FinesColumnStickyMonodisperse/parCFDDEMrun.sh +++ b/tutorials/cfdemSolverRhoPimple/FinesColumnStickyMonodisperse/parCFDDEMrun.sh @@ -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 diff --git a/tutorials/cfdemSolverRhoPimple/FinesColumnStickyMonodisperse/parDEMrun.sh b/tutorials/cfdemSolverRhoPimple/FinesColumnStickyMonodisperse/parDEMrun.sh index 43f0cadf..85e05ab2 100755 --- a/tutorials/cfdemSolverRhoPimple/FinesColumnStickyMonodisperse/parDEMrun.sh +++ b/tutorials/cfdemSolverRhoPimple/FinesColumnStickyMonodisperse/parDEMrun.sh @@ -3,7 +3,7 @@ #===================================================================# # DEM run script for FinesColumn testcase # init FinesColumn -# Thomas Lichtenegger - January 2017 +# Thomas Lichtenegger - June 2023 #===================================================================# #- source CFDEM env vars diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/chemistryProperties b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/chemistryProperties new file mode 100644 index 00000000..fe5b84f2 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/chemistryProperties @@ -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; +} + + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/combustionProperties b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/combustionProperties new file mode 100644 index 00000000..ecec7ca7 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/combustionProperties @@ -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; //PaSR; //PaSR; +combustionModel none; // OF6 + +active false; //true; + +//noCombustionCoeffs +noneCoeffs +{ +} + + +PaSRCoeffs +{ + Cmix 1.0; + turbulentReaction off; +} + + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/couplingProperties b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/couplingProperties new file mode 100755 index 00000000..6e3562df --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/couplingProperties @@ -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 +{ +} + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/foam.dat b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/foam.dat new file mode 100644 index 00000000..631bbb1d --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/foam.dat @@ -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; + } +} + + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/foam.inp b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/foam.inp new file mode 100644 index 00000000..77a83ba0 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/foam.inp @@ -0,0 +1,12 @@ +species +( + CO + CO2 + N2 +); + +reactions +{ +} + + diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/g b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/g new file mode 100755 index 00000000..317bdd50 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/g @@ -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 ); + + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/liggghtsCommands b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/liggghtsCommands new file mode 100755 index 00000000..148b23b3 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/liggghtsCommands @@ -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; +}*/ diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/recProperties b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/recProperties new file mode 100755 index 00000000..9c5cccbd --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/recProperties @@ -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 +{ +} + + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/thermophysicalProperties b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/thermophysicalProperties new file mode 100644 index 00000000..6628cbac --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/thermophysicalProperties @@ -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; + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/transportProperties b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/transportProperties new file mode 100644 index 00000000..b0eb45e4 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/transportProperties @@ -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; +} + + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/turbulenceProperties b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/turbulenceProperties new file mode 100644 index 00000000..b2fc0a43 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/constant/turbulenceProperties @@ -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; + + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/decomposeParDB.sh b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/decomposeParDB.sh new file mode 100755 index 00000000..740146a8 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/decomposeParDB.sh @@ -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 diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/CO b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/CO new file mode 100755 index 00000000..f25b94a8 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/CO @@ -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; + } +} + + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/CO2 b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/CO2 new file mode 100755 index 00000000..8638fd72 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/CO2 @@ -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; + } +} + + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/Ksl b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/Ksl new file mode 100755 index 00000000..c6c06fb4 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/Ksl @@ -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; + } + +} + + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/N2 b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/N2 new file mode 100755 index 00000000..aa13ddf4 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/N2 @@ -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; + } +} + + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/T b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/T new file mode 100755 index 00000000..e485064c --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/T @@ -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; + } +} + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/U b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/U new file mode 100644 index 00000000..e920040a --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/U @@ -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; + } +} + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/Us b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/Us new file mode 100755 index 00000000..c44672e0 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/Us @@ -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; + } +} + + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/addSource b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/addSource new file mode 100644 index 00000000..c7b35d43 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/addSource @@ -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; + } +} + + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/p b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/p new file mode 100644 index 00000000..17b97e3b --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/p @@ -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; + } +} + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/particleForces b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/particleForces new file mode 100755 index 00000000..37ded0ae --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/particleForces @@ -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; + } +} + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/voidfraction b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/voidfraction new file mode 100755 index 00000000..66037164 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/orig.0/voidfraction @@ -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; + } +} + + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/system/blockMeshDict b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/system/blockMeshDict new file mode 100644 index 00000000..1b763758 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/system/blockMeshDict @@ -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 +( +); diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/system/blockMeshDict.m4 b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/system/blockMeshDict.m4 new file mode 100644 index 00000000..b09d8535 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/system/blockMeshDict.m4 @@ -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 +( +); diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/system/controlDict b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/system/controlDict new file mode 100644 index 00000000..1f8e15b2 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/system/controlDict @@ -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; + } +} + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/system/decomposeParDict b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/system/decomposeParDict new file mode 100644 index 00000000..e39a11c3 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/system/decomposeParDict @@ -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; +} + +// ************************************************************************* // + diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/system/fvOptions b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/system/fvOptions new file mode 100644 index 00000000..cc95b6ea --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/system/fvOptions @@ -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; + } +} + diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/system/fvSchemes b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/system/fvSchemes new file mode 100644 index 00000000..6f001b8e --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/system/fvSchemes @@ -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 ; +} + + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/system/fvSolution b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/system/fvSolution new file mode 100644 index 00000000..09d7d97a --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/CFD/system/fvSolution @@ -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; + } +} + + +// ************************************************************************* // diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/DEMparameters b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/DEMparameters new file mode 100644 index 00000000..20d276c3 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/DEMparameters @@ -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}) diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.chemistry b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.chemistry new file mode 100644 index 00000000..10aa3855 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.chemistry @@ -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 diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.coupling_integration b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.coupling_integration new file mode 100644 index 00000000..4f7143c1 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.coupling_integration @@ -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 diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.liggghts_dataDrivenCFD b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.liggghts_dataDrivenCFD new file mode 100755 index 00000000..ba710515 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.liggghts_dataDrivenCFD @@ -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 diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.liggghts_init b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.liggghts_init new file mode 100755 index 00000000..74798cbd --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.liggghts_init @@ -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 diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.material_parameters b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.material_parameters new file mode 100644 index 00000000..462dff86 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.material_parameters @@ -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 diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.monitor b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.monitor new file mode 100644 index 00000000..e62c3a65 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.monitor @@ -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 + diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.output b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.output new file mode 100644 index 00000000..4352da35 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/in.output @@ -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] diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/liggghts.restart b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/liggghts.restart new file mode 100644 index 00000000..80bec21b Binary files /dev/null and b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/liggghts.restart differ diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/monitor/.gitkeep b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/DEM/monitor/.gitkeep new file mode 100644 index 00000000..e69de29b diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/README b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/README new file mode 100644 index 00000000..0e538c60 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/README @@ -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. diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/mass_layer_opt.dat b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/mass_layer_opt.dat new file mode 100644 index 00000000..fe8f6efa --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/mass_layer_opt.dat @@ -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 diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/mesh.sh b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/mesh.sh new file mode 100755 index 00000000..a3543858 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/mesh.sh @@ -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 diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/parCFDDEMrun.sh b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/parCFDDEMrun.sh new file mode 100755 index 00000000..53707ef1 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/parCFDDEMrun.sh @@ -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 diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/parDEMrun.sh b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/parDEMrun.sh new file mode 100755 index 00000000..d19fe3b8 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/parDEMrun.sh @@ -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 + diff --git a/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/target_mass_loss_pellets.txt b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/target_mass_loss_pellets.txt new file mode 100644 index 00000000..5af0a8d9 --- /dev/null +++ b/tutorials/rcfdemSolverRhoSteadyPimpleChem/reductionTest/target_mass_loss_pellets.txt @@ -0,0 +1,1668 @@ +# t [sec] mass loss [g] +5 0 +10 0.13 +15 0.48 +20 1.31 +25 1.71 +30 2.05 +35 2.58 +40 3 +45 3.39 +50 4.07 +55 4.4 +60 4.78 +65 5.31 +70 5.67 +75 6.21 +80 6.54 +85 7.02 +90 7.36 +95 7.67 +100 7.98 +105 8.25 +110 8.61 +115 8.89 +120 9.41 +125 9.78 +130 10.21 +135 10.52 +140 10.74 +145 11.11 +150 11.39 +155 11.74 +160 12.09 +165 12.71 +170 12.81 +175 13.1 +180 13.38 +185 13.66 +190 14.14 +195 14.31 +200 14.61 +205 14.88 +210 15.23 +215 15.55 +220 15.81 +225 16.14 +230 16.46 +235 16.8 +240 17.33 +245 17.43 +250 17.72 +255 17.92 +260 18.29 +265 18.62 +270 18.88 +275 19.03 +280 19.43 +285 19.67 +290 19.92 +295 20.23 +300 20.47 +305 20.77 +310 20.96 +315 21.26 +320 21.56 +325 21.8 +330 22.12 +335 22.23 +340 22.49 +345 22.79 +350 23.19 +355 23.34 +360 23.5 +365 23.74 +370 24 +375 24.34 +380 24.58 +385 24.79 +390 25.14 +395 25.29 +400 25.51 +405 25.73 +410 25.89 +415 26.21 +420 26.44 +425 26.7 +430 27.03 +435 27.17 +440 27.4 +445 27.67 +450 27.82 +455 28.03 +460 28.23 +465 28.5 +470 28.81 +475 29.03 +480 29.31 +485 29.54 +490 29.71 +495 29.99 +500 30.28 +505 30.39 +510 30.65 +515 30.83 +520 31 +525 31.15 +530 31.34 +535 31.53 +540 31.87 +545 32.02 +550 32.22 +555 32.3 +560 32.55 +565 32.67 +570 32.89 +575 33.04 +580 33.22 +585 33.38 +590 33.7 +595 33.85 +600 33.93 +605 34.15 +610 34.39 +615 34.51 +620 34.68 +625 34.85 +630 35.13 +635 35.34 +640 35.47 +645 35.65 +650 35.8 +655 35.97 +660 36.2 +665 36.41 +670 36.6 +675 36.82 +680 37.02 +685 37.13 +690 37.2 +695 37.4 +700 37.54 +705 37.71 +710 37.9 +715 38.16 +720 38.34 +725 38.46 +730 38.66 +735 38.78 +740 38.89 +745 39.09 +750 39.25 +755 39.4 +760 39.5 +765 39.71 +770 39.83 +775 39.99 +780 40.17 +785 40.26 +790 40.43 +795 40.65 +800 40.83 +805 41.01 +810 41.16 +815 41.29 +820 41.45 +825 41.63 +830 41.76 +835 41.87 +840 41.99 +845 42.22 +850 42.41 +855 42.56 +860 42.75 +865 42.82 +870 43.15 +875 43.2 +880 43.38 +885 43.46 +890 43.67 +895 43.73 +900 43.85 +905 43.92 +910 44.14 +915 44.2 +920 44.33 +925 44.45 +930 44.59 +935 44.77 +940 44.84 +945 44.97 +950 45.26 +955 45.36 +960 45.46 +965 45.71 +970 45.83 +975 45.92 +980 46.06 +985 46.14 +990 46.31 +995 46.39 +1000 46.48 +1005 46.7 +1010 46.76 +1015 46.96 +1020 46.97 +1025 47.11 +1030 47.24 +1035 47.39 +1040 47.47 +1045 47.53 +1050 47.66 +1055 47.81 +1060 47.98 +1065 48.15 +1070 48.29 +1075 48.39 +1080 48.55 +1085 48.62 +1090 48.82 +1095 48.95 +1100 49.07 +1105 49.11 +1110 49.3 +1115 49.36 +1120 49.53 +1125 49.57 +1130 49.69 +1135 49.77 +1140 49.9 +1145 50.07 +1150 50.21 +1155 50.32 +1160 50.41 +1165 50.49 +1170 50.57 +1175 50.71 +1180 50.71 +1185 50.84 +1190 50.99 +1195 51.13 +1200 51.28 +1205 51.35 +1210 51.43 +1215 51.56 +1220 51.69 +1225 51.76 +1230 51.78 +1235 51.79 +1240 52 +1245 52.02 +1250 52.35 +1255 52.43 +1260 52.54 +1265 52.63 +1270 52.68 +1275 52.78 +1280 52.89 +1285 52.99 +1290 53.14 +1295 53.2 +1300 53.36 +1305 53.55 +1310 53.55 +1315 53.65 +1320 53.81 +1325 53.85 +1330 53.95 +1335 54.08 +1340 54.16 +1345 54.28 +1350 54.45 +1355 54.42 +1360 54.54 +1365 54.54 +1370 54.74 +1375 54.76 +1380 54.88 +1385 54.97 +1390 55.02 +1395 55.16 +1400 55.19 +1405 55.38 +1410 55.37 +1415 55.56 +1420 55.56 +1425 55.58 +1430 55.68 +1435 55.9 +1440 55.92 +1445 56.01 +1450 56.09 +1455 56.24 +1460 56.32 +1465 56.46 +1470 56.59 +1475 56.56 +1480 56.66 +1485 56.79 +1490 56.92 +1495 56.86 +1500 56.98 +1505 57.26 +1510 57.27 +1515 57.44 +1520 57.52 +1525 57.43 +1530 57.39 +1535 57.74 +1540 57.81 +1545 57.88 +1550 57.73 +1555 58.05 +1560 58.2 +1565 58.39 +1570 58.46 +1575 58.59 +1580 58.53 +1585 58.68 +1590 58.76 +1595 58.87 +1600 59 +1605 59.07 +1610 59.23 +1615 59.32 +1620 59.5 +1625 59.6 +1630 59.67 +1635 59.83 +1640 59.89 +1645 59.91 +1650 60.04 +1655 60.09 +1660 60.16 +1665 60.27 +1670 60.35 +1675 60.47 +1680 60.58 +1685 60.64 +1690 60.78 +1695 60.84 +1700 60.99 +1705 61.11 +1710 61.17 +1715 61.32 +1720 61.46 +1725 61.58 +1730 61.61 +1735 61.62 +1740 61.75 +1745 61.84 +1750 61.97 +1755 62.04 +1760 62.09 +1765 62.09 +1770 62.2 +1775 62.18 +1780 62.27 +1785 62.34 +1790 62.37 +1795 62.43 +1800 62.5 +1805 62.58 +1810 62.72 +1815 62.75 +1820 62.68 +1825 62.79 +1830 62.73 +1835 62.83 +1840 62.97 +1845 63.07 +1850 63.18 +1855 63.27 +1860 63.3 +1865 63.39 +1870 63.53 +1875 63.54 +1880 63.53 +1885 63.69 +1890 63.81 +1895 63.89 +1900 63.76 +1905 63.72 +1910 63.76 +1915 63.67 +1920 63.79 +1925 63.99 +1930 64.05 +1935 64.17 +1940 64.26 +1945 64.33 +1950 64.49 +1955 64.58 +1960 64.63 +1965 64.82 +1970 64.93 +1975 65.07 +1980 65.17 +1985 65.16 +1990 65.25 +1995 65.29 +2000 65.33 +2005 65.17 +2010 65.16 +2015 65.06 +2020 65.03 +2025 65.31 +2030 65.32 +2035 65.31 +2040 65.27 +2045 65.59 +2050 65.71 +2055 65.6 +2060 65.57 +2065 65.86 +2070 65.72 +2075 65.8 +2080 65.85 +2085 65.94 +2090 66.19 +2095 66.38 +2100 66.39 +2105 66.3 +2110 66.31 +2115 66.34 +2120 66.42 +2125 66.47 +2130 66.7 +2135 66.73 +2140 67.06 +2145 67.1 +2150 67.28 +2155 67.16 +2160 67.08 +2165 67.38 +2170 67.55 +2175 67.75 +2180 67.9 +2185 67.8 +2190 67.88 +2195 68.11 +2200 68.05 +2205 68.13 +2210 68.05 +2215 68 +2220 68.14 +2225 68.08 +2230 68.35 +2235 68.3 +2240 68.55 +2245 68.67 +2250 68.78 +2255 68.86 +2260 68.83 +2265 69.03 +2270 68.99 +2275 69.11 +2280 69.14 +2285 69.11 +2290 69.29 +2295 69.22 +2300 69.11 +2305 69.39 +2310 69.35 +2315 69.5 +2320 69.52 +2325 69.68 +2330 69.63 +2335 69.96 +2340 69.96 +2345 69.98 +2350 70.23 +2355 70.27 +2360 70.28 +2365 70.38 +2370 70.51 +2375 70.63 +2380 70.55 +2385 70.54 +2390 70.66 +2395 70.81 +2400 70.72 +2405 70.9 +2410 71.02 +2415 71.15 +2420 71.21 +2425 71.31 +2430 71.51 +2435 71.6 +2440 71.6 +2445 71.63 +2450 71.6 +2455 71.59 +2460 71.54 +2465 71.79 +2470 71.64 +2475 71.74 +2480 71.63 +2485 71.59 +2490 71.63 +2495 71.46 +2500 71.59 +2505 71.69 +2510 71.92 +2515 72.1 +2520 71.9 +2525 72.11 +2530 72.23 +2535 72.46 +2540 72.28 +2545 72.51 +2550 72.63 +2555 72.77 +2560 72.88 +2565 72.91 +2570 72.97 +2575 72.99 +2580 73.13 +2585 73.2 +2590 73.27 +2595 73.16 +2600 72.97 +2605 73.36 +2610 73.12 +2615 73.39 +2620 73.23 +2625 73.21 +2630 73.22 +2635 73.52 +2640 73.54 +2645 73.57 +2650 73.6 +2655 73.58 +2660 73.81 +2665 73.93 +2670 74.03 +2675 73.99 +2680 74.17 +2685 74.18 +2690 74.36 +2695 74.52 +2700 74.54 +2705 74.64 +2710 74.69 +2715 74.79 +2720 74.9 +2725 74.98 +2730 75.08 +2735 75.07 +2740 75.13 +2745 75.18 +2750 75.24 +2755 75.27 +2760 75.04 +2765 75.16 +2770 75.05 +2775 74.95 +2780 74.87 +2785 75.22 +2790 75.35 +2795 75.5 +2800 75.66 +2805 75.72 +2810 75.74 +2815 75.77 +2820 75.91 +2825 75.95 +2830 75.96 +2835 76.09 +2840 76.2 +2845 76.29 +2850 76.21 +2855 76.31 +2860 76.4 +2865 76.54 +2870 76.5 +2875 76.77 +2880 76.74 +2885 76.85 +2890 76.88 +2895 76.96 +2900 77 +2905 76.99 +2910 77.02 +2915 77.16 +2920 77.23 +2925 77.29 +2930 77.43 +2935 77.11 +2940 77.21 +2945 77.45 +2950 77.55 +2955 77.57 +2960 77.68 +2965 77.74 +2970 77.74 +2975 77.81 +2980 77.9 +2985 77.97 +2990 77.99 +2995 78.08 +3000 78.17 +3005 78.28 +3010 78.27 +3015 78.26 +3020 78.36 +3025 78.4 +3030 78.32 +3035 78.51 +3040 78.55 +3045 78.68 +3050 78.64 +3055 78.69 +3060 78.67 +3065 78.82 +3070 78.54 +3075 78.82 +3080 78.68 +3085 78.71 +3090 78.81 +3095 78.79 +3100 78.94 +3105 79.04 +3110 79.1 +3115 79.18 +3120 79.07 +3125 79.31 +3130 79.5 +3135 79.62 +3140 79.62 +3145 79.66 +3150 79.56 +3155 79.65 +3160 79.73 +3165 79.78 +3170 79.92 +3175 79.93 +3180 80.04 +3185 80.05 +3190 80.03 +3195 80.1 +3200 80.26 +3205 80.41 +3210 80.5 +3215 80.44 +3220 80.5 +3225 80.67 +3230 80.63 +3235 80.67 +3240 80.66 +3245 80.74 +3250 80.76 +3255 80.83 +3260 80.86 +3265 80.93 +3270 81.01 +3275 81 +3280 81.09 +3285 81.11 +3290 81.15 +3295 81.16 +3300 81.3 +3305 81.19 +3310 81.32 +3315 81.37 +3320 81.34 +3325 81.33 +3330 81.24 +3335 81.33 +3340 81.52 +3345 81.51 +3350 81.55 +3355 81.68 +3360 81.72 +3365 81.8 +3370 81.8 +3375 81.86 +3380 81.96 +3385 81.96 +3390 82.02 +3395 82.23 +3400 82.28 +3405 82.19 +3410 82.24 +3415 82.35 +3420 82.4 +3425 82.56 +3430 82.55 +3435 82.55 +3440 82.62 +3445 82.65 +3450 82.56 +3455 82.69 +3460 82.78 +3465 82.89 +3470 82.88 +3475 82.9 +3480 82.9 +3485 83.03 +3490 83.1 +3495 83 +3500 83.12 +3505 83.13 +3510 83.17 +3515 83.07 +3520 83.25 +3525 83.28 +3530 83.3 +3535 83.35 +3540 83.2 +3545 83.33 +3550 83.62 +3555 83.6 +3560 83.48 +3565 83.55 +3570 83.53 +3575 83.69 +3580 83.74 +3585 83.77 +3590 83.76 +3595 83.92 +3600 83.87 +3605 83.9 +3610 83.9 +3615 84.06 +3620 84.13 +3625 84.19 +3630 84.21 +3635 84.33 +3640 84.35 +3645 84.48 +3650 84.49 +3655 84.54 +3660 84.54 +3665 84.6 +3670 84.43 +3675 84.59 +3680 84.57 +3685 84.65 +3690 84.77 +3695 84.89 +3700 84.9 +3705 84.96 +3710 84.87 +3715 84.98 +3720 85.06 +3725 85.11 +3730 85.28 +3735 85.31 +3740 85.38 +3745 85.34 +3750 85.11 +3755 85.33 +3760 85.38 +3765 85.42 +3770 85.35 +3775 85.41 +3780 85.66 +3785 85.6 +3790 85.62 +3795 85.63 +3800 85.82 +3805 85.83 +3810 85.73 +3815 85.81 +3820 85.92 +3825 85.88 +3830 85.95 +3835 85.93 +3840 85.94 +3845 86.02 +3850 86.11 +3855 86.12 +3860 86.13 +3865 86.16 +3870 86.2 +3875 86.17 +3880 86.4 +3885 86.45 +3890 86.42 +3895 86.43 +3900 86.62 +3905 86.63 +3910 86.7 +3915 86.66 +3920 86.89 +3925 86.86 +3930 86.92 +3935 86.92 +3940 86.92 +3945 87.07 +3950 87.04 +3955 87.04 +3960 87.02 +3965 86.89 +3970 87 +3975 87.09 +3980 87.25 +3985 87.29 +3990 87.29 +3995 87.35 +4000 87.44 +4005 87.46 +4010 87.53 +4015 87.56 +4020 87.65 +4025 87.7 +4030 87.82 +4035 87.94 +4040 87.88 +4045 87.95 +4050 87.95 +4055 88 +4060 88.07 +4065 88.02 +4070 88.09 +4075 88.14 +4080 88.13 +4085 88.23 +4090 88.34 +4095 88.25 +4100 88.31 +4105 88.38 +4110 88.26 +4115 88.4 +4120 88.46 +4125 88.46 +4130 88.51 +4135 88.59 +4140 88.67 +4145 88.69 +4150 88.65 +4155 88.62 +4160 88.63 +4165 88.76 +4170 88.76 +4175 88.85 +4180 89.02 +4185 89.06 +4190 89.14 +4195 88.99 +4200 89.15 +4205 89.06 +4210 89.22 +4215 89.27 +4220 89.07 +4225 89.14 +4230 89.23 +4235 89.34 +4240 89.27 +4245 89.43 +4250 89.49 +4255 89.62 +4260 89.71 +4265 89.59 +4270 89.56 +4275 89.69 +4280 89.61 +4285 89.72 +4290 89.63 +4295 89.76 +4300 89.73 +4305 89.82 +4310 89.87 +4315 89.96 +4320 89.93 +4325 90.17 +4330 90.24 +4335 90.22 +4340 90.39 +4345 90.44 +4350 90.46 +4355 90.6 +4360 90.63 +4365 90.72 +4370 90.77 +4375 90.76 +4380 90.81 +4385 90.94 +4390 91.01 +4395 90.86 +4400 90.82 +4405 90.85 +4410 90.98 +4415 91.12 +4420 91.15 +4425 91.22 +4430 91.27 +4435 91.36 +4440 91.43 +4445 91.45 +4450 91.55 +4455 91.56 +4460 91.6 +4465 91.65 +4470 91.71 +4475 91.76 +4480 91.79 +4485 91.89 +4490 91.92 +4495 92.05 +4500 91.97 +4505 91.91 +4510 92.01 +4515 92.06 +4520 92.18 +4525 92.22 +4530 92.23 +4535 92.26 +4540 92.31 +4545 92.41 +4550 92.44 +4555 92.42 +4560 92.27 +4565 92.34 +4570 92.38 +4575 92.54 +4580 92.63 +4585 92.66 +4590 92.71 +4595 92.82 +4600 92.84 +4605 93 +4610 92.93 +4615 93.07 +4620 93.08 +4625 93.06 +4630 92.96 +4635 93.17 +4640 93.23 +4645 93.14 +4650 93.12 +4655 92.89 +4660 93.11 +4665 93.33 +4670 93.32 +4675 93.38 +4680 93.34 +4685 93.39 +4690 93.5 +4695 93.55 +4700 93.59 +4705 93.55 +4710 93.47 +4715 93.44 +4720 93.61 +4725 93.78 +4730 93.76 +4735 93.68 +4740 93.75 +4745 93.82 +4750 93.72 +4755 93.84 +4760 93.95 +4765 94.02 +4770 94.05 +4775 94.02 +4780 94.18 +4785 94.12 +4790 94.14 +4795 94.2 +4800 94.16 +4805 94.07 +4810 94.13 +4815 94.26 +4820 94.3 +4825 94.33 +4830 94.39 +4835 94.27 +4840 94.3 +4845 94.15 +4850 94.23 +4855 94.39 +4860 94.35 +4865 94.42 +4870 94.56 +4875 94.63 +4880 94.75 +4885 94.81 +4890 94.8 +4895 94.84 +4900 94.87 +4905 94.87 +4910 94.98 +4915 95.03 +4920 94.99 +4925 95.09 +4930 95.11 +4935 95.12 +4940 95.17 +4945 95.28 +4950 95.35 +4955 95.27 +4960 95.36 +4965 95.3 +4970 95.23 +4975 95.07 +4980 95.11 +4985 95.2 +4990 95.2 +4995 95.29 +5000 95.48 +5005 95.57 +5010 95.67 +5015 95.73 +5020 95.69 +5025 95.86 +5030 95.94 +5035 95.92 +5040 95.92 +5045 96.04 +5050 96.1 +5055 96.24 +5060 96.27 +5065 96.27 +5070 96.27 +5075 96.31 +5080 96.54 +5085 96.56 +5090 96.53 +5095 96.5 +5100 96.29 +5105 96.41 +5110 96.62 +5115 96.65 +5120 96.69 +5125 96.8 +5130 96.85 +5135 96.93 +5140 96.94 +5145 97 +5150 97.13 +5155 97.23 +5160 97.32 +5165 97.27 +5170 97.25 +5175 97.3 +5180 97.37 +5185 97.36 +5190 97.51 +5195 97.53 +5200 97.51 +5205 97.56 +5210 97.63 +5215 97.62 +5220 97.72 +5225 97.72 +5230 97.8 +5235 97.86 +5240 97.83 +5245 97.79 +5250 97.76 +5255 97.78 +5260 97.89 +5265 97.88 +5270 97.91 +5275 97.99 +5280 97.93 +5285 97.99 +5290 98.06 +5295 98.12 +5300 97.98 +5305 98.05 +5310 98.22 +5315 98.23 +5320 98.25 +5325 98.25 +5330 98.32 +5335 98.18 +5340 98.17 +5345 98.26 +5350 98.42 +5355 98.39 +5360 98.5 +5365 98.49 +5370 98.49 +5375 98.43 +5380 98.5 +5385 98.49 +5390 98.58 +5395 98.58 +5400 98.68 +5405 98.71 +5410 98.79 +5415 98.89 +5420 98.86 +5425 98.85 +5430 98.86 +5435 98.98 +5440 98.99 +5445 98.99 +5450 99.06 +5455 99.12 +5460 99.15 +5465 99.22 +5470 99.26 +5475 99.26 +5480 99.36 +5485 99.31 +5490 99.34 +5495 99.29 +5500 99.33 +5505 99.36 +5510 99.42 +5515 99.38 +5520 99.38 +5525 99.48 +5530 99.57 +5535 99.61 +5540 99.64 +5545 99.71 +5550 99.8 +5555 99.8 +5560 99.88 +5565 99.91 +5570 99.89 +5575 99.96 +5580 100.02 +5585 100.06 +5590 100.1 +5595 100.12 +5600 99.94 +5605 99.81 +5610 99.88 +5615 99.95 +5620 100 +5625 100.06 +5630 100.06 +5635 100.09 +5640 100.09 +5645 100.13 +5650 100.13 +5655 100.27 +5660 100.27 +5665 100.28 +5670 100.33 +5675 100.44 +5680 100.43 +5685 100.23 +5690 100.29 +5695 100.18 +5700 100.33 +5705 100.48 +5710 100.54 +5715 100.53 +5720 100.61 +5725 100.62 +5730 100.63 +5735 100.73 +5740 100.69 +5745 100.76 +5750 100.81 +5755 100.73 +5760 100.86 +5765 100.95 +5770 100.99 +5775 101.02 +5780 101.04 +5785 100.99 +5790 101.12 +5795 101.14 +5800 101.21 +5805 101.2 +5810 101.17 +5815 101.27 +5820 101.28 +5825 101.35 +5830 101.33 +5835 101.45 +5840 101.53 +5845 101.58 +5850 101.57 +5855 101.64 +5860 101.66 +5865 101.68 +5870 101.76 +5875 101.78 +5880 101.82 +5885 101.81 +5890 101.83 +5895 101.81 +5900 101.84 +5905 101.86 +5910 102.02 +5915 102.03 +5920 102.16 +5925 102.22 +5930 102.07 +5935 102.22 +5940 101.93 +5945 101.93 +5950 101.86 +5955 102.03 +5960 102.05 +5965 102.07 +5970 101.91 +5975 101.79 +5980 102.04 +5985 102.29 +5990 102.38 +5995 102.34 +6000 102.41 +6005 102.47 +6010 102.54 +6015 102.62 +6020 102.78 +6025 102.81 +6030 102.79 +6035 102.83 +6040 102.9 +6045 102.94 +6050 103.1 +6055 103.14 +6060 103.09 +6065 103.15 +6070 103.26 +6075 103.29 +6080 103.36 +6085 103.27 +6090 103.19 +6095 103.28 +6100 103.3 +6105 103.43 +6110 103.44 +6115 103.48 +6120 103.51 +6125 103.72 +6130 103.8 +6135 103.7 +6140 103.79 +6145 103.84 +6150 103.93 +6155 104.14 +6160 104.1 +6165 103.94 +6170 103.92 +6175 104.01 +6180 104.21 +6185 104.14 +6190 104.05 +6195 104.03 +6200 104.1 +6205 104.24 +6210 104.2 +6215 104.16 +6220 104.16 +6225 104.24 +6230 104.32 +6235 104.36 +6240 104.31 +6245 104.37 +6250 104.26 +6255 104.29 +6260 104.47 +6265 104.46 +6270 104.48 +6275 104.32 +6280 104.32 +6285 104.31 +6290 104.07 +6295 104.33 +6300 104.39 +6305 104.33 +6310 104.63 +6315 104.51 +6320 104.63 +6325 104.6 +6330 104.59 +6335 104.65 +6340 104.65 +6345 104.74 +6350 104.77 +6355 104.82 +6360 104.96 +6365 104.99 +6370 105.17 +6375 105.14 +6380 105.1 +6385 105.2 +6390 105.29 +6395 105.31 +6400 105.33 +6405 105.45 +6410 105.49 +6415 105.45 +6420 105.33 +6425 105.55 +6430 105.45 +6435 105.48 +6440 105.49 +6445 105.31 +6450 105.38 +6455 105.37 +6460 105.18 +6465 105.12 +6470 105.02 +6475 104.94 +6480 104.92 +6485 104.94 +6490 104.94 +6495 104.85 +6500 105 +6505 105.03 +6510 105.09 +6515 105.19 +6520 105.25 +6525 105.07 +6530 105.25 +6535 105.29 +6540 105.33 +6545 105.32 +6550 105.36 +6555 105.28 +6560 105.37 +6565 105.43 +6570 105.5 +6575 105.47 +6580 105.43 +6585 105.47 +6590 105.45 +6595 105.47 +6600 105.52 +6605 105.69 +6610 105.76 +6615 105.59 +6620 105.74 +6625 105.72 +6630 105.72 +6635 105.63 +6640 105.75 +6645 105.75 +6650 105.77 +6655 105.85 +6660 106 +6665 105.81 +6670 105.89 +6675 105.83 +6680 105.96 +6685 105.97 +6690 105.93 +6695 105.93 +6700 106.43 +6705 105.96 +6710 105.92 +6715 106.41 +6720 106.27 +6725 106.34 +6730 106.34 +6735 106.36 +6740 106.34 +6745 106.37 +6750 106.65 +6755 106.61 +6760 106.59 +6765 106.27 +6770 106.55 +6775 106.53 +6780 106.93 +6785 106.98 +6790 106.55 +6795 106.58 +6800 106.47 +6805 106.62 +6810 106.66 +6815 106.76 +6820 106.71 +6825 106.81 +6830 106.77 +6835 106.81 +6840 106.8 +6845 106.95 +6850 107.02 +6855 107.08 +6860 107.36 +6865 107.41 +6870 107.53 +6875 107.49 +6880 107.6 +6885 107.66 +6890 107.71 +6895 107.76 +6900 107.88 +6905 107.8 +6910 107.62 +6915 107.75 +6920 107.72 +6925 107.55 +6930 107.59 +6935 107.75 +6940 107.76 +6945 107.75 +6950 107.75 +6955 107.79 +6960 108.11 +6965 108.03 +6970 108.02 +6975 108.09 +6980 108 +6985 108.11 +6990 108.15 +6995 108.29 +7000 108.12 +7005 108.28 +7010 108.27 +7015 108.16 +7020 108.27 +7025 108.19 +7030 108.23 +7035 108.17 +7040 108.29 +7045 108.23 +7050 108.26 +7055 108.2 +7060 108.19 +7065 108.48 +7070 108.57 +7075 108.61 +7080 108.43 +7085 108.33 +7090 108.3 +7095 108.38 +7100 108.25 +7105 108.36 +7110 108.47 +7115 108.74 +7120 108.71 +7125 108.83 +7130 108.94 +7135 109.05 +7140 109.12 +7145 109.29 +7150 109.4 +7155 109.48 +7160 109.63 +7165 109.58 +7170 109.66 +7175 109.57 +7180 109.31 +7185 109.47 +7190 109.42 +7195 109.46 +7200 109.52 +7205 109.46 +7210 109.64 +7215 109.81 +7220 109.79 +7225 109.87 +7230 109.81 +7235 109.92 +7240 109.86 +7245 109.97 +7250 109.91 +7255 109.99 +7260 110.07 +7265 110.14 +7270 110 +7275 110.24 +7280 110.25 +7285 110.24 +7290 110.35 +7295 110.32 +7300 110.28 +7305 110.2 +7310 110.26 +7315 110.34 +7320 110.51 +7325 110.43 +7330 110.49 +7335 110.51 +7340 110.45 +7345 110.62 +7350 110.63 +7355 110.58 +7360 110.47 +7365 110.45 +7370 110.55 +7375 110.53 +7380 110.53 +7385 110.43 +7390 110.47 +7395 110.56 +7400 110.74 +7405 110.69 +7410 110.73 +7415 110.74 +7420 110.5 +7425 110.29 +7430 110.37 +7435 110.42 +7440 110.51 +7445 110.58 +7450 110.7 +7455 110.69 +7460 110.71 +7465 110.6 +7470 110.56 +7475 110.75 +7480 110.75 +7485 110.72 +7490 110.76 +7495 111.07 +7500 111.12 +7505 111.2 +7510 111.08 +7515 111.07 +7520 111.14 +7525 111.27 +7530 111.37 +7535 111.48 +7540 111.59 +7545 111.62 +7550 111.71 +7555 111.75 +7560 111.8 +7565 111.92 +7570 111.89 +7575 111.83 +7580 111.64 +7585 111.63 +7590 111.56 +7595 111.49 +7600 111.5 +7605 111.49 +7610 111.4 +7615 111.4 +7620 111.51 +7625 111.37 +7630 111.53 +7635 111.5 +7640 111.54 +7645 111.74 +7650 111.66 +7655 111.69 +7660 111.65 +7665 111.63 +7670 111.74 +7675 111.94 +7680 111.87 +7685 111.74 +7690 111.81 +7695 111.79 +7700 111.88 +7705 111.66 +7710 111.67 +7715 111.7 +7720 111.58 +7725 111.87 +7730 111.84 +7735 111.95 +7740 111.8 +7745 111.69 +7750 111.73 +7755 111.82 +7760 111.92 +7765 111.87 +7770 112.11 +7775 112.26 +7780 112.31 +7785 112.42 +7790 112.51 +7795 112.58 +7800 112.64 +7805 112.71 +7810 112.73 +7815 112.73 +7820 112.77 +7825 112.83 +7830 112.8 +7835 112.72 +7840 112.82 +7845 112.86 +7850 112.92 +7855 112.97 +7860 113.04 +7865 113.03 +7870 113.03 +7875 113.05 +7880 113.17 +7885 113.17 +7890 113.17 +7895 113.17 +7900 113.18 +7905 113.19 +7910 113.21 +7915 113.04 +7920 113 +7925 113.2 +7930 113.05 +7935 112.91 +7940 113.3 +7945 113.3 +7950 113.38 +7955 113.41 +7960 113.54 +7965 113.49 +7970 113.65 +7975 113.63 +7980 113.57 +7985 113.69 +7990 113.72 +7995 114.06 +8000 113.96 +8005 113.97 +8010 113.95 +8015 114.03 +8020 113.83 +8025 113.93 +8030 114.05 +8035 114.09 +8040 114.07 +8045 113.71 +8050 113.72 +8055 113.94 +8060 114.03 +8065 114.06 +8070 114.05 +8075 114.07 +8080 113.8 +8085 113.94 +8090 114.02 +8095 114.22 +8100 114.24 +8105 114.28 +8110 114.29 +8115 114.28 +8120 114.3 +8125 114.39 +8130 114.47 +8135 114.58 +8140 114.51 +8145 114.49 +8150 114.37 +8155 114.45 +8160 114.37 +8165 114.38 +8170 114.42 +8175 114.48 +8180 114.49 +8185 114.53 +8190 114.48 +8195 114.56 +8200 114.63 +8205 114.61 +8210 114.48 +8215 114.58 +8220 114.56 +8225 114.55 +8230 114.68 +8235 114.68 +8240 114.62 +8245 114.65 +8250 114.65 +8255 114.76 +8260 114.84 +8265 114.81 +8270 114.66 +8275 114.46 +8280 114.63 +8285 114.86 +8290 114.91 +8295 114.89 +8300 114.76 +8305 115 +8310 114.98 +8315 115.08 +8320 115.19 +8325 115.17 +8330 115.06 +8335 115.14