Updated tetFem motion tutorial

This commit is contained in:
Hrvoje Jasak 2010-10-13 14:13:10 +01:00
parent 325350fbe5
commit 0c8297c859
9 changed files with 227 additions and 71 deletions

View file

@ -23,6 +23,7 @@ dynamicFvMesh dynamicMotionSolverFvMesh;
solver laplaceFaceDecomposition;
diffusivity file;
// diffusivity uniform;
frozenDiffusion off;

View file

@ -1,4 +1,4 @@
/*---------------------------------------------------------------------------*\
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.4 |
@ -10,12 +10,6 @@ FoamFile
{
version 2.0;
format ascii;
root "";
case "";
instance "";
local "";
class dictionary;
object blockMeshDict;
}
@ -50,11 +44,11 @@ vertices
blocks
(
hex (0 1 5 4 0 1 13 12) (15 15 1) simpleGrading (1 1 1)
hex (2 3 7 6 2 3 15 14) (20 20 1) simpleGrading (2 0.25 1)
hex (4 5 9 8 12 13 17 16) (15 15 1) simpleGrading (1 1 1)
hex (5 6 10 9 13 14 18 17) (50 15 1) simpleGrading (1 1 1)
hex (6 7 11 10 14 15 19 18) (20 15 1) simpleGrading (2 1 1)
hex (0 1 5 4 0 1 13 12) (5 5 1) simpleGrading (1 1 1)
hex (2 3 7 6 2 3 15 14) (7 7 1) simpleGrading (2 0.25 1)
hex (4 5 9 8 12 13 17 16) (5 5 1) simpleGrading (1 1 1)
hex (5 6 10 9 13 14 18 17) (15 5 1) simpleGrading (1 1 1)
hex (6 7 11 10 14 15 19 18) (7 5 1) simpleGrading (2 1 1)
);
edges

View file

@ -0,0 +1,112 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / 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;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
convertToMeters 0.001;
vertices
(
(-7.5 0 0)
(-7 0 0)
(-3.5 0 0)
(0 0 0)
(-7.5 0.75 -0.0327457)
(-7 0.75 -0.0327457)
(-3.5 2 -0.0873219)
(0 2 -0.0873219)
(-7.5 2.5 -0.109152)
(-7 2.5 -0.109152)
(-3.5 2.5 -0.109152)
(0 2.5 -0.109152)
(-7.5 0.75 0.0327457)
(-7 0.75 0.0327457)
(-3.5 2 0.0873219)
(0 2 0.0873219)
(-7.5 2.5 0.109152)
(-7 2.5 0.109152)
(-3.5 2.5 0.109152)
(0 2.5 0.109152)
);
blocks
(
hex (0 1 5 4 0 1 13 12) (15 15 1) simpleGrading (1 1 1)
hex (2 3 7 6 2 3 15 14) (20 20 1) simpleGrading (2 0.25 1)
hex (4 5 9 8 12 13 17 16) (15 15 1) simpleGrading (1 1 1)
hex (5 6 10 9 13 14 18 17) (50 15 1) simpleGrading (1 1 1)
hex (6 7 11 10 14 15 19 18) (20 15 1) simpleGrading (2 1 1)
);
edges
(
);
patches
(
wall movingWall
(
(1 5 13 1)
(5 6 14 13)
(2 2 14 6)
)
patch farFieldMoving
(
(9 17 18 10)
)
wall fixedWall
(
(3 7 15 3)
(7 11 19 15)
)
empty axis
(
(0 1 1 0)
(2 3 3 2)
)
patch left
(
(0 0 12 4)
(4 12 16 8)
)
patch farField
(
(8 16 17 9)
(10 18 19 11)
)
wedge back
(
(0 4 5 1)
(2 6 7 3)
(4 8 9 5)
(5 9 10 6)
(6 10 11 7)
)
wedge front
(
(0 1 13 12)
(2 3 15 14)
(12 13 17 16)
(13 14 18 17)
(14 15 19 18)
)
);
mergePatchPairs
(
);
// ************************************************************************* //

View file

@ -20,50 +20,50 @@ FoamFile
movingWall
{
type wall;
nFaces 85;
startFace 3665;
nFaces 27;
startFace 374;
}
farFieldMoving
{
type patch;
nFaces 50;
startFace 3750;
nFaces 15;
startFace 401;
}
fixedWall
{
type wall;
nFaces 35;
startFace 3800;
nFaces 12;
startFace 416;
}
axis
{
type empty;
nFaces 0;
startFace 3835;
startFace 428;
}
left
{
type patch;
nFaces 30;
startFace 3835;
nFaces 10;
startFace 428;
}
farField
{
type patch;
nFaces 35;
startFace 3865;
nFaces 12;
startFace 438;
}
back
{
type wedge;
nFaces 1900;
startFace 3900;
nFaces 209;
startFace 450;
}
front
{
type wedge;
nFaces 1900;
startFace 5800;
nFaces 209;
startFace 659;
}
)

View file

@ -1,4 +1,4 @@
/*---------------------------------------------------------------------------*\
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.4 |
@ -10,12 +10,6 @@ FoamFile
{
version 2.0;
format ascii;
root "";
case "";
instance "";
local "";
class dictionary;
object transportProperties;
}

View file

@ -26,11 +26,11 @@ stopAt endTime;
endTime 0.003;
deltaT 5e-06;
deltaT 2e-05;
writeControl timeStep;
writeInterval 40;
writeInterval 1;
purgeWrite 0;
@ -38,7 +38,7 @@ writeFormat ascii;
writePrecision 6;
writeCompression uncompressed;
writeCompression compressed;
timeFormat general;
@ -50,4 +50,5 @@ adjustTimeStep no;
maxCo 0.2;
// ************************************************************************* //

View file

@ -0,0 +1,50 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object decomposeParDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
numberOfSubdomains 4;
method simple;
simpleCoeffs
{
n (2 2 1);
delta 0.001;
}
hierarchicalCoeffs
{
n (1 1 1);
delta 0.001;
order xyz;
}
metisCoeffs
{
processorWeights 4{1};
}
manualCoeffs
{
dataFile "";
}
distributed no;
roots
(
);
// ************************************************************************* //

View file

@ -10,12 +10,6 @@ FoamFile
{
version 2.0;
format ascii;
root "";
case "";
instance "";
local "";
class dictionary;
object fvSchemes;
}

View file

@ -20,22 +20,32 @@ solvers
{
motionU
{
solver amgSolver;
cycle W-cycle;
policy AAMG;
nPreSweeps 0;
nPostSweeps 2;
groupSize 4;
minCoarseEqns 30;
nMaxLevels 100;
scale on;
smoother GaussSeidel;
solver PCG;
preconditioner DIC;
minIter 0;
maxIter 500;
tolerance 1e-8;
relTol 0.0;
maxIter 1000;
tolerance 1e-6;
relTol 0.01;
};
// motionU
// {
// solver amgSolver;
// cycle W-cycle;
// policy AAMG;
// nPreSweeps 0;
// nPostSweeps 2;
// groupSize 4;
// minCoarseEqns 30;
// nMaxLevels 100;
// scale on;
// smoother GaussSeidel;
// minIter 0;
// maxIter 500;
// tolerance 1e-8;
// relTol 0.0;
// };
}