zhyang-dev/CFDEMcoupling-PFM
1
0
Fork 0
mirror of https://github.com/ParticulateFlow/CFDEMcoupling-PFM.git synced 2025-12-08 06:37:44 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

Added mesh generation for tutorial FinesColumn.

Browse Source
This commit is contained in:
Thomas Lichtenegger
2017-01-23 09:36:43 +01:00
parent cccb78b4d3
commit 9ce3eb5d24
40 changed files with 36 additions and 106180 deletions
Download Patch File Download Diff File Expand all files Collapse all files

41
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/Ksl
View File

@ -1,41 +0,0 @@
/*--------------------------------*- 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;
}
}
// ************************************************************************* //

39
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/KslNext
View File

@ -1,39 +0,0 @@
/*--------------------------------*- 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 KslNext;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -3 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
walls
{
type zeroGradient;
}
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
}
// ************************************************************************* //

40
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/KslPrev
View File

@ -1,40 +0,0 @@
/*--------------------------------*- 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 KslPrev;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -3 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
walls
{
type zeroGradient;
}
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
}
// ************************************************************************* //

43
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/T
View File

@ -1,43 +0,0 @@
/*--------------------------------*- 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 293;
boundaryField
{
walls
{
type fixedValue;
value uniform 293;
}
inlet
{
type fixedValue;
value uniform 293;
}
outlet
{
type zeroGradient;
}
}
// ************************************************************************* //

44
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/U
View File

@ -1,44 +0,0 @@
/*--------------------------------*- 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 uniformFixedValue;
uniformValue constant (0 0 0);
}
inlet
{
type interstitialInletVelocity;
inletVelocity uniform (0 0 1.24);//0.88);
alpha voidfraction;
value $internalField;
}
outlet
{
type zeroGradient;
}
}
// ************************************************************************* //

40
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/Us
View File

@ -1,40 +0,0 @@
/*--------------------------------*- 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;
}
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
}
// ************************************************************************* //

40
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/UsNext
View File

@ -1,40 +0,0 @@
/*--------------------------------*- 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 UsNext;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
walls
{
type zeroGradient;
}
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
}
// ************************************************************************* //

39
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/UsPrev
View File

@ -1,39 +0,0 @@
/*--------------------------------*- 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 UsPrev;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
walls
{
type zeroGradient;
}
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
}
// ************************************************************************* //

42
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/addSource
View File

@ -1,42 +0,0 @@
/*--------------------------------*- 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;
}
}
// ************************************************************************* //

48
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/alphaDyn
View File

@ -1,48 +0,0 @@
/*--------------------------------*- 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 alphaDyn;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 0 0 0 0];
internalField uniform 0.0;
boundaryField
{
walls
{
type zeroGradient;
}
inlet
{
type fixedValue;
value uniform 0.01;
/* type groovyBC;
uDyn uDyn;
variables "gF=0.03;rhoF=1396;";
valueExpression "gF/(rhoF*mag(uDyn)+0.01)";
fractionExpression "1";
value uniform 0.0;
*/
}
outlet
{
type zeroGradient;
}
}
// ************************************************************************* //

40
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/alphaSt
View File

@ -1,40 +0,0 @@
/*--------------------------------*- 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 alphaSt;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 0 0 0 0];
internalField uniform 0;
boundaryField
{
walls
{
type zeroGradient;
}
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
}
// ************************************************************************* //

38
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/f
View File

@ -1,38 +0,0 @@
/*--------------------------------*- 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 f;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -2 -2 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
walls
{
type zeroGradient;
}
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
}
// ************************************************************************* //

44
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/p
View File

@ -1,44 +0,0 @@
/*--------------------------------*- 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;
}
}
// ************************************************************************* //

38
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/particleForces
View File

@ -1,38 +0,0 @@
/*--------------------------------*- 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;
}
}
// ************************************************************************* //

39
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/rho
View File

@ -1,39 +0,0 @@
/*--------------------------------*- 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 rho;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -3 0 0 0 0 0];
internalField uniform 1.1885;
boundaryField
{
walls
{
type zeroGradient;
}
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
}
// ************************************************************************* //

40
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/sSmoothField
View File

@ -1,40 +0,0 @@
/*--------------------------------*- 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 sSmoothField;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 0 0 0 0];
internalField uniform 1;
boundaryField
{
walls
{
type zeroGradient;
}
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
}
// ************************************************************************* //

54
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/turb/epsilon
View File

@ -1,54 +0,0 @@
/*--------------------------------*- 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;
location "0";
object epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -3 0 0 0 0];
internalField uniform 0.000765;
boundaryField
{
top
{
type epsilonWallFunction;
value uniform 0.000765;
}
bottom
{
type epsilonWallFunction;
value uniform 0.000765;
}
side
{
type epsilonWallFunction;
value uniform 0.000765;
}
inlet
{
type empty;
}
outlet_up
{
type zeroGradient;
}
outlet_down
{
type zeroGradient;
}
}
// ************************************************************************* //

54
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/turb/k
View File

@ -1,54 +0,0 @@
/*--------------------------------*- 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;
location "0";
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0.0;//0325;
boundaryField
{
top
{
type kqRWallFunction;
value uniform 0.00325;
}
bottom
{
type kqRWallFunction;
value uniform 0.00325;
}
side
{
type kqRWallFunction;
value uniform 0.00325;
}
inlet
{
type empty;
}
outlet_up
{
type zeroGradient;
}
outlet_down
{
type zeroGradient;
}
}
// ************************************************************************* //

49
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/turb/nuTilda
View File

@ -1,49 +0,0 @@
/*--------------------------------*- 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 nuTilda;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
top
{
type zeroGradient;
}
bottom
{
type zeroGradient;
}
side
{
type zeroGradient;
}
inlet
{
type empty;
}
outlet_up
{
type zeroGradient;
}
outlet_down
{
type zeroGradient;
}
}
// ************************************************************************* //

54
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/turb/nut
View File

@ -1,54 +0,0 @@
/*--------------------------------*- 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;
location "0";
object nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
top
{
type nutkWallFunction;
value uniform 0;
}
bottom
{
type nutkWallFunction;
value uniform 0;
}
side
{
type nutkWallFunction;
value uniform 0;
}
inlet
{
type empty;
}
outlet_up
{
type zeroGradient;
}
outlet_down
{
type zeroGradient;
}
}
// ************************************************************************* //

47
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/uDyn
View File

@ -1,47 +0,0 @@
/*--------------------------------*- 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 uDyn;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
walls
{
type slip;
}
defaultFaces
{
type empty;
}
inlet
{
type calculated;
value $internalField;
// type uniformFixedValue;
// uniformValue constant (0 0 1);
}
outlet
{
type zeroGradient;
}
}
// ************************************************************************* //

40
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/vSmoothField
View File

@ -1,40 +0,0 @@
/*--------------------------------*- 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 vSmoothField;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
walls
{
type zeroGradient;
}
inlet
{
type zeroGradient;
}
outlet
{
type zeroGradient;
}
}
// ************************************************************************* //

42
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/voidfraction
View File

@ -1,42 +0,0 @@
/*--------------------------------*- 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;
type zeroGradient;
}
outlet
{
type zeroGradient;
}
}
// ************************************************************************* //

42
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/voidfractionNext
View File

@ -1,42 +0,0 @@
/*--------------------------------*- 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 voidfractionNext;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 0 0 0 0];
internalField uniform 1;
boundaryField
{
walls
{
type zeroGradient;
}
inlet
{
//type fixedValue;
//value uniform 1.0;
type zeroGradient;
}
outlet
{
type zeroGradient;
}
}
// ************************************************************************* //

42
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/voidfractionPrev
View File

@ -1,42 +0,0 @@
/*--------------------------------*- 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 voidfractionPrev;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 0 0 0 0];
internalField uniform 1;
boundaryField
{
walls
{
type zeroGradient;
}
inlet
{
//type fixedValue;
//value uniform 1.0;
type zeroGradient;
}
outlet
{
type zeroGradient;
}
}
// ************************************************************************* //

43
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/constant/polyMesh/boundary
View File

@ -1,43 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/polyMesh";
object boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
3
(
inlet
{
type wall;
inGroups 1(wall);
nFaces 192;
startFace 29488;
}
outlet
{
type wall;
inGroups 1(wall);
nFaces 192;
startFace 29680;
}
walls
{
type wall;
inGroups 1(wall);
nFaces 1696;
startFace 29872;
}
)
// ************************************************************************* //

224
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/constant/polyMesh/faceZones
View File

@ -1,224 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class regIOobject;
location "constant/polyMesh";
object faceZones;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
1
(
plane1
{
type faceZone;
faceLabels List<label>
192
(
3954
3959
3963
3967
3971
3975
3979
3982
3987
3991
3994
3997
4000
4003
4006
4008
4012
4016
4019
4022
4025
4028
4031
4033
4037
4041
4044
4047
4050
4053
4056
4058
4062
4066
4069
4072
4075
4078
4081
4083
4087
4091
4094
4097
4100
4103
4106
4108
4112
4116
4119
4122
4125
4128
4131
4133
4136
4140
4143
4146
4149
4152
4155
4157
12710
12714
12717
12720
12723
12726
12729
12731
12735
12739
12742
12745
12748
12751
12754
12756
12760
12764
12767
12770
12773
12776
12779
12781
12784
12787
12789
12791
12793
12795
12797
12798
17478
17481
17484
17486
17489
17492
17495
17497
17500
17503
17506
17508
17511
17514
17517
17519
17522
17525
17528
17530
17533
17536
17539
17541
17544
17547
17550
17552
17554
17557
17560
17562
22110
22113
22116
22118
22121
22124
22127
22129
22132
22135
22138
22140
22143
22146
22149
22151
22154
22157
22160
22162
22165
22168
22171
22173
22176
22179
22182
22184
22186
22189
22192
22194
26666
26669
26672
26674
26677
26680
26683
26685
26688
26691
26694
26696
26699
26702
26705
26707
26710
26713
26716
26718
26721
26724
26727
26729
26732
26735
26738
26740
26742
26744
26746
26747
)
;
flipMap List<bool> 192{0};
}
)
// ************************************************************************* //

31592
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/constant/polyMesh/faces
View File

File diff suppressed because it is too large Load Diff

29513
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/constant/polyMesh/neighbour
View File

File diff suppressed because it is too large Load Diff

31593
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/constant/polyMesh/owner
View File

File diff suppressed because it is too large Load Diff

11310
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/constant/polyMesh/points
View File

File diff suppressed because it is too large Load Diff

215
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/constant/polyMesh/sets/plane1
View File

@ -1,215 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class faceSet;
location "constant/polyMesh/sets";
object plane1;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
192
(
4097
4100
4103
4106
4108
4112
4116
4119
4122
4125
4128
4131
4133
4136
26666
4140
26669
4143
26672
26674
4146
26677
4149
26680
4152
26683
4155
26685
4157
26688
26691
26694
17478
26696
17481
26699
17484
26702
17486
26705
17489
26707
17492
26710
17495
26713
17497
26716
17500
26718
22110
17503
26721
22113
17506
26724
22116
17508
22118
26727
17511
26729
22121
17514
26732
22124
17517
26735
22127
17519
22129
26738
17522
26740
22132
17525
26742
22135
26744
17528
26746
22138
17530
26747
22140
17533
22143
17536
22146
17539
22149
17541
22151
17544
22154
17547
22157
17550
22160
17552
22162
17554
22165
17557
22168
17560
17562
22171
22173
22176
22179
22182
22184
22186
22189
22192
22194
3954
3959
3963
3967
3971
3975
3979
3982
3987
3991
3994
3997
4000
4003
12710
4006
4008
12714
4012
12717
12720
4016
12723
4019
12726
4022
12729
4025
12731
4028
12735
4031
4033
12739
4037
12742
12745
4041
12748
4044
12751
4047
12754
4050
12756
4053
12760
4056
4058
12764
4062
12767
12770
4066
12773
4069
12776
4072
12779
4075
12781
4078
12784
4081
12787
4083
12789
12791
4087
12793
12795
4091
12797
12798
4094
)
// ************************************************************************* //

215
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/constant/polyMesh/sets/plane1FaceSet
View File

@ -1,215 +0,0 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class faceSet;
location "constant/polyMesh/sets";
object plane1FaceSet;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
192
(
17478
17481
17484
17486
17489
17492
17495
17497
17500
17503
17506
17508
17511
17514
17517
17519
17522
17525
17528
17530
17533
17536
17539
17541
17544
17547
17550
17552
17554
17557
17560
17562
3954
3959
3963
3967
3971
3975
3979
3982
3987
3991
3994
3997
4000
4003
4006
4008
4012
4016
4019
4022
4025
4028
4031
4033
4037
4041
4044
4047
4050
4053
4056
4058
4062
4066
4069
4072
4075
4078
4081
4083
4087
4091
4094
4097
4100
4103
4106
4108
4112
4116
4119
4122
4125
4128
4131
4133
4136
4140
4143
4146
4149
4152
4155
4157
22110
22113
22116
22118
22121
22124
22127
22129
22132
22135
22138
22140
22143
22146
22149
22151
22154
22157
22160
22162
22165
22168
22171
22173
22176
22179
22182
22184
22186
22189
22192
22194
26666
26669
26672
26674
26677
26680
26683
26685
26688
26691
26694
26696
26699
26702
26705
26707
26710
26713
26716
26718
26721
26724
26727
26729
26732
26735
26738
26740
26742
26744
26746
26747
12710
12714
12717
12720
12723
12726
12729
12731
12735
12739
12742
12745
12748
12751
12754
12756
12760
12764
12767
12770
12773
12776
12779
12781
12784
12787
12789
12791
12793
12795
12797
12798
)
// ************************************************************************* //

8
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/geometry Normal file
View File

@ -0,0 +1,8 @@
hsqrt2 0.7071068;
r1 0.05;
r2 0.025;
mr1 -0.05;
mr2 -0.025;
zmin 0.0;
zmax 0.53;
#inputMode merge

96
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/system/controlDict
View File

@ -106,100 +106,4 @@ functions
}
}
/*functions
{
difference
{
// Load the library containing the 'coded' functionObject
functionObjectLibs ("libutilityFunctionObjects.so");
type coded;
// Name of on-the-fly generated functionObject
redirectType error;
code
#{
// Lookup U
Info<< "Looking up field U\n" << endl;
const volVectorField& U = mesh().lookupObject<volVectorField>("U");
Info<< "Reading inlet velocity uInfX\n" << endl;
scalar ULeft = 0.0;
label leftI = mesh().boundaryMesh().findPatchID("left");
const fvPatchVectorField& fvp = U.boundaryField()[leftI];
if (fvp.size())
{
ULeft = fvp[0].x();
}
reduce(ULeft, maxOp<scalar>());
dimensionedScalar uInfX
(
"uInfx",
dimensionSet(0, 1, -1, 0, 0),
ULeft
);
Info << "U at inlet = " << uInfX.value() << " m/s" << endl;
scalar magCylinder = 0.0;
label cylI = mesh().boundaryMesh().findPatchID("cylinder");
const fvPatchVectorField& cylFvp = mesh().C().boundaryField()[cylI];
if (cylFvp.size())
{
magCylinder = mag(cylFvp[0]);
}
reduce(magCylinder, maxOp<scalar>());
dimensionedScalar radius
(
"radius",
dimensionSet(0, 1, 0, 0, 0),
magCylinder
);
Info << "Cylinder radius = " << radius.value() << " m" << endl;
volVectorField UA
(
IOobject
(
"UA",
mesh().time().timeName(),
U.mesh(),
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
U
);
Info<< "\nEvaluating analytical solution" << endl;
const volVectorField& centres = UA.mesh().C();
volScalarField magCentres(mag(centres));
volScalarField theta(acos((centres & vector(1,0,0))/magCentres));
volVectorField cs2theta
(
cos(2*theta)*vector(1,0,0)
+ sin(2*theta)*vector(0,1,0)
);
UA = uInfX*(dimensionedVector(vector(1,0,0))
- pow((radius/magCentres),2)*cs2theta);
// Force writing of UA (since time has not changed)
UA.write();
volScalarField error("error", mag(U-UA)/mag(UA));
Info<<"Writing relative error in U to " << error.objectPath()
<< endl;
error.write();
#};
}
}
*/
// ************************************************************************* //

193
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/system/controlDict.save
View File

@ -1,193 +0,0 @@
/*--------------------------------*- 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 10;
deltaT 0.00125;
writeControl timeStep;
writeInterval 400;
purgeWrite 0;
writeFormat ascii;
writePrecision 6;
writeCompression off;
timeFormat general;
timePrecision 6;
runTimeModifiable true;
functions
{
probes1
{
type probes;
functionObjectLibs ("libsampling.so");
dictionary probesDict;
}
volInt
{
type cellSource;
functionObjectLibs ("libfieldFunctionObjects.so");
enabled true;
outputControl timeStep;
log true;
valueOutput false;
source all;
operation volIntegrate;
fields
(
alphaSt
);
}
inflow
{
type faceSource;
functionObjectLibs ("libfieldFunctionObjects.so");
outputControl timeStep;
log true;
// Output field values as well
valueOutput false;
source in;
sourceName outlet; // replace patch name!!!
operation sum;
fields
(
rho*phi*alpha // maybe fix eqn. for your purpose!!!
);
}
}
/*functions
{
difference
{
// Load the library containing the 'coded' functionObject
functionObjectLibs ("libutilityFunctionObjects.so");
type coded;
// Name of on-the-fly generated functionObject
redirectType error;
code
#{
// Lookup U
Info<< "Looking up field U\n" << endl;
const volVectorField& U = mesh().lookupObject<volVectorField>("U");
Info<< "Reading inlet velocity uInfX\n" << endl;
scalar ULeft = 0.0;
label leftI = mesh().boundaryMesh().findPatchID("left");
const fvPatchVectorField& fvp = U.boundaryField()[leftI];
if (fvp.size())
{
ULeft = fvp[0].x();
}
reduce(ULeft, maxOp<scalar>());
dimensionedScalar uInfX
(
"uInfx",
dimensionSet(0, 1, -1, 0, 0),
ULeft
);
Info << "U at inlet = " << uInfX.value() << " m/s" << endl;
scalar magCylinder = 0.0;
label cylI = mesh().boundaryMesh().findPatchID("cylinder");
const fvPatchVectorField& cylFvp = mesh().C().boundaryField()[cylI];
if (cylFvp.size())
{
magCylinder = mag(cylFvp[0]);
}
reduce(magCylinder, maxOp<scalar>());
dimensionedScalar radius
(
"radius",
dimensionSet(0, 1, 0, 0, 0),
magCylinder
);
Info << "Cylinder radius = " << radius.value() << " m" << endl;
volVectorField UA
(
IOobject
(
"UA",
mesh().time().timeName(),
U.mesh(),
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
U
);
Info<< "\nEvaluating analytical solution" << endl;
const volVectorField& centres = UA.mesh().C();
volScalarField magCentres(mag(centres));
volScalarField theta(acos((centres & vector(1,0,0))/magCentres));
volVectorField cs2theta
(
cos(2*theta)*vector(1,0,0)
+ sin(2*theta)*vector(0,1,0)
);
UA = uInfX*(dimensionedVector(vector(1,0,0))
- pow((radius/magCentres),2)*cs2theta);
// Force writing of UA (since time has not changed)
UA.write();
volScalarField error("error", mag(U-UA)/mag(UA));
Info<<"Writing relative error in U to " << error.objectPath()
<< endl;
error.write();
#};
}
}
*/
// ************************************************************************* //

82
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/system/createPatchDict
View File

@ -1,82 +0,0 @@
/*--------------------------------*- 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;
object createPatchDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
// This application/dictionary controls:
// - optional: create new patches from boundary faces (either given as
// a set of patches or as a faceSet)
// - always: order faces on coupled patches such that they are opposite. This
// is done for all coupled faces, not just for any patches created.
// - optional: synchronise points on coupled patches.
// 1. Create cyclic:
// - specify where the faces should come from
// - specify the type of cyclic. If a rotational specify the rotationAxis
// and centre to make matching easier
// - always create both halves in one invocation with correct 'neighbourPatch'
// setting.
// - optionally pointSync true to guarantee points to line up.
// 2. Correct incorrect cyclic:
// This will usually fail upon loading:
// "face 0 area does not match neighbour 2 by 0.0100005%"
// " -- possible face ordering problem."
// - in polyMesh/boundary file:
// - loosen matchTolerance of all cyclics to get case to load
// - or change patch type from 'cyclic' to 'patch'
// and regenerate cyclic as above
// Do a synchronisation of coupled points after creation of any patches.
// Note: this does not work with points that are on multiple coupled patches
// with transformations (i.e. cyclics).
pointSync false;
// Patches to create.
patches
(
{
name inlet;
patchInfo
{
type patch;
}
constructFrom set;
set inlet;
}
{
name outlet_up;
patchInfo
{
type patch;
}
constructFrom set;
set outlet_up;
}
{
name outlet_down;
patchInfo
{
type patch;
}
constructFrom set;
set outlet_down;
}
);
// ************************************************************************* //

46
tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/0/p_kin → tutorials/cfdemSolverRhoPimple/FinesColumn/CFD/system/topoSetDict
View File

@ -9,37 +9,35 @@ FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
class dictionary;
object topoSetDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 1e5;
boundaryField
{
walls
actions
(
{
name plane1FaceSet;
type faceSet;
action new;
source boxToFace;
sourceInfo
{
type zeroGradient;
box (-0.055 -0.055 0.199)(0.055 0.055 0.201);
}
}
defaultFaces
{
name plane1;
type faceZoneSet;
action new;
source setToFaceZone;
sourceInfo
{
type empty;
faceSet plane1FaceSet;
}
inlet
{
type zeroGradient;
}
outlet
{
type uniformFixedValue;
uniformValue constant 1e5;
}
}
}
);
// ************************************************************************* //

3
tutorials/cfdemSolverRhoPimple/FinesColumn/README
View File

@ -1,7 +1,8 @@
Test case for fines transport through packed particle column
1) fill a cylinder with sperical particles (parDEMrun.sh)
2) introduce gas flow with fines (parCFDDEMrun.sh)
2) create CFD mesh (mesh.sh)
3) introduce gas flow with fines (parCFDDEMrun.sh)
Comments:

4
tutorials/cfdemSolverRhoPimple/FinesColumn/mesh.sh Executable file
View File

@ -0,0 +1,4 @@
cd CFD
blockMesh
topoSet -dict system/topoSetDict
cp -r orig.0/ 0