From 74eefa88c6399bb1527ce1470c7dac1377233d08 Mon Sep 17 00:00:00 2001 From: Hrvoje Jasak Date: Tue, 15 Oct 2019 09:00:39 +0100 Subject: [PATCH] Immersed boundary tutorial: two interfering IB patches --- tutorials/immersedBoundary/twoIbPatches/.foam | 0 .../immersedBoundary/twoIbPatches/.gitignore | 52 + .../immersedBoundary/twoIbPatches/0_org/U | 67 + .../twoIbPatches/0_org/epsilon | 74 + .../immersedBoundary/twoIbPatches/0_org/k | 66 + .../immersedBoundary/twoIbPatches/0_org/nut | 67 + .../immersedBoundary/twoIbPatches/0_org/omega | 70 + .../immersedBoundary/twoIbPatches/0_org/p | 69 + .../immersedBoundary/twoIbPatches/Allclean | 10 + .../immersedBoundary/twoIbPatches/Allrun | 16 + .../twoIbPatches/constant/RASProperties | 25 + .../twoIbPatches/constant/dynamicMeshDict | 56 + .../twoIbPatches/constant/transportProperties | 21 + .../constant/triSurface/ibCylinder.ftr | 205 ++ .../constant/triSurface/ibCylinder.stl | 674 ++++++ .../constant/triSurface/ibCylinder2.ftr | 205 ++ .../constant/triSurface/ibCylinder2.stl | 674 ++++++ .../constant/turbulenceProperties | 21 + .../twoIbPatches/save/blockMeshDict | 70 + .../twoIbPatches/save/boundary | 71 + .../immersedBoundary/twoIbPatches/save/k | 56 + .../immersedBoundary/twoIbPatches/save/nut | 57 + .../immersedBoundary/twoIbPatches/save/omega | 60 + .../immersedBoundary/twoIbPatches/save/p | 1938 +++++++++++++++++ .../twoIbPatches/system/controlDict | 152 ++ .../twoIbPatches/system/decomposeParDict | 34 + .../twoIbPatches/system/fvSchemes | 54 + .../twoIbPatches/system/fvSolution | 78 + .../twoIbPatches/system/mapFieldsDict | 31 + .../twoIbPatches/twoIbPatches.foam | 0 30 files changed, 4973 insertions(+) create mode 100644 tutorials/immersedBoundary/twoIbPatches/.foam create mode 100644 tutorials/immersedBoundary/twoIbPatches/.gitignore create mode 100644 tutorials/immersedBoundary/twoIbPatches/0_org/U create mode 100644 tutorials/immersedBoundary/twoIbPatches/0_org/epsilon create mode 100644 tutorials/immersedBoundary/twoIbPatches/0_org/k create mode 100644 tutorials/immersedBoundary/twoIbPatches/0_org/nut create mode 100644 tutorials/immersedBoundary/twoIbPatches/0_org/omega create mode 100644 tutorials/immersedBoundary/twoIbPatches/0_org/p create mode 100755 tutorials/immersedBoundary/twoIbPatches/Allclean create mode 100755 tutorials/immersedBoundary/twoIbPatches/Allrun create mode 100644 tutorials/immersedBoundary/twoIbPatches/constant/RASProperties create mode 100644 tutorials/immersedBoundary/twoIbPatches/constant/dynamicMeshDict create mode 100644 tutorials/immersedBoundary/twoIbPatches/constant/transportProperties create mode 100644 tutorials/immersedBoundary/twoIbPatches/constant/triSurface/ibCylinder.ftr create mode 100644 tutorials/immersedBoundary/twoIbPatches/constant/triSurface/ibCylinder.stl create mode 100644 tutorials/immersedBoundary/twoIbPatches/constant/triSurface/ibCylinder2.ftr create mode 100644 tutorials/immersedBoundary/twoIbPatches/constant/triSurface/ibCylinder2.stl create mode 100644 tutorials/immersedBoundary/twoIbPatches/constant/turbulenceProperties create mode 100644 tutorials/immersedBoundary/twoIbPatches/save/blockMeshDict create mode 100644 tutorials/immersedBoundary/twoIbPatches/save/boundary create mode 100644 tutorials/immersedBoundary/twoIbPatches/save/k create mode 100644 tutorials/immersedBoundary/twoIbPatches/save/nut create mode 100644 tutorials/immersedBoundary/twoIbPatches/save/omega create mode 100644 tutorials/immersedBoundary/twoIbPatches/save/p create mode 100644 tutorials/immersedBoundary/twoIbPatches/system/controlDict create mode 100644 tutorials/immersedBoundary/twoIbPatches/system/decomposeParDict create mode 100644 tutorials/immersedBoundary/twoIbPatches/system/fvSchemes create mode 100644 tutorials/immersedBoundary/twoIbPatches/system/fvSolution create mode 100644 tutorials/immersedBoundary/twoIbPatches/system/mapFieldsDict create mode 100644 tutorials/immersedBoundary/twoIbPatches/twoIbPatches.foam diff --git a/tutorials/immersedBoundary/twoIbPatches/.foam b/tutorials/immersedBoundary/twoIbPatches/.foam new file mode 100644 index 000000000..e69de29bb diff --git a/tutorials/immersedBoundary/twoIbPatches/.gitignore b/tutorials/immersedBoundary/twoIbPatches/.gitignore new file mode 100644 index 000000000..f42669c4e --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/.gitignore @@ -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 diff --git a/tutorials/immersedBoundary/twoIbPatches/0_org/U b/tutorials/immersedBoundary/twoIbPatches/0_org/U new file mode 100644 index 000000000..47c6575eb --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/0_org/U @@ -0,0 +1,67 @@ +/*--------------------------------*- 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); + } + ibCylinder2 + { + 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; + } +} + + +// ************************************************************************* // diff --git a/tutorials/immersedBoundary/twoIbPatches/0_org/epsilon b/tutorials/immersedBoundary/twoIbPatches/0_org/epsilon new file mode 100644 index 000000000..459a61353 --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/0_org/epsilon @@ -0,0 +1,74 @@ +/*--------------------------------*- 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; + } + ibCylinder2 + { + 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; + } +} + + +// ************************************************************************* // diff --git a/tutorials/immersedBoundary/twoIbPatches/0_org/k b/tutorials/immersedBoundary/twoIbPatches/0_org/k new file mode 100644 index 000000000..7f82f1b74 --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/0_org/k @@ -0,0 +1,66 @@ +/*--------------------------------*- 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 yes; + deadValue 1.5e-4; + value uniform 1.5e-4; + } + ibCylinder2 + { + type immersedBoundaryKqRWallFunction; + patchType immersedBoundary; + setDeadValue yes; + 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; + } +} + + +// ************************************************************************* // diff --git a/tutorials/immersedBoundary/twoIbPatches/0_org/nut b/tutorials/immersedBoundary/twoIbPatches/0_org/nut new file mode 100644 index 000000000..09b6680a9 --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/0_org/nut @@ -0,0 +1,67 @@ +/*--------------------------------*- 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; + } + ibCylinder2 + { + 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; + } +} + + +// ************************************************************************* // diff --git a/tutorials/immersedBoundary/twoIbPatches/0_org/omega b/tutorials/immersedBoundary/twoIbPatches/0_org/omega new file mode 100644 index 000000000..824189067 --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/0_org/omega @@ -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 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; + } + ibCylinder2 + { + 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; + } +} + + +// ************************************************************************* // diff --git a/tutorials/immersedBoundary/twoIbPatches/0_org/p b/tutorials/immersedBoundary/twoIbPatches/0_org/p new file mode 100644 index 000000000..0d8e58174 --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/0_org/p @@ -0,0 +1,69 @@ +/*--------------------------------*- 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 yes; + deadValue 0; + value uniform 0; + } + ibCylinder2 + { + type mixedIb; + patchType immersedBoundary; + triValue uniform 0; + triGradient uniform 0; + triValueFraction uniform 0; + setDeadValue yes; + deadValue 0; + value uniform 0; + } + in + { + type zeroGradient; + } + out + { + type fixedValue; + value uniform 0; + } + top + { + type zeroGradient; + } + bottom + { + type zeroGradient; + } + frontAndBack + { + type empty; + } +} + + +// ************************************************************************* // diff --git a/tutorials/immersedBoundary/twoIbPatches/Allclean b/tutorials/immersedBoundary/twoIbPatches/Allclean new file mode 100755 index 000000000..1ec2089d8 --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/Allclean @@ -0,0 +1,10 @@ +#!/bin/sh + +# Source tutorial clean functions +. $WM_PROJECT_DIR/bin/tools/CleanFunctions + +cleanCase +\rm -f *.dat +\rm -rf 0 +\rm -rf constant/polyMesh +\rm -rf postProcessing diff --git a/tutorials/immersedBoundary/twoIbPatches/Allrun b/tutorials/immersedBoundary/twoIbPatches/Allrun new file mode 100755 index 000000000..dea08785a --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/Allrun @@ -0,0 +1,16 @@ +#!/bin/sh +# Source tutorial run functions +. $WM_PROJECT_DIR/bin/tools/RunFunctions + +# Get application name +application="pimpleDyMIbFoam" + +runApplication blockMesh -dict save/blockMeshDict +cp save/boundary constant/polyMesh/ + +mkdir 0 +\cp -f 0_org/* 0/ + +runApplication potentialFoam +runApplication writeIbMasks +runApplication $application diff --git a/tutorials/immersedBoundary/twoIbPatches/constant/RASProperties b/tutorials/immersedBoundary/twoIbPatches/constant/RASProperties new file mode 100644 index 000000000..b0012fd4b --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/constant/RASProperties @@ -0,0 +1,25 @@ +/*--------------------------------*- 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; + + +// ************************************************************************* // diff --git a/tutorials/immersedBoundary/twoIbPatches/constant/dynamicMeshDict b/tutorials/immersedBoundary/twoIbPatches/constant/dynamicMeshDict new file mode 100644 index 000000000..e6f7b10e0 --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/constant/dynamicMeshDict @@ -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 dictionary; + object dynamicMeshDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +// dynamicFvMesh staticFvMesh; +dynamicFvMesh immersedBoundarySolidBodyMotionFvMesh; + +immersedBoundarySolidBodyMotionFvMeshCoeffs +{ + motionFunctions + ( + + ibCylinder + { +// solidBodyMotionFunction translation; +// translationCoeffs +// { +// velocity (0.1 0 0); +// } + + solidBodyMotionFunction linearOscillation; + linearOscillationCoeffs + { + amplitude (0.5 0 0); + period 2.5; + } + } + + ibCylinder2 + { + solidBodyMotionFunction linearOscillation; + linearOscillationCoeffs + { + amplitude (0 0.05 0); + period 2.5; + } + } + + + ); +} + + +// ************************************************************************* // diff --git a/tutorials/immersedBoundary/twoIbPatches/constant/transportProperties b/tutorials/immersedBoundary/twoIbPatches/constant/transportProperties new file mode 100644 index 000000000..d29c3feb3 --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/constant/transportProperties @@ -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; + + +// ************************************************************************* // diff --git a/tutorials/immersedBoundary/twoIbPatches/constant/triSurface/ibCylinder.ftr b/tutorials/immersedBoundary/twoIbPatches/constant/triSurface/ibCylinder.ftr new file mode 100644 index 000000000..44af86500 --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/constant/triSurface/ibCylinder.ftr @@ -0,0 +1,205 @@ +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) +) + diff --git a/tutorials/immersedBoundary/twoIbPatches/constant/triSurface/ibCylinder.stl b/tutorials/immersedBoundary/twoIbPatches/constant/triSurface/ibCylinder.stl new file mode 100644 index 000000000..efdf1e681 --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/constant/triSurface/ibCylinder.stl @@ -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.11896e-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.11896e-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.11896e-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.11896e-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.11896e-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.11896e-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.11896e-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.11896e-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 diff --git a/tutorials/immersedBoundary/twoIbPatches/constant/triSurface/ibCylinder2.ftr b/tutorials/immersedBoundary/twoIbPatches/constant/triSurface/ibCylinder2.ftr new file mode 100644 index 000000000..5c0d9a881 --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/constant/triSurface/ibCylinder2.ftr @@ -0,0 +1,205 @@ +1( +cylinder +patch) + +96 +( +(0.45 8.32665e-17 0.05) +(0.45 -8.32665e-17 -0.05) +(0.452139 0.0326315 -0.05) +(0.452139 0.0326315 0.05) +(0.452139 -0.0326315 -0.05) +(0.452139 -0.0326315 0.05) +(0.458519 0.064705 -0.05) +(0.458519 0.064705 0.05) +(0.458519 -0.064705 -0.05) +(0.458519 -0.064705 0.05) +(0.46903 0.095671 -0.05) +(0.46903 0.095671 0.05) +(0.46903 -0.095671 -0.05) +(0.46903 -0.095671 0.05) +(0.483494 0.125 -0.05) +(0.483494 0.125 0.05) +(0.483494 -0.125 -0.05) +(0.483494 -0.125 0.05) +(0.501662 0.152191 -0.05) +(0.501662 0.152191 0.05) +(0.501662 -0.152191 -0.05) +(0.501662 -0.152191 0.05) +(0.523224 0.176776 -0.05) +(0.523224 0.176776 0.05) +(0.523224 -0.176776 -0.05) +(0.523224 -0.176776 0.05) +(0.547809 0.198338 -0.05) +(0.547809 0.198338 0.05) +(0.547809 -0.198338 0.05) +(0.547809 -0.198338 -0.05) +(0.575 0.216506 -0.05) +(0.575 0.216506 0.05) +(0.575 -0.216506 0.05) +(0.575 -0.216506 -0.05) +(0.604329 0.23097 -0.05) +(0.604329 0.23097 0.05) +(0.604329 -0.23097 0.05) +(0.604329 -0.23097 -0.05) +(0.635295 0.241481 -0.05) +(0.635295 0.241481 0.05) +(0.635295 -0.241481 0.05) +(0.635295 -0.241481 -0.05) +(0.667368 0.247861 -0.05) +(0.667368 0.247861 0.05) +(0.667368 -0.247861 0.05) +(0.667368 -0.247861 -0.05) +(0.7 0.25 0.05) +(0.7 0.25 -0.05) +(0.7 -0.25 -0.05) +(0.7 -0.25 0.05) +(0.732632 0.247861 -0.05) +(0.732632 0.247861 0.05) +(0.732632 -0.247861 0.05) +(0.732632 -0.247861 -0.05) +(0.764705 0.241481 -0.05) +(0.764705 0.241481 0.05) +(0.764705 -0.241481 0.05) +(0.764705 -0.241481 -0.05) +(0.795671 0.23097 -0.05) +(0.795671 0.23097 0.05) +(0.795671 -0.23097 0.05) +(0.795671 -0.23097 -0.05) +(0.825 0.216506 -0.05) +(0.825 0.216506 0.05) +(0.825 -0.216506 0.05) +(0.825 -0.216506 -0.05) +(0.852191 0.198338 -0.05) +(0.852191 0.198338 0.05) +(0.852191 -0.198338 0.05) +(0.852191 -0.198338 -0.05) +(0.876776 -0.176776 -0.05) +(0.876776 -0.176776 0.05) +(0.876776 0.176776 -0.05) +(0.876776 0.176776 0.05) +(0.898338 -0.152191 -0.05) +(0.898338 -0.152191 0.05) +(0.898338 0.152191 -0.05) +(0.898338 0.152191 0.05) +(0.916506 -0.125 -0.05) +(0.916506 -0.125 0.05) +(0.916506 0.125 -0.05) +(0.916506 0.125 0.05) +(0.93097 -0.095671 -0.05) +(0.93097 -0.095671 0.05) +(0.93097 0.095671 -0.05) +(0.93097 0.095671 0.05) +(0.941481 -0.064705 -0.05) +(0.941481 -0.064705 0.05) +(0.941481 0.064705 -0.05) +(0.941481 0.064705 0.05) +(0.947861 -0.0326315 -0.05) +(0.947861 -0.0326315 0.05) +(0.947861 0.0326315 -0.05) +(0.947861 0.0326315 0.05) +(0.95 -8.32665e-17 0.05) +(0.95 8.32665e-17 -0.05) +) + + +96 +( +((70 74 75) 0) +((75 71 70) 0) +((74 78 79) 0) +((79 75 74) 0) +((78 82 83) 0) +((83 79 78) 0) +((82 86 87) 0) +((87 83 82) 0) +((86 90 91) 0) +((91 87 86) 0) +((94 91 90) 0) +((90 95 94) 0) +((95 92 93) 0) +((93 94 95) 0) +((92 88 89) 0) +((89 93 92) 0) +((88 84 85) 0) +((85 89 88) 0) +((84 80 81) 0) +((81 85 84) 0) +((80 76 77) 0) +((77 81 80) 0) +((76 72 73) 0) +((73 77 76) 0) +((72 66 67) 0) +((67 73 72) 0) +((66 62 63) 0) +((63 67 66) 0) +((62 58 59) 0) +((59 63 62) 0) +((58 54 55) 0) +((55 59 58) 0) +((54 50 51) 0) +((51 55 54) 0) +((46 51 50) 0) +((50 47 46) 0) +((47 42 43) 0) +((43 46 47) 0) +((42 38 39) 0) +((39 43 42) 0) +((38 34 35) 0) +((35 39 38) 0) +((34 30 31) 0) +((31 35 34) 0) +((30 26 27) 0) +((27 31 30) 0) +((26 22 23) 0) +((23 27 26) 0) +((22 18 19) 0) +((19 23 22) 0) +((18 14 15) 0) +((15 19 18) 0) +((14 10 11) 0) +((11 15 14) 0) +((10 6 7) 0) +((7 11 10) 0) +((6 2 3) 0) +((3 7 6) 0) +((0 3 2) 0) +((2 1 0) 0) +((1 4 5) 0) +((5 0 1) 0) +((4 8 9) 0) +((9 5 4) 0) +((8 12 13) 0) +((13 9 8) 0) +((12 16 17) 0) +((17 13 12) 0) +((16 20 21) 0) +((21 17 16) 0) +((20 24 25) 0) +((25 21 20) 0) +((70 71 68) 0) +((68 69 70) 0) +((69 68 64) 0) +((64 65 69) 0) +((65 64 60) 0) +((60 61 65) 0) +((61 60 56) 0) +((56 57 61) 0) +((57 56 52) 0) +((52 53 57) 0) +((48 53 52) 0) +((52 49 48) 0) +((49 44 45) 0) +((45 48 49) 0) +((45 44 40) 0) +((40 41 45) 0) +((41 40 36) 0) +((36 37 41) 0) +((37 36 32) 0) +((32 33 37) 0) +((33 32 28) 0) +((28 29 33) 0) +((29 28 25) 0) +((25 24 29) 0) +) + diff --git a/tutorials/immersedBoundary/twoIbPatches/constant/triSurface/ibCylinder2.stl b/tutorials/immersedBoundary/twoIbPatches/constant/triSurface/ibCylinder2.stl new file mode 100644 index 000000000..53603fd5f --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/constant/triSurface/ibCylinder2.stl @@ -0,0 +1,674 @@ +solid cylinder + facet normal 0.751817 -0.659372 0 + outer loop + vertex 0.876776 -0.176776 -0.05 + vertex 0.898338 -0.152191 -0.05 + vertex 0.898338 -0.152191 0.05 + endloop + endfacet + facet normal 0.751817 -0.659372 0 + outer loop + vertex 0.898338 -0.152191 0.05 + vertex 0.876776 -0.176776 0.05 + vertex 0.876776 -0.176776 -0.05 + endloop + endfacet + facet normal 0.831476 -0.555561 0 + outer loop + vertex 0.898338 -0.152191 -0.05 + vertex 0.916506 -0.125 -0.05 + vertex 0.916506 -0.125 0.05 + endloop + endfacet + facet normal 0.831476 -0.555561 0 + outer loop + vertex 0.916506 -0.125 0.05 + vertex 0.898338 -0.152191 0.05 + vertex 0.898338 -0.152191 -0.05 + endloop + endfacet + facet normal 0.896866 -0.442302 0 + outer loop + vertex 0.916506 -0.125 -0.05 + vertex 0.93097 -0.095671 -0.05 + vertex 0.93097 -0.095671 0.05 + endloop + endfacet + facet normal 0.896866 -0.442302 0 + outer loop + vertex 0.93097 -0.095671 0.05 + vertex 0.916506 -0.125 0.05 + vertex 0.916506 -0.125 -0.05 + endloop + endfacet + facet normal 0.946935 -0.321425 0 + outer loop + vertex 0.93097 -0.095671 -0.05 + vertex 0.941481 -0.064705 -0.05 + vertex 0.941481 -0.064705 0.05 + endloop + endfacet + facet normal 0.946935 -0.321425 0 + outer loop + vertex 0.941481 -0.064705 0.05 + vertex 0.93097 -0.095671 0.05 + vertex 0.93097 -0.095671 -0.05 + endloop + endfacet + facet normal 0.980784 -0.195096 0 + outer loop + vertex 0.941481 -0.064705 -0.05 + vertex 0.947861 -0.0326315 -0.05 + vertex 0.947861 -0.0326315 0.05 + endloop + endfacet + facet normal 0.980784 -0.195096 0 + outer loop + vertex 0.947861 -0.0326315 0.05 + vertex 0.941481 -0.064705 0.05 + vertex 0.941481 -0.064705 -0.05 + endloop + endfacet + facet normal 0.997858 -0.0654098 0 + outer loop + vertex 0.95 -8.32665e-17 0.05 + vertex 0.947861 -0.0326315 0.05 + vertex 0.947861 -0.0326315 -0.05 + endloop + endfacet + facet normal 0.997858 -0.0654098 -1.07752e-16 + outer loop + vertex 0.947861 -0.0326315 -0.05 + vertex 0.95 8.32665e-17 -0.05 + vertex 0.95 -8.32665e-17 0.05 + endloop + endfacet + facet normal 0.997858 0.0654098 0 + outer loop + vertex 0.95 8.32665e-17 -0.05 + vertex 0.947861 0.0326315 -0.05 + vertex 0.947861 0.0326315 0.05 + endloop + endfacet + facet normal 0.997858 0.0654098 1.07752e-16 + outer loop + vertex 0.947861 0.0326315 0.05 + vertex 0.95 -8.32665e-17 0.05 + vertex 0.95 8.32665e-17 -0.05 + endloop + endfacet + facet normal 0.980784 0.195096 0 + outer loop + vertex 0.947861 0.0326315 -0.05 + vertex 0.941481 0.064705 -0.05 + vertex 0.941481 0.064705 0.05 + endloop + endfacet + facet normal 0.980784 0.195096 0 + outer loop + vertex 0.941481 0.064705 0.05 + vertex 0.947861 0.0326315 0.05 + vertex 0.947861 0.0326315 -0.05 + endloop + endfacet + facet normal 0.946935 0.321425 0 + outer loop + vertex 0.941481 0.064705 -0.05 + vertex 0.93097 0.095671 -0.05 + vertex 0.93097 0.095671 0.05 + endloop + endfacet + facet normal 0.946935 0.321425 0 + outer loop + vertex 0.93097 0.095671 0.05 + vertex 0.941481 0.064705 0.05 + vertex 0.941481 0.064705 -0.05 + endloop + endfacet + facet normal 0.896866 0.442302 0 + outer loop + vertex 0.93097 0.095671 -0.05 + vertex 0.916506 0.125 -0.05 + vertex 0.916506 0.125 0.05 + endloop + endfacet + facet normal 0.896866 0.442302 0 + outer loop + vertex 0.916506 0.125 0.05 + vertex 0.93097 0.095671 0.05 + vertex 0.93097 0.095671 -0.05 + endloop + endfacet + facet normal 0.831476 0.555561 0 + outer loop + vertex 0.916506 0.125 -0.05 + vertex 0.898338 0.152191 -0.05 + vertex 0.898338 0.152191 0.05 + endloop + endfacet + facet normal 0.831476 0.555561 0 + outer loop + vertex 0.898338 0.152191 0.05 + vertex 0.916506 0.125 0.05 + vertex 0.916506 0.125 -0.05 + endloop + endfacet + facet normal 0.751817 0.659372 0 + outer loop + vertex 0.898338 0.152191 -0.05 + vertex 0.876776 0.176776 -0.05 + vertex 0.876776 0.176776 0.05 + endloop + endfacet + facet normal 0.751817 0.659372 0 + outer loop + vertex 0.876776 0.176776 0.05 + vertex 0.898338 0.152191 0.05 + vertex 0.898338 0.152191 -0.05 + endloop + endfacet + facet normal 0.659372 0.751817 0 + outer loop + vertex 0.876776 0.176776 -0.05 + vertex 0.852191 0.198338 -0.05 + vertex 0.852191 0.198338 0.05 + endloop + endfacet + facet normal 0.659372 0.751817 0 + outer loop + vertex 0.852191 0.198338 0.05 + vertex 0.876776 0.176776 0.05 + vertex 0.876776 0.176776 -0.05 + endloop + endfacet + facet normal 0.555561 0.831476 0 + outer loop + vertex 0.852191 0.198338 -0.05 + vertex 0.825 0.216506 -0.05 + vertex 0.825 0.216506 0.05 + endloop + endfacet + facet normal 0.555561 0.831476 0 + outer loop + vertex 0.825 0.216506 0.05 + vertex 0.852191 0.198338 0.05 + vertex 0.852191 0.198338 -0.05 + endloop + endfacet + facet normal 0.442302 0.896866 0 + outer loop + vertex 0.825 0.216506 -0.05 + vertex 0.795671 0.23097 -0.05 + vertex 0.795671 0.23097 0.05 + endloop + endfacet + facet normal 0.442302 0.896866 0 + outer loop + vertex 0.795671 0.23097 0.05 + vertex 0.825 0.216506 0.05 + vertex 0.825 0.216506 -0.05 + endloop + endfacet + facet normal 0.321425 0.946935 0 + outer loop + vertex 0.795671 0.23097 -0.05 + vertex 0.764705 0.241481 -0.05 + vertex 0.764705 0.241481 0.05 + endloop + endfacet + facet normal 0.321425 0.946935 0 + outer loop + vertex 0.764705 0.241481 0.05 + vertex 0.795671 0.23097 0.05 + vertex 0.795671 0.23097 -0.05 + endloop + endfacet + facet normal 0.195096 0.980784 0 + outer loop + vertex 0.764705 0.241481 -0.05 + vertex 0.732632 0.247861 -0.05 + vertex 0.732632 0.247861 0.05 + endloop + endfacet + facet normal 0.195096 0.980784 0 + outer loop + vertex 0.732632 0.247861 0.05 + vertex 0.764705 0.241481 0.05 + vertex 0.764705 0.241481 -0.05 + endloop + endfacet + facet normal 0.0654098 0.997858 0 + outer loop + vertex 0.7 0.25 0.05 + vertex 0.732632 0.247861 0.05 + vertex 0.732632 0.247861 -0.05 + endloop + endfacet + facet normal 0.0654098 0.997858 -1.49195e-16 + outer loop + vertex 0.732632 0.247861 -0.05 + vertex 0.7 0.25 -0.05 + vertex 0.7 0.25 0.05 + endloop + endfacet + facet normal -0.0654098 0.997858 0 + outer loop + vertex 0.7 0.25 -0.05 + vertex 0.667368 0.247861 -0.05 + vertex 0.667368 0.247861 0.05 + endloop + endfacet + facet normal -0.0654098 0.997858 1.49195e-16 + outer loop + vertex 0.667368 0.247861 0.05 + vertex 0.7 0.25 0.05 + vertex 0.7 0.25 -0.05 + endloop + endfacet + facet normal -0.195096 0.980784 0 + outer loop + vertex 0.667368 0.247861 -0.05 + vertex 0.635295 0.241481 -0.05 + vertex 0.635295 0.241481 0.05 + endloop + endfacet + facet normal -0.195096 0.980784 0 + outer loop + vertex 0.635295 0.241481 0.05 + vertex 0.667368 0.247861 0.05 + vertex 0.667368 0.247861 -0.05 + endloop + endfacet + facet normal -0.321425 0.946935 0 + outer loop + vertex 0.635295 0.241481 -0.05 + vertex 0.604329 0.23097 -0.05 + vertex 0.604329 0.23097 0.05 + endloop + endfacet + facet normal -0.321425 0.946935 0 + outer loop + vertex 0.604329 0.23097 0.05 + vertex 0.635295 0.241481 0.05 + vertex 0.635295 0.241481 -0.05 + endloop + endfacet + facet normal -0.442302 0.896866 0 + outer loop + vertex 0.604329 0.23097 -0.05 + vertex 0.575 0.216506 -0.05 + vertex 0.575 0.216506 0.05 + endloop + endfacet + facet normal -0.442302 0.896866 0 + outer loop + vertex 0.575 0.216506 0.05 + vertex 0.604329 0.23097 0.05 + vertex 0.604329 0.23097 -0.05 + endloop + endfacet + facet normal -0.555561 0.831476 0 + outer loop + vertex 0.575 0.216506 -0.05 + vertex 0.547809 0.198338 -0.05 + vertex 0.547809 0.198338 0.05 + endloop + endfacet + facet normal -0.555561 0.831476 0 + outer loop + vertex 0.547809 0.198338 0.05 + vertex 0.575 0.216506 0.05 + vertex 0.575 0.216506 -0.05 + endloop + endfacet + facet normal -0.659372 0.751817 0 + outer loop + vertex 0.547809 0.198338 -0.05 + vertex 0.523224 0.176776 -0.05 + vertex 0.523224 0.176776 0.05 + endloop + endfacet + facet normal -0.659372 0.751817 0 + outer loop + vertex 0.523224 0.176776 0.05 + vertex 0.547809 0.198338 0.05 + vertex 0.547809 0.198338 -0.05 + endloop + endfacet + facet normal -0.751817 0.659372 0 + outer loop + vertex 0.523224 0.176776 -0.05 + vertex 0.501662 0.152191 -0.05 + vertex 0.501662 0.152191 0.05 + endloop + endfacet + facet normal -0.751817 0.659372 0 + outer loop + vertex 0.501662 0.152191 0.05 + vertex 0.523224 0.176776 0.05 + vertex 0.523224 0.176776 -0.05 + endloop + endfacet + facet normal -0.831476 0.555561 0 + outer loop + vertex 0.501662 0.152191 -0.05 + vertex 0.483494 0.125 -0.05 + vertex 0.483494 0.125 0.05 + endloop + endfacet + facet normal -0.831476 0.555561 0 + outer loop + vertex 0.483494 0.125 0.05 + vertex 0.501662 0.152191 0.05 + vertex 0.501662 0.152191 -0.05 + endloop + endfacet + facet normal -0.896866 0.442302 0 + outer loop + vertex 0.483494 0.125 -0.05 + vertex 0.46903 0.095671 -0.05 + vertex 0.46903 0.095671 0.05 + endloop + endfacet + facet normal -0.896866 0.442302 0 + outer loop + vertex 0.46903 0.095671 0.05 + vertex 0.483494 0.125 0.05 + vertex 0.483494 0.125 -0.05 + endloop + endfacet + facet normal -0.946935 0.321425 0 + outer loop + vertex 0.46903 0.095671 -0.05 + vertex 0.458519 0.064705 -0.05 + vertex 0.458519 0.064705 0.05 + endloop + endfacet + facet normal -0.946935 0.321425 0 + outer loop + vertex 0.458519 0.064705 0.05 + vertex 0.46903 0.095671 0.05 + vertex 0.46903 0.095671 -0.05 + endloop + endfacet + facet normal -0.980784 0.195096 0 + outer loop + vertex 0.458519 0.064705 -0.05 + vertex 0.452139 0.0326315 -0.05 + vertex 0.452139 0.0326315 0.05 + endloop + endfacet + facet normal -0.980784 0.195096 0 + outer loop + vertex 0.452139 0.0326315 0.05 + vertex 0.458519 0.064705 0.05 + vertex 0.458519 0.064705 -0.05 + endloop + endfacet + facet normal -0.997858 0.0654098 0 + outer loop + vertex 0.45 8.32665e-17 0.05 + vertex 0.452139 0.0326315 0.05 + vertex 0.452139 0.0326315 -0.05 + endloop + endfacet + facet normal -0.997858 0.0654098 -1.07752e-16 + outer loop + vertex 0.452139 0.0326315 -0.05 + vertex 0.45 -8.32665e-17 -0.05 + vertex 0.45 8.32665e-17 0.05 + endloop + endfacet + facet normal -0.997858 -0.0654098 0 + outer loop + vertex 0.45 -8.32665e-17 -0.05 + vertex 0.452139 -0.0326315 -0.05 + vertex 0.452139 -0.0326315 0.05 + endloop + endfacet + facet normal -0.997858 -0.0654098 1.07752e-16 + outer loop + vertex 0.452139 -0.0326315 0.05 + vertex 0.45 8.32665e-17 0.05 + vertex 0.45 -8.32665e-17 -0.05 + endloop + endfacet + facet normal -0.980784 -0.195096 0 + outer loop + vertex 0.452139 -0.0326315 -0.05 + vertex 0.458519 -0.064705 -0.05 + vertex 0.458519 -0.064705 0.05 + endloop + endfacet + facet normal -0.980784 -0.195096 0 + outer loop + vertex 0.458519 -0.064705 0.05 + vertex 0.452139 -0.0326315 0.05 + vertex 0.452139 -0.0326315 -0.05 + endloop + endfacet + facet normal -0.946935 -0.321425 0 + outer loop + vertex 0.458519 -0.064705 -0.05 + vertex 0.46903 -0.095671 -0.05 + vertex 0.46903 -0.095671 0.05 + endloop + endfacet + facet normal -0.946935 -0.321425 0 + outer loop + vertex 0.46903 -0.095671 0.05 + vertex 0.458519 -0.064705 0.05 + vertex 0.458519 -0.064705 -0.05 + endloop + endfacet + facet normal -0.896866 -0.442302 0 + outer loop + vertex 0.46903 -0.095671 -0.05 + vertex 0.483494 -0.125 -0.05 + vertex 0.483494 -0.125 0.05 + endloop + endfacet + facet normal -0.896866 -0.442302 0 + outer loop + vertex 0.483494 -0.125 0.05 + vertex 0.46903 -0.095671 0.05 + vertex 0.46903 -0.095671 -0.05 + endloop + endfacet + facet normal -0.831476 -0.555561 0 + outer loop + vertex 0.483494 -0.125 -0.05 + vertex 0.501662 -0.152191 -0.05 + vertex 0.501662 -0.152191 0.05 + endloop + endfacet + facet normal -0.831476 -0.555561 0 + outer loop + vertex 0.501662 -0.152191 0.05 + vertex 0.483494 -0.125 0.05 + vertex 0.483494 -0.125 -0.05 + endloop + endfacet + facet normal -0.751817 -0.659372 0 + outer loop + vertex 0.501662 -0.152191 -0.05 + vertex 0.523224 -0.176776 -0.05 + vertex 0.523224 -0.176776 0.05 + endloop + endfacet + facet normal -0.751817 -0.659372 0 + outer loop + vertex 0.523224 -0.176776 0.05 + vertex 0.501662 -0.152191 0.05 + vertex 0.501662 -0.152191 -0.05 + endloop + endfacet + facet normal 0.659372 -0.751817 0 + outer loop + vertex 0.876776 -0.176776 -0.05 + vertex 0.876776 -0.176776 0.05 + vertex 0.852191 -0.198338 0.05 + endloop + endfacet + facet normal 0.659372 -0.751817 0 + outer loop + vertex 0.852191 -0.198338 0.05 + vertex 0.852191 -0.198338 -0.05 + vertex 0.876776 -0.176776 -0.05 + endloop + endfacet + facet normal 0.555561 -0.831476 0 + outer loop + vertex 0.852191 -0.198338 -0.05 + vertex 0.852191 -0.198338 0.05 + vertex 0.825 -0.216506 0.05 + endloop + endfacet + facet normal 0.555561 -0.831476 0 + outer loop + vertex 0.825 -0.216506 0.05 + vertex 0.825 -0.216506 -0.05 + vertex 0.852191 -0.198338 -0.05 + endloop + endfacet + facet normal 0.442302 -0.896866 0 + outer loop + vertex 0.825 -0.216506 -0.05 + vertex 0.825 -0.216506 0.05 + vertex 0.795671 -0.23097 0.05 + endloop + endfacet + facet normal 0.442302 -0.896866 0 + outer loop + vertex 0.795671 -0.23097 0.05 + vertex 0.795671 -0.23097 -0.05 + vertex 0.825 -0.216506 -0.05 + endloop + endfacet + facet normal 0.321425 -0.946935 0 + outer loop + vertex 0.795671 -0.23097 -0.05 + vertex 0.795671 -0.23097 0.05 + vertex 0.764705 -0.241481 0.05 + endloop + endfacet + facet normal 0.321425 -0.946935 0 + outer loop + vertex 0.764705 -0.241481 0.05 + vertex 0.764705 -0.241481 -0.05 + vertex 0.795671 -0.23097 -0.05 + endloop + endfacet + facet normal 0.195096 -0.980784 0 + outer loop + vertex 0.764705 -0.241481 -0.05 + vertex 0.764705 -0.241481 0.05 + vertex 0.732632 -0.247861 0.05 + endloop + endfacet + facet normal 0.195096 -0.980784 0 + outer loop + vertex 0.732632 -0.247861 0.05 + vertex 0.732632 -0.247861 -0.05 + vertex 0.764705 -0.241481 -0.05 + endloop + endfacet + facet normal 0.0654098 -0.997858 0 + outer loop + vertex 0.7 -0.25 -0.05 + vertex 0.732632 -0.247861 -0.05 + vertex 0.732632 -0.247861 0.05 + endloop + endfacet + facet normal 0.0654098 -0.997858 1.49195e-16 + outer loop + vertex 0.732632 -0.247861 0.05 + vertex 0.7 -0.25 0.05 + vertex 0.7 -0.25 -0.05 + endloop + endfacet + facet normal -0.0654098 -0.997858 0 + outer loop + vertex 0.7 -0.25 0.05 + vertex 0.667368 -0.247861 0.05 + vertex 0.667368 -0.247861 -0.05 + endloop + endfacet + facet normal -0.0654098 -0.997858 -1.49195e-16 + outer loop + vertex 0.667368 -0.247861 -0.05 + vertex 0.7 -0.25 -0.05 + vertex 0.7 -0.25 0.05 + endloop + endfacet + facet normal -0.195096 -0.980784 0 + outer loop + vertex 0.667368 -0.247861 -0.05 + vertex 0.667368 -0.247861 0.05 + vertex 0.635295 -0.241481 0.05 + endloop + endfacet + facet normal -0.195096 -0.980784 -0 + outer loop + vertex 0.635295 -0.241481 0.05 + vertex 0.635295 -0.241481 -0.05 + vertex 0.667368 -0.247861 -0.05 + endloop + endfacet + facet normal -0.321425 -0.946935 0 + outer loop + vertex 0.635295 -0.241481 -0.05 + vertex 0.635295 -0.241481 0.05 + vertex 0.604329 -0.23097 0.05 + endloop + endfacet + facet normal -0.321425 -0.946935 -0 + outer loop + vertex 0.604329 -0.23097 0.05 + vertex 0.604329 -0.23097 -0.05 + vertex 0.635295 -0.241481 -0.05 + endloop + endfacet + facet normal -0.442302 -0.896866 0 + outer loop + vertex 0.604329 -0.23097 -0.05 + vertex 0.604329 -0.23097 0.05 + vertex 0.575 -0.216506 0.05 + endloop + endfacet + facet normal -0.442302 -0.896866 -0 + outer loop + vertex 0.575 -0.216506 0.05 + vertex 0.575 -0.216506 -0.05 + vertex 0.604329 -0.23097 -0.05 + endloop + endfacet + facet normal -0.555561 -0.831476 0 + outer loop + vertex 0.575 -0.216506 -0.05 + vertex 0.575 -0.216506 0.05 + vertex 0.547809 -0.198338 0.05 + endloop + endfacet + facet normal -0.555561 -0.831476 -0 + outer loop + vertex 0.547809 -0.198338 0.05 + vertex 0.547809 -0.198338 -0.05 + vertex 0.575 -0.216506 -0.05 + endloop + endfacet + facet normal -0.659372 -0.751817 0 + outer loop + vertex 0.547809 -0.198338 -0.05 + vertex 0.547809 -0.198338 0.05 + vertex 0.523224 -0.176776 0.05 + endloop + endfacet + facet normal -0.659372 -0.751817 -0 + outer loop + vertex 0.523224 -0.176776 0.05 + vertex 0.523224 -0.176776 -0.05 + vertex 0.547809 -0.198338 -0.05 + endloop + endfacet +endsolid cylinder diff --git a/tutorials/immersedBoundary/twoIbPatches/constant/turbulenceProperties b/tutorials/immersedBoundary/twoIbPatches/constant/turbulenceProperties new file mode 100644 index 000000000..6a2ea973c --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/constant/turbulenceProperties @@ -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; + + +// ************************************************************************* // diff --git a/tutorials/immersedBoundary/twoIbPatches/save/blockMeshDict b/tutorials/immersedBoundary/twoIbPatches/save/blockMeshDict new file mode 100644 index 000000000..5ed7f91cc --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/save/blockMeshDict @@ -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 +( +); + + +// ************************************************************************* // diff --git a/tutorials/immersedBoundary/twoIbPatches/save/boundary b/tutorials/immersedBoundary/twoIbPatches/save/boundary new file mode 100644 index 000000000..24a03b663 --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/save/boundary @@ -0,0 +1,71 @@ +/*--------------------------------*- 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; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +7 +( + ibCylinder + { + type immersedBoundary; + nFaces 0; + startFace 3650; + + internalFlow no; + isWall yes; + } + ibCylinder2 + { + 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; + } +) + + +// ************************************************************************* // diff --git a/tutorials/immersedBoundary/twoIbPatches/save/k b/tutorials/immersedBoundary/twoIbPatches/save/k new file mode 100644 index 000000000..54a010cc2 --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/save/k @@ -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; + } +} + + +// ************************************************************************* // diff --git a/tutorials/immersedBoundary/twoIbPatches/save/nut b/tutorials/immersedBoundary/twoIbPatches/save/nut new file mode 100644 index 000000000..b262b36c2 --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/save/nut @@ -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; + } +} + + +// ************************************************************************* // diff --git a/tutorials/immersedBoundary/twoIbPatches/save/omega b/tutorials/immersedBoundary/twoIbPatches/save/omega new file mode 100644 index 000000000..1840db118 --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/save/omega @@ -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; + } +} + + +// ************************************************************************* // diff --git a/tutorials/immersedBoundary/twoIbPatches/save/p b/tutorials/immersedBoundary/twoIbPatches/save/p new file mode 100644 index 000000000..7d3f2768d --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/save/p @@ -0,0 +1,1938 @@ +/*--------------------------------*- 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 volScalarField; + location "0"; + object p; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 2 -2 0 0 0 0]; + +internalField nonuniform List +1875 +( +-0.00144758 +-0.00443607 +-0.00770641 +-0.011461 +-0.0159321 +-0.0213957 +-0.0281872 +-0.0367211 +-0.0475127 +-0.0612051 +-0.0785984 +-0.100682 +-0.12866 +-0.163972 +-0.208273 +-0.263358 +-0.330973 +-0.412451 +-0.508106 +-0.616377 +-0.732832 +-0.849441 +-0.954733 +-1.03542 +-1.07943 +-1.07944 +-1.03543 +-0.954751 +-0.849468 +-0.732871 +-0.616429 +-0.508175 +-0.412539 +-0.331085 +-0.2635 +-0.208452 +-0.164197 +-0.128943 +-0.101038 +-0.0790489 +-0.0617759 +-0.0482376 +-0.0376434 +-0.0293626 +-0.0228954 +-0.017848 +-0.0139104 +-0.0108398 +-0.00844596 +-0.00658015 +-0.00512613 +-0.00399316 +-0.00311045 +-0.00242279 +-0.00188711 +-0.00146983 +-0.00114481 +-0.000891677 +-0.000694535 +-0.000541002 +-0.000421442 +-0.000328356 +-0.000255904 +-0.000199535 +-0.000155701 +-0.000121647 +-9.52346e-05 +-7.48123e-05 +-5.91028e-05 +-4.71203e-05 +-3.81083e-05 +-3.14977e-05 +-2.68736e-05 +-2.39471e-05 +-2.25331e-05 +-0.00136684 +-0.00417312 +-0.00724599 +-0.0107771 +-0.014987 +-0.0201391 +-0.0265555 +-0.0346367 +-0.0448849 +-0.0579333 +-0.07458 +-0.0958261 +-0.122915 +-0.157364 +-0.20097 +-0.25576 +-0.323822 +-0.406947 +-0.505981 +-0.619827 +-0.744195 +-0.870531 +-0.985994 +-1.07534 +-1.1244 +-1.1244 +-1.07535 +-0.98601 +-0.870554 +-0.744229 +-0.619873 +-0.506042 +-0.407026 +-0.323924 +-0.255889 +-0.201133 +-0.157569 +-0.123175 +-0.0961548 +-0.0749969 +-0.0584629 +-0.0455587 +-0.0354955 +-0.0276515 +-0.0215391 +-0.016777 +-0.0130673 +-0.0101777 +-0.00792689 +-0.00617381 +-0.00480838 +-0.00374491 +-0.00291664 +-0.00227156 +-0.00176914 +-0.00137783 +-0.0010731 +-0.000835791 +-0.000650982 +-0.000507057 +-0.000394985 +-0.000307737 +-0.000239835 +-0.000187005 +-0.000145919 +-0.000113998 +-8.9242e-05 +-7.0104e-05 +-5.5385e-05 +-4.41573e-05 +-3.57112e-05 +-2.95153e-05 +-2.51825e-05 +-2.24401e-05 +-2.11121e-05 +-0.00120851 +-0.00366062 +-0.00634931 +-0.00944493 +-0.0131452 +-0.0176881 +-0.0233685 +-0.0305577 +-0.03973 +-0.0514952 +-0.0666408 +-0.0861834 +-0.11143 +-0.14404 +-0.186082 +-0.240034 +-0.308688 +-0.394819 +-0.500491 +-0.625799 +-0.767024 +-0.914642 +-1.05273 +-1.16165 +-1.22233 +-1.22233 +-1.16165 +-1.05274 +-0.914658 +-0.767048 +-0.625833 +-0.500537 +-0.394881 +-0.308768 +-0.240138 +-0.186215 +-0.14421 +-0.111646 +-0.08646 +-0.0669937 +-0.051946 +-0.0403063 +-0.0312948 +-0.0243119 +-0.0188962 +-0.0146928 +-0.0114282 +-0.00889132 +-0.00691908 +-0.00538524 +-0.00419199 +-0.0032635 +-0.00254087 +-0.00197837 +-0.00154047 +-0.00119956 +-0.000934149 +-0.000727501 +-0.000566586 +-0.000441285 +-0.00034373 +-0.000267797 +-0.000208703 +-0.000162723 +-0.000126963 +-9.91813e-05 +-7.76409e-05 +-6.09925e-05 +-4.81884e-05 +-3.84195e-05 +-3.10694e-05 +-2.56773e-05 +-2.19058e-05 +-1.95169e-05 +-1.83581e-05 +-0.000979084 +-0.00292516 +-0.00506423 +-0.00753563 +-0.010503 +-0.0141668 +-0.0187796 +-0.0246672 +-0.0322566 +-0.0421143 +-0.0549975 +-0.0719244 +-0.094265 +-0.123854 +-0.163118 +-0.215187 +-0.283927 +-0.373742 +-0.488915 +-0.63208 +-0.801349 +-0.986126 +-1.16464 +-1.30958 +-1.39241 +-1.39241 +-1.30958 +-1.16464 +-0.98613 +-0.801358 +-0.632095 +-0.48894 +-0.373779 +-0.283978 +-0.215257 +-0.16321 +-0.123975 +-0.0944228 +-0.0721291 +-0.0552622 +-0.0424559 +-0.032697 +-0.0252342 +-0.0195092 +-0.015105 +-0.0117091 +-0.00908541 +-0.00705507 +-0.00548184 +-0.00426153 +-0.00331418 +-0.00257822 +-0.00200616 +-0.00156131 +-0.00121529 +-0.000946084 +-0.000736602 +-0.000573551 +-0.000446619 +-0.000347807 +-0.000270897 +-0.000211042 +-0.00016446 +-0.000128212 +-0.000100024 +-7.81313e-05 +-6.11612e-05 +-4.80468e-05 +-3.79601e-05 +-3.02636e-05 +-2.44718e-05 +-2.02216e-05 +-1.7247e-05 +-1.53618e-05 +-1.4447e-05 +-0.000688848 +-0.00200626 +-0.00346141 +-0.00515388 +-0.00720334 +-0.00976023 +-0.0130203 +-0.0172453 +-0.0227918 +-0.030154 +-0.0400252 +-0.0533864 +-0.0716339 +-0.0967539 +-0.131554 +-0.179941 +-0.247192 +-0.340082 +-0.466552 +-0.634267 +-0.846631 +-1.09328 +-1.33985 +-1.54889 +-1.67389 +-1.67388 +-1.54888 +-1.33984 +-1.09327 +-0.846619 +-0.63426 +-0.466551 +-0.340089 +-0.24721 +-0.179971 +-0.1316 +-0.0968185 +-0.0717229 +-0.0535067 +-0.0401854 +-0.0303657 +-0.0230694 +-0.0176074 +-0.0134911 +-0.0103705 +-0.00799283 +-0.00617348 +-0.00477648 +-0.0037007 +-0.00287037 +-0.00222829 +-0.00173103 +-0.00134544 +-0.00104618 +-0.000813767 +-0.000633172 +-0.000492769 +-0.000383561 +-0.000298595 +-0.000232487 +-0.000181051 +-0.000141026 +-0.000109876 +-8.5641e-05 +-6.67998e-05 +-5.21713e-05 +-4.08351e-05 +-3.20761e-05 +-2.53403e-05 +-2.02004e-05 +-1.63308e-05 +-1.34886e-05 +-1.14982e-05 +-1.02375e-05 +-9.62736e-06 +-0.000352233 +-0.000955516 +-0.00163221 +-0.00243519 +-0.00343163 +-0.0047113 +-0.00639889 +-0.00867253 +-0.0117922 +-0.0161437 +-0.0223068 +-0.03116 +-0.0440424 +-0.0629961 +-0.0911234 +-0.133088 +-0.195765 +-0.288975 +-0.426064 +-0.623851 +-0.900271 +-1.25372 +-1.61013 +-1.93646 +-2.14909 +-2.14909 +-1.93644 +-1.6101 +-1.25369 +-0.900234 +-0.623816 +-0.426035 +-0.288952 +-0.195749 +-0.133078 +-0.0911208 +-0.0630024 +-0.0440597 +-0.0311912 +-0.0223557 +-0.0162153 +-0.0118929 +-0.00881053 +-0.0065848 +-0.0049588 +-0.00375829 +-0.00286364 +-0.00219154 +-0.00168317 +-0.00129646 +-0.00100091 +-0.000774154 +-0.000599625 +-0.000464974 +-0.000360901 +-0.000280339 +-0.000217889 +-0.000169425 +-0.000131787 +-0.000102545 +-7.9814e-05 +-6.21347e-05 +-4.83811e-05 +-3.76857e-05 +-2.9376e-05 +-2.29284e-05 +-1.79359e-05 +-1.40826e-05 +-1.11218e-05 +-8.8615e-06 +-7.15611e-06 +-5.90059e-06 +-5.02165e-06 +-4.46802e-06 +-4.20356e-06 +1.21647e-05 +0.000165199 +0.000314747 +0.000459229 +0.000589983 +0.000686231 +0.000706462 +0.00057386 +0.000151815 +-0.00079702 +-0.00267868 +-0.00618014 +-0.0124538 +-0.0234181 +-0.0422405 +-0.0740982 +-0.127328 +-0.215053 +-0.35728 +-0.584132 +-0.951306 +-1.522 +-1.86108 +-2.29373 +-2.66051 +-2.6605 +-2.2937 +-1.86103 +-1.52194 +-0.95124 +-0.58407 +-0.357225 +-0.215003 +-0.127281 +-0.0740525 +-0.0421945 +-0.0233704 +-0.0124033 +-0.00612551 +-0.00261849 +-0.000729524 +0.00022883 +0.000663197 +0.000811643 +0.000811744 +0.000741545 +0.000644098 +0.000542167 +0.000446927 +0.000363093 +0.000291927 +0.000232939 +0.000184835 +0.000146041 +0.000115005 +9.03281e-05 +7.08091e-05 +5.54312e-05 +4.33487e-05 +3.38754e-05 +2.64637e-05 +2.06762e-05 +1.61623e-05 +1.26434e-05 +9.90086e-06 +7.76338e-06 +6.09735e-06 +4.80133e-06 +3.80117e-06 +3.04148e-06 +2.47621e-06 +2.0645e-06 +1.77348e-06 +1.58269e-06 +1.48425e-06 +0.0003821 +0.00128652 +0.00225876 +0.00335008 +0.00461293 +0.00610007 +0.00786096 +0.00993332 +0.0123261 +0.0149874 +0.0177458 +0.0202059 +0.0215616 +0.0202694 +0.0134729 +-0.00399936 +-0.0413509 +-0.114598 +-0.250117 +-0.482076 +-0.916583 +-1.23485 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +-1.23477 +-0.916492 +-0.481988 +-0.250042 +-0.11453 +-0.0412827 +-0.00392715 +0.0135531 +0.0203616 +0.0216703 +0.0203365 +0.017905 +0.0151835 +0.0125698 +0.0102382 +0.00824458 +0.00658483 +0.00522755 +0.00413141 +0.00325407 +0.00255645 +0.00200443 +0.00156924 +0.00122712 +0.000958735 +0.000748525 +0.000584077 +0.000455558 +0.000355205 +0.000276898 +0.000215822 +0.000168207 +0.000131103 +0.000102203 +7.97015e-05 +6.2189e-05 +4.85651e-05 +3.79737e-05 +2.97511e-05 +2.33871e-05 +1.84908e-05 +1.47604e-05 +1.19591e-05 +9.9033e-06 +8.45827e-06 +7.53482e-06 +7.08117e-06 +0.000732885 +0.00233566 +0.00407417 +0.00605045 +0.0083765 +0.011178 +0.0145967 +0.0187898 +0.0239239 +0.0301566 +0.0375982 +0.0462307 +0.055751 +0.0652649 +0.0727004 +0.0736636 +0.0591751 +0.0109813 +-0.109255 +-0.368743 +-0.742767 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +-0.742681 +-0.368646 +-0.109175 +0.0110554 +0.059252 +0.0737502 +0.0728027 +0.0653893 +0.0559048 +0.0464232 +0.0378405 +0.0304633 +0.0243134 +0.0192858 +0.0152295 +0.0119865 +0.0094108 +0.00737478 +0.00577105 +0.00451117 +0.00352338 +0.00275012 +0.00214552 +0.0016732 +0.00130446 +0.00101675 +0.00079235 +0.000617402 +0.000481042 +0.000374782 +0.000291998 +0.000227523 +0.000177323 +0.000138247 +0.000107838 +8.41834e-05 +6.57999e-05 +5.1537e-05 +4.0507e-05 +3.20241e-05 +2.55569e-05 +2.06933e-05 +1.71207e-05 +1.46125e-05 +1.30161e-05 +1.22379e-05 +0.00103926 +0.00324147 +0.00563897 +0.00837869 +0.0116258 +0.0155723 +0.0204453 +0.0265152 +0.0341017 +0.0435744 +0.0553399 +0.0698006 +0.0872543 +0.107667 +0.130186 +0.152091 +0.166554 +0.157874 +0.0917091 +-0.108466 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +-0.108375 +0.0917777 +0.157939 +0.166627 +0.15218 +0.130299 +0.107811 +0.0874393 +0.0700389 +0.0556466 +0.0439691 +0.0346094 +0.027168 +0.0212846 +0.0166511 +0.0130123 +0.0101605 +0.00792885 +0.0061845 +0.0048222 +0.00375898 +0.00292957 +0.00228281 +0.00177861 +0.00138564 +0.00107943 +0.00084085 +0.000654985 +0.000510207 +0.000397454 +0.000309662 +0.000241316 +0.000188115 +0.000146714 +0.000114512 +8.94921e-05 +7.00875e-05 +5.50863e-05 +4.35494e-05 +3.47493e-05 +2.81272e-05 +2.32625e-05 +1.98507e-05 +1.76838e-05 +1.66309e-05 +0.00127776 +0.00394021 +0.00684447 +0.0101727 +0.0141325 +0.0189687 +0.0249782 +0.032525 +0.042059 +0.0541334 +0.0694204 +0.0887137 +0.112897 +0.142827 +0.179011 +0.22078 +0.264199 +0.296613 +0.28465 +0.0821147 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.0821964 +0.284695 +0.296659 +0.264259 +0.220863 +0.179125 +0.14298 +0.113101 +0.0889823 +0.0697719 +0.0545911 +0.0426529 +0.0332938 +0.0259713 +0.0202501 +0.0157842 +0.0123004 +0.00958394 +0.00746648 +0.00581634 +0.0045306 +0.00352892 +0.0027486 +0.00214078 +0.00166734 +0.0012986 +0.00101141 +0.000787737 +0.000613556 +0.000477933 +0.000372346 +0.000290146 +0.000226161 +0.000176369 +0.000137645 +0.000107564 +8.42392e-05 +6.62087e-05 +5.23398e-05 +4.17578e-05 +3.37934e-05 +2.79443e-05 +2.38454e-05 +2.12447e-05 +1.99824e-05 +0.00142925 +0.00438127 +0.00760469 +0.0113043 +0.0157148 +0.0211156 +0.027849 +0.0363414 +0.04713 +0.0608937 +0.0784899 +0.100992 +0.129715 +0.166196 +0.212035 +0.26837 +0.334335 +0.402657 +0.431935 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.431951 +0.402681 +0.334381 +0.268445 +0.212146 +0.166352 +0.129928 +0.101277 +0.0788675 +0.061389 +0.0477762 +0.0371811 +0.0289369 +0.0225224 +0.0175312 +0.0136472 +0.0106245 +0.0082718 +0.00644043 +0.00501476 +0.00390482 +0.00304065 +0.0023678 +0.00184389 +0.00143594 +0.00111827 +0.000870906 +0.000678307 +0.000528357 +0.000411616 +0.000320731 +0.000249986 +0.000194937 +0.000152132 +0.000118883 +9.31032e-05 +7.31734e-05 +5.78415e-05 +4.61423e-05 +3.73379e-05 +3.08739e-05 +2.63459e-05 +2.34733e-05 +2.20784e-05 +0.0014812 +0.00453203 +0.00786443 +0.011691 +0.0162557 +0.02185 +0.028832 +0.0376501 +0.0488721 +0.0632217 +0.081623 +0.105251 +0.135578 +0.174389 +0.223677 +0.285146 +0.358574 +0.436929 +0.492123 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.492125 +0.436945 +0.358615 +0.285217 +0.223786 +0.174546 +0.135793 +0.105541 +0.0820091 +0.0637295 +0.0495357 +0.0385136 +0.029952 +0.0232992 +0.018128 +0.014107 +0.0109795 +0.00854642 +0.00665318 +0.00517976 +0.00403291 +0.00314015 +0.00244513 +0.00190402 +0.0014827 +0.00115466 +0.000899232 +0.000700364 +0.000545534 +0.000424988 +0.000331141 +0.000258092 +0.000201257 +0.000157065 +0.000122738 +9.612e-05 +7.5541e-05 +5.97094e-05 +4.76298e-05 +3.85404e-05 +3.18679e-05 +2.71936e-05 +2.42268e-05 +2.27849e-05 +0.00142925 +0.00438126 +0.00760465 +0.0113043 +0.0157148 +0.0211156 +0.0278489 +0.0363414 +0.04713 +0.0608936 +0.0784899 +0.100992 +0.129715 +0.166196 +0.212035 +0.26837 +0.334335 +0.402657 +0.431935 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.431951 +0.402681 +0.334381 +0.268445 +0.212146 +0.166352 +0.129928 +0.101278 +0.0788675 +0.061389 +0.0477762 +0.0371811 +0.0289369 +0.0225224 +0.0175312 +0.0136472 +0.0106245 +0.00827181 +0.00644045 +0.00501477 +0.00390482 +0.00304065 +0.0023678 +0.00184388 +0.00143593 +0.00111826 +0.000870917 +0.000678327 +0.000528373 +0.000411618 +0.000320723 +0.000249976 +0.000194933 +0.000152132 +0.000118882 +9.30978e-05 +7.31631e-05 +5.78288e-05 +4.61297e-05 +3.73264e-05 +3.08627e-05 +2.63327e-05 +2.3456e-05 +2.20565e-05 +0.00127776 +0.00394019 +0.00684441 +0.0101727 +0.0141324 +0.0189687 +0.0249781 +0.032525 +0.042059 +0.0541334 +0.0694204 +0.0887137 +0.112897 +0.142827 +0.179011 +0.22078 +0.264199 +0.296613 +0.28465 +0.0821145 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.0821964 +0.284695 +0.296659 +0.264259 +0.220863 +0.179125 +0.14298 +0.113101 +0.0889823 +0.0697719 +0.0545911 +0.0426529 +0.0332938 +0.0259713 +0.0202502 +0.0157843 +0.0123004 +0.00958396 +0.00746651 +0.00581636 +0.00453061 +0.00352891 +0.00274859 +0.00214076 +0.00166733 +0.00129859 +0.0010114 +0.000787751 +0.000613585 +0.000477961 +0.000372357 +0.000290142 +0.000226151 +0.000176361 +0.00013764 +0.000107556 +8.42262e-05 +6.61916e-05 +5.23202e-05 +4.17368e-05 +3.37707e-05 +2.79189e-05 +2.38156e-05 +2.12089e-05 +1.99403e-05 +0.00103925 +0.00324143 +0.00563891 +0.00837863 +0.0116258 +0.0155723 +0.0204452 +0.0265151 +0.0341016 +0.0435744 +0.0553399 +0.0698006 +0.0872543 +0.107667 +0.130186 +0.152091 +0.166554 +0.157874 +0.091709 +-0.108466 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +-0.108375 +0.0917777 +0.157939 +0.166627 +0.15218 +0.130299 +0.107811 +0.0874393 +0.0700389 +0.0556466 +0.0439691 +0.0346094 +0.027168 +0.0212846 +0.0166511 +0.0130123 +0.0101605 +0.00792888 +0.00618454 +0.00482223 +0.00375898 +0.00292956 +0.00228278 +0.00177859 +0.00138563 +0.00107942 +0.000840848 +0.000654998 +0.000510235 +0.000397487 +0.000309686 +0.000241326 +0.000188114 +0.000146704 +0.000114495 +8.94703e-05 +7.00654e-05 +5.50666e-05 +4.353e-05 +3.4725e-05 +2.80937e-05 +2.32192e-05 +1.97999e-05 +1.76278e-05 +1.6571e-05 +0.000732871 +0.00233561 +0.00407411 +0.00605039 +0.00837643 +0.0111779 +0.0145966 +0.0187898 +0.0239238 +0.0301565 +0.0375981 +0.0462307 +0.0557509 +0.0652649 +0.0727004 +0.0736636 +0.0591751 +0.0109812 +-0.109255 +-0.368744 +-0.742767 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +-0.742681 +-0.368646 +-0.109175 +0.0110554 +0.0592521 +0.0737502 +0.0728027 +0.0653893 +0.0559048 +0.0464232 +0.0378406 +0.0304634 +0.0243134 +0.0192859 +0.0152295 +0.0119865 +0.00941083 +0.00737481 +0.00577109 +0.0045112 +0.00352341 +0.00275013 +0.0021455 +0.00167316 +0.00130443 +0.00101673 +0.000792345 +0.000617409 +0.000481061 +0.000374813 +0.000292035 +0.000227558 +0.000177348 +0.000138255 +0.000107826 +8.41558e-05 +6.57656e-05 +5.15077e-05 +4.04873e-05 +3.20078e-05 +2.55314e-05 +2.06495e-05 +1.70585e-05 +1.45391e-05 +1.29394e-05 +1.21618e-05 +0.00038209 +0.00128649 +0.00225871 +0.00335001 +0.00461285 +0.0061 +0.00786087 +0.00993324 +0.012326 +0.0149874 +0.0177458 +0.0202058 +0.0215616 +0.0202694 +0.0134729 +-0.00399938 +-0.0413509 +-0.114598 +-0.250117 +-0.482076 +-0.916583 +-1.23485 +0.5 +0.5 +0.5 +0.5 +0.5 +0.5 +-1.23477 +-0.916492 +-0.481988 +-0.250042 +-0.11453 +-0.0412827 +-0.00392713 +0.0135531 +0.0203616 +0.0216704 +0.0203365 +0.0179051 +0.0151836 +0.0125698 +0.0102383 +0.00824462 +0.00658485 +0.00522757 +0.00413144 +0.00325411 +0.00255648 +0.00200446 +0.00156925 +0.00122711 +0.000958701 +0.000748482 +0.000584046 +0.000455552 +0.000355223 +0.000276931 +0.000215864 +0.000168251 +0.000131142 +0.00010223 +7.9709e-05 +6.21742e-05 +4.85325e-05 +3.79353e-05 +2.97198e-05 +2.33674e-05 +1.84743e-05 +1.47311e-05 +1.19054e-05 +9.82506e-06 +8.36538e-06 +7.43918e-06 +6.98929e-06 +1.21534e-05 +0.000165157 +0.000314672 +0.000459138 +0.000589901 +0.000686146 +0.000706352 +0.000573759 +0.000151749 +-0.000797056 +-0.0026787 +-0.00618014 +-0.0124538 +-0.0234181 +-0.0422406 +-0.0740983 +-0.127328 +-0.215053 +-0.357281 +-0.584132 +-0.951306 +-1.522 +-1.86108 +-2.29373 +-2.66051 +-2.6605 +-2.2937 +-1.86103 +-1.52194 +-0.95124 +-0.58407 +-0.357225 +-0.215003 +-0.127281 +-0.0740525 +-0.0421945 +-0.0233704 +-0.0124032 +-0.00612544 +-0.00261841 +-0.000729425 +0.000228927 +0.000663262 +0.000811684 +0.000811773 +0.000741566 +0.00064412 +0.00054219 +0.000446946 +0.000363109 +0.000291939 +0.000232937 +0.000184808 +0.000146001 +0.000114976 +9.03287e-05 +7.08393e-05 +5.54777e-05 +4.33993e-05 +3.39225e-05 +2.65007e-05 +2.06963e-05 +1.61615e-05 +1.26235e-05 +9.86909e-06 +7.72962e-06 +6.06927e-06 +4.77962e-06 +3.77758e-06 +3.0028e-06 +2.41245e-06 +1.97532e-06 +1.66776e-06 +1.47238e-06 +1.37686e-06 +-0.000352246 +-0.000955574 +-0.00163231 +-0.0024353 +-0.00343173 +-0.00471139 +-0.00639898 +-0.00867264 +-0.0117923 +-0.0161438 +-0.0223068 +-0.03116 +-0.0440424 +-0.0629961 +-0.0911234 +-0.133088 +-0.195765 +-0.288975 +-0.426064 +-0.623851 +-0.900271 +-1.25372 +-1.61013 +-1.93646 +-2.14909 +-2.14909 +-1.93644 +-1.6101 +-1.25369 +-0.900234 +-0.623816 +-0.426035 +-0.288952 +-0.195749 +-0.133078 +-0.0911207 +-0.0630023 +-0.0440597 +-0.0311911 +-0.0223556 +-0.0162152 +-0.0118928 +-0.00881045 +-0.00658476 +-0.00495877 +-0.00375828 +-0.00286363 +-0.00219152 +-0.00168316 +-0.00129645 +-0.0010009 +-0.000774144 +-0.000599637 +-0.000465001 +-0.000360916 +-0.000280321 +-0.000217844 +-0.000169371 +-0.000131738 +-0.000102505 +-7.97859e-05 +-6.21229e-05 +-4.83883e-05 +-3.77074e-05 +-2.94033e-05 +-2.29547e-05 +-1.79605e-05 +-1.41097e-05 +-1.11584e-05 +-8.91453e-06 +-7.23034e-06 +-5.99603e-06 +-5.13321e-06 +-4.58771e-06 +-4.3251e-06 +-0.000688859 +-0.00200633 +-0.00346153 +-0.00515401 +-0.00720346 +-0.00976034 +-0.0130204 +-0.0172454 +-0.0227919 +-0.0301541 +-0.0400252 +-0.0533864 +-0.0716339 +-0.0967539 +-0.131554 +-0.179941 +-0.247192 +-0.340082 +-0.466552 +-0.634267 +-0.846631 +-1.09328 +-1.33985 +-1.54889 +-1.67389 +-1.67388 +-1.54888 +-1.33984 +-1.09327 +-0.846619 +-0.63426 +-0.466551 +-0.340089 +-0.24721 +-0.179971 +-0.1316 +-0.0968184 +-0.0717228 +-0.0535066 +-0.0401853 +-0.0303655 +-0.0230692 +-0.0176074 +-0.0134911 +-0.0103705 +-0.00799283 +-0.00617349 +-0.00477649 +-0.00370071 +-0.00287037 +-0.00222828 +-0.001731 +-0.00134543 +-0.00104619 +-0.000813762 +-0.000633135 +-0.000492706 +-0.000383498 +-0.000298551 +-0.000232462 +-0.000181038 +-0.000141021 +-0.000109882 +-8.56555e-05 +-6.68187e-05 +-5.21922e-05 +-4.08611e-05 +-3.21132e-05 +-2.53927e-05 +-2.02692e-05 +-1.64154e-05 +-1.35882e-05 +-1.16117e-05 +-1.03623e-05 +-9.75997e-06 +-0.000979096 +-0.00292521 +-0.00506434 +-0.00753576 +-0.0105032 +-0.0141669 +-0.0187797 +-0.0246673 +-0.0322567 +-0.0421143 +-0.0549975 +-0.0719244 +-0.094265 +-0.123854 +-0.163118 +-0.215187 +-0.283927 +-0.373742 +-0.488915 +-0.63208 +-0.801349 +-0.986126 +-1.16464 +-1.30958 +-1.39241 +-1.3924 +-1.30958 +-1.16464 +-0.98613 +-0.801358 +-0.632095 +-0.48894 +-0.373779 +-0.283978 +-0.215257 +-0.16321 +-0.123975 +-0.0944227 +-0.072129 +-0.055262 +-0.0424557 +-0.0326968 +-0.0252341 +-0.0195091 +-0.015105 +-0.0117091 +-0.00908543 +-0.0070551 +-0.00548187 +-0.00426154 +-0.00331416 +-0.00257818 +-0.00200613 +-0.00156129 +-0.00121527 +-0.000946041 +-0.000736535 +-0.000573484 +-0.000446576 +-0.000347792 +-0.000270898 +-0.000211044 +-0.000164458 +-0.000128212 +-0.00010003 +-7.81484e-05 +-6.11938e-05 +-4.80975e-05 +-3.80282e-05 +-3.0346e-05 +-2.45656e-05 +-2.03257e-05 +-1.73626e-05 +-1.54898e-05 +-1.45864e-05 +-0.00120855 +-0.00366067 +-0.00634938 +-0.00944505 +-0.0131453 +-0.0176883 +-0.0233686 +-0.0305579 +-0.0397301 +-0.0514953 +-0.0666409 +-0.0861835 +-0.11143 +-0.14404 +-0.186082 +-0.240034 +-0.308688 +-0.394819 +-0.500491 +-0.625799 +-0.767024 +-0.914642 +-1.05273 +-1.16165 +-1.22233 +-1.22233 +-1.16165 +-1.05274 +-0.914658 +-0.767048 +-0.625833 +-0.500537 +-0.394881 +-0.308768 +-0.240138 +-0.186215 +-0.14421 +-0.111646 +-0.0864598 +-0.0669936 +-0.0519459 +-0.0403061 +-0.0312947 +-0.0243118 +-0.0188962 +-0.0146928 +-0.0114282 +-0.00889137 +-0.00691911 +-0.00538524 +-0.00419197 +-0.00326346 +-0.00254082 +-0.00197832 +-0.00154043 +-0.00119952 +-0.000934094 +-0.000727439 +-0.000566539 +-0.000441265 +-0.000343731 +-0.000267799 +-0.000208696 +-0.000162709 +-0.000126955 +-9.91938e-05 +-7.76802e-05 +-6.10565e-05 +-4.827e-05 +-3.85123e-05 +-3.11707e-05 +-2.57874e-05 +-2.20261e-05 +-1.96482e-05 +-1.85e-05 +-0.00136688 +-0.00417321 +-0.0072461 +-0.0107772 +-0.0149871 +-0.0201392 +-0.0265557 +-0.0346368 +-0.044885 +-0.0579334 +-0.0745801 +-0.0958261 +-0.122915 +-0.157364 +-0.20097 +-0.255761 +-0.323822 +-0.406947 +-0.505981 +-0.619827 +-0.744195 +-0.870531 +-0.985994 +-1.07534 +-1.1244 +-1.1244 +-1.07535 +-0.98601 +-0.870554 +-0.744229 +-0.619873 +-0.506042 +-0.407026 +-0.323924 +-0.255889 +-0.201133 +-0.157569 +-0.123174 +-0.0961547 +-0.0749967 +-0.0584627 +-0.0455586 +-0.0354954 +-0.0276514 +-0.0215391 +-0.0167771 +-0.0130674 +-0.0101777 +-0.00792692 +-0.0061738 +-0.00480836 +-0.00374488 +-0.0029166 +-0.0022715 +-0.00176908 +-0.00137779 +-0.00107306 +-0.000835741 +-0.000650936 +-0.000507028 +-0.000394975 +-0.000307733 +-0.000239823 +-0.000186984 +-0.000145905 +-0.000114007 +-8.92835e-05 +-7.01754e-05 +-5.54752e-05 +-4.4257e-05 +-3.58182e-05 +-2.96318e-05 +-2.53094e-05 +-2.25751e-05 +-2.12527e-05 +-0.0014476 +-0.00443613 +-0.00770654 +-0.0114612 +-0.0159322 +-0.0213958 +-0.0281874 +-0.0367213 +-0.0475129 +-0.0612052 +-0.0785984 +-0.100682 +-0.12866 +-0.163972 +-0.208273 +-0.263358 +-0.330973 +-0.412451 +-0.508107 +-0.616377 +-0.732832 +-0.849441 +-0.954733 +-1.03542 +-1.07943 +-1.07944 +-1.03543 +-0.954751 +-0.849468 +-0.732871 +-0.616429 +-0.508175 +-0.412539 +-0.331085 +-0.2635 +-0.208452 +-0.164197 +-0.128942 +-0.101038 +-0.0790487 +-0.0617758 +-0.0482375 +-0.0376433 +-0.0293625 +-0.0228954 +-0.017848 +-0.0139105 +-0.0108398 +-0.008446 +-0.00658015 +-0.00512611 +-0.00399313 +-0.00311042 +-0.00242274 +-0.00188705 +-0.00146978 +-0.00114477 +-0.000891633 +-0.000694493 +-0.00054097 +-0.000421422 +-0.000328343 +-0.00025589 +-0.000199517 +-0.000155689 +-0.000121656 +-9.52744e-05 +-7.4882e-05 +-5.9193e-05 +-4.72216e-05 +-3.82179e-05 +-3.1618e-05 +-2.70054e-05 +-2.40853e-05 +-2.26701e-05 +) +; + +boundaryField +{ + ibCylinder + { + type mixedIb; + patchType immersedBoundary; + triValue uniform 0; + triGradient uniform 0; + triValueFraction uniform 0; + setDeadValue yes; + deadValue 0; + value nonuniform 0(); + } + in + { + type zeroGradient; + } + out + { + type fixedValue; + value uniform 0; + } + top + { + type zeroGradient; + } + bottom + { + type zeroGradient; + } + frontAndBack + { + type empty; + } +} + + +// ************************************************************************* // diff --git a/tutorials/immersedBoundary/twoIbPatches/system/controlDict b/tutorials/immersedBoundary/twoIbPatches/system/controlDict new file mode 100644 index 000000000..1c558907a --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/system/controlDict @@ -0,0 +1,152 @@ +/*--------------------------------*- 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; + // immersedBoundaryFvPatch 10; + // polyMesh 2; + // fvMesh 2; + /* + immersedBoundary 2; + immersedBoundaryCorrectedMeshFields 2; + immersedBoundaryEpsilonWallFunction 2; + immersedBoundaryKqRWallFunction 2; + immersedBoundaryNutWallFunction 2; + immersedBoundaryOmegaWallFunction 2; + foo 2; + */ +} + +application pimpleDyMIbFoam; + +startFrom startTime; + +startTime 0; + +stopAt endTime; + +endTime 5; + +deltaT 0.01; + +// writeControl runTime; +// writeInterval 0.2; + +writeControl timeStep; +writeInterval 5; + +purgeWrite 0; + +writeFormat binary; + +writePrecision 6; + +writeCompression uncompressed; + +timeFormat general; + +timePrecision 6; + +runTimeModifiable yes; + +adjustTimeStep no; + +maxCo 0.2; + +libs +( + "libimmersedBoundary.so" + "libimmersedBoundaryDynamicMesh.so" + // "libimmersedBoundaryTurbulence.so" +"libincompressibleImmersedBoundaryTurbulence.so" +); + +functions +( + forces + { + type forces; + functionObjectLibs ("libforces.so"); + + outputControl timeStep; + outputInterval 1; + patches ( ibCylinder2); + + 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; + } + */ +); + + +// ************************************************************************* // diff --git a/tutorials/immersedBoundary/twoIbPatches/system/decomposeParDict b/tutorials/immersedBoundary/twoIbPatches/system/decomposeParDict new file mode 100644 index 000000000..4eda43bf8 --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/system/decomposeParDict @@ -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 +( +); + + +// ************************************************************************* // diff --git a/tutorials/immersedBoundary/twoIbPatches/system/fvSchemes b/tutorials/immersedBoundary/twoIbPatches/system/fvSchemes new file mode 100644 index 000000000..cba11a7d9 --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/system/fvSchemes @@ -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 linearUpwindV grad(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; +} + + +// ************************************************************************* // diff --git a/tutorials/immersedBoundary/twoIbPatches/system/fvSolution b/tutorials/immersedBoundary/twoIbPatches/system/fvSolution new file mode 100644 index 000000000..14ff72bdc --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/system/fvSolution @@ -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 + { + } +} + + +// ************************************************************************* // diff --git a/tutorials/immersedBoundary/twoIbPatches/system/mapFieldsDict b/tutorials/immersedBoundary/twoIbPatches/system/mapFieldsDict new file mode 100644 index 000000000..5373b1b22 --- /dev/null +++ b/tutorials/immersedBoundary/twoIbPatches/system/mapFieldsDict @@ -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 +); + + +// ************************************************************************* // diff --git a/tutorials/immersedBoundary/twoIbPatches/twoIbPatches.foam b/tutorials/immersedBoundary/twoIbPatches/twoIbPatches.foam new file mode 100644 index 000000000..e69de29bb