Studying the Viscous Flow Around a Cylinder Using OpenFoam Part 2

File: /home/kmarc/Documents/ME702/OpenfoamScripts/blockMeshDict Page 1 of 5 /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: v1606+ | | \\ / A nd | Web: www.OpenFOAM.com | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; object blockMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // convertToMeters 0.1; vertices ( (0.5 0 -0.5) (1 0 -0.5) (10 0 -0.5) (10 0.707107 -0.5) (0.707107 0.707107 -0.5) (0.353553 0.353553 -0.5) (10 2 -0.5) (0.707107 2 -0.5) (0 2 -0.5) (0 1 -0.5) (0 0.5 -0.5) (-0.5 0 -0.5) (-1 0 -0.5) (-2 0 -0.5) (-2 0.707107 -0.5) (-0.707107 0.707107 -0.5) (-0.353553 0.353553 -0.5) (-2 2 -0.5) (-0.707107 2 -0.5) (0.5 0 0.5) (1 0 0.5) (10 0 0.5) (10 0.707107 0.5) (0.707107 0.707107 0.5) (0.353553 0.353553 0.5) (10 2 0.5) (0.707107 2 0.5) (0 2 0.5) (0 1 0.5) (0 0.5 0.5) (-0.5 0 0.5) (-1 0 0.5) (-2 0 0.5) (-2 0.707107 0.5) (-0.707107 0.707107 0.5) (-0.353553 0.353553 0.5) (-2 2 0.5) (-0.707107 2 0.5) File: /home/kmarc/Documents/ME702/OpenfoamScripts/blockMeshDict (10 -0.707107 -0.5) (0.707107 -0.707107 -0.5) (0.353553 -0.353553 -0.5) (10 -2 -0.5) (0.707107 -2 -0.5) (0 -2 -0.5) (0 -1 -0.5) (0 -0.5 -0.5) (-2 -0.707107 -0.5) (-0.707107 -0.707107 -0.5) (-0.353553 -0.353553 -0.5) (-2 -2 -0.5) (-0.707107 -2 -0.5) (10 -0.707107 0.5) (0.707107 -0.707107 0.5) (0.353553 -0.353553 0.5) (10 -2 0.5) (0.707107 -2 0.5) (0 -2 0.5) (0 -1 0.5) (0 -0.5 0.5) (-2 -0.707107 0.5) (-0.707107 -0.707107 0.5) (-0.353553 -0.353553 0.5) (-2 -2 0.5) (-0.707107 -2 0.5) ); blocks ( hex hex hex hex hex hex hex hex hex hex (5 4 9 10 24 23 28 29) (80 20 1) simpleGrading (10 1 1) (0 1 4 5 19 20 23 24) (80 20 1) simpleGrading (10 1 1) (1 2 3 4 20 21 22 23) (200 20 1) simpleGrading (1 1 1) (4 3 6 7 23 22 25 26) (200 40 1) simpleGrading (1 1 1) (9 4 7 8 28 23 26 27) (20 40 1) simpleGrading (1 1 1) (16 10 9 15 35 29 28 34) (20 80 1) simpleGrading (1 10 1) (11 16 15 12 30 35 34 31) (20 80 1) simpleGrading (1 10 1) (12 15 14 13 31 34 33 32) (20 20 1) simpleGrading (1 1 1) (15 18 17 14 34 37 36 33) (40 20 1) simpleGrading (1 1 1) (9 8 18 15 28 27 37 34) (40 20 1) simpleGrading (1 1 1) hex hex hex hex hex hex hex hex hex hex (40 45 44 39 53 58 57 52) (20 80 1) simpleGrading (1 10 1) (0 40 39 1 19 53 52 20) (20 80 1) simpleGrading (1 10 1) (1 39 38 2 20 52 51 21) (20 200 1) simpleGrading (1 1 1) (39 42 41 38 52 55 54 51) (40 200 1) simpleGrading (1 1 1) (44 43 42 39 57 56 55 52) (40 20 1) simpleGrading (1 1 1) (48 47 44 45 61 60 57 58) (80 20 1) simpleGrading (10 1 1) (11 12 47 48 30 31 60 61) (80 20 1) simpleGrading (10 1 1) (12 13 46 47 31 32 59 60) (20 20 1) simpleGrading (1 1 1) (47 46 49 50 60 59 62 63) (20 40 1) simpleGrading (1 1 1) (44 47 50 43 57 60 63 56) (20 40 1) simpleGrading (1 1 1) ); edges ( arc 0 5 (0.469846 0.17101 -0.5) arc 5 10 (0.17101 0.469846 -0.5) Page 2 of 5 File: /home/kmarc/Documents/ME702/OpenfoamScripts/blockMeshDict arc arc arc arc arc arc arc arc arc arc arc arc arc arc 1 4 (0.939693 0.34202 -0.5) 4 9 (0.34202 0.939693 -0.5) 19 24 (0.469846 0.17101 0.5) 24 29 (0.17101 0.469846 0.5) 20 23 (0.939693 0.34202 0.5) 23 28 (0.34202 0.939693 0.5) 11 16 (-0.469846 0.17101 -0.5) 16 10 (-0.17101 0.469846 -0.5) 12 15 (-0.939693 0.34202 -0.5) 15 9 (-0.34202 0.939693 -0.5) 30 35 (-0.469846 0.17101 0.5) 35 29 (-0.17101 0.469846 0.5) 31 34 (-0.939693 0.34202 0.5) 34 28 (-0.34202 0.939693 0.5) arc arc arc arc arc arc arc arc arc arc arc arc arc arc arc arc 0 40 (0.469846 -0.17101 -0.5) 40 45 (0.17101 -0.469846 -0.5) 1 39 (0.939693 -0.34202 -0.5) 39 44 (0.34202 -0.939693 -0.5) 19 53 (0.469846 -0.17101 0.5) 53 58 (0.17101 -0.469846 0.5) 20 52 (0.939693 -0.34202 0.5) 52 57 (0.34202 -0.939693 0.5) 11 48 (-0.469846 -0.17101 -0.5) 48 45 (-0.17101 -0.469846 -0.5) 12 47 (-0.939693 -0.34202 -0.5) 47 44 (-0.34202 -0.939693 -0.5) 30 61 (-0.469846 -0.17101 0.5) 61 58 (-0.17101 -0.469846 0.5) 31 60 (-0.939693 -0.34202 0.5) 60 57 (-0.34202 -0.939693 0.5) ); boundary ( top { type symmetryPlane; faces ( (7 8 27 26) (6 7 26 25) (8 18 37 27) (18 17 36 37) ); } bottom { type symmetryPlane; faces ( (49 50 63 62) (50 43 56 63) (43 42 55 56) (42 41 54 55) ); } inlet Page 3 of 5 File: /home/kmarc/Documents/ME702/OpenfoamScripts/blockMeshDict { type patch; faces ( (14 13 32 33) (17 14 33 36) (46 13 32 59) (46 49 62 59) ); } outlet { type patch; faces ( (2 3 22 21) (3 6 25 22) (38 51 21 2) (41 54 51 38) ); } cylinder { type wall; faces ( (10 5 24 29) (5 0 19 24) (16 10 29 35) (11 16 35 30) (48 11 30 61) (45 48 61 58) (40 45 58 53) (0 40 53 19) ); } frontAndBack { type empty; faces ( (5 10 9 4) (24 23 28 29) (0 5 4 1) (19 20 23 24) (1 4 3 2) (20 21 22 23) (4 7 6 3) (23 22 25 26) (4 9 8 7) (28 23 26 27) (16 15 9 10) (35 29 28 34) (12 15 16 11) (31 30 35 34) (13 14 15 12) Page 4 of 5 File: /home/kmarc/Documents/ME702/OpenfoamScripts/blockMeshDict (32 (14 (33 (15 (34 31 17 34 18 28 (45 (58 (40 (53 (39 (52 (39 (52 (44 (57 (47 (60 (12 (31 (13 (32 (49 (62 (50 (63 40 39 44) 57 52 53) 0 1 39) 52 20 19) 1 2 38) 51 21 20) 38 41 42) 55 54 51) 39 42 43) 56 55 52) 48 45 44) 57 58 61) 11 48 47) 60 61 30) 12 47 46) 59 60 31) 46 47 50) 63 60 59) 47 44 43) 56 57 60) Page 5 of 5 34 33) 18 15) 37 36) 8 9) 27 37) ); } ); mergePatchPairs ( ); // ************************************************************************* // B. U File: /home/kmarc/Documents/ME702/OpenfoamScripts/U Page 1 of 1 /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.4.0 | | \\ / A nd | Web: 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 (5 0 0); boundaryField { top { type } symmetryPlane; bottom { type } symmetryPlane; inlet { type freestreamValue } freestream; uniform (5 0 0); outlet { type freestreamValue } freestream; uniform (5 0 0); cylinder { type value } frontAndBack { type } fixedValue; uniform (0 0 0); empty; } // ************************************************************************* // C. p File: /home/kmarc/Documents/ME702/OpenfoamScripts/p Page 1 of 1 /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.4.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class volScalarField; object p; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -2 0 0 0 0]; internalField uniform 0; boundaryField { top { type } bottom { type } inlet { type } outlet { type } frontAndBack { type } cylinder { type } symmetryPlane; symmetryPlane; freestreamPressure; freestreamPressure; empty; zeroGradient; } // ************************************************************************* // D. transportProperties File: /home/kmarc/Documents/ME702/O…oamScripts/transportProperties Page 1 of 1 /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.4.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "constant"; object transportProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // transportModel Newtonian; nu nu [0 2 -1 0 0 0 0] 0.003106; // ************************************************************************* // E. RASProperties File: /home/kmarc/Documents/ME702/OpenfoamScripts/RASProperties Page 1 of 1 /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.4.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "constant"; object RASProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // RASModel laminar; turbulence off; printCoeffs off; // ************************************************************************* // F. controlDict File: /home/kmarc/Documents/ME702/OpenfoamScripts/controlDict Page 1 of 1 /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.4.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object controlDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // application icoFoam; startFrom startTime; startTime 0; stopAt endTime; endTime 1.75; deltaT 0.00001; writeControl timeStep; writeInterval 1000; purgeWrite 0; writeFormat ascii; writePrecision 6; writeCompression off; timeFormat general; timePrecision 6; runTimeModifiable true; // ************************************************************************* // G. fvSchemes File: /home/kmarc/Documents/ME702/OpenfoamScripts/fvSchemes Page 1 of 1 /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.4.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object fvSchemes; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // ddtSchemes { default } gradSchemes { default grad(p) } divSchemes { default div(phi,U) } Euler; Gauss linear; Gauss linear; none; Gauss linear; laplacianSchemes { default none; laplacian(nu,U) Gauss linear orthogonal; laplacian((1|A(U)),p) Gauss linear orthogonal; } interpolationSchemes { default interpolate(HbyA) } snGradSchemes { default } fluxRequired { default p } linear; linear; orthogonal; no; ; // ************************************************************************* // H. fvSolution File: /home/kmarc/Documents/ME702/OpenfoamScripts/fvSolution Page 1 of 1 /*--------------------------------*- C++ -*----------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.4.0 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ FoamFile { version 2.0; format ascii; class dictionary; location "system"; object fvSolution; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // solvers { p { solver preconditioner tolerance relTol PCG; DIC; 1e-06; 0; solver preconditioner tolerance relTol PBiCG; DILU; 1e-05; 0; } U { } } PISO { nCorrectors 2; nNonOrthogonalCorrectors 3; pRefCell 0; pRefValue 0; } // ************************************************************************* //