ENH: added column, Goldschmidt and injectionChannel tutorials for MPPICFoam

This commit is contained in:
william
2013-12-20 10:31:06 +00:00
parent bdd50222b4
commit 005ff02085
34 changed files with 29687 additions and 0 deletions

View File

@ -0,0 +1,57 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format binary;
class volVectorField;
location "0";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0 0 0);
boundaryField
{
walls
{
type fixedValue;
value uniform (0 0 0);
}
lowerInlet
{
type fixedValue;
value uniform (18.7939 6.8404 0);
}
upperInlet
{
type fixedValue;
value uniform (18.7939 6.8404 0);
}
outlet
{
type inletOutlet;
phi phi.air;
inletValue uniform (0 0 0);
value uniform (0 0 0);
}
frontAndBack
{
type symmetry;
}
}
// ************************************************************************* //

View File

@ -0,0 +1,57 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 0;
boundaryField
{
walls
{
type fixedFluxPressure;
phi phi.air;
value $internalField;
}
lowerInlet
{
type fixedFluxPressure;
phi phi.air;
value $internalField;
}
upperInlet
{
type fixedFluxPressure;
phi phi.air;
value $internalField;
}
outlet
{
type fixedValue;
phi phi.air;
value uniform 0;
}
frontAndBack
{
type symmetry;
}
}
// ************************************************************************* //

View File

@ -0,0 +1,22 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / 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 0 );
// ************************************************************************* //

View File

@ -0,0 +1,193 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object particleProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solution
{
active true;
coupled true;
transient yes;
cellValueSourceCorrection on;
interpolationSchemes
{
rho.air cell;
U.air cellPoint;
mu.air cell;
}
averagingMethod basic;
integrationSchemes
{
U Euler;
}
sourceTerms
{
schemes
{
U semiImplicit 1;
}
}
}
constantProperties
{
rho0 1000;
alphaMax 0.9;
}
subModels
{
particleForces
{
WenYuDrag
{
alphac alpha.air;
}
}
injectionModels
{
/*
% matlab/octave calculation of flow rate
alpha=0.1
U=20
A=0.01^2
Q=U*A
d=650e-6
v=(4/3)*pi*(d/2)^3
n=1
rate=alpha*Q/(v*n)
*/
model1
{
type patchInjection;
parcelBasisType fixed;
patchName lowerInlet;
U0 (18.7939 6.8404 0);
nParticle 1;
parcelsPerSecond 1390885;
sizeDistribution
{
type normal;
normalDistribution
{
expectation 650e-6;
variance 25e-6;
minValue 500e-6;
maxValue 800e-6;
}
}
flowRateProfile constant 1;
massTotal 0;
SOI 0;
duration 60;
}
model2
{
type patchInjection;
parcelBasisType fixed;
patchName upperInlet;
U0 (18.7939 6.8404 0);
nParticle 1;
parcelsPerSecond 1390885;
sizeDistribution
{
type normal;
normalDistribution
{
expectation 650e-6;
variance 25e-6;
minValue 500e-6;
maxValue 800e-6;
}
}
flowRateProfile constant 1;
massTotal 0;
SOI 0;
duration 60;
}
}
dispersionModel none;
patchInteractionModel localInteraction;
localInteractionCoeffs
{
patches
(
frontAndBack
{
type rebound;
e 1;
mu 0;
}
walls
{
type rebound;
e 1;
mu 0;
}
);
}
heatTransferModel none;
surfaceFilmModel none;
packingModel none;
dampingModel relaxation;
relaxationCoeffs
{
timeScaleModel
{
type nonEquilibrium;
alphaPacked 0.58;
e 0.9;
}
}
isotropyModel stochastic;
stochasticCoeffs
{
timeScaleModel
{
type isotropic;
alphaPacked 0.58;
e 0.9;
}
}
stochasticCollisionModel none;
radiation off;
}
cloudFunctions
{}
// ************************************************************************* //

View File

@ -0,0 +1,181 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
convertToMeters 0.01;
vertices
(
( 0 -3 -0.5)
( 2 -3 -0.5)
( 0 -3.727940468532405 -0.5)
( 0.34202014332567 -4.667633089318313 -0.5)
( 2.34202014332567 -3.939692620785910 -0.5)
( 4.92380440016309 -3 -0.5)
(41 -3 -0.5)
(41 -6 -0.5)
(53 -6 -0.5)
(53 -3 -0.5)
(53 3 -0.5)
(53 6 -0.5)
(41 6 -0.5)
(41 3 -0.5)
( 4.92380440016309 3 -0.5)
( 2.34202014332567 3.939692620785910 -0.5)
( 0.34202014332567 4.667633089318313 -0.5)
( 0 3.727940468532405 -0.5)
( 2 3 -0.5)
( 0 3 -0.5)
( 3.46190220008154 -2.3 -0.5)
( 3.46190220008154 2.3 -0.5)
( 0 -3 0.5)
( 2 -3 0.5)
( 0 -3.727940468532405 0.5)
( 0.34202014332567 -4.667633089318313 0.5)
( 2.34202014332567 -3.939692620785910 0.5)
( 4.92380440016309 -3 0.5)
(41 -3 0.5)
(41 -6 0.5)
(53 -6 0.5)
(53 -3 0.5)
(53 3 0.5)
(53 6 0.5)
(41 6 0.5)
(41 3 0.5)
( 4.92380440016309 3 0.5)
( 2.34202014332567 3.939692620785910 0.5)
( 0.34202014332567 4.667633089318313 0.5)
( 0 3.727940468532405 0.5)
( 2 3 0.5)
( 0 3 0.5)
( 3.46190220008154 -2.3 0.5)
( 3.46190220008154 2.3 0.5)
);
blocks
(
hex ( 0 1 18 19 22 23 40 41) ( 4 34 2) simpleGrading (1 1 1)
hex ( 1 20 21 18 23 42 43 40) ( 7 34 2) simpleGrading (1 1 1)
hex (20 5 14 21 42 27 36 43) ( 7 34 2) simpleGrading (1 1 1)
hex ( 5 6 13 14 27 28 35 36) (72 34 2) simpleGrading (1 1 1)
hex ( 6 9 10 13 28 31 32 35) (24 34 2) simpleGrading (1 1 1)
hex ( 7 8 9 6 29 30 31 28) (24 17 2) simpleGrading (1 1 1)
hex (13 10 11 12 35 32 33 34) (24 17 2) simpleGrading (1 1 1)
hex ( 3 4 1 2 25 26 23 24) ( 7 7 2) simpleGrading (1 1 1)
hex ( 4 5 20 1 26 27 42 23) ( 7 7 2) simpleGrading (1 1 1)
hex (18 21 14 15 40 43 36 37) ( 7 7 2) simpleGrading (1 1 1)
hex (17 18 15 16 39 40 37 38) ( 7 7 2) simpleGrading (1 1 1)
);
edges
(
);
boundary
(
walls
{
type wall;
faces
(
( 0 1 23 22)
( 1 2 24 23)
( 3 4 26 25)
( 4 5 27 26)
( 5 6 28 27)
( 6 7 29 28)
( 7 8 30 29)
(11 12 34 33)
(12 13 35 34)
(13 14 36 35)
(14 15 37 36)
(15 16 38 37)
(17 18 40 39)
(18 19 41 40)
(19 0 22 41)
);
}
lowerInlet
{
type patch;
faces
(
(2 3 25 24)
);
}
upperInlet
{
type patch;
faces
(
(16 17 39 38)
);
}
outlet
{
type patch;
faces
(
(8 9 31 30)
(9 10 32 31)
(10 11 33 32)
);
}
frontAndBack
{
type symmetry;
faces
(
( 0 1 18 19)
( 1 20 21 18)
(20 5 14 21)
( 5 6 13 14)
( 6 9 10 13)
( 7 8 9 6)
(13 10 11 12)
( 3 4 1 2)
( 4 5 20 1)
(18 21 14 15)
(17 18 15 16)
(22 23 40 41)
(23 42 43 40)
(42 27 36 43)
(27 28 35 36)
(28 31 32 35)
(29 30 31 28)
(35 32 33 34)
(25 26 23 24)
(26 27 42 23)
(40 43 36 37)
(39 40 37 38)
);
}
);
// ************************************************************************* //

View File

@ -0,0 +1,25 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
contiuousPhaseName air;
rho.air 1.2;
transportModel Newtonian;
nu 1e-05;
// ************************************************************************* //

View File

@ -0,0 +1,20 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object turbulenceProperties.air;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
simulationType laminar;
// ************************************************************************* //

View File

@ -0,0 +1,48 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
application MPPICFoam;
startFrom startTime;
startTime 0.0;
stopAt endTime;
endTime 0.1;
deltaT 1e-4;
writeControl runTime;
writeInterval 0.001;
purgeWrite 0;
writeFormat ascii;
writePrecision 6;
writeCompression uncompressed;
timeFormat general;
timePrecision 6;
runTimeModifiable yes;
// ************************************************************************* //

View File

@ -0,0 +1,56 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
default Euler;
}
gradSchemes
{
default Gauss linear;
}
divSchemes
{
default none;
div(alphaPhic,U.air) Gauss linearUpwindV unlimited;
div(((alpha.air*nuEff)*dev2(T(grad(U.air))))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear corrected;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default corrected;
}
fluxRequired
{
default no;
p;
}
// ************************************************************************* //

View File

@ -0,0 +1,75 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: dev |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
p
{
solver GAMG;
tolerance 1e-06;
relTol 0.01;
smoother GaussSeidel;
cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}
pFinal
{
solver GAMG;
tolerance 1e-06;
relTol 0;
smoother GaussSeidel;
cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}
"(U|k|omega).*"
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0.1;
}
"(U|k|omega).*Final"
{
$U;
tolerance 1e-05;
relTol 0;
}
}
PIMPLE
{
nOuterCorrectors 1;
nCorrectors 2;
momentumPredictor yes;
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefValue 0;
}
relaxationFactors
{
}
// ************************************************************************* //