exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

New tutorial: moving mesh turbulent immersed boundary with adaptive mesh refinement, 2-D

Browse source
This commit is contained in:
Hrvoje Jasak 2018-09-06 18:27:10 +01:00
parent e38c7e68ec
commit 007a935a81
27 changed files with 4118 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

52
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/.gitignore vendored Normal file
View file

@ -0,0 +1,52 @@
# git-ls-files --others --exclude-from=.git/info/exclude
# Lines that start with '#' are comments.
# editor and misc backup files - anywhere
*~
.*~
*.bak
*.bak[0-9][0-9]
*.orig
*.orig[0-9][0-9]
\#*\#
# file-browser settings - anywhere
.directory
# CVS recovered versions - anywhere
.#*
# SVN directories - anywhere
.svn/
# OpenFOAM results
[1-9]*/
!/0/
processor*
*/polyMesh/*
!*/polyMesh/blockMeshDict
cellToRegion*
log*
# packages - anywhere
*.tar.bz2
*.tar.gz
*.tar
*.tgz
*.gtgz
# Pictures and movies
*.png
*.jpg
*.jpeg
*.bmp
*.png
*.avi
*.mp4
*.mpg
#end-of-file

59
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/0_org/U Normal file
View file

@ -0,0 +1,59 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (1 0 0);
boundaryField
{
ibCylinder
{
type movingImmersedBoundaryVelocity;
patchType immersedBoundary;
setDeadValue yes;
deadValue (0 0 0);
value uniform (0 0 0);
}
in
{
type fixedValue;
value uniform (1 0 0);
}
out
{
type inletOutlet;
inletValue uniform (0 0 0);
value uniform (1 0 0);
}
top
{
type fixedValue;
value uniform (0 0 0);
}
bottom
{
type fixedValue;
value uniform (0 0 0);
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //

62
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/0_org/epsilon Normal file
View file

@ -0,0 +1,62 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -3 0 0 0 0];
internalField uniform 2e-4;
boundaryField
{
ibCylinder
{
type immersedBoundaryEpsilonWallFunction;
patchType immersedBoundary;
Cmu 0.09;
kappa 0.41;
E 9.8;
setDeadValue true;
deadValue 2e-4;
value uniform 2e-4;
}
top
{
type zeroGradient;
}
in
{
type fixedValue;
value uniform 2e-4;
}
out
{
type inletOutlet;
inletValue uniform 2e-4;
value uniform 2e-4;
}
bottom
{
type zeroGradient;
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //

58
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/0_org/k Normal file
View file

@ -0,0 +1,58 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: http://www.foam-extend.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 1.5e-4;
boundaryField
{
ibCylinder
{
type immersedBoundaryKqRWallFunction;
patchType immersedBoundary;
setDeadValue true;
deadValue 1.5e-4;
value uniform 1.5e-4;
}
top
{
type zeroGradient;
}
in
{
type fixedValue;
value uniform 1.5e-4;
}
out
{
type inletOutlet;
inletValue uniform 1.5e-4;
value uniform 1.5e-4;
}
bottom
{
type zeroGradient;
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //

59
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/0_org/nut Normal file
View file

@ -0,0 +1,59 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: http://www.foam-extend.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
{
ibCylinder
{
type immersedBoundaryNutWallFunction;
patchType immersedBoundary;
setDeadValue true;
deadValue 1e-5;
value uniform 0;
}
top
{
type zeroGradient;
value uniform 0;
}
in
{
type zeroGradient;
value uniform 0;
}
out
{
type zeroGradient;
value uniform 0;
}
bottom
{
type zeroGradient;
value uniform 0;
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //

60
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/0_org/omega Normal file
View file

@ -0,0 +1,60 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 -1 0 0 0 0];
internalField uniform 500;
boundaryField
{
ibCylinder
{
type immersedBoundaryOmegaWallFunction;
patchType immersedBoundary;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 500;
}
top
{
type zeroGradient;
}
in
{
type fixedValue;
value uniform 500;
}
out
{
type inletOutlet;
inletValue uniform 500;
value uniform 500;
}
bottom
{
type zeroGradient;
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //

58
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/0_org/p Normal file
View file

@ -0,0 +1,58 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: http://www.foam-extend.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
{
ibCylinder
{
type mixedIb;
patchType immersedBoundary;
triValue uniform 0;
triGradient uniform 0;
triValueFraction uniform 0;
setDeadValue no;
deadValue 0;
value uniform 0;
}
in
{
type zeroGradient;
}
out
{
type fixedValue;
value uniform 0;
}
top
{
type zeroGradient;
}
bottom
{
type zeroGradient;
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //

11
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/Allclean Executable file
View file

@ -0,0 +1,11 @@
#!/bin/sh
# Source tutorial clean functions
. $WM_PROJECT_DIR/bin/tools/CleanFunctions
cleanCase
\rm -f *.dat
\rm -rf 0
\rm -f constant/polyMesh/boundary
\rm -rf postProcessing

16
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/Allrun Executable file
View file

@ -0,0 +1,16 @@
#!/bin/sh
# Source tutorial run functions
. $WM_PROJECT_DIR/bin/tools/RunFunctions
# Get application name
application="pimpleDyMIbFoam"
runApplication blockMesh
\cp save/boundary constant/polyMesh/
mkdir 0
\cp -f 0_org/* 0/
runApplication potentialFoam
runApplication writeIbMasks
runApplication $application

26
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/constant/RASProperties Normal file
View file

@ -0,0 +1,26 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object RASProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
RASModel kEpsilon;
// RASModel kOmegaSST;
turbulence on;
printCoeffs on;
nuRatio 1000;
// ************************************************************************* //

83
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/constant/dynamicMeshDict Normal file
View file

@ -0,0 +1,83 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object dynamicMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
// dynamicFvMesh staticFvMesh;
// dynamicFvMesh immersedBoundarySolidBodyMotionFvMesh;
dynamicFvMesh immersedBoundaryDynamicRefineSolidBodyMotionFvMesh;
immersedBoundaryDynamicRefineSolidBodyMotionFvMeshCoeffs
{
// Immersed boundary controls
motionFunctions
(
ibCylinder
{
solidBodyMotionFunction linearOscillation;
linearOscillationCoeffs
{
amplitude (0.5 0 0);
period 2.5;
}
}
);
// Dynamic mesh procedure controls
// Refine every refineInterval step
refineInterval 2;
// Unrefine every unrefineInterval step
unrefineInterval 2;
// Separate refinement/unrefinement steps. In case this is switched on,
// if both refinement and unrefinement should have been performed in a
// single step, unrefinement is skipped. Switched off by default, meaning
// that it should be safe to perform both at the same time
separateUpdates false;
// Refinement selection criteria
refinementSelection
{
type immersedBoundaryRefinement;
refinementDistance 0.08;
unrefinementDistance 0.08;
}
// Polyhedral refinement engine controls
active yes;
// Maximum number of cells to allow (not strictly controlled)
maxCells 20000;
// Maximum refinement level
maxRefinementLevel 2;
// Number of buffer layers between refinement levels
nRefinementBufferLayers 1;
// Number of buffer layers for unrefinement in order to run away from the
// region that is getting refined at the same time in order to avoid point
// level inconsistencies
nUnrefinementBufferLayers 3;
// Whether to use edge based consistency check. Needed when one allows more
// than 2 refinement levels (automatically switched on)
edgeBasedConsistency true;
}
// ************************************************************************* //

70
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/constant/polyMesh/blockMeshDict Normal file
View file

@ -0,0 +1,70 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
convertToMeters 1;
vertices
(
(-1 -0.5 -0.04)
( 2 -0.5 -0.04)
( 2 0.5 -0.04)
(-1 0.5 -0.04)
(-1 -0.5 0.04)
( 2 -0.5 0.04)
( 2 0.5 0.04)
(-1 0.5 0.04)
);
blocks
(
hex (0 1 2 3 4 5 6 7) (75 25 1) simpleGrading (1 1 1)
);
edges
(
);
patches
(
patch in
(
(0 4 7 3)
)
patch out
(
(2 6 5 1)
)
patch top
(
(3 7 6 2)
)
patch bottom
(
(1 5 4 0)
)
empty frontAndBack
(
(0 3 2 1)
(4 5 6 7)
)
);
mergePatchPairs
(
);
// ************************************************************************* //

21
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/constant/transportProperties Normal file
View file

@ -0,0 +1,21 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
transportModel Newtonian;
nu nu [0 2 -1 0 0 0 0] 1e-6;
// ************************************************************************* //

210
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/constant/triSurface/ibCylinder.ftr Normal file
View file

@ -0,0 +1,210 @@
1
(
cylinder
patch
)
96
(
(0.176776 -0.176776 -0.05)
(0.176776 -0.176776 0.05)
(0.176776 0.176776 -0.05)
(0.176776 0.176776 0.05)
(-0.176776 0.176776 -0.05)
(-0.176776 0.176776 0.05)
(-0.176776 -0.176776 -0.05)
(-0.176776 -0.176776 0.05)
(0.241481 -0.064705 -0.05)
(0.241481 -0.064705 0.05)
(0.241481 0.064705 -0.05)
(0.241481 0.064705 0.05)
(0.064705 0.241481 -0.05)
(0.064705 0.241481 0.05)
(-0.064705 0.241481 -0.05)
(-0.064705 0.241481 0.05)
(-0.241481 0.064705 -0.05)
(-0.241481 0.064705 0.05)
(-0.241481 -0.064705 -0.05)
(-0.241481 -0.064705 0.05)
(0.064705 -0.241481 0.05)
(0.064705 -0.241481 -0.05)
(-0.064705 -0.241481 0.05)
(-0.064705 -0.241481 -0.05)
(0.216506 -0.125 -0.05)
(0.216506 -0.125 0.05)
(0.216506 0.125 -0.05)
(0.216506 0.125 0.05)
(0.125 0.216506 -0.05)
(0.125 0.216506 0.05)
(-0.125 0.216506 -0.05)
(-0.125 0.216506 0.05)
(-0.216506 0.125 -0.05)
(-0.216506 0.125 0.05)
(-0.216506 -0.125 -0.05)
(-0.216506 -0.125 0.05)
(0.125 -0.216506 0.05)
(0.125 -0.216506 -0.05)
(-0.125 -0.216506 0.05)
(-0.125 -0.216506 -0.05)
(0.247861 -0.0326315 -0.05)
(0.247861 -0.0326315 0.05)
(0.247861 0.0326315 -0.05)
(0.247861 0.0326315 0.05)
(0.0326315 0.247861 -0.05)
(0.0326315 0.247861 0.05)
(-0.0326315 0.247861 -0.05)
(-0.0326315 0.247861 0.05)
(-0.247861 0.0326315 -0.05)
(-0.247861 0.0326315 0.05)
(-0.247861 -0.0326315 -0.05)
(-0.247861 -0.0326315 0.05)
(0.0326315 -0.247861 0.05)
(0.0326315 -0.247861 -0.05)
(-0.0326315 -0.247861 0.05)
(-0.0326315 -0.247861 -0.05)
(0.25 -8.32665e-17 0.05)
(0.25 8.32665e-17 -0.05)
(8.32665e-17 0.25 0.05)
(-8.32665e-17 0.25 -0.05)
(-0.25 8.32665e-17 0.05)
(-0.25 -8.32665e-17 -0.05)
(8.32665e-17 -0.25 -0.05)
(-8.32665e-17 -0.25 0.05)
(0.198338 -0.152191 -0.05)
(0.198338 -0.152191 0.05)
(0.198338 0.152191 -0.05)
(0.198338 0.152191 0.05)
(0.152191 0.198338 -0.05)
(0.152191 0.198338 0.05)
(-0.152191 0.198338 -0.05)
(-0.152191 0.198338 0.05)
(-0.198338 0.152191 -0.05)
(-0.198338 0.152191 0.05)
(-0.198338 -0.152191 -0.05)
(-0.198338 -0.152191 0.05)
(0.152191 -0.198338 0.05)
(0.152191 -0.198338 -0.05)
(-0.152191 -0.198338 0.05)
(-0.152191 -0.198338 -0.05)
(0.23097 -0.095671 -0.05)
(0.23097 -0.095671 0.05)
(0.23097 0.095671 -0.05)
(0.23097 0.095671 0.05)
(0.095671 0.23097 -0.05)
(0.095671 0.23097 0.05)
(-0.095671 0.23097 -0.05)
(-0.095671 0.23097 0.05)
(-0.23097 0.095671 -0.05)
(-0.23097 0.095671 0.05)
(-0.23097 -0.095671 -0.05)
(-0.23097 -0.095671 0.05)
(0.095671 -0.23097 0.05)
(0.095671 -0.23097 -0.05)
(-0.095671 -0.23097 0.05)
(-0.095671 -0.23097 -0.05)
)
96
(
((0 64 65) 0)
((65 1 0) 0)
((64 24 25) 0)
((25 65 64) 0)
((24 80 81) 0)
((81 25 24) 0)
((80 8 9) 0)
((9 81 80) 0)
((8 40 41) 0)
((41 9 8) 0)
((56 41 40) 0)
((40 57 56) 0)
((57 42 43) 0)
((43 56 57) 0)
((42 10 11) 0)
((11 43 42) 0)
((10 82 83) 0)
((83 11 10) 0)
((82 26 27) 0)
((27 83 82) 0)
((26 66 67) 0)
((67 27 26) 0)
((66 2 3) 0)
((3 67 66) 0)
((2 68 69) 0)
((69 3 2) 0)
((68 28 29) 0)
((29 69 68) 0)
((28 84 85) 0)
((85 29 28) 0)
((84 12 13) 0)
((13 85 84) 0)
((12 44 45) 0)
((45 13 12) 0)
((58 45 44) 0)
((44 59 58) 0)
((59 46 47) 0)
((47 58 59) 0)
((46 14 15) 0)
((15 47 46) 0)
((14 86 87) 0)
((87 15 14) 0)
((86 30 31) 0)
((31 87 86) 0)
((30 70 71) 0)
((71 31 30) 0)
((70 4 5) 0)
((5 71 70) 0)
((4 72 73) 0)
((73 5 4) 0)
((72 32 33) 0)
((33 73 72) 0)
((32 88 89) 0)
((89 33 32) 0)
((88 16 17) 0)
((17 89 88) 0)
((16 48 49) 0)
((49 17 16) 0)
((60 49 48) 0)
((48 61 60) 0)
((61 50 51) 0)
((51 60 61) 0)
((50 18 19) 0)
((19 51 50) 0)
((18 90 91) 0)
((91 19 18) 0)
((90 34 35) 0)
((35 91 90) 0)
((34 74 75) 0)
((75 35 34) 0)
((74 6 7) 0)
((7 75 74) 0)
((0 1 76) 0)
((76 77 0) 0)
((77 76 36) 0)
((36 37 77) 0)
((37 36 92) 0)
((92 93 37) 0)
((93 92 20) 0)
((20 21 93) 0)
((21 20 52) 0)
((52 53 21) 0)
((62 53 52) 0)
((52 63 62) 0)
((63 54 55) 0)
((55 62 63) 0)
((55 54 22) 0)
((22 23 55) 0)
((23 22 94) 0)
((94 95 23) 0)
((95 94 38) 0)
((38 39 95) 0)
((39 38 78) 0)
((78 79 39) 0)
((79 78 7) 0)
((7 6 79) 0)
)

674
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/constant/triSurface/ibCylinder.stl Normal file
View file

@ -0,0 +1,674 @@
solid cylinder
facet normal 0.751817 -0.659372 0
outer loop
vertex 0.176776 -0.176776 -0.05
vertex 0.198338 -0.152191 -0.05
vertex 0.198338 -0.152191 0.05
endloop
endfacet
facet normal 0.751817 -0.659372 0
outer loop
vertex 0.198338 -0.152191 0.05
vertex 0.176776 -0.176776 0.05
vertex 0.176776 -0.176776 -0.05
endloop
endfacet
facet normal 0.831476 -0.555561 0
outer loop
vertex 0.198338 -0.152191 -0.05
vertex 0.216506 -0.125 -0.05
vertex 0.216506 -0.125 0.05
endloop
endfacet
facet normal 0.831476 -0.555561 0
outer loop
vertex 0.216506 -0.125 0.05
vertex 0.198338 -0.152191 0.05
vertex 0.198338 -0.152191 -0.05
endloop
endfacet
facet normal 0.896866 -0.442302 0
outer loop
vertex 0.216506 -0.125 -0.05
vertex 0.23097 -0.095671 -0.05
vertex 0.23097 -0.095671 0.05
endloop
endfacet
facet normal 0.896866 -0.442302 0
outer loop
vertex 0.23097 -0.095671 0.05
vertex 0.216506 -0.125 0.05
vertex 0.216506 -0.125 -0.05
endloop
endfacet
facet normal 0.946935 -0.321425 0
outer loop
vertex 0.23097 -0.095671 -0.05
vertex 0.241481 -0.064705 -0.05
vertex 0.241481 -0.064705 0.05
endloop
endfacet
facet normal 0.946935 -0.321425 0
outer loop
vertex 0.241481 -0.064705 0.05
vertex 0.23097 -0.095671 0.05
vertex 0.23097 -0.095671 -0.05
endloop
endfacet
facet normal 0.980784 -0.195096 0
outer loop
vertex 0.241481 -0.064705 -0.05
vertex 0.247861 -0.0326315 -0.05
vertex 0.247861 -0.0326315 0.05
endloop
endfacet
facet normal 0.980784 -0.195096 0
outer loop
vertex 0.247861 -0.0326315 0.05
vertex 0.241481 -0.064705 0.05
vertex 0.241481 -0.064705 -0.05
endloop
endfacet
facet normal 0.997858 -0.0654098 0
outer loop
vertex 0.25 -8.32665e-17 0.05
vertex 0.247861 -0.0326315 0.05
vertex 0.247861 -0.0326315 -0.05
endloop
endfacet
facet normal 0.997858 -0.0654098 -1.07752e-16
outer loop
vertex 0.247861 -0.0326315 -0.05
vertex 0.25 8.32665e-17 -0.05
vertex 0.25 -8.32665e-17 0.05
endloop
endfacet
facet normal 0.997858 0.0654098 0
outer loop
vertex 0.25 8.32665e-17 -0.05
vertex 0.247861 0.0326315 -0.05
vertex 0.247861 0.0326315 0.05
endloop
endfacet
facet normal 0.997858 0.0654098 1.07752e-16
outer loop
vertex 0.247861 0.0326315 0.05
vertex 0.25 -8.32665e-17 0.05
vertex 0.25 8.32665e-17 -0.05
endloop
endfacet
facet normal 0.980784 0.195096 0
outer loop
vertex 0.247861 0.0326315 -0.05
vertex 0.241481 0.064705 -0.05
vertex 0.241481 0.064705 0.05
endloop
endfacet
facet normal 0.980784 0.195096 0
outer loop
vertex 0.241481 0.064705 0.05
vertex 0.247861 0.0326315 0.05
vertex 0.247861 0.0326315 -0.05
endloop
endfacet
facet normal 0.946935 0.321425 0
outer loop
vertex 0.241481 0.064705 -0.05
vertex 0.23097 0.095671 -0.05
vertex 0.23097 0.095671 0.05
endloop
endfacet
facet normal 0.946935 0.321425 0
outer loop
vertex 0.23097 0.095671 0.05
vertex 0.241481 0.064705 0.05
vertex 0.241481 0.064705 -0.05
endloop
endfacet
facet normal 0.896866 0.442302 0
outer loop
vertex 0.23097 0.095671 -0.05
vertex 0.216506 0.125 -0.05
vertex 0.216506 0.125 0.05
endloop
endfacet
facet normal 0.896866 0.442302 0
outer loop
vertex 0.216506 0.125 0.05
vertex 0.23097 0.095671 0.05
vertex 0.23097 0.095671 -0.05
endloop
endfacet
facet normal 0.831476 0.555561 0
outer loop
vertex 0.216506 0.125 -0.05
vertex 0.198338 0.152191 -0.05
vertex 0.198338 0.152191 0.05
endloop
endfacet
facet normal 0.831476 0.555561 0
outer loop
vertex 0.198338 0.152191 0.05
vertex 0.216506 0.125 0.05
vertex 0.216506 0.125 -0.05
endloop
endfacet
facet normal 0.751817 0.659372 0
outer loop
vertex 0.198338 0.152191 -0.05
vertex 0.176776 0.176776 -0.05
vertex 0.176776 0.176776 0.05
endloop
endfacet
facet normal 0.751817 0.659372 0
outer loop
vertex 0.176776 0.176776 0.05
vertex 0.198338 0.152191 0.05
vertex 0.198338 0.152191 -0.05
endloop
endfacet
facet normal 0.659372 0.751817 0
outer loop
vertex 0.176776 0.176776 -0.05
vertex 0.152191 0.198338 -0.05
vertex 0.152191 0.198338 0.05
endloop
endfacet
facet normal 0.659372 0.751817 0
outer loop
vertex 0.152191 0.198338 0.05
vertex 0.176776 0.176776 0.05
vertex 0.176776 0.176776 -0.05
endloop
endfacet
facet normal 0.555561 0.831476 0
outer loop
vertex 0.152191 0.198338 -0.05
vertex 0.125 0.216506 -0.05
vertex 0.125 0.216506 0.05
endloop
endfacet
facet normal 0.555561 0.831476 0
outer loop
vertex 0.125 0.216506 0.05
vertex 0.152191 0.198338 0.05
vertex 0.152191 0.198338 -0.05
endloop
endfacet
facet normal 0.442302 0.896866 0
outer loop
vertex 0.125 0.216506 -0.05
vertex 0.095671 0.23097 -0.05
vertex 0.095671 0.23097 0.05
endloop
endfacet
facet normal 0.442302 0.896866 0
outer loop
vertex 0.095671 0.23097 0.05
vertex 0.125 0.216506 0.05
vertex 0.125 0.216506 -0.05
endloop
endfacet
facet normal 0.321425 0.946935 0
outer loop
vertex 0.095671 0.23097 -0.05
vertex 0.064705 0.241481 -0.05
vertex 0.064705 0.241481 0.05
endloop
endfacet
facet normal 0.321425 0.946935 0
outer loop
vertex 0.064705 0.241481 0.05
vertex 0.095671 0.23097 0.05
vertex 0.095671 0.23097 -0.05
endloop
endfacet
facet normal 0.195096 0.980784 0
outer loop
vertex 0.064705 0.241481 -0.05
vertex 0.0326315 0.247861 -0.05
vertex 0.0326315 0.247861 0.05
endloop
endfacet
facet normal 0.195096 0.980784 0
outer loop
vertex 0.0326315 0.247861 0.05
vertex 0.064705 0.241481 0.05
vertex 0.064705 0.241481 -0.05
endloop
endfacet
facet normal 0.0654098 0.997858 0
outer loop
vertex 8.32665e-17 0.25 0.05
vertex 0.0326315 0.247861 0.05
vertex 0.0326315 0.247861 -0.05
endloop
endfacet
facet normal 0.0654098 0.997858 -1.07752e-16
outer loop
vertex 0.0326315 0.247861 -0.05
vertex -8.32665e-17 0.25 -0.05
vertex 8.32665e-17 0.25 0.05
endloop
endfacet
facet normal -0.0654098 0.997858 0
outer loop
vertex -8.32665e-17 0.25 -0.05
vertex -0.0326315 0.247861 -0.05
vertex -0.0326315 0.247861 0.05
endloop
endfacet
facet normal -0.0654098 0.997858 1.07752e-16
outer loop
vertex -0.0326315 0.247861 0.05
vertex 8.32665e-17 0.25 0.05
vertex -8.32665e-17 0.25 -0.05
endloop
endfacet
facet normal -0.195096 0.980784 0
outer loop
vertex -0.0326315 0.247861 -0.05
vertex -0.064705 0.241481 -0.05
vertex -0.064705 0.241481 0.05
endloop
endfacet
facet normal -0.195096 0.980784 0
outer loop
vertex -0.064705 0.241481 0.05
vertex -0.0326315 0.247861 0.05
vertex -0.0326315 0.247861 -0.05
endloop
endfacet
facet normal -0.321425 0.946935 0
outer loop
vertex -0.064705 0.241481 -0.05
vertex -0.095671 0.23097 -0.05
vertex -0.095671 0.23097 0.05
endloop
endfacet
facet normal -0.321425 0.946935 0
outer loop
vertex -0.095671 0.23097 0.05
vertex -0.064705 0.241481 0.05
vertex -0.064705 0.241481 -0.05
endloop
endfacet
facet normal -0.442302 0.896866 0
outer loop
vertex -0.095671 0.23097 -0.05
vertex -0.125 0.216506 -0.05
vertex -0.125 0.216506 0.05
endloop
endfacet
facet normal -0.442302 0.896866 0
outer loop
vertex -0.125 0.216506 0.05
vertex -0.095671 0.23097 0.05
vertex -0.095671 0.23097 -0.05
endloop
endfacet
facet normal -0.555561 0.831476 0
outer loop
vertex -0.125 0.216506 -0.05
vertex -0.152191 0.198338 -0.05
vertex -0.152191 0.198338 0.05
endloop
endfacet
facet normal -0.555561 0.831476 0
outer loop
vertex -0.152191 0.198338 0.05
vertex -0.125 0.216506 0.05
vertex -0.125 0.216506 -0.05
endloop
endfacet
facet normal -0.659372 0.751817 0
outer loop
vertex -0.152191 0.198338 -0.05
vertex -0.176776 0.176776 -0.05
vertex -0.176776 0.176776 0.05
endloop
endfacet
facet normal -0.659372 0.751817 0
outer loop
vertex -0.176776 0.176776 0.05
vertex -0.152191 0.198338 0.05
vertex -0.152191 0.198338 -0.05
endloop
endfacet
facet normal -0.751817 0.659372 0
outer loop
vertex -0.176776 0.176776 -0.05
vertex -0.198338 0.152191 -0.05
vertex -0.198338 0.152191 0.05
endloop
endfacet
facet normal -0.751817 0.659372 0
outer loop
vertex -0.198338 0.152191 0.05
vertex -0.176776 0.176776 0.05
vertex -0.176776 0.176776 -0.05
endloop
endfacet
facet normal -0.831476 0.555561 0
outer loop
vertex -0.198338 0.152191 -0.05
vertex -0.216506 0.125 -0.05
vertex -0.216506 0.125 0.05
endloop
endfacet
facet normal -0.831476 0.555561 0
outer loop
vertex -0.216506 0.125 0.05
vertex -0.198338 0.152191 0.05
vertex -0.198338 0.152191 -0.05
endloop
endfacet
facet normal -0.896866 0.442302 0
outer loop
vertex -0.216506 0.125 -0.05
vertex -0.23097 0.095671 -0.05
vertex -0.23097 0.095671 0.05
endloop
endfacet
facet normal -0.896866 0.442302 0
outer loop
vertex -0.23097 0.095671 0.05
vertex -0.216506 0.125 0.05
vertex -0.216506 0.125 -0.05
endloop
endfacet
facet normal -0.946935 0.321425 0
outer loop
vertex -0.23097 0.095671 -0.05
vertex -0.241481 0.064705 -0.05
vertex -0.241481 0.064705 0.05
endloop
endfacet
facet normal -0.946935 0.321425 0
outer loop
vertex -0.241481 0.064705 0.05
vertex -0.23097 0.095671 0.05
vertex -0.23097 0.095671 -0.05
endloop
endfacet
facet normal -0.980784 0.195096 0
outer loop
vertex -0.241481 0.064705 -0.05
vertex -0.247861 0.0326315 -0.05
vertex -0.247861 0.0326315 0.05
endloop
endfacet
facet normal -0.980784 0.195096 0
outer loop
vertex -0.247861 0.0326315 0.05
vertex -0.241481 0.064705 0.05
vertex -0.241481 0.064705 -0.05
endloop
endfacet
facet normal -0.997858 0.0654098 0
outer loop
vertex -0.25 8.32665e-17 0.05
vertex -0.247861 0.0326315 0.05
vertex -0.247861 0.0326315 -0.05
endloop
endfacet
facet normal -0.997858 0.0654098 -1.07752e-16
outer loop
vertex -0.247861 0.0326315 -0.05
vertex -0.25 -8.32665e-17 -0.05
vertex -0.25 8.32665e-17 0.05
endloop
endfacet
facet normal -0.997858 -0.0654098 0
outer loop
vertex -0.25 -8.32665e-17 -0.05
vertex -0.247861 -0.0326315 -0.05
vertex -0.247861 -0.0326315 0.05
endloop
endfacet
facet normal -0.997858 -0.0654098 1.07752e-16
outer loop
vertex -0.247861 -0.0326315 0.05
vertex -0.25 8.32665e-17 0.05
vertex -0.25 -8.32665e-17 -0.05
endloop
endfacet
facet normal -0.980784 -0.195096 0
outer loop
vertex -0.247861 -0.0326315 -0.05
vertex -0.241481 -0.064705 -0.05
vertex -0.241481 -0.064705 0.05
endloop
endfacet
facet normal -0.980784 -0.195096 0
outer loop
vertex -0.241481 -0.064705 0.05
vertex -0.247861 -0.0326315 0.05
vertex -0.247861 -0.0326315 -0.05
endloop
endfacet
facet normal -0.946935 -0.321425 0
outer loop
vertex -0.241481 -0.064705 -0.05
vertex -0.23097 -0.095671 -0.05
vertex -0.23097 -0.095671 0.05
endloop
endfacet
facet normal -0.946935 -0.321425 0
outer loop
vertex -0.23097 -0.095671 0.05
vertex -0.241481 -0.064705 0.05
vertex -0.241481 -0.064705 -0.05
endloop
endfacet
facet normal -0.896866 -0.442302 0
outer loop
vertex -0.23097 -0.095671 -0.05
vertex -0.216506 -0.125 -0.05
vertex -0.216506 -0.125 0.05
endloop
endfacet
facet normal -0.896866 -0.442302 0
outer loop
vertex -0.216506 -0.125 0.05
vertex -0.23097 -0.095671 0.05
vertex -0.23097 -0.095671 -0.05
endloop
endfacet
facet normal -0.831476 -0.555561 0
outer loop
vertex -0.216506 -0.125 -0.05
vertex -0.198338 -0.152191 -0.05
vertex -0.198338 -0.152191 0.05
endloop
endfacet
facet normal -0.831476 -0.555561 0
outer loop
vertex -0.198338 -0.152191 0.05
vertex -0.216506 -0.125 0.05
vertex -0.216506 -0.125 -0.05
endloop
endfacet
facet normal -0.751817 -0.659372 0
outer loop
vertex -0.198338 -0.152191 -0.05
vertex -0.176776 -0.176776 -0.05
vertex -0.176776 -0.176776 0.05
endloop
endfacet
facet normal -0.751817 -0.659372 0
outer loop
vertex -0.176776 -0.176776 0.05
vertex -0.198338 -0.152191 0.05
vertex -0.198338 -0.152191 -0.05
endloop
endfacet
facet normal 0.659372 -0.751817 0
outer loop
vertex 0.176776 -0.176776 -0.05
vertex 0.176776 -0.176776 0.05
vertex 0.152191 -0.198338 0.05
endloop
endfacet
facet normal 0.659372 -0.751817 0
outer loop
vertex 0.152191 -0.198338 0.05
vertex 0.152191 -0.198338 -0.05
vertex 0.176776 -0.176776 -0.05
endloop
endfacet
facet normal 0.555561 -0.831476 0
outer loop
vertex 0.152191 -0.198338 -0.05
vertex 0.152191 -0.198338 0.05
vertex 0.125 -0.216506 0.05
endloop
endfacet
facet normal 0.555561 -0.831476 0
outer loop
vertex 0.125 -0.216506 0.05
vertex 0.125 -0.216506 -0.05
vertex 0.152191 -0.198338 -0.05
endloop
endfacet
facet normal 0.442302 -0.896866 0
outer loop
vertex 0.125 -0.216506 -0.05
vertex 0.125 -0.216506 0.05
vertex 0.095671 -0.23097 0.05
endloop
endfacet
facet normal 0.442302 -0.896866 0
outer loop
vertex 0.095671 -0.23097 0.05
vertex 0.095671 -0.23097 -0.05
vertex 0.125 -0.216506 -0.05
endloop
endfacet
facet normal 0.321425 -0.946935 0
outer loop
vertex 0.095671 -0.23097 -0.05
vertex 0.095671 -0.23097 0.05
vertex 0.064705 -0.241481 0.05
endloop
endfacet
facet normal 0.321425 -0.946935 0
outer loop
vertex 0.064705 -0.241481 0.05
vertex 0.064705 -0.241481 -0.05
vertex 0.095671 -0.23097 -0.05
endloop
endfacet
facet normal 0.195096 -0.980784 0
outer loop
vertex 0.064705 -0.241481 -0.05
vertex 0.064705 -0.241481 0.05
vertex 0.0326315 -0.247861 0.05
endloop
endfacet
facet normal 0.195096 -0.980784 0
outer loop
vertex 0.0326315 -0.247861 0.05
vertex 0.0326315 -0.247861 -0.05
vertex 0.064705 -0.241481 -0.05
endloop
endfacet
facet normal 0.0654098 -0.997858 0
outer loop
vertex 8.32665e-17 -0.25 -0.05
vertex 0.0326315 -0.247861 -0.05
vertex 0.0326315 -0.247861 0.05
endloop
endfacet
facet normal 0.0654098 -0.997858 1.07752e-16
outer loop
vertex 0.0326315 -0.247861 0.05
vertex -8.32665e-17 -0.25 0.05
vertex 8.32665e-17 -0.25 -0.05
endloop
endfacet
facet normal -0.0654098 -0.997858 0
outer loop
vertex -8.32665e-17 -0.25 0.05
vertex -0.0326315 -0.247861 0.05
vertex -0.0326315 -0.247861 -0.05
endloop
endfacet
facet normal -0.0654098 -0.997858 -1.07752e-16
outer loop
vertex -0.0326315 -0.247861 -0.05
vertex 8.32665e-17 -0.25 -0.05
vertex -8.32665e-17 -0.25 0.05
endloop
endfacet
facet normal -0.195096 -0.980784 0
outer loop
vertex -0.0326315 -0.247861 -0.05
vertex -0.0326315 -0.247861 0.05
vertex -0.064705 -0.241481 0.05
endloop
endfacet
facet normal -0.195096 -0.980784 -0
outer loop
vertex -0.064705 -0.241481 0.05
vertex -0.064705 -0.241481 -0.05
vertex -0.0326315 -0.247861 -0.05
endloop
endfacet
facet normal -0.321425 -0.946935 0
outer loop
vertex -0.064705 -0.241481 -0.05
vertex -0.064705 -0.241481 0.05
vertex -0.095671 -0.23097 0.05
endloop
endfacet
facet normal -0.321425 -0.946935 -0
outer loop
vertex -0.095671 -0.23097 0.05
vertex -0.095671 -0.23097 -0.05
vertex -0.064705 -0.241481 -0.05
endloop
endfacet
facet normal -0.442302 -0.896866 0
outer loop
vertex -0.095671 -0.23097 -0.05
vertex -0.095671 -0.23097 0.05
vertex -0.125 -0.216506 0.05
endloop
endfacet
facet normal -0.442302 -0.896866 -0
outer loop
vertex -0.125 -0.216506 0.05
vertex -0.125 -0.216506 -0.05
vertex -0.095671 -0.23097 -0.05
endloop
endfacet
facet normal -0.555561 -0.831476 0
outer loop
vertex -0.125 -0.216506 -0.05
vertex -0.125 -0.216506 0.05
vertex -0.152191 -0.198338 0.05
endloop
endfacet
facet normal -0.555561 -0.831476 -0
outer loop
vertex -0.152191 -0.198338 0.05
vertex -0.152191 -0.198338 -0.05
vertex -0.125 -0.216506 -0.05
endloop
endfacet
facet normal -0.659372 -0.751817 0
outer loop
vertex -0.152191 -0.198338 -0.05
vertex -0.152191 -0.198338 0.05
vertex -0.176776 -0.176776 0.05
endloop
endfacet
facet normal -0.659372 -0.751817 -0
outer loop
vertex -0.176776 -0.176776 0.05
vertex -0.176776 -0.176776 -0.05
vertex -0.152191 -0.198338 -0.05
endloop
endfacet
endsolid cylinder

21
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/constant/turbulenceProperties Normal file
View file

@ -0,0 +1,21 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object turbulenceProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
simulationType RASModel;
// ************************************************************************* //

70
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/save/blockMeshDict Normal file
View file

@ -0,0 +1,70 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
convertToMeters 1;
vertices
(
(-1 -0.5 -0.04)
( 2 -0.5 -0.04)
( 2 0.5 -0.04)
(-1 0.5 -0.04)
(-1 -0.5 0.04)
( 2 -0.5 0.04)
( 2 0.5 0.04)
(-1 0.5 0.04)
);
blocks
(
hex (0 1 2 3 4 5 6 7) (75 25 1) simpleGrading (1 1 1)
);
edges
(
);
patches
(
patch in
(
(0 4 7 3)
)
patch out
(
(2 6 5 1)
)
patch top
(
(3 7 6 2)
)
patch bottom
(
(1 5 4 0)
)
empty frontAndBack
(
(0 3 2 1)
(4 5 6 7)
)
);
mergePatchPairs
(
);
// ************************************************************************* //

62
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/save/boundary Normal file
View file

@ -0,0 +1,62 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/polyMesh";
object boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
6
(
ibCylinder
{
type immersedBoundary;
nFaces 0;
startFace 3650;
internalFlow no;
isWall yes;
}
in
{
type patch;
nFaces 25;
startFace 3650;
}
out
{
type patch;
nFaces 25;
startFace 3675;
}
top
{
type patch;
nFaces 75;
startFace 3700;
}
bottom
{
type patch;
nFaces 75;
startFace 3775;
}
frontAndBack
{
type empty;
nFaces 3750;
startFace 3850;
}
)
// ************************************************************************* //

56
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/save/k Normal file
View file

@ -0,0 +1,56 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: http://www.foam-extend.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.375;
boundaryField
{
ibCylinder
{
type immersedBoundaryKqRWallFunction;
patchType immersedBoundary;
value uniform 0.375;
}
top
{
type zeroGradient;
}
in
{
type fixedValue;
value uniform 0.375;
}
out
{
type inletOutlet;
inletValue uniform 0.375;
value uniform 0.375;
}
bottom
{
type zeroGradient;
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //

57
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/save/nut Normal file
View file

@ -0,0 +1,57 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: http://www.foam-extend.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
{
ibCylinder
{
type immersedBoundaryNutWallFunction;
patchType immersedBoundary;
value uniform 0;
}
top
{
type zeroGradient;
value uniform 0;
}
in
{
type zeroGradient;
value uniform 0;
}
out
{
type zeroGradient;
value uniform 0;
}
bottom
{
type zeroGradient;
value uniform 0;
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //

60
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/save/omega Normal file
View file

@ -0,0 +1,60 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 -1 0 0 0 0];
internalField uniform 500;
boundaryField
{
ibCylinder
{
type immersedBoundaryOmegaWallFunction;
patchType immersedBoundary;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 500;
}
top
{
type zeroGradient;
}
in
{
type fixedValue;
value uniform 500;
}
out
{
type inletOutlet;
inletValue uniform 500;
value uniform 500;
}
bottom
{
type zeroGradient;
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //

1938
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/save/p Normal file
View file

File diff suppressed because it is too large Load diff

138
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/system/controlDict Normal file
View file

@ -0,0 +1,138 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
DebugSwitches
{
// immersedBoundary 1;
// polyMesh 2;
// fvMesh 2;
}
application pimpleDyMFoam;
startFrom startTime;
startTime 0;
stopAt endTime;
endTime 5;
deltaT 0.02;
// writeControl runTime;
// writeInterval 0.2;
writeControl timeStep;
writeInterval 1;
purgeWrite 0;
writeFormat ascii;
writePrecision 6;
writeCompression uncompressed;
timeFormat general;
timePrecision 6;
runTimeModifiable yes;
adjustTimeStep no;
maxCo 0.2;
libs
(
"libimmersedBoundary.so"
"libimmersedBoundaryDynamicMesh.so"
"libimmersedBoundaryTurbulence.so"
);
functions
(
forces
{
type forces;
functionObjectLibs ("libforces.so");
outputControl timeStep;
outputInterval 1;
patches ( ibCylinder );
pName p;
UName U;
rhoName rhoInf;
rhoInf 1;
log true;
CofR ( 0 0 0 );
Aref 0.05;
Uref 1;
}
probes
{
type probes;
functionObjectLibs ("libsampling.so");
outputControl timeStep;
outputInterval 1;
probeLocations
(
( 1.5 0 0)
( 1.75 0 0)
);
fields
(
p
);
}
minMaxP
{
type minMaxField;
functionObjectLibs ("libfieldFunctionObjects.so");
name p;
}
minMaxU
{
type minMaxField;
functionObjectLibs ("libfieldFunctionObjects.so");
name U;
}
minMaxNut
{
type minMaxField;
functionObjectLibs ("libfieldFunctionObjects.so");
name nut;
}
minMaxK
{
type minMaxField;
functionObjectLibs ("libfieldFunctionObjects.so");
name k;
}
minMaxNut
{
type minMaxField;
functionObjectLibs ("libfieldFunctionObjects.so");
name epsilon;
}
);
// ************************************************************************* //

34
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/system/decomposeParDict Normal file
View file

@ -0,0 +1,34 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object decomposeParDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
numberOfSubdomains 3;
method scotch;
simpleCoeffs
{
n (5 2 1);
delta 0.001;
}
distributed no;
roots
(
);
// ************************************************************************* //

54
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/system/fvSchemes Normal file
View file

@ -0,0 +1,54 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
default Euler;
}
gradSchemes
{
default cellLimited Gauss linear 1;
}
divSchemes
{
default none;
div(phi,U) Gauss upwind;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(phi,omega) Gauss upwind;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}
laplacianSchemes
{
default Gauss linear limited 0.5;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default limited 0.5;
}
// ************************************************************************* //

78
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/system/fvSolution Normal file
View file

@ -0,0 +1,78 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
"p|pFinal|pcorr"
{
solver CG;
preconditioner Cholesky;
minIter 1;
maxIter 1000;
tolerance 1e-7;
relTol 0.0;
}
"U|UFinal|k|epsilon|omega"
{
solver BiCGStab;
preconditioner ILU0;
minIter 1;
maxIter 1000;
tolerance 1e-08;
relTol 0;
}
}
SIMPLE
{
nNonOrthogonalCorrectors 1;
pRefPoint (0 -0.45 0.05);
pRefValue 0;
}
PIMPLE
{
nOuterCorrectors 4;
nCorrectors 2;
nNonOrthogonalCorrectors 0;
pRefPoint (0 -0.45 0.05);
pRefValue 0;
checkMeshCourantNo no;
correctPhi yes;
}
PISO
{
nNonOrthogonalCorrectors 1;
}
relaxationFactors
{
equations
{
U 1;
UFinal 1;
}
fields
{
}
}
// ************************************************************************* //

31
tutorials/immersedBoundary/refiningMovingCylinderInChannelTurbulent/system/mapFieldsDict Normal file
View file

@ -0,0 +1,31 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | foam-extend: Open Source CFD |
| \\ / O peration | Version: 4.0 |
| \\ / A nd | Web: http://www.foam-extend.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object mapFieldsDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
// List of pairs of source/target patches for mapping
patchMap
(
// lid movingWall
);
// List of target patches cutting the source domain (these need to be
// handled specially e.g. interpolated from internal values)
cuttingPatches
(
// fixedWalls
);
// ************************************************************************* //