From 7b74fd7c9590fa757cd28c124a4d7cc05df57f35 Mon Sep 17 00:00:00 2001 From: Hrvoje Jasak Date: Mon, 18 Oct 2010 22:30:28 +0100 Subject: [PATCH] Partial overlap GGI turbo passage --- .../icoDyMFoam/turboPassageRotating/0/U | 61 +++++++++ .../icoDyMFoam/turboPassageRotating/0/p | 58 ++++++++ .../icoDyMFoam/turboPassageRotating/Allclean | 9 ++ .../icoDyMFoam/turboPassageRotating/Allrun | 11 ++ .../constant/dynamicMeshDict | 39 ++++++ .../constant/polyMesh/blockMeshDict | 125 ++++++++++++++++++ .../constant/polyMesh/boundary | 78 +++++++++++ .../constant/transportProperties | 19 +++ .../turboPassageRotating/system/controlDict | 53 ++++++++ .../turboPassageRotating/system/fvSchemes | 66 +++++++++ .../turboPassageRotating/system/fvSolution | 75 +++++++++++ 11 files changed, 594 insertions(+) create mode 100644 tutorials/incompressible/icoDyMFoam/turboPassageRotating/0/U create mode 100644 tutorials/incompressible/icoDyMFoam/turboPassageRotating/0/p create mode 100755 tutorials/incompressible/icoDyMFoam/turboPassageRotating/Allclean create mode 100755 tutorials/incompressible/icoDyMFoam/turboPassageRotating/Allrun create mode 100644 tutorials/incompressible/icoDyMFoam/turboPassageRotating/constant/dynamicMeshDict create mode 100644 tutorials/incompressible/icoDyMFoam/turboPassageRotating/constant/polyMesh/blockMeshDict create mode 100644 tutorials/incompressible/icoDyMFoam/turboPassageRotating/constant/polyMesh/boundary create mode 100644 tutorials/incompressible/icoDyMFoam/turboPassageRotating/constant/transportProperties create mode 100644 tutorials/incompressible/icoDyMFoam/turboPassageRotating/system/controlDict create mode 100644 tutorials/incompressible/icoDyMFoam/turboPassageRotating/system/fvSchemes create mode 100644 tutorials/incompressible/icoDyMFoam/turboPassageRotating/system/fvSolution diff --git a/tutorials/incompressible/icoDyMFoam/turboPassageRotating/0/U b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/0/U new file mode 100644 index 000000000..8e1d66bf1 --- /dev/null +++ b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/0/U @@ -0,0 +1,61 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volVectorField; + object U; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 1 -1 0 0 0 0]; + +internalField uniform (9.6592582628906829 2.5881904510252074 0); + +boundaryField +{ + inlet + { + type surfaceNormalFixedValue; + refValue uniform -10; + value uniform (9.6592582628906829 2.5881904510252074 0); + } + outlet + { + type zeroGradient; + } + fixedWalls + { + type fixedValue; + value uniform (0 0 0); + } + movingwalls + { + type movingWallVelocity; + value uniform (0 0 0); + } + rotor_cyclics + { + type cyclic; + } + stator_cyclics + { + type cyclic; + } + interface1 + { + type overlapGgi; + } + interface2 + { + type overlapGgi; + } +} + +// ************************************************************************* // diff --git a/tutorials/incompressible/icoDyMFoam/turboPassageRotating/0/p b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/0/p new file mode 100644 index 000000000..7935fd878 --- /dev/null +++ b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/0/p @@ -0,0 +1,58 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class volScalarField; + object p; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dimensions [0 2 -2 0 0 0 0]; + +internalField uniform 0; + +boundaryField +{ + inlet + { + type zeroGradient; + } + outlet + { + type fixedValue; + value uniform 0; + } + fixedWalls + { + type zeroGradient; + } + movingwalls + { + type zeroGradient; + } + rotor_cyclics + { + type cyclic; + } + stator_cyclics + { + type cyclic; + } + interface1 + { + type overlapGgi; + } + interface2 + { + type overlapGgi; + } +} + +// ************************************************************************* // diff --git a/tutorials/incompressible/icoDyMFoam/turboPassageRotating/Allclean b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/Allclean new file mode 100755 index 000000000..0a256acee --- /dev/null +++ b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/Allclean @@ -0,0 +1,9 @@ +#!/bin/sh + +# Source tutorial clean functions +. $WM_PROJECT_DIR/bin/tools/CleanFunctions + +cleanCase +rm -rf {VTK,log*} +rm constant/polyMesh/{cellZones*,faces*,neighbour*,owner*,points*,sets} +find -iname "*~" | xargs rm -f diff --git a/tutorials/incompressible/icoDyMFoam/turboPassageRotating/Allrun b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/Allrun new file mode 100755 index 000000000..41229747c --- /dev/null +++ b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/Allrun @@ -0,0 +1,11 @@ +#!/bin/sh +# Source tutorial run functions +. $WM_PROJECT_DIR/bin/tools/RunFunctions + +application="icoDyMFoam" + +runApplication blockMesh +cp constant/polyMesh/boundary.orig constant/polyMesh/boundary +runApplication regionCellSets +runApplication setsToZones -noFlipMap +runApplication $application diff --git a/tutorials/incompressible/icoDyMFoam/turboPassageRotating/constant/dynamicMeshDict b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/constant/dynamicMeshDict new file mode 100644 index 000000000..9dd0c3494 --- /dev/null +++ b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/constant/dynamicMeshDict @@ -0,0 +1,39 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.openfoam.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ + +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object dynamicMeshDict; +} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +dynamicFvMesh turboFvMesh; + +turboFvMeshCoeffs +{ + coordinateSystem + { + type cylindrical; + origin (0 0 0); + axis (0 0 1); + direction (1 0 0); + } + + rpm + { + cellRegion0 60; + cellRegion1 0; + } +} + + +// ************************************************************************* // diff --git a/tutorials/incompressible/icoDyMFoam/turboPassageRotating/constant/polyMesh/blockMeshDict b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/constant/polyMesh/blockMeshDict new file mode 100644 index 000000000..f860a222d --- /dev/null +++ b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/constant/polyMesh/blockMeshDict @@ -0,0 +1,125 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object blockMeshDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +convertToMeters 0.01; + +vertices +( + ( 10. 0. 0. ) + ( 10. 0. 10. ) + ( 18. 0. 0. ) + ( 18. 0. 10. ) + ( 20. 0. 0. ) + ( 20. 0. 10. ) + ( 19.69615506 3.47296355 0. ) + ( 19.69615506 3.47296355 10. ) + ( 21.66577057 3.82025991 0. ) + ( 21.66577057 3.82025991 10. ) + ( 29.54423259 5.20944533 0. ) + ( 29.54423259 5.20944533 10. ) + ( 8.66025404 5. 0. ) + ( 8.66025404 5. 10. ) + ( 15.58845727 9. 0. ) + ( 15.58845727 9. 10. ) + ( 17.32050808 10. 0. ) + ( 17.32050808 10. 10. ) + ( 15.32088886 12.85575219 0. ) + ( 15.32088886 12.85575219 10. ) + ( 16.85297775 14.14132741 0. ) + ( 16.85297775 14.14132741 10. ) + ( 22.98133329 19.28362829 0. ) + ( 22.98133329 19.28362829 10. ) +); + +blocks +( + hex (0 12 13 1 2 14 15 3 ) (20 10 20) simpleGrading (1 1 1) + hex (2 14 15 3 4 16 17 5 ) (20 10 5) simpleGrading (1 1 1) + hex (6 18 19 7 8 20 21 9 ) (20 10 5) simpleGrading (1 1 1) + hex (8 20 21 9 10 22 23 11 ) (20 10 20) simpleGrading (1 1 1) +); + +edges +( +arc 0 12 ( 9.65925826 2.58819045 0. ) +arc 1 13 ( 9.65925826 2.58819045 10. ) +arc 2 14 ( 17.38666487 4.65874281 0. ) +arc 3 15 ( 17.38666487 4.65874281 10. ) +arc 4 16 ( 19.31851653 5.1763809 0. ) +arc 5 17 ( 19.31851653 5.1763809 10. ) +arc 6 18 ( 19.31851653 5.1763809 0. ) +arc 7 19 ( 19.31851653 5.1763809 10. ) +arc 8 20 ( 21.25036818 5.69401899 0. ) +arc 9 21 ( 21.25036818 5.69401899 10. ) +arc 10 22 ( 28.97777479 7.76457135 0. ) +arc 11 23 ( 28.97777479 7.76457135 10. ) +); + +patches +( + patch inlet + ( + (0 12 13 1) + ) + patch outlet + ( + (10 22 23 11) + ) + wall fixedWalls + ( + (6 18 20 8) + (7 19 21 9) + (8 20 22 10) + (9 21 23 11) + (20 22 23 21) + (8 10 11 9) + ) + patch movingwalls + ( + (2 14 16 4) + (3 15 17 5) + (0 12 14 2) + (1 13 15 3) + (12 14 15 13) + (0 2 3 1) + ) + cyclic rotor_cyclics + ( + (2 4 5 3) + (14 16 17 15) + ) + cyclic stator_cyclics + ( + (6 8 9 7) + (18 20 21 19) + ) + overlapGgi interface1 + ( + (4 5 17 16) + ) + overlapGgi interface2 + ( + (6 7 19 18) + ) +); + + +mergePatchPairs +( +); + +// ************************************************************************* // + diff --git a/tutorials/incompressible/icoDyMFoam/turboPassageRotating/constant/polyMesh/boundary b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/constant/polyMesh/boundary new file mode 100644 index 000000000..d97c073f7 --- /dev/null +++ b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/constant/polyMesh/boundary @@ -0,0 +1,78 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5-dev | +| \\ / A nd | Revision: exported | +| \\/ M anipulation | Web: http://www.OpenFOAM.org | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class polyBoundaryMesh; + location "constant/polyMesh"; + object boundary; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +8 +( + inlet + { + type patch; + nFaces 200; + startFace 28100; + } + outlet + { + type patch; + nFaces 200; + startFace 28300; + } + fixedWalls + { + type wall; + nFaces 1400; + startFace 28500; + } + movingwalls + { + type patch; + nFaces 1400; + startFace 29900; + } + rotor_cyclics + { + type cyclic; + nFaces 100; + startFace 31300; + featureCos 0.9; + } + stator_cyclics + { + type cyclic; + nFaces 100; + startFace 31400; + featureCos 0.9; + } + interface1 + { + type overlapGgi; + nFaces 200; + startFace 31500; + rotationAxis (0 0 1); + nCopies 12; + shadowPatch interface2; + } + interface2 + { + type overlapGgi; + nFaces 200; + startFace 31700; + rotationAxis (0 0 1); + nCopies 12; + shadowPatch interface1; + } +) + +// ************************************************************************* // diff --git a/tutorials/incompressible/icoDyMFoam/turboPassageRotating/constant/transportProperties b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/constant/transportProperties new file mode 100644 index 000000000..2f1250ed4 --- /dev/null +++ b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/constant/transportProperties @@ -0,0 +1,19 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object transportProperties; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // +transportModel Newtonian; +nu nu [0 2 -1 0 0 0 0] 0.02; + +// ************************************************************************* // diff --git a/tutorials/incompressible/icoDyMFoam/turboPassageRotating/system/controlDict b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/system/controlDict new file mode 100644 index 000000000..d6a1b3313 --- /dev/null +++ b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/system/controlDict @@ -0,0 +1,53 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.OpenFOAM.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object controlDict; +} +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +application icoDyMFoam; + +startFrom latestTime; + +startTime 0; + +stopAt endTime; + +endTime 1; + +deltaT 0.0001; + +writeControl adjustableRunTime; + +writeInterval 0.01; + +purgeWrite 0; + +writeFormat ascii; + +writePrecision 6; + +writeCompression compressed; + +timeFormat general; + +timePrecision 4; + +runTimeModifiable yes; +adjustTimeStep yes; + +maxCo 0.8; + +maxDeltaT 0.01; + + +// ************************************************************************* // diff --git a/tutorials/incompressible/icoDyMFoam/turboPassageRotating/system/fvSchemes b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/system/fvSchemes new file mode 100644 index 000000000..21fb55235 --- /dev/null +++ b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/system/fvSchemes @@ -0,0 +1,66 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.openfoam.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ + +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object fvSchemes; +} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +ddtSchemes +{ + default Euler; +} + +gradSchemes +{ + default Gauss linear; + grad(p) Gauss linear; +} + +divSchemes +{ + default none; + div(phi,U) Gauss upwind; +} + +laplacianSchemes +{ + default none; + laplacian(nu,U) Gauss linear corrected; + laplacian(rAU,pcorr) Gauss linear corrected; + laplacian(rAU,p) Gauss linear corrected; + + laplacian((1|A(U)),p) Gauss linear corrected; +} + +interpolationSchemes +{ + default linear; + interpolate(HbyA) linear; + interpolate(1|A) linear; +} + +snGradSchemes +{ + default corrected; +} + +fluxRequired +{ + default no; + pcorr; + p; +} + + +// ************************************************************************* // diff --git a/tutorials/incompressible/icoDyMFoam/turboPassageRotating/system/fvSolution b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/system/fvSolution new file mode 100644 index 000000000..8b80bdb05 --- /dev/null +++ b/tutorials/incompressible/icoDyMFoam/turboPassageRotating/system/fvSolution @@ -0,0 +1,75 @@ +/*--------------------------------*- C++ -*----------------------------------*\ +| ========= | | +| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | +| \\ / O peration | Version: 1.5 | +| \\ / A nd | Web: http://www.openfoam.org | +| \\/ M anipulation | | +\*---------------------------------------------------------------------------*/ + +FoamFile +{ + version 2.0; + format ascii; + class dictionary; + object fvSolution; +} + +// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // + +solvers +{ + pcorr + { + solver PCG; + preconditioner DIC; + + minIter 0; + maxIter 1000; + tolerance 0.01; + relTol 0; + } + + p + { + solver PCG; + preconditioner DIC; + + minIter 0; + maxIter 1000; + tolerance 1e-07; + relTol 0; + } + + pFinal + { + solver PCG; + preconditioner DIC; + + minIter 0; + maxIter 1000; + tolerance 1e-07; + relTol 0; + } + + U + { + solver BiCGStab; + preconditioner DILU; + + minIter 0; + maxIter 1000; + tolerance 1e-07; + relTol 0; + } +} + +PISO +{ + nCorrectors 4; + nNonOrthogonalCorrectors 0; + pRefCell 0; + pRefValue 0; +} + + +// ************************************************************************* //