├── Benchmark ├── casebase │ ├── 0 │ │ ├── JxB │ │ ├── p_rgh │ │ ├── nut │ │ ├── U │ │ ├── alpha.melt │ │ ├── alpha.melt.orig │ │ ├── k │ │ ├── epsilon │ │ └── omega │ ├── meshElmer │ │ ├── mesh.header │ │ ├── mesh.names │ │ └── mesh.boundary │ ├── runTravisTest.sh │ ├── constant │ │ ├── polyMesh │ │ │ └── boundary │ │ ├── dynamicMeshDict │ │ ├── g │ │ ├── turbulenceProperties │ │ └── transportProperties │ ├── system │ │ ├── mapFieldsDict │ │ ├── decomposeParDict │ │ ├── extrudeMeshDict │ │ ├── setFieldsDict │ │ ├── controlDict │ │ ├── blockMeshDict │ │ ├── fvSchemes │ │ ├── fvSolution │ │ └── snappyHexMeshDict │ └── case.sif ├── templatedir │ ├── test.sh │ ├── sim.template │ └── case.template ├── README.md └── Dakota files │ ├── EOF_LPS.sh │ └── EOF_LPS.in ├── Examples ├── UQ │ ├── casebase │ │ ├── 0 │ │ │ ├── JxB │ │ │ ├── p_rgh │ │ │ ├── nut │ │ │ ├── U │ │ │ ├── alpha.melt │ │ │ ├── alpha.melt.orig │ │ │ ├── k │ │ │ ├── omega │ │ │ └── epsilon │ │ ├── meshElmer │ │ │ ├── mesh.header │ │ │ ├── mesh.names │ │ │ └── mesh.boundary │ │ ├── runTravisTest.sh │ │ ├── constant │ │ │ ├── polyMesh │ │ │ │ └── boundary │ │ │ ├── dynamicMeshDict │ │ │ ├── g │ │ │ ├── turbulenceProperties │ │ │ └── transportProperties │ │ ├── system │ │ │ ├── mapFieldsDict │ │ │ ├── decomposeParDict │ │ │ ├── extrudeMeshDict │ │ │ ├── setFieldsDict │ │ │ ├── controlDict │ │ │ ├── blockMeshDict │ │ │ ├── fvSchemes │ │ │ ├── fvSolution │ │ │ └── snappyHexMeshDict │ │ └── case.sif │ ├── README.md │ ├── templatedir │ │ ├── test.sh │ │ ├── sim.template │ │ └── case.template │ └── Dakota files │ │ ├── EOF_LHS.in │ │ └── EOFD_driver.sh └── Surrogate Modeling │ ├── casebase │ ├── 0 │ │ ├── JxB │ │ ├── p_rgh │ │ ├── nut │ │ ├── U │ │ ├── alpha.melt │ │ ├── alpha.melt.orig │ │ ├── k │ │ ├── epsilon │ │ └── omega │ ├── meshElmer │ │ ├── mesh.header │ │ ├── mesh.names │ │ └── mesh.boundary │ ├── runTravisTest.sh │ ├── constant │ │ ├── polyMesh │ │ │ └── boundary │ │ ├── dynamicMeshDict │ │ ├── g │ │ ├── turbulenceProperties │ │ └── transportProperties │ ├── system │ │ ├── mapFieldsDict │ │ ├── decomposeParDict │ │ ├── extrudeMeshDict │ │ ├── setFieldsDict │ │ ├── controlDict │ │ ├── blockMeshDict │ │ ├── fvSchemes │ │ ├── fvSolution │ │ └── snappyHexMeshDict │ └── case.sif │ ├── README.md │ ├── templatedir │ ├── test.sh │ ├── sim.template │ └── case.template │ └── Dakota files │ ├── EOFD_driver.sh │ └── EOF_LHS_SM.in ├── README.md └── src ├── test.sh └── Interface_program.sh /Benchmark/casebase/meshElmer/mesh.header: -------------------------------------------------------------------------------- 1 | 2867 5628 515 2 | 2 3 | 202 515 4 | 303 5628 5 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/meshElmer/mesh.header: -------------------------------------------------------------------------------- 1 | 2867 5628 515 2 | 2 3 | 202 515 4 | 303 5628 5 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/meshElmer/mesh.header: -------------------------------------------------------------------------------- 1 | 2867 5628 515 2 | 2 3 | 202 515 4 | 303 5628 5 | -------------------------------------------------------------------------------- /Benchmark/casebase/meshElmer/mesh.names: -------------------------------------------------------------------------------- 1 | ! ----- names for bodies ----- 2 | $ melt = 1 3 | $ air = 2 4 | $ indTop = 3 5 | $ indBtm = 4 6 | ! ----- names for boundaries ----- 7 | $ out = 5 8 | $ empty = 6 9 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/meshElmer/mesh.names: -------------------------------------------------------------------------------- 1 | ! ----- names for bodies ----- 2 | $ melt = 1 3 | $ air = 2 4 | $ indTop = 3 5 | $ indBtm = 4 6 | ! ----- names for boundaries ----- 7 | $ out = 5 8 | $ empty = 6 9 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/meshElmer/mesh.names: -------------------------------------------------------------------------------- 1 | ! ----- names for bodies ----- 2 | $ melt = 1 3 | $ air = 2 4 | $ indTop = 3 5 | $ indBtm = 4 6 | ! ----- names for boundaries ----- 7 | $ out = 5 8 | $ empty = 6 9 | -------------------------------------------------------------------------------- /Benchmark/casebase/runTravisTest.sh: -------------------------------------------------------------------------------- 1 | #!/bin/bash 2 | cd "$( cd "$( dirname "${BASH_SOURCE[0]}" )" && pwd )" 3 | 4 | nFails=0 5 | 6 | wmake $EOF_HOME/solvers/mhdInterFoam6 7 | 8 | ElmerGrid 2 2 meshElmer -metis 2 9 | setFields 10 | decomposePar 11 | 12 | for case in *.sif; do 13 | rm TEST.PASSED* 14 | mpirun -n 2 mhdInterFoam -parallel : -n 2 ElmerSolver_mpi $case 15 | if [[ $(< TEST.PASSED*) != "1" ]]; then 16 | nFails=$((nFails+1)) 17 | fi 18 | done 19 | 20 | exit $nFails 21 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/runTravisTest.sh: -------------------------------------------------------------------------------- 1 | #!/bin/bash 2 | cd "$( cd "$( dirname "${BASH_SOURCE[0]}" )" && pwd )" 3 | 4 | nFails=0 5 | 6 | wmake $EOF_HOME/solvers/mhdInterFoam6 7 | 8 | ElmerGrid 2 2 meshElmer -metis 2 9 | setFields 10 | decomposePar 11 | 12 | for case in *.sif; do 13 | rm TEST.PASSED* 14 | mpirun -n 2 mhdInterFoam -parallel : -n 2 ElmerSolver_mpi $case 15 | if [[ $(< TEST.PASSED*) != "1" ]]; then 16 | nFails=$((nFails+1)) 17 | fi 18 | done 19 | 20 | exit $nFails 21 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/README.md: -------------------------------------------------------------------------------- 1 | The UQ study was extended to surrogate modeling. 2 | Design and analysis of computer experiments (DACE) technique was employed to generate the training data, a sample size of 20 was chosen.  3 | A linear regression fit was applied to the training data to construct the global surrogate model.  4 | Latin hypercube sampling study, with a sample size of 100, was then performed on the global surrogate model.  5 | The DACE technique was shown to reduce the total run time by 45%. 6 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/runTravisTest.sh: -------------------------------------------------------------------------------- 1 | #!/bin/bash 2 | cd "$( cd "$( dirname "${BASH_SOURCE[0]}" )" && pwd )" 3 | 4 | nFails=0 5 | 6 | wmake $EOF_HOME/solvers/mhdInterFoam6 7 | 8 | ElmerGrid 2 2 meshElmer -metis 2 9 | setFields 10 | decomposePar 11 | 12 | for case in *.sif; do 13 | rm TEST.PASSED* 14 | mpirun -n 2 mhdInterFoam -parallel : -n 2 ElmerSolver_mpi $case 15 | if [[ $(< TEST.PASSED*) != "1" ]]; then 16 | nFails=$((nFails+1)) 17 | fi 18 | done 19 | 20 | exit $nFails 21 | -------------------------------------------------------------------------------- /Benchmark/casebase/constant/polyMesh/boundary: -------------------------------------------------------------------------------- 1 | /*--------------------------------------------------------------------*\ 2 | | | 3 | | File exported from Salome Platform using SalomeToFoamExporter | 4 | | | 5 | \*--------------------------------------------------------------------*/ 6 | FoamFile 7 | { 8 | version 2.0; 9 | format ascii; 10 | class polyBoundaryMesh; 11 | location "constant/polyMesh"; 12 | object boundary; 13 | } 14 | 15 | 3 16 | ( 17 | right 18 | { 19 | type wedge; 20 | nFaces 1925; 21 | startFace 3760; 22 | } 23 | out 24 | { 25 | type patch; 26 | nFaces 125; 27 | startFace 5685; 28 | } 29 | left 30 | { 31 | type wedge; 32 | nFaces 1925; 33 | startFace 5810; 34 | } 35 | ) 36 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/constant/polyMesh/boundary: -------------------------------------------------------------------------------- 1 | /*--------------------------------------------------------------------*\ 2 | | | 3 | | File exported from Salome Platform using SalomeToFoamExporter | 4 | | | 5 | \*--------------------------------------------------------------------*/ 6 | FoamFile 7 | { 8 | version 2.0; 9 | format ascii; 10 | class polyBoundaryMesh; 11 | location "constant/polyMesh"; 12 | object boundary; 13 | } 14 | 15 | 3 16 | ( 17 | right 18 | { 19 | type wedge; 20 | nFaces 1925; 21 | startFace 3760; 22 | } 23 | out 24 | { 25 | type patch; 26 | nFaces 125; 27 | startFace 5685; 28 | } 29 | left 30 | { 31 | type wedge; 32 | nFaces 1925; 33 | startFace 5810; 34 | } 35 | ) 36 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/constant/polyMesh/boundary: -------------------------------------------------------------------------------- 1 | /*--------------------------------------------------------------------*\ 2 | | | 3 | | File exported from Salome Platform using SalomeToFoamExporter | 4 | | | 5 | \*--------------------------------------------------------------------*/ 6 | FoamFile 7 | { 8 | version 2.0; 9 | format ascii; 10 | class polyBoundaryMesh; 11 | location "constant/polyMesh"; 12 | object boundary; 13 | } 14 | 15 | 3 16 | ( 17 | right 18 | { 19 | type wedge; 20 | nFaces 1925; 21 | startFace 3760; 22 | } 23 | out 24 | { 25 | type patch; 26 | nFaces 125; 27 | startFace 5685; 28 | } 29 | left 30 | { 31 | type wedge; 32 | nFaces 1925; 33 | startFace 5810; 34 | } 35 | ) 36 | -------------------------------------------------------------------------------- /Benchmark/casebase/constant/dynamicMeshDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "constant"; 14 | object dynamicMeshDict; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dynamicFvMesh staticFvMesh; 19 | 20 | 21 | // ************************************************************************* // 22 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/constant/dynamicMeshDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "constant"; 14 | object dynamicMeshDict; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dynamicFvMesh staticFvMesh; 19 | 20 | 21 | // ************************************************************************* // 22 | -------------------------------------------------------------------------------- /Benchmark/casebase/system/mapFieldsDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "system"; 14 | object mapFieldsDict; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | patchMap ( ); 19 | 20 | cuttingPatches ( ); 21 | 22 | 23 | // ************************************************************************* // 24 | -------------------------------------------------------------------------------- /Benchmark/casebase/system/decomposeParDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "system"; 14 | object decomposeParDict; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | numberOfSubdomains 2; 19 | 20 | method scotch; 21 | 22 | // ************************************************************************* // 23 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/constant/dynamicMeshDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "constant"; 14 | object dynamicMeshDict; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dynamicFvMesh staticFvMesh; 19 | 20 | 21 | // ************************************************************************* // 22 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/system/mapFieldsDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "system"; 14 | object mapFieldsDict; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | patchMap ( ); 19 | 20 | cuttingPatches ( ); 21 | 22 | 23 | // ************************************************************************* // 24 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/system/decomposeParDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "system"; 14 | object decomposeParDict; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | numberOfSubdomains 2; 19 | 20 | method scotch; 21 | 22 | // ************************************************************************* // 23 | -------------------------------------------------------------------------------- /Benchmark/casebase/constant/g: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class uniformDimensionedVectorField; 13 | location "constant"; 14 | object g; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dimensions [0 1 -2 0 0 0 0]; 19 | value (0 -9.81 0); 20 | 21 | 22 | // ************************************************************************* // 23 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/system/mapFieldsDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "system"; 14 | object mapFieldsDict; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | patchMap ( ); 19 | 20 | cuttingPatches ( ); 21 | 22 | 23 | // ************************************************************************* // 24 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/constant/g: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class uniformDimensionedVectorField; 13 | location "constant"; 14 | object g; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dimensions [0 1 -2 0 0 0 0]; 19 | value (0 -9.81 0); 20 | 21 | 22 | // ************************************************************************* // 23 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/system/decomposeParDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "system"; 14 | object decomposeParDict; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | numberOfSubdomains 2; 19 | 20 | method scotch; 21 | 22 | // ************************************************************************* // 23 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/constant/g: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class uniformDimensionedVectorField; 13 | location "constant"; 14 | object g; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dimensions [0 1 -2 0 0 0 0]; 19 | value (0 -9.81 0); 20 | 21 | 22 | // ************************************************************************* // 23 | -------------------------------------------------------------------------------- /Benchmark/casebase/constant/turbulenceProperties: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "constant"; 14 | object turbulenceProperties; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | simulationType RAS; 19 | 20 | RAS 21 | { 22 | RASModel kOmegaSST; 23 | 24 | turbulence on; 25 | 26 | printCoeffs on; 27 | } 28 | 29 | 30 | // ************************************************************************* // 31 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/constant/turbulenceProperties: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "constant"; 14 | object turbulenceProperties; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | simulationType RAS; 19 | 20 | RAS 21 | { 22 | RASModel kOmegaSST; 23 | 24 | turbulence on; 25 | 26 | printCoeffs on; 27 | } 28 | 29 | 30 | // ************************************************************************* // 31 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/constant/turbulenceProperties: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "constant"; 14 | object turbulenceProperties; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | simulationType RAS; 19 | 20 | RAS 21 | { 22 | RASModel kOmegaSST; 23 | 24 | turbulence on; 25 | 26 | printCoeffs on; 27 | } 28 | 29 | 30 | // ************************************************************************* // 31 | -------------------------------------------------------------------------------- /Benchmark/casebase/0/JxB: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volVectorField; 13 | location "0"; 14 | object JxB; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dimensions [1 -2 -2 0 0 0 0]; 19 | 20 | internalField uniform (0 0 0); 21 | 22 | boundaryField 23 | { 24 | left 25 | { 26 | type wedge; 27 | } 28 | right 29 | { 30 | type wedge; 31 | } 32 | ".*" 33 | { 34 | type zeroGradient; 35 | } 36 | } 37 | 38 | // ************************************************************************* // 39 | -------------------------------------------------------------------------------- /Benchmark/casebase/system/extrudeMeshDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | object extrudeProperties; 14 | } 15 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 16 | 17 | constructFrom patch; 18 | sourceCase "$FOAM_CASE"; 19 | 20 | sourcePatches (front); 21 | exposedPatchName back; 22 | 23 | extrudeModel wedge; 24 | 25 | sectorCoeffs 26 | { 27 | axisPt (0 0 0); 28 | axis (1 0 0); 29 | angle 1; 30 | } 31 | 32 | flipNormals false; 33 | mergeFaces false; 34 | 35 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 36 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/0/JxB: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volVectorField; 13 | location "0"; 14 | object JxB; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dimensions [1 -2 -2 0 0 0 0]; 19 | 20 | internalField uniform (0 0 0); 21 | 22 | boundaryField 23 | { 24 | left 25 | { 26 | type wedge; 27 | } 28 | right 29 | { 30 | type wedge; 31 | } 32 | ".*" 33 | { 34 | type zeroGradient; 35 | } 36 | } 37 | 38 | // ************************************************************************* // 39 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/system/extrudeMeshDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | object extrudeProperties; 14 | } 15 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 16 | 17 | constructFrom patch; 18 | sourceCase "$FOAM_CASE"; 19 | 20 | sourcePatches (front); 21 | exposedPatchName back; 22 | 23 | extrudeModel wedge; 24 | 25 | sectorCoeffs 26 | { 27 | axisPt (0 0 0); 28 | axis (1 0 0); 29 | angle 1; 30 | } 31 | 32 | flipNormals false; 33 | mergeFaces false; 34 | 35 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 36 | -------------------------------------------------------------------------------- /Benchmark/casebase/0/p_rgh: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volScalarField; 13 | object p_rgh; 14 | } 15 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 16 | 17 | dimensions [1 -1 -2 0 0 0 0]; 18 | 19 | internalField uniform 0; 20 | 21 | boundaryField 22 | { 23 | 24 | left 25 | { 26 | type wedge; 27 | } 28 | right 29 | { 30 | type wedge; 31 | } 32 | ".*" 33 | { 34 | type totalPressure; 35 | p0 uniform 0; 36 | } 37 | } 38 | 39 | // ************************************************************************* // 40 | -------------------------------------------------------------------------------- /Benchmark/casebase/system/setFieldsDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "system"; 14 | object setFieldsDict; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | defaultFieldValues 19 | ( 20 | volScalarFieldValue alpha.melt 0 21 | ); 22 | 23 | regions 24 | ( 25 | sphereToCell 26 | { 27 | centre (0.0 0.005 0.0); 28 | radius 0.0128; 29 | fieldValues 30 | ( 31 | volScalarFieldValue alpha.melt 1 32 | ); 33 | } 34 | ); 35 | 36 | 37 | // ************************************************************************* // 38 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/0/JxB: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volVectorField; 13 | location "0"; 14 | object JxB; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dimensions [1 -2 -2 0 0 0 0]; 19 | 20 | internalField uniform (0 0 0); 21 | 22 | boundaryField 23 | { 24 | left 25 | { 26 | type wedge; 27 | } 28 | right 29 | { 30 | type wedge; 31 | } 32 | ".*" 33 | { 34 | type zeroGradient; 35 | } 36 | } 37 | 38 | // ************************************************************************* // 39 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/system/extrudeMeshDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | object extrudeProperties; 14 | } 15 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 16 | 17 | constructFrom patch; 18 | sourceCase "$FOAM_CASE"; 19 | 20 | sourcePatches (front); 21 | exposedPatchName back; 22 | 23 | extrudeModel wedge; 24 | 25 | sectorCoeffs 26 | { 27 | axisPt (0 0 0); 28 | axis (1 0 0); 29 | angle 1; 30 | } 31 | 32 | flipNormals false; 33 | mergeFaces false; 34 | 35 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 36 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/system/setFieldsDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "system"; 14 | object setFieldsDict; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | defaultFieldValues 19 | ( 20 | volScalarFieldValue alpha.melt 0 21 | ); 22 | 23 | regions 24 | ( 25 | sphereToCell 26 | { 27 | centre (0.0 0.005 0.0); 28 | radius 0.0128; 29 | fieldValues 30 | ( 31 | volScalarFieldValue alpha.melt 1 32 | ); 33 | } 34 | ); 35 | 36 | 37 | // ************************************************************************* // 38 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/0/p_rgh: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volScalarField; 13 | object p_rgh; 14 | } 15 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 16 | 17 | dimensions [1 -1 -2 0 0 0 0]; 18 | 19 | internalField uniform 0; 20 | 21 | boundaryField 22 | { 23 | 24 | left 25 | { 26 | type wedge; 27 | } 28 | right 29 | { 30 | type wedge; 31 | } 32 | ".*" 33 | { 34 | type totalPressure; 35 | p0 uniform 0; 36 | } 37 | } 38 | 39 | // ************************************************************************* // 40 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/system/setFieldsDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "system"; 14 | object setFieldsDict; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | defaultFieldValues 19 | ( 20 | volScalarFieldValue alpha.melt 0 21 | ); 22 | 23 | regions 24 | ( 25 | sphereToCell 26 | { 27 | centre (0.0 0.005 0.0); 28 | radius 0.0128; 29 | fieldValues 30 | ( 31 | volScalarFieldValue alpha.melt 1 32 | ); 33 | } 34 | ); 35 | 36 | 37 | // ************************************************************************* // 38 | -------------------------------------------------------------------------------- /Benchmark/casebase/0/nut: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volScalarField; 13 | location "0"; 14 | object nut; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dimensions [0 2 -1 0 0 0 0]; 19 | 20 | internalField uniform 0; 21 | 22 | boundaryField 23 | { 24 | left 25 | { 26 | type wedge; 27 | } 28 | right 29 | { 30 | type wedge; 31 | } 32 | ".*" 33 | { 34 | type calculated; 35 | value uniform 0; 36 | } 37 | } 38 | 39 | 40 | // ************************************************************************* // 41 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/0/p_rgh: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volScalarField; 13 | object p_rgh; 14 | } 15 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 16 | 17 | dimensions [1 -1 -2 0 0 0 0]; 18 | 19 | internalField uniform 0; 20 | 21 | boundaryField 22 | { 23 | 24 | left 25 | { 26 | type wedge; 27 | } 28 | right 29 | { 30 | type wedge; 31 | } 32 | ".*" 33 | { 34 | type totalPressure; 35 | p0 uniform 0; 36 | } 37 | } 38 | 39 | // ************************************************************************* // 40 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/0/nut: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volScalarField; 13 | location "0"; 14 | object nut; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dimensions [0 2 -1 0 0 0 0]; 19 | 20 | internalField uniform 0; 21 | 22 | boundaryField 23 | { 24 | left 25 | { 26 | type wedge; 27 | } 28 | right 29 | { 30 | type wedge; 31 | } 32 | ".*" 33 | { 34 | type calculated; 35 | value uniform 0; 36 | } 37 | } 38 | 39 | 40 | // ************************************************************************* // 41 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/0/nut: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volScalarField; 13 | location "0"; 14 | object nut; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dimensions [0 2 -1 0 0 0 0]; 19 | 20 | internalField uniform 0; 21 | 22 | boundaryField 23 | { 24 | left 25 | { 26 | type wedge; 27 | } 28 | right 29 | { 30 | type wedge; 31 | } 32 | ".*" 33 | { 34 | type calculated; 35 | value uniform 0; 36 | } 37 | } 38 | 39 | 40 | // ************************************************************************* // 41 | -------------------------------------------------------------------------------- /Benchmark/casebase/0/U: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volVectorField; 13 | location "0"; 14 | object U; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dimensions [0 1 -1 0 0 0 0]; 19 | 20 | internalField uniform (0 0 0); 21 | 22 | boundaryField 23 | { 24 | left 25 | { 26 | type wedge; 27 | } 28 | right 29 | { 30 | type wedge; 31 | } 32 | ".*" 33 | { 34 | type pressureInletOutletVelocity; 35 | value uniform (0 0 0); 36 | } 37 | } 38 | 39 | 40 | // ************************************************************************* // 41 | -------------------------------------------------------------------------------- /Benchmark/casebase/0/alpha.melt: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volScalarField; 13 | object alpha.melt; 14 | } 15 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 16 | 17 | dimensions [0 0 0 0 0 0 0]; 18 | 19 | internalField uniform 0; 20 | 21 | boundaryField 22 | { 23 | left 24 | { 25 | type wedge; 26 | } 27 | right 28 | { 29 | type wedge; 30 | } 31 | ".*" 32 | { 33 | type inletOutlet; 34 | inletValue uniform 0; 35 | value uniform 0; 36 | } 37 | } 38 | 39 | 40 | // ************************************************************************* // 41 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/0/U: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volVectorField; 13 | location "0"; 14 | object U; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dimensions [0 1 -1 0 0 0 0]; 19 | 20 | internalField uniform (0 0 0); 21 | 22 | boundaryField 23 | { 24 | left 25 | { 26 | type wedge; 27 | } 28 | right 29 | { 30 | type wedge; 31 | } 32 | ".*" 33 | { 34 | type pressureInletOutletVelocity; 35 | value uniform (0 0 0); 36 | } 37 | } 38 | 39 | 40 | // ************************************************************************* // 41 | -------------------------------------------------------------------------------- /Benchmark/casebase/0/alpha.melt.orig: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volScalarField; 13 | object alpha.melt; 14 | } 15 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 16 | 17 | dimensions [0 0 0 0 0 0 0]; 18 | 19 | internalField uniform 0; 20 | 21 | boundaryField 22 | { 23 | left 24 | { 25 | type wedge; 26 | } 27 | right 28 | { 29 | type wedge; 30 | } 31 | ".*" 32 | { 33 | type inletOutlet; 34 | inletValue uniform 0; 35 | value uniform 0; 36 | } 37 | } 38 | 39 | 40 | // ************************************************************************* // 41 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/0/alpha.melt: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volScalarField; 13 | object alpha.melt; 14 | } 15 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 16 | 17 | dimensions [0 0 0 0 0 0 0]; 18 | 19 | internalField uniform 0; 20 | 21 | boundaryField 22 | { 23 | left 24 | { 25 | type wedge; 26 | } 27 | right 28 | { 29 | type wedge; 30 | } 31 | ".*" 32 | { 33 | type inletOutlet; 34 | inletValue uniform 0; 35 | value uniform 0; 36 | } 37 | } 38 | 39 | 40 | // ************************************************************************* // 41 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/0/alpha.melt.orig: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volScalarField; 13 | object alpha.melt; 14 | } 15 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 16 | 17 | dimensions [0 0 0 0 0 0 0]; 18 | 19 | internalField uniform 0; 20 | 21 | boundaryField 22 | { 23 | left 24 | { 25 | type wedge; 26 | } 27 | right 28 | { 29 | type wedge; 30 | } 31 | ".*" 32 | { 33 | type inletOutlet; 34 | inletValue uniform 0; 35 | value uniform 0; 36 | } 37 | } 38 | 39 | 40 | // ************************************************************************* // 41 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/0/U: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volVectorField; 13 | location "0"; 14 | object U; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dimensions [0 1 -1 0 0 0 0]; 19 | 20 | internalField uniform (0 0 0); 21 | 22 | boundaryField 23 | { 24 | left 25 | { 26 | type wedge; 27 | } 28 | right 29 | { 30 | type wedge; 31 | } 32 | ".*" 33 | { 34 | type pressureInletOutletVelocity; 35 | value uniform (0 0 0); 36 | } 37 | } 38 | 39 | 40 | // ************************************************************************* // 41 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/0/alpha.melt: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volScalarField; 13 | object alpha.melt; 14 | } 15 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 16 | 17 | dimensions [0 0 0 0 0 0 0]; 18 | 19 | internalField uniform 0; 20 | 21 | boundaryField 22 | { 23 | left 24 | { 25 | type wedge; 26 | } 27 | right 28 | { 29 | type wedge; 30 | } 31 | ".*" 32 | { 33 | type inletOutlet; 34 | inletValue uniform 0; 35 | value uniform 0; 36 | } 37 | } 38 | 39 | 40 | // ************************************************************************* // 41 | -------------------------------------------------------------------------------- /Benchmark/casebase/0/k: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volScalarField; 13 | location "0"; 14 | object k; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dimensions [0 2 -2 0 0 0 0]; 19 | 20 | internalField uniform 1e-5; 21 | 22 | boundaryField 23 | { 24 | left 25 | { 26 | type wedge; 27 | } 28 | right 29 | { 30 | type wedge; 31 | } 32 | ".*" 33 | { 34 | type inletOutlet; 35 | inletValue $internalField; 36 | value $internalField; 37 | } 38 | } 39 | 40 | 41 | // ************************************************************************* // 42 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/0/alpha.melt.orig: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volScalarField; 13 | object alpha.melt; 14 | } 15 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 16 | 17 | dimensions [0 0 0 0 0 0 0]; 18 | 19 | internalField uniform 0; 20 | 21 | boundaryField 22 | { 23 | left 24 | { 25 | type wedge; 26 | } 27 | right 28 | { 29 | type wedge; 30 | } 31 | ".*" 32 | { 33 | type inletOutlet; 34 | inletValue uniform 0; 35 | value uniform 0; 36 | } 37 | } 38 | 39 | 40 | // ************************************************************************* // 41 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/0/k: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volScalarField; 13 | location "0"; 14 | object k; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dimensions [0 2 -2 0 0 0 0]; 19 | 20 | internalField uniform 1e-5; 21 | 22 | boundaryField 23 | { 24 | left 25 | { 26 | type wedge; 27 | } 28 | right 29 | { 30 | type wedge; 31 | } 32 | ".*" 33 | { 34 | type inletOutlet; 35 | inletValue $internalField; 36 | value $internalField; 37 | } 38 | } 39 | 40 | 41 | // ************************************************************************* // 42 | -------------------------------------------------------------------------------- /Benchmark/casebase/0/epsilon: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volScalarField; 13 | location "0"; 14 | object epsilon; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dimensions [0 2 -3 0 0 0 0]; 19 | 20 | internalField uniform 1e-6; 21 | 22 | boundaryField 23 | { 24 | left 25 | { 26 | type wedge; 27 | } 28 | right 29 | { 30 | type wedge; 31 | } 32 | ".*" 33 | { 34 | type inletOutlet; 35 | inletValue $internalField; 36 | value $internalField; 37 | } 38 | } 39 | 40 | 41 | // ************************************************************************* // 42 | -------------------------------------------------------------------------------- /Benchmark/casebase/0/omega: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volScalarField; 13 | location "0"; 14 | object omega; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dimensions [0 0 -1 0 0 0 0]; 19 | 20 | internalField uniform 2; 21 | 22 | boundaryField 23 | { 24 | left 25 | { 26 | type wedge; 27 | } 28 | right 29 | { 30 | type wedge; 31 | } 32 | ".*" 33 | { 34 | type inletOutlet; 35 | inletValue $internalField; 36 | value $internalField; 37 | } 38 | } 39 | 40 | 41 | 42 | // ************************************************************************* // 43 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/0/omega: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volScalarField; 13 | location "0"; 14 | object omega; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dimensions [0 0 -1 0 0 0 0]; 19 | 20 | internalField uniform 2; 21 | 22 | boundaryField 23 | { 24 | left 25 | { 26 | type wedge; 27 | } 28 | right 29 | { 30 | type wedge; 31 | } 32 | ".*" 33 | { 34 | type inletOutlet; 35 | inletValue $internalField; 36 | value $internalField; 37 | } 38 | } 39 | 40 | 41 | 42 | // ************************************************************************* // 43 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/0/k: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volScalarField; 13 | location "0"; 14 | object k; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dimensions [0 2 -2 0 0 0 0]; 19 | 20 | internalField uniform 1e-5; 21 | 22 | boundaryField 23 | { 24 | left 25 | { 26 | type wedge; 27 | } 28 | right 29 | { 30 | type wedge; 31 | } 32 | ".*" 33 | { 34 | type inletOutlet; 35 | inletValue $internalField; 36 | value $internalField; 37 | } 38 | } 39 | 40 | 41 | // ************************************************************************* // 42 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/0/epsilon: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volScalarField; 13 | location "0"; 14 | object epsilon; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dimensions [0 2 -3 0 0 0 0]; 19 | 20 | internalField uniform 1e-6; 21 | 22 | boundaryField 23 | { 24 | left 25 | { 26 | type wedge; 27 | } 28 | right 29 | { 30 | type wedge; 31 | } 32 | ".*" 33 | { 34 | type inletOutlet; 35 | inletValue $internalField; 36 | value $internalField; 37 | } 38 | } 39 | 40 | 41 | // ************************************************************************* // 42 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/0/epsilon: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volScalarField; 13 | location "0"; 14 | object epsilon; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dimensions [0 2 -3 0 0 0 0]; 19 | 20 | internalField uniform 1e-6; 21 | 22 | boundaryField 23 | { 24 | left 25 | { 26 | type wedge; 27 | } 28 | right 29 | { 30 | type wedge; 31 | } 32 | ".*" 33 | { 34 | type inletOutlet; 35 | inletValue $internalField; 36 | value $internalField; 37 | } 38 | } 39 | 40 | 41 | // ************************************************************************* // 42 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/0/omega: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class volScalarField; 13 | location "0"; 14 | object omega; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | dimensions [0 0 -1 0 0 0 0]; 19 | 20 | internalField uniform 2; 21 | 22 | boundaryField 23 | { 24 | left 25 | { 26 | type wedge; 27 | } 28 | right 29 | { 30 | type wedge; 31 | } 32 | ".*" 33 | { 34 | type inletOutlet; 35 | inletValue $internalField; 36 | value $internalField; 37 | } 38 | } 39 | 40 | 41 | 42 | // ************************************************************************* // 43 | -------------------------------------------------------------------------------- /README.md: -------------------------------------------------------------------------------- 1 | This program, Elmer FEM-OpenFoam-Dakota interface, is released under license GNU GENERAL PUBLIC LICENSE Version 3, 29 June 2007. 2 | 3 | A program for 3 way coupling between open-source software Elmer FEM, OpenFoam and Dakota has been developed. Furthermore, it was succesfully applied to benchmarks. The program along with examples have been provided. 4 | 5 | Capabilities: 6 | The interface developed provides an open-source framework that integrates data analytics, finite element and finite volume methods. This framework enables a richer understanding of multi-physics/phase predictions in a time sensitive manner. Its key features include uncertainty quantification, surrogate modeling, parallel programming and robust post-processing capabilities. 7 | 8 | If you use this program, please cite this repository -- Sandip, A. 2020, Elmer FEM-OpenFOAM-Dakota, v1.0, Zenodo, DOI: 10.5281/zenodo.3934329 9 | 10 | 11 | About: 12 | This program is mantained by Anjali Sandip. 13 | 14 | The software versions are listed below. 15 | The interface program can be extended to other versions of Elmer FEM, OpenFoam and Dakota with minor modifications. 16 | 17 | - Sandia Dakota version 6.10 18 | - Elmer version 8.4 19 | - OpenFoam version 6.0 20 | - EOF Library version 1.0 21 | - Operating system: Linux 22 | 23 | Support email for questions: anjali.sandip@und.edu 24 | 25 | 26 | 27 | 28 | 29 | 30 | 31 | 32 | -------------------------------------------------------------------------------- /Examples/UQ/README.md: -------------------------------------------------------------------------------- 1 | 2 | The benchmark case [1] was extended to a UQ study. 3 | 4 | - Problem objective: Determine the effect of melt density, surface tension, ac frequency and current density on the melt shape 5 | - Data analysis: Uncertainty quantification/Latin hypercube sampling study was performed, with a sample size of 100. 6 | The input parameters were distributed normally. 7 | Asynchronous local parallelism was implemented, 4 evaluations performed concurrently 8 | - Input parameter/s: melt density (2500 +/- 2000 kg/m3), surface tension (5.1 +/- 4.9 N/m), ac frequency (21250 +/- 18750 Hz), current density (5e7 +/- 3e7 A/m2) 9 | - Output parameter/s: melt shape (Lorentz force and flow velocity) at the final timestep 10 | - The turbulence model chosen for this study was k − ω SST. 11 | - More details on the test case (boundary conditions, etc.), can be found here [2] 12 | 13 | References: 14 | [1]. Spitans, Sergejs, Egbert Baake, Bernard Nacke, and Andris Jakovics. "Numerical modeling of free surface dynamics of melt in an alternate electromagnetic field. Part II: Conventional electromagnetic levitation." Metallurgical and Materials Transactions B 47, no. 1 (2016): 522-536. 15 | [2]. Vencels, J., A. Jakovics, and V. Geza. "Simulation of 3D MHD with free surface using open-source EOF-Library: levitating liquid metal in an alternating electromagnetic field." Magnetohydrodynamics 53, no. 4 (2017): 643-652. 16 | 17 | 18 | 19 | 20 | 21 | 22 | 23 | -------------------------------------------------------------------------------- /src/test.sh: -------------------------------------------------------------------------------- 1 | # copyright (c) 2020 ASandip 2 | #This program is free software: you can redistribute it and/or 3 | # modify it under the terms of the GNU General Public License as 4 | # published by the Free Software Foundation, either version 3 of the License, or 5 | # (at your option) any later version. 6 | # This program is distributed in the hope that it will be useful, 7 | # but WITHOUT ANY WARRANTY; without even the implied warranty of 8 | # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 9 | # GNU General Public License for more details. 10 | 11 | # You should have received a copy of the GNU General Public License 12 | # along with this program. If not, see . 13 | 14 | 15 | 16 | #!/bin/bash 17 | 18 | #Source EOF library environment 19 | . $HOME/EOF-Library/etc/bashrc 20 | 21 | #Source OpenFoam environment 22 | source ~/openfoam6/etc/bashrc 23 | 24 | 25 | #Compile EOF Library 26 | eofCompile >> log.out 2>&1 27 | 28 | #Compile the chosen solver for the test; for example, mhdInterFoam6 29 | wmake ../Solvers/mhdInterFoam6 >> log.out 2>&1 30 | 31 | 32 | #Copy Elmer and OpenFoam case files to Dakota Gui workspace 33 | cp -r ../casebase/* . 34 | 35 | #Prepare test 36 | setFields >> log.out 2>&1 37 | decomposePar >> log.out 2>&1 38 | ElmerGrid 2 2 meshElmer -metis 2 39 | 40 | 41 | #Run openFoam on n processes and Elmer on n processes, define n 42 | mpirun --allow-run-as-root -np n mhdInterFoam -parallel : -np n ElmerSolver_mpi case.sif >> log.out 2>&1 43 | -------------------------------------------------------------------------------- /Benchmark/casebase/system/controlDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "system"; 14 | object controlDict; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | application mhdInterFoam; 19 | 20 | startFrom startTime; 21 | 22 | startTime 0; 23 | 24 | stopAt endTime; 25 | 26 | endTime 0.5; 27 | 28 | deltaT 1e-6; 29 | 30 | writeControl adjustableRunTime; 31 | 32 | writeInterval 0.01; 33 | 34 | purgeWrite 0; 35 | 36 | writeFormat ascii; 37 | 38 | writePrecision 5; 39 | 40 | writeCompression uncompressed; 41 | 42 | timeFormat general; 43 | 44 | timePrecision 6; 45 | 46 | runTimeModifiable yes; 47 | 48 | adjustTimeStep yes; 49 | 50 | maxCo 0.8; 51 | maxAlphaCo 0.4; 52 | 53 | maxDeltaT 1; 54 | 55 | 56 | // ************************************************************************* // 57 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/system/controlDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "system"; 14 | object controlDict; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | application mhdInterFoam; 19 | 20 | startFrom startTime; 21 | 22 | startTime 0; 23 | 24 | stopAt endTime; 25 | 26 | endTime 0.5; 27 | 28 | deltaT 1e-6; 29 | 30 | writeControl adjustableRunTime; 31 | 32 | writeInterval 0.01; 33 | 34 | purgeWrite 0; 35 | 36 | writeFormat ascii; 37 | 38 | writePrecision 5; 39 | 40 | writeCompression uncompressed; 41 | 42 | timeFormat general; 43 | 44 | timePrecision 6; 45 | 46 | runTimeModifiable yes; 47 | 48 | adjustTimeStep yes; 49 | 50 | maxCo 0.8; 51 | maxAlphaCo 0.4; 52 | 53 | maxDeltaT 1; 54 | 55 | 56 | // ************************************************************************* // 57 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/system/controlDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "system"; 14 | object controlDict; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | application mhdInterFoam; 19 | 20 | startFrom startTime; 21 | 22 | startTime 0; 23 | 24 | stopAt endTime; 25 | 26 | endTime 0.5; 27 | 28 | deltaT 1e-6; 29 | 30 | writeControl adjustableRunTime; 31 | 32 | writeInterval 0.01; 33 | 34 | purgeWrite 0; 35 | 36 | writeFormat ascii; 37 | 38 | writePrecision 5; 39 | 40 | writeCompression uncompressed; 41 | 42 | timeFormat general; 43 | 44 | timePrecision 6; 45 | 46 | runTimeModifiable yes; 47 | 48 | adjustTimeStep yes; 49 | 50 | maxCo 0.8; 51 | maxAlphaCo 0.4; 52 | 53 | maxDeltaT 1; 54 | 55 | 56 | // ************************************************************************* // 57 | -------------------------------------------------------------------------------- /Benchmark/casebase/system/blockMeshDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | o | | 3 | | o o | HELYX-OS | 4 | | o O o | Version: v2.3.1 | 5 | | o o | Web: http://www.engys.com | 6 | | o | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location system; 14 | object blockMeshDict; 15 | } 16 | 17 | convertToMeters 1; 18 | vertices 19 | ( 20 | ( 0 -0.026 -0.001) 21 | ( 0.021 -0.026 -0.001) 22 | ( 0.021 0.026 -0.001) 23 | ( 0 0.026 -0.001) 24 | ( 0 -0.026 0.001) 25 | ( 0.021 -0.026 0.001) 26 | ( 0.021 0.026 0.001) 27 | ( 0 0.026 0.001) 28 | ); 29 | blocks 30 | ( hex 31 | ( 0 1 2 3 4 5 6 7) 32 | ( 63 156 1) simpleGrading 33 | ( 1 1 1) 34 | ); 35 | edges 36 | ( 37 | ); 38 | patches 39 | ( wall ffminx 40 | ( 41 | ( 0 4 7 3)) wall ffmaxx 42 | ( 43 | ( 1 2 6 5)) wall ffminy 44 | ( 45 | ( 0 1 5 4)) wall ffmaxy 46 | ( 47 | ( 3 7 6 2)) wall ffminz 48 | ( 49 | ( 0 3 2 1)) wall ffmaxz 50 | ( 51 | ( 4 5 6 7)) 52 | ); 53 | mergePatchPairs 54 | ( 55 | ); 56 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/system/blockMeshDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | o | | 3 | | o o | HELYX-OS | 4 | | o O o | Version: v2.3.1 | 5 | | o o | Web: http://www.engys.com | 6 | | o | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location system; 14 | object blockMeshDict; 15 | } 16 | 17 | convertToMeters 1; 18 | vertices 19 | ( 20 | ( 0 -0.026 -0.001) 21 | ( 0.021 -0.026 -0.001) 22 | ( 0.021 0.026 -0.001) 23 | ( 0 0.026 -0.001) 24 | ( 0 -0.026 0.001) 25 | ( 0.021 -0.026 0.001) 26 | ( 0.021 0.026 0.001) 27 | ( 0 0.026 0.001) 28 | ); 29 | blocks 30 | ( hex 31 | ( 0 1 2 3 4 5 6 7) 32 | ( 63 156 1) simpleGrading 33 | ( 1 1 1) 34 | ); 35 | edges 36 | ( 37 | ); 38 | patches 39 | ( wall ffminx 40 | ( 41 | ( 0 4 7 3)) wall ffmaxx 42 | ( 43 | ( 1 2 6 5)) wall ffminy 44 | ( 45 | ( 0 1 5 4)) wall ffmaxy 46 | ( 47 | ( 3 7 6 2)) wall ffminz 48 | ( 49 | ( 0 3 2 1)) wall ffmaxz 50 | ( 51 | ( 4 5 6 7)) 52 | ); 53 | mergePatchPairs 54 | ( 55 | ); 56 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/system/blockMeshDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | o | | 3 | | o o | HELYX-OS | 4 | | o O o | Version: v2.3.1 | 5 | | o o | Web: http://www.engys.com | 6 | | o | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location system; 14 | object blockMeshDict; 15 | } 16 | 17 | convertToMeters 1; 18 | vertices 19 | ( 20 | ( 0 -0.026 -0.001) 21 | ( 0.021 -0.026 -0.001) 22 | ( 0.021 0.026 -0.001) 23 | ( 0 0.026 -0.001) 24 | ( 0 -0.026 0.001) 25 | ( 0.021 -0.026 0.001) 26 | ( 0.021 0.026 0.001) 27 | ( 0 0.026 0.001) 28 | ); 29 | blocks 30 | ( hex 31 | ( 0 1 2 3 4 5 6 7) 32 | ( 63 156 1) simpleGrading 33 | ( 1 1 1) 34 | ); 35 | edges 36 | ( 37 | ); 38 | patches 39 | ( wall ffminx 40 | ( 41 | ( 0 4 7 3)) wall ffmaxx 42 | ( 43 | ( 1 2 6 5)) wall ffminy 44 | ( 45 | ( 0 1 5 4)) wall ffmaxy 46 | ( 47 | ( 3 7 6 2)) wall ffminz 48 | ( 49 | ( 0 3 2 1)) wall ffmaxz 50 | ( 51 | ( 4 5 6 7)) 52 | ); 53 | mergePatchPairs 54 | ( 55 | ); 56 | -------------------------------------------------------------------------------- /Benchmark/templatedir/test.sh: -------------------------------------------------------------------------------- 1 | # copyright (c) 2020 ASandip 2 | #This program is free software: you can redistribute it and/or 3 | # modify it under the terms of the GNU General Public License as 4 | # published by the Free Software Foundation, either version 3 of the License, or 5 | # (at your option) any later version. 6 | # This program is distributed in the hope that it will be useful, 7 | # but WITHOUT ANY WARRANTY; without even the implied warranty of 8 | # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 9 | # GNU General Public License for more details. 10 | 11 | # You should have received a copy of the GNU General Public License 12 | # along with this program. If not, see . 13 | 14 | 15 | 16 | #!/bin/bash 17 | 18 | #Source EOF library environment 19 | . $HOME/EOF-Library/etc/bashrc 20 | 21 | #Source OpenFoam environment 22 | source ~/openfoam6/etc/bashrc 23 | 24 | 25 | #Compile EOF Library 26 | eofCompile >> log.out 2>&1 27 | 28 | #Compile the chosen solver for the test, mhdInterFoam6 29 | wmake ../Solvers/mhdInterFoam6 >> log.out 2>&1 30 | 31 | 32 | #Copy Elmer and OpenFoam case files to Dakota Gui workspace 33 | cp -r ../casebase/* . 34 | 35 | #Transfer new set of values for the chosen input parameters to the case files 36 | cp transportProperties constant/transportProperties 37 | cp case_d.sif case.sif 38 | 39 | #Prepare test 40 | setFields >> log.out 2>&1 41 | decomposePar >> log.out 2>&1 42 | ElmerGrid 2 2 meshElmer -metis 2 43 | 44 | 45 | #Run openFoam on 2 processes and Elmer on 2 processes 46 | mpirun --allow-run-as-root -np 2 mhdInterFoam -parallel : -np 2 ElmerSolver_mpi case.sif >> log.out 2>&1 47 | -------------------------------------------------------------------------------- /Examples/UQ/templatedir/test.sh: -------------------------------------------------------------------------------- 1 | # copyright (c) 2020 ASandip 2 | #This program is free software: you can redistribute it and/or 3 | # modify it under the terms of the GNU General Public License as 4 | # published by the Free Software Foundation, either version 3 of the License, or 5 | # (at your option) any later version. 6 | # This program is distributed in the hope that it will be useful, 7 | # but WITHOUT ANY WARRANTY; without even the implied warranty of 8 | # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 9 | # GNU General Public License for more details. 10 | 11 | # You should have received a copy of the GNU General Public License 12 | # along with this program. If not, see . 13 | 14 | 15 | 16 | #!/bin/bash 17 | 18 | #Source EOF library environment 19 | . $HOME/EOF-Library/etc/bashrc 20 | 21 | #Source OpenFoam environment 22 | source ~/openfoam6/etc/bashrc 23 | 24 | 25 | #Compile EOF Library 26 | eofCompile >> log.out 2>&1 27 | 28 | #Compile the chosen solver for the test, mhdInterFoam6 29 | wmake ../Solvers/mhdInterFoam6 >> log.out 2>&1 30 | 31 | 32 | #Copy Elmer and OpenFoam case files to Dakota Gui workspace 33 | cp -r ../casebase/* . 34 | 35 | #Transfer new set of values for the chosen input parameters to the case files 36 | cp transportProperties constant/transportProperties 37 | cp case_d.sif case.sif 38 | 39 | #Prepare test 40 | setFields >> log.out 2>&1 41 | decomposePar >> log.out 2>&1 42 | ElmerGrid 2 2 meshElmer -metis 2 43 | 44 | 45 | #Run openFoam on 2 processes and Elmer on 2 processes 46 | mpirun --allow-run-as-root -np 2 mhdInterFoam -parallel : -np 1 ElmerSolver_mpi case.sif >> log.out 2>&1 47 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/templatedir/test.sh: -------------------------------------------------------------------------------- 1 | # copyright (c) 2020 ASandip 2 | #This program is free software: you can redistribute it and/or 3 | # modify it under the terms of the GNU General Public License as 4 | # published by the Free Software Foundation, either version 3 of the License, or 5 | # (at your option) any later version. 6 | # This program is distributed in the hope that it will be useful, 7 | # but WITHOUT ANY WARRANTY; without even the implied warranty of 8 | # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 9 | # GNU General Public License for more details. 10 | 11 | # You should have received a copy of the GNU General Public License 12 | # along with this program. If not, see . 13 | 14 | 15 | 16 | #!/bin/bash 17 | 18 | #Source EOF library environment 19 | . $HOME/EOF-Library/etc/bashrc 20 | 21 | #Source OpenFoam environment 22 | source ~/openfoam6/etc/bashrc 23 | 24 | 25 | #Compile EOF Library 26 | eofCompile >> log.out 2>&1 27 | 28 | #Compile the chosen solver for the test, mhdInterFoam6 29 | wmake ../Solvers/mhdInterFoam6 >> log.out 2>&1 30 | 31 | 32 | #Copy Elmer and OpenFoam case files to Dakota Gui workspace 33 | cp -r ../casebase/* . 34 | 35 | #Transfer new set of values for the chosen input parameters to the case files 36 | cp transportProperties constant/transportProperties 37 | cp case_d.sif case.sif 38 | 39 | #Prepare test 40 | setFields >> log.out 2>&1 41 | decomposePar >> log.out 2>&1 42 | ElmerGrid 2 2 meshElmer -metis 2 43 | 44 | 45 | #Run openFoam on 2 processes and Elmer on 2 processes 46 | mpirun --allow-run-as-root -np 2 mhdInterFoam -parallel : -np 1 ElmerSolver_mpi case.sif >> log.out 2>&1 47 | -------------------------------------------------------------------------------- /Benchmark/templatedir/sim.template: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "constant"; 14 | object transportProperties; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | phases (melt air); 19 | 20 | melt 21 | { 22 | transportModel Newtonian; 23 | nu [0 2 -1 0 0 0 0] 1e-4; 24 | rho [1 -3 0 0 0 0 0] {r}; 25 | } 26 | 27 | air 28 | { 29 | transportModel Newtonian; 30 | nu [0 2 -1 0 0 0 0] 1e-4; //-5 31 | rho [1 -3 0 0 0 0 0] 5; 32 | } 33 | 34 | sigma [ 1 0 -2 0 0 0 0] {s}; 35 | elcond_melt [-1 -3 3 0 0 2 0] 3850000; 36 | 37 | // Maximum allowable relative field difference in any cell comparing 38 | // to last Elmer update. This option controls frequency Elmer is called. 39 | // 1 - once 40 | // (0,1) - when relative difference in any cell exceeds given value 41 | // 0 - every iteration 42 | maxRelDiff 0.5; 43 | 44 | // ************************************************************************* // 45 | -------------------------------------------------------------------------------- /Examples/UQ/templatedir/sim.template: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "constant"; 14 | object transportProperties; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | phases (melt air); 19 | 20 | melt 21 | { 22 | transportModel Newtonian; 23 | nu [0 2 -1 0 0 0 0] 1e-4; 24 | rho [1 -3 0 0 0 0 0] {r}; 25 | } 26 | 27 | air 28 | { 29 | transportModel Newtonian; 30 | nu [0 2 -1 0 0 0 0] 1e-4; //-5 31 | rho [1 -3 0 0 0 0 0] 5; 32 | } 33 | 34 | sigma [ 1 0 -2 0 0 0 0] {s}; 35 | elcond_melt [-1 -3 3 0 0 2 0] 3850000; 36 | 37 | // Maximum allowable relative field difference in any cell comparing 38 | // to last Elmer update. This option controls frequency Elmer is called. 39 | // 1 - once 40 | // (0,1) - when relative difference in any cell exceeds given value 41 | // 0 - every iteration 42 | maxRelDiff 0.5; 43 | 44 | // ************************************************************************* // 45 | -------------------------------------------------------------------------------- /Benchmark/casebase/constant/transportProperties: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "constant"; 14 | object transportProperties; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | phases (melt air); 19 | 20 | melt 21 | { 22 | transportModel Newtonian; 23 | nu [0 2 -1 0 0 0 0] 1e-4; 24 | rho [1 -3 0 0 0 0 0] 1380; 25 | } 26 | 27 | air 28 | { 29 | transportModel Newtonian; 30 | nu [0 2 -1 0 0 0 0] 1e-4; //-5 31 | rho [1 -3 0 0 0 0 0] 5; 32 | } 33 | 34 | sigma [ 1 0 -2 0 0 0 0] 0.94; 35 | elcond_melt [-1 -3 3 0 0 2 0] 3850000; 36 | 37 | // Maximum allowable relative field difference in any cell comparing 38 | // to last Elmer update. This option controls frequency Elmer is called. 39 | // 1 - once 40 | // (0,1) - when relative difference in any cell exceeds given value 41 | // 0 - every iteration 42 | maxRelDiff 0.5; 43 | 44 | // ************************************************************************* // 45 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/templatedir/sim.template: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "constant"; 14 | object transportProperties; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | phases (melt air); 19 | 20 | melt 21 | { 22 | transportModel Newtonian; 23 | nu [0 2 -1 0 0 0 0] 1e-4; 24 | rho [1 -3 0 0 0 0 0] {r}; 25 | } 26 | 27 | air 28 | { 29 | transportModel Newtonian; 30 | nu [0 2 -1 0 0 0 0] 1e-4; //-5 31 | rho [1 -3 0 0 0 0 0] 5; 32 | } 33 | 34 | sigma [ 1 0 -2 0 0 0 0] {s}; 35 | elcond_melt [-1 -3 3 0 0 2 0] 3850000; 36 | 37 | // Maximum allowable relative field difference in any cell comparing 38 | // to last Elmer update. This option controls frequency Elmer is called. 39 | // 1 - once 40 | // (0,1) - when relative difference in any cell exceeds given value 41 | // 0 - every iteration 42 | maxRelDiff 0.5; 43 | 44 | // ************************************************************************* // 45 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/constant/transportProperties: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "constant"; 14 | object transportProperties; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | phases (melt air); 19 | 20 | melt 21 | { 22 | transportModel Newtonian; 23 | nu [0 2 -1 0 0 0 0] 1e-4; 24 | rho [1 -3 0 0 0 0 0] 1380; 25 | } 26 | 27 | air 28 | { 29 | transportModel Newtonian; 30 | nu [0 2 -1 0 0 0 0] 1e-4; //-5 31 | rho [1 -3 0 0 0 0 0] 5; 32 | } 33 | 34 | sigma [ 1 0 -2 0 0 0 0] 0.94; 35 | elcond_melt [-1 -3 3 0 0 2 0] 3850000; 36 | 37 | // Maximum allowable relative field difference in any cell comparing 38 | // to last Elmer update. This option controls frequency Elmer is called. 39 | // 1 - once 40 | // (0,1) - when relative difference in any cell exceeds given value 41 | // 0 - every iteration 42 | maxRelDiff 0.5; 43 | 44 | // ************************************************************************* // 45 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/constant/transportProperties: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "constant"; 14 | object transportProperties; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | phases (melt air); 19 | 20 | melt 21 | { 22 | transportModel Newtonian; 23 | nu [0 2 -1 0 0 0 0] 1e-4; 24 | rho [1 -3 0 0 0 0 0] 1380; 25 | } 26 | 27 | air 28 | { 29 | transportModel Newtonian; 30 | nu [0 2 -1 0 0 0 0] 1e-4; //-5 31 | rho [1 -3 0 0 0 0 0] 5; 32 | } 33 | 34 | sigma [ 1 0 -2 0 0 0 0] 0.94; 35 | elcond_melt [-1 -3 3 0 0 2 0] 3850000; 36 | 37 | // Maximum allowable relative field difference in any cell comparing 38 | // to last Elmer update. This option controls frequency Elmer is called. 39 | // 1 - once 40 | // (0,1) - when relative difference in any cell exceeds given value 41 | // 0 - every iteration 42 | maxRelDiff 0.5; 43 | 44 | // ************************************************************************* // 45 | -------------------------------------------------------------------------------- /Benchmark/README.md: -------------------------------------------------------------------------------- 1 | 2 | 3 | 4 | Free surface dynamics of melt in an alternated electromagnetic field (2D, Transient) 5 | 6 | The numerical predictions were compared against experimental data [1]. 7 | 8 | - Data analysis: Listed parameter study Combinations of values for input parameters used in the experimental study were chosen. https://dakota.sandia.gov/sites/default/files/docs/6.10/html-ref/method-list_parameter_study.html 9 | - Input parameter/s: melt density (kg/m3), surface tension (N/m), ac frequency (Hz), current density (A/m2) 10 | - Output parameter/s: Melt shape (Lorentz force and flow velocity) at the final timestep 11 | - The turbulence model chosen for this study was k − ω SST. 12 | - Elmer and OpenFoam files courtesy https://github.com/jvencels/EOF-Library 13 | - EOF Library [2] was employed to invoke communication between Elmer and OpenFoam during run time 14 | - More details on the test case (boundary conditions, etc.), can be found here [3] 15 | - The melt shapes obtained for the combinations of input parameter values were in good agreement with experimental data [1]. 16 | 17 | References: 18 | [1]. Spitans, Sergejs, Egbert Baake, Bernard Nacke, and Andris Jakovics. "Numerical modeling of free surface dynamics of melt in an alternate electromagnetic field. Part II: Conventional electromagnetic levitation." Metallurgical and Materials Transactions B 47, no. 1 (2016): 522-536. 19 | [2]. Vencels, Juris, Andris Jakovics, Vadims Geza, and Mihails Scepanskis. "EOF Library: Open-Source Elmer and OpenFOAM Coupler for simulation of MHD with free surface." arXiv preprint arXiv:1707.04080 (2017). 20 | [3]. Vencels, J., A. Jakovics, and V. Geza. "Simulation of 3D MHD with free surface using open-source EOF-Library: levitating liquid metal in an alternating electromagnetic field." Magnetohydrodynamics 53, no. 4 (2017): 643-652. 21 | -------------------------------------------------------------------------------- /Benchmark/casebase/system/fvSchemes: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "system"; 14 | object fvSchemes; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | ddtSchemes 19 | { 20 | default Euler; 21 | } 22 | 23 | gradSchemes 24 | { 25 | default Gauss linear; 26 | } 27 | 28 | divSchemes 29 | { 30 | div(rhoPhi,U) Gauss linearUpwindV grad(U); 31 | div(phi,alpha) Gauss vanLeer; 32 | div(phirb,alpha) Gauss linear; 33 | div(phi,k) Gauss upwind; 34 | div(phi,epsilon) Gauss upwind; 35 | div(phi,R) Gauss upwind; 36 | div(R) Gauss linear; 37 | div(phi,nuTilda) Gauss upwind; 38 | div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; 39 | "div\(phi,(k|omega)\)" Gauss upwind; 40 | } 41 | 42 | laplacianSchemes 43 | { 44 | default Gauss linear corrected; 45 | } 46 | 47 | interpolationSchemes 48 | { 49 | default linear; 50 | elcond cellPoint; 51 | } 52 | 53 | snGradSchemes 54 | { 55 | default corrected; 56 | } 57 | 58 | wallDist 59 | { 60 | method meshWave; 61 | } 62 | 63 | // ************************************************************************* // 64 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/system/fvSchemes: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "system"; 14 | object fvSchemes; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | ddtSchemes 19 | { 20 | default Euler; 21 | } 22 | 23 | gradSchemes 24 | { 25 | default Gauss linear; 26 | } 27 | 28 | divSchemes 29 | { 30 | div(rhoPhi,U) Gauss linearUpwindV grad(U); 31 | div(phi,alpha) Gauss vanLeer; 32 | div(phirb,alpha) Gauss linear; 33 | div(phi,k) Gauss upwind; 34 | div(phi,epsilon) Gauss upwind; 35 | div(phi,R) Gauss upwind; 36 | div(R) Gauss linear; 37 | div(phi,nuTilda) Gauss upwind; 38 | div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; 39 | "div\(phi,(k|omega)\)" Gauss upwind; 40 | } 41 | 42 | laplacianSchemes 43 | { 44 | default Gauss linear corrected; 45 | } 46 | 47 | interpolationSchemes 48 | { 49 | default linear; 50 | elcond cellPoint; 51 | } 52 | 53 | snGradSchemes 54 | { 55 | default corrected; 56 | } 57 | 58 | wallDist 59 | { 60 | method meshWave; 61 | } 62 | 63 | // ************************************************************************* // 64 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/system/fvSchemes: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location "system"; 14 | object fvSchemes; 15 | } 16 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 17 | 18 | ddtSchemes 19 | { 20 | default Euler; 21 | } 22 | 23 | gradSchemes 24 | { 25 | default Gauss linear; 26 | } 27 | 28 | divSchemes 29 | { 30 | div(rhoPhi,U) Gauss linearUpwindV grad(U); 31 | div(phi,alpha) Gauss vanLeer; 32 | div(phirb,alpha) Gauss linear; 33 | div(phi,k) Gauss upwind; 34 | div(phi,epsilon) Gauss upwind; 35 | div(phi,R) Gauss upwind; 36 | div(R) Gauss linear; 37 | div(phi,nuTilda) Gauss upwind; 38 | div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear; 39 | "div\(phi,(k|omega)\)" Gauss upwind; 40 | } 41 | 42 | laplacianSchemes 43 | { 44 | default Gauss linear corrected; 45 | } 46 | 47 | interpolationSchemes 48 | { 49 | default linear; 50 | elcond cellPoint; 51 | } 52 | 53 | snGradSchemes 54 | { 55 | default corrected; 56 | } 57 | 58 | wallDist 59 | { 60 | method meshWave; 61 | } 62 | 63 | // ************************************************************************* // 64 | -------------------------------------------------------------------------------- /Benchmark/casebase/system/fvSolution: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | object fvSolution; 14 | } 15 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 16 | 17 | solvers 18 | { 19 | alpha.melt 20 | { 21 | nAlphaCorr 2; 22 | nAlphaSubCycles 2; 23 | cAlpha 2; 24 | } 25 | 26 | pcorr 27 | { 28 | solver PCG; 29 | preconditioner DIC; 30 | tolerance 1e-7; 31 | relTol 0; 32 | } 33 | 34 | pcorrFinal 35 | { 36 | $pcorr; 37 | relTol 0; 38 | } 39 | 40 | p_rgh 41 | { 42 | solver PCG; 43 | preconditioner DIC; 44 | tolerance 1e-07; 45 | relTol 0.05; 46 | } 47 | 48 | p_rghFinal 49 | { 50 | $p_rgh; 51 | relTol 0; 52 | } 53 | 54 | "(U|k|omega)" 55 | { 56 | solver smoothSolver; 57 | smoother symGaussSeidel; 58 | tolerance 1e-7; 59 | relTol 0.1; 60 | } 61 | 62 | "(U|k|omega)Final" 63 | { 64 | $U; 65 | relTol 0; 66 | } 67 | } 68 | 69 | PIMPLE 70 | { 71 | momentumPredictor yes; 72 | nOuterCorrectors 1; 73 | nCorrectors 2; 74 | nNonOrthogonalCorrectors 1; 75 | 76 | pRefPoint (0.0 0.0 0.0); 77 | pRefValue 0; 78 | } 79 | 80 | 81 | // ************************************************************************* // 82 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/system/fvSolution: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | object fvSolution; 14 | } 15 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 16 | 17 | solvers 18 | { 19 | alpha.melt 20 | { 21 | nAlphaCorr 2; 22 | nAlphaSubCycles 2; 23 | cAlpha 2; 24 | } 25 | 26 | pcorr 27 | { 28 | solver PCG; 29 | preconditioner DIC; 30 | tolerance 1e-7; 31 | relTol 0; 32 | } 33 | 34 | pcorrFinal 35 | { 36 | $pcorr; 37 | relTol 0; 38 | } 39 | 40 | p_rgh 41 | { 42 | solver PCG; 43 | preconditioner DIC; 44 | tolerance 1e-07; 45 | relTol 0.05; 46 | } 47 | 48 | p_rghFinal 49 | { 50 | $p_rgh; 51 | relTol 0; 52 | } 53 | 54 | "(U|k|omega)" 55 | { 56 | solver smoothSolver; 57 | smoother symGaussSeidel; 58 | tolerance 1e-7; 59 | relTol 0.1; 60 | } 61 | 62 | "(U|k|omega)Final" 63 | { 64 | $U; 65 | relTol 0; 66 | } 67 | } 68 | 69 | PIMPLE 70 | { 71 | momentumPredictor yes; 72 | nOuterCorrectors 1; 73 | nCorrectors 2; 74 | nNonOrthogonalCorrectors 1; 75 | 76 | pRefPoint (0.0 0.0 0.0); 77 | pRefValue 0; 78 | } 79 | 80 | 81 | // ************************************************************************* // 82 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/system/fvSolution: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | ========= | | 3 | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | 4 | | \\ / O peration | Version: dev | 5 | | \\ / A nd | Web: www.OpenFOAM.org | 6 | | \\/ M anipulation | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | object fvSolution; 14 | } 15 | // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // 16 | 17 | solvers 18 | { 19 | alpha.melt 20 | { 21 | nAlphaCorr 2; 22 | nAlphaSubCycles 2; 23 | cAlpha 2; 24 | } 25 | 26 | pcorr 27 | { 28 | solver PCG; 29 | preconditioner DIC; 30 | tolerance 1e-7; 31 | relTol 0; 32 | } 33 | 34 | pcorrFinal 35 | { 36 | $pcorr; 37 | relTol 0; 38 | } 39 | 40 | p_rgh 41 | { 42 | solver PCG; 43 | preconditioner DIC; 44 | tolerance 1e-07; 45 | relTol 0.05; 46 | } 47 | 48 | p_rghFinal 49 | { 50 | $p_rgh; 51 | relTol 0; 52 | } 53 | 54 | "(U|k|omega)" 55 | { 56 | solver smoothSolver; 57 | smoother symGaussSeidel; 58 | tolerance 1e-7; 59 | relTol 0.1; 60 | } 61 | 62 | "(U|k|omega)Final" 63 | { 64 | $U; 65 | relTol 0; 66 | } 67 | } 68 | 69 | PIMPLE 70 | { 71 | momentumPredictor yes; 72 | nOuterCorrectors 1; 73 | nCorrectors 2; 74 | nNonOrthogonalCorrectors 1; 75 | 76 | pRefPoint (0.0 0.0 0.0); 77 | pRefValue 0; 78 | } 79 | 80 | 81 | // ************************************************************************* // 82 | -------------------------------------------------------------------------------- /src/Interface_program.sh: -------------------------------------------------------------------------------- 1 | # copyright (c) 2020 ASandip 2 | #This program is free software: you can redistribute it and/or 3 | # modify it under the terms of the GNU General Public License as 4 | # published by the Free Software Foundation, either version 3 of the License, or 5 | # (at your option) any later version. 6 | # This program is distributed in the hope that it will be useful, 7 | # but WITHOUT ANY WARRANTY; without even the implied warranty of 8 | # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 9 | # GNU General Public License for more details. 10 | 11 | # You should have received a copy of the GNU General Public License 12 | # along with this program. If not, see . 13 | 14 | 15 | 16 | #!/bin/bash 17 | 18 | # $1 and $2 are special variables in bash that contain the 1st and 2nd 19 | # command line arguments to the script, which are the names of the 20 | # Dakota parameters and results files, respectively. 21 | 22 | params=$1 23 | results=$2 24 | 25 | ############################################################################### 26 | ## 27 | ##Pre-processing Phase -- Transfer new set of values for the chosen input parameters 28 | ##from Dakota to Elmer and OpenFoam (OF) 29 | ## 30 | ############################################################################### 31 | 32 | dprepro $params Elmer.template case.sif 33 | 34 | dprepro $params OF.template OF_file 35 | 36 | 37 | ############################################################################### 38 | ## 39 | ## Execution Phase -- Run your simulation 40 | ## 41 | ############################################################################### 42 | 43 | #Invoke bash script, test.sh, which in turn invokes EOF library that begins 44 | ##the simulation (Elmer and OpenFoam, run in parallel) 45 | ## and enables communication between Elmer FEM and OpenFOAM during runtime 46 | 47 | sudo bash test.sh 48 | 49 | ############################################################################### 50 | ## 51 | ## Post-processing Phase -- Extract (or calculate) quantities of interest 52 | ## from your simulation's output (Elmer and OpenFoam) and write them to a properly-formatted 53 | ## Dakota results file. 54 | ## 55 | ############################################################################## 56 | 57 | 58 | simulation_results > $results 59 | 60 | 61 | 62 | 63 | 64 | -------------------------------------------------------------------------------- /Examples/UQ/Dakota files/EOF_LHS.in: -------------------------------------------------------------------------------- 1 | environment 2 | #Store results in tabular data file 3 | tabular_data 4 | tabular_data_file = 'UQ_results.dat' 5 | 6 | method 7 | #Define data analysis -- UQ study 8 | sampling 9 | sample_type lhs #lhs -- latin hypercube sampling 10 | samples 100 11 | seed 1987 12 | 13 | variables 14 | #Define input parameter/s 15 | normal_uncertain 4 #Input parameters vary normally 16 | descriptors 'r' 's' 'f' 'I' #Input parameters, 'r' melt density, 's' surface tension, 'f' AC frequency, 'I' current density 17 | means 2500 5.1 21250 5e7 18 | std_deviations 2000 4.9 18750 3e7 19 | 20 | 21 | responses 22 | #Define output parameter/s 23 | response_functions 1 24 | num_field_responses 1 lengths = 2886 num_coordinates_per_field = 3 #length of lorentz force vector field = 2886, coordinates = x,y and z 25 | descriptors 'lorentz_force' 26 | no_gradients 27 | no_hessians 28 | 29 | 30 | interface 31 | #Define Dakota - EOF Library interface 32 | fork #External simulation interface type 33 | asynchronous #Asynchronous local parallelism applied 34 | evaluation_concurrency = 4 #4 iterations performed simultaneously 35 | analysis_driver ='EOFD_driver.sh' #Calls the program that interfaces Dakota with EOF Library for this test case 36 | parameters_file='params.in' #Dakota parameters file 37 | results_file='/Dakota/Gui/Workspace/Project/results.out' #Dakota results file 38 | work_directory directory_tag #Creates and tags work directory for every iteration, the work directory stores case files along with the results for that iteration 39 | 40 | 41 | copy_files = 'templatedir/*' #Files stored in 'templatedir' folder would be copied to work directory corresponding every iteration 42 | named 'workdir' file_save directory_save #Work directory, for every iteration, named 'workdir' saved in the Dakota workspace 43 | 44 | 45 | 46 | 47 | 48 | 49 | 50 | -------------------------------------------------------------------------------- /Benchmark/Dakota files/EOF_LPS.sh: -------------------------------------------------------------------------------- 1 | # copyright (c) 2020 ASandip 2 | #This program is free software: you can redistribute it and/or 3 | # modify it under the terms of the GNU General Public License as 4 | # published by the Free Software Foundation, either version 3 of the License, or 5 | # (at your option) any later version. 6 | # This program is distributed in the hope that it will be useful, 7 | # but WITHOUT ANY WARRANTY; without even the implied warranty of 8 | # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 9 | # GNU General Public License for more details. 10 | 11 | # You should have received a copy of the GNU General Public License 12 | # along with this program. If not, see . 13 | 14 | 15 | 16 | #!/bin/bash 17 | 18 | # $1 and $2 are special variables in bash that contain the 1st and 2nd 19 | # command line arguments to the script, which are the names of the 20 | # Dakota parameters and results files, respectively. 21 | 22 | params=$1 23 | results=$2 24 | 25 | ############################################################################### 26 | ## 27 | ## Pre-processing Phase -- Generate/configure an input file for your simulation 28 | ## by substiting in parameter values from the Dakota paramters file. 29 | ## 30 | ############################################################################### 31 | 32 | #Transfer new set of values for the chosen input parameters -- melt density, 33 | #surface tension, ac frequency and current density -- from Dakota to Elmer and OpenFoam 34 | 35 | dprepro $params sim.template transportProperties 36 | 37 | dprepro --inline "@ @" $params case.template case_d.sif 38 | 39 | 40 | ############################################################################### 41 | ## 42 | ## Execution Phase -- Run your simulation 43 | ## 44 | ############################################################################### 45 | 46 | #Execute bash script, test.sh: Invokes EOF library that begins the simulation (Elmer and OpenFoam, run in parallel) 47 | ## and enables communication between Elmer FEM and OpenFOAM during runtime 48 | 49 | sudo bash test.sh 50 | 51 | ############################################################################### 52 | ## 53 | ## Post-processing Phase -- Extract (or calculate) quantities of interest 54 | ## from your simulation's output and write them to a properly-formatted 55 | ## Dakota results file. 56 | ## 57 | ############################################################################## 58 | 59 | # Quantity of interest -- lorentz force at the final time step 60 | 61 | awk 'NR == 23, NR==984 { print $1, $2, $3 }' processor1/0.49/JxB > lorentz_force.txt 62 | sed 's/[()]//g' lorentz_force.txt > lf.txt 63 | cp lf.txt $results 64 | 65 | 66 | 67 | 68 | -------------------------------------------------------------------------------- /Examples/UQ/Dakota files/EOFD_driver.sh: -------------------------------------------------------------------------------- 1 | # copyright (c) 2020 ASandip 2 | #This program is free software: you can redistribute it and/or 3 | # modify it under the terms of the GNU General Public License as 4 | # published by the Free Software Foundation, either version 3 of the License, or 5 | # (at your option) any later version. 6 | # This program is distributed in the hope that it will be useful, 7 | # but WITHOUT ANY WARRANTY; without even the implied warranty of 8 | # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 9 | # GNU General Public License for more details. 10 | 11 | # You should have received a copy of the GNU General Public License 12 | # along with this program. If not, see . 13 | 14 | 15 | 16 | #!/bin/bash 17 | 18 | # $1 and $2 are special variables in bash that contain the 1st and 2nd 19 | # command line arguments to the script, which are the names of the 20 | # Dakota parameters and results files, respectively. 21 | 22 | params=$1 23 | results=$2 24 | 25 | ############################################################################### 26 | ## 27 | ## Pre-processing Phase -- Generate/configure an input file for your simulation 28 | ## by substiting in parameter values from the Dakota paramters file. 29 | ## 30 | ############################################################################### 31 | 32 | #Transfer new set of values for the chosen input parameters -- melt density, 33 | #surface tension, ac frequency and current density -- from Dakota to Elmer and OpenFoam 34 | 35 | dprepro $params sim.template transportProperties 36 | 37 | dprepro --inline "@ @" $params case.template case_d.sif 38 | 39 | 40 | ############################################################################### 41 | ## 42 | ## Execution Phase -- Run your simulation 43 | ## 44 | ############################################################################### 45 | 46 | #Execute bash script, test.sh: Invokes EOF library that begins the simulation (Elmer and OpenFoam, run in parallel) 47 | ## and enables communication between Elmer FEM and OpenFOAM during runtime 48 | 49 | sudo bash test.sh 50 | 51 | ############################################################################### 52 | ## 53 | ## Post-processing Phase -- Extract (or calculate) quantities of interest 54 | ## from your simulation's output and write them to a properly-formatted 55 | ## Dakota results file. 56 | ## 57 | ############################################################################## 58 | 59 | # Quantity of interest -- lorentz force at the final time step 60 | 61 | awk 'NR == 23, NR==984 { print $1, $2, $3 }' processor1/0.49/JxB > lorentz_force.txt 62 | sed 's/[()]//g' lorentz_force.txt > lf.txt 63 | cp lf.txt $results 64 | 65 | 66 | 67 | 68 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/Dakota files/EOFD_driver.sh: -------------------------------------------------------------------------------- 1 | # copyright (c) 2020 ASandip 2 | #This program is free software: you can redistribute it and/or 3 | # modify it under the terms of the GNU General Public License as 4 | # published by the Free Software Foundation, either version 3 of the License, or 5 | # (at your option) any later version. 6 | # This program is distributed in the hope that it will be useful, 7 | # but WITHOUT ANY WARRANTY; without even the implied warranty of 8 | # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 9 | # GNU General Public License for more details. 10 | 11 | # You should have received a copy of the GNU General Public License 12 | # along with this program. If not, see . 13 | 14 | 15 | 16 | #!/bin/bash 17 | 18 | # $1 and $2 are special variables in bash that contain the 1st and 2nd 19 | # command line arguments to the script, which are the names of the 20 | # Dakota parameters and results files, respectively. 21 | 22 | params=$1 23 | results=$2 24 | 25 | ############################################################################### 26 | ## 27 | ## Pre-processing Phase -- Generate/configure an input file for your simulation 28 | ## by substiting in parameter values from the Dakota paramters file. 29 | ## 30 | ############################################################################### 31 | 32 | #Transfer new set of values for the chosen input parameters -- melt density, 33 | #surface tension, ac frequency and current density -- from Dakota to Elmer and OpenFoam 34 | 35 | dprepro $params sim.template transportProperties 36 | 37 | dprepro --inline "@ @" $params case.template case_d.sif 38 | 39 | 40 | ############################################################################### 41 | ## 42 | ## Execution Phase -- Run your simulation 43 | ## 44 | ############################################################################### 45 | 46 | #Execute bash script, test.sh: Invokes EOF library that begins the simulation (Elmer and OpenFoam, run in parallel) 47 | ## and enables communication between Elmer FEM and OpenFOAM during runtime 48 | 49 | sudo bash test.sh 50 | 51 | ############################################################################### 52 | ## 53 | ## Post-processing Phase -- Extract (or calculate) quantities of interest 54 | ## from your simulation's output and write them to a properly-formatted 55 | ## Dakota results file. 56 | ## 57 | ############################################################################## 58 | 59 | # Quantity of interest -- lorentz force at the final time step 60 | 61 | awk 'NR == 23, NR==984 { print $1, $2, $3 }' processor1/0.49/JxB > lorentz_force.txt 62 | sed 's/[()]//g' lorentz_force.txt > lf.txt 63 | cp lf.txt $results 64 | 65 | 66 | 67 | 68 | -------------------------------------------------------------------------------- /Benchmark/Dakota files/EOF_LPS.in: -------------------------------------------------------------------------------- 1 | environment 2 | #Store results in tabular data file 3 | tabular_data 4 | tabular_data_file = 'LPS_results.dat' 5 | 6 | method 7 | #Define data analysis -- parameter study 8 | list_parameter_study 9 | list_of_points = 500 0.94 9800 3e7 10 | 1380 0.94 9800 3e7 11 | 2380 0.94 9800 3e7 12 | 3380 0.94 9800 3e7 13 | 4500 0.94 9800 3e7 14 | 2380 0.2 9800 3e7 15 | 2380 0.5 9800 3e7 16 | 2380 3.0 9800 3e7 17 | 2380 10.0 9800 3e7 18 | 2380 0.94 2500 3e7 19 | 2380 0.94 4500 3e7 20 | 2380 0.94 18000 3e7 21 | 2380 0.94 40000 3e7 22 | 2380 0.94 9800 2e7 23 | 2380 0.94 9800 2.5e7 24 | 2380 0.94 9800 5e7 25 | 2380 0.94 9800 8e7 26 | 27 | 28 | variables 29 | #Define input parameter/s 30 | continuous_design = 4 #Input parameters vary continuously 31 | descriptors 'r' 's' 'f' 'I' #Input parameters, 'r' melt density, 's' surface tension, 'f' AC frequency, 'I' current density 32 | 33 | 34 | responses 35 | #Define output parameter/s 36 | response_functions 1 37 | num_field_responses 1 lengths = 2886 num_coordinates_per_field = 3 #length of lorentz force vector field = 2886, coordinates = x,y and z 38 | descriptors 'lorentz_force' 39 | no_gradients 40 | no_hessians 41 | 42 | 43 | interface 44 | #Define Dakota - EOF Library interface 45 | fork #External simulation interface type 46 | analysis_driver ='EOF_LPS.sh' #Calls the program that interfaces Dakota with EOF Library for this test case 47 | parameters_file='params.in' #Dakota parameters file 48 | results_file='/Dakota/Gui/Workspace/Project/results.out' #Dakota results file 49 | work_directory directory_tag #Creates and tags work directory for every iteration, the work directory stores case files along with the results for that iteration 50 | 51 | copy_files = 'templatedir/*' #Files stored in 'templatedir' folder would be copied to work directory corresponding every iteration 52 | named 'workdir' file_save directory #Work directory, for every iteration, named 'workdir' saved in the Dakota workspace 53 | 54 | 55 | 56 | 57 | 58 | 59 | 60 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/Dakota files/EOF_LHS_SM.in: -------------------------------------------------------------------------------- 1 | environment 2 | #Store results in tabular data file 3 | tabular_data 4 | tabular_data_file = 'UQ_SM_results.dat' 5 | 6 | method 7 | id_method = 'UQ_on_surrogate' 8 | #Define data analysis -- UQ study 9 | sampling #latin hypercube sampling study was performed on the global surrogate model 10 | sample_type lhs #lhs -- latin hypercube sampling 11 | samples 100 12 | seed 1987 13 | 14 | 15 | model_pointer = 'surr_model' 16 | output verbose 17 | 18 | 19 | model 20 | id_model = 'surr_model' 21 | surrogate global 22 | dace_method_pointer = 'DACE' 23 | polynomial linear #linear regression fit was applied to the training data to construct the global surrogate model 24 | 25 | method 26 | id_method = 'DACE' #Design and analysis of computer experiments 'DACE' was employed to generate training data 27 | sampling 28 | samples = 20 #A sample size of 20 was chosen to generate the training data 29 | model_pointer = 'MHD_model' 30 | 31 | model 32 | id_model = 'MHD_model' 33 | single 34 | interface_pointer = 'MHD_interface' 35 | 36 | variables 37 | #Define input parameter/s 38 | normal_uncertain 4 #Input parameters vary normally 39 | descriptors 'r' 's' 'f' 'I' #Input parameters, 'r' melt density, 's' surface tension, 'f' AC frequency, 'I' current density 40 | means 2500 5.1 21250 5e7 41 | std_deviations 2000 4.9 18750 3e7 42 | 43 | 44 | responses 45 | #Define output parameter/s 46 | response_functions 1 47 | num_field_responses 1 lengths = 2886 num_coordinates_per_field = 3 #length of lorentz force vector field = 2886, coordinates = x,y and z 48 | descriptors 'lorentz_force' 49 | no_gradients 50 | no_hessians 51 | 52 | 53 | interface 54 | #Define Dakota - EOF Library interface 55 | 56 | id_interface = 'MHD_interface' 57 | 58 | fork #External simulation interface type 59 | asynchronous #Asynchronous local parallelism applied 60 | evaluation_concurrency = 4 #4 iterations performed simultaneously 61 | analysis_driver ='EOFD_driver.sh' #Calls the program that interfaces Dakota with EOF Library for this test case 62 | parameters_file='params.in' #Dakota parameters file 63 | results_file='/Dakota/Gui/Workspace/Project/results.out' #Dakota results file 64 | work_directory directory_tag #Creates and tags work directory for every iteration, the work directory stores case files along with the results for that iteration 65 | 66 | 67 | copy_files = 'templatedir/*' #Files stored in 'templatedir' folder would be copied to work directory corresponding every iteration 68 | named 'workdir' file_save directory_save #Work directory, for every iteration, named 'workdir' saved in the Dakota workspace 69 | 70 | 71 | 72 | 73 | 74 | -------------------------------------------------------------------------------- /Benchmark/templatedir/case.template: -------------------------------------------------------------------------------- 1 | Header 2 | CHECK KEYWORDS Warn 3 | Mesh DB "." "meshElmer" 4 | Include Path "" 5 | Results Directory "" 6 | End 7 | 8 | $frequency = @f@ 9 | $current = @I@ 10 | 11 | Simulation 12 | Max Output Level = 3 13 | Coordinate System = Axi Symmetric 14 | Angular Frequency = $2*pi*@f@ 15 | 16 | Simulation Type = Steady 17 | Steady State Max Iterations = 10000000 18 | Steady State Min Iterations = 10000000 19 | 20 | Solver Input File = case.sif 21 | Post File = case.vtu 22 | Output Intervals = 2 23 | End 24 | 25 | Constants 26 | Permittivity of Vacuum = 8.8542e-12 27 | Permeability of Vacuum = 1.257e-6 28 | End 29 | 30 | Body 1 31 | Target Bodies(1) = 3 32 | Name = "inductorsTop" 33 | Equation = 1 34 | Material = 1 35 | Body Force = 1 36 | End 37 | 38 | Body 2 39 | Target Bodies(1) = 4 40 | Name = "inductorsBtm" 41 | Equation = 1 42 | Material = 1 43 | Body Force = 2 44 | End 45 | 46 | Body 3 47 | Target Bodies(1) = 1 48 | Name = "melt" 49 | Equation = 1 50 | Material = 3 51 | End 52 | 53 | Body 4 54 | Target Bodies(1) = 2 55 | Name = "air" 56 | Equation = 1 57 | Material = 2 58 | End 59 | 60 | Solver 1 61 | Equation = OpenFOAM2Elmer 62 | Procedure = "OpenFOAM2Elmer" "OpenFOAM2ElmerSolver" 63 | 64 | Target Variable 1 = String "OF conductivity" 65 | End 66 | 67 | Solver 2 68 | Equation = "MDynamics" 69 | Variable = AV[AV re:1 AV im:1] 70 | Procedure = "MagnetoDynamics2D" "MagnetoDynamics2DHarmonic" 71 | 72 | Linear System Symmetric = True 73 | 74 | Linear System Solver = "Iterative" 75 | Linear System Preconditioning = ILU0 76 | Linear System Residual Output = 10 77 | Linear System Max Iterations = 500 78 | Linear System Iterative Method = BiCGStabl 79 | Linear System Convergence Tolerance = 1.0e-8 80 | 81 | Steady State Convergence Tolerance = 1e-08 82 | Nonlinear System Max Iterations = 3 83 | Nonlinear System Convergence Tolerance = 1.0e-6 84 | Nonlinear System Relaxation Factor = 1 85 | 86 | Exported Variable 1 = -dg "OF conductivity" 87 | Exported Variable 1 Mask = String "material: mapping" 88 | End 89 | 90 | Solver 3 91 | Equation = "Postprocess" 92 | 93 | Potential Variable = String "AV" 94 | Procedure = "MagnetoDynamics" "MagnetoDynamicsCalcFields" 95 | 96 | Linear System Solver = Iterative 97 | Linear System Iterative Method = CG 98 | Linear System Symmetric = True 99 | Linear System Max Iterations = 100 100 | Linear System Preconditioning = None 101 | Linear System Convergence Tolerance = 1.0e-8 102 | Steady State Convergence Tolerance = 1e-08 103 | 104 | Nonlinear System Max Iterations = 3 105 | Nonlinear System Convergence Tolerance = 1.0e-8 106 | 107 | Calculate Current Density = Logical True 108 | Calculate Magnetic Field Strength = Logical True 109 | Calculate JxB = Logical True 110 | End 111 | 112 | Solver 4 113 | Equation = ElmerToOpenFOAM 114 | Procedure = "Elmer2OpenFOAM" "Elmer2OpenFOAMSolver" 115 | 116 | Target Variable 1 = jxb re 1 117 | Target Variable 2 = jxb re 2 118 | Target Variable 3 = jxb re 3 119 | End 120 | 121 | Equation 1 122 | Name = "MGDyn" 123 | Active Solvers(4) = 1 2 3 4 124 | End 125 | 126 | Material 1 127 | Name = "Copper" 128 | Electric Conductivity = 0 129 | Relative Permittivity = 1 130 | Relative Permeability = 1 131 | End 132 | 133 | Material 2 134 | Name = "Air" 135 | Electric Conductivity = 0 136 | Relative Permittivity = 1 137 | Relative Permeability = 1 138 | End 139 | 140 | Material 3 141 | Name = "Melt" 142 | Electric Conductivity = Equals "OF conductivity" 143 | Relative Permittivity = 1 144 | Relative Permeability = 1 145 | mapping = Logical True 146 | End 147 | 148 | Body Force 1 149 | Name = "Current pos" 150 | Current Density = Real @I@ 151 | End 152 | 153 | Body Force 2 154 | Name = "Current neg" 155 | Current Density = Real -@I@ 156 | End 157 | 158 | Boundary Condition 1 159 | Target Boundaries(1) = 5 160 | AV re {e} = Real 0.0 161 | AV im {e} = Real 0.0 162 | End 163 | 164 | Solver 2 165 | Reference Norm = Real 6.16456343E-04 166 | Reference Norm Tolerance = Real 1E-8 167 | End 168 | -------------------------------------------------------------------------------- /Examples/UQ/templatedir/case.template: -------------------------------------------------------------------------------- 1 | Header 2 | CHECK KEYWORDS Warn 3 | Mesh DB "." "meshElmer" 4 | Include Path "" 5 | Results Directory "" 6 | End 7 | 8 | $frequency = @f@ 9 | $current = @I@ 10 | 11 | Simulation 12 | Max Output Level = 3 13 | Coordinate System = Axi Symmetric 14 | Angular Frequency = $2*pi*@f@ 15 | 16 | Simulation Type = Steady 17 | Steady State Max Iterations = 10000000 18 | Steady State Min Iterations = 10000000 19 | 20 | Solver Input File = case.sif 21 | Post File = case.vtu 22 | Output Intervals = 2 23 | End 24 | 25 | Constants 26 | Permittivity of Vacuum = 8.8542e-12 27 | Permeability of Vacuum = 1.257e-6 28 | End 29 | 30 | Body 1 31 | Target Bodies(1) = 3 32 | Name = "inductorsTop" 33 | Equation = 1 34 | Material = 1 35 | Body Force = 1 36 | End 37 | 38 | Body 2 39 | Target Bodies(1) = 4 40 | Name = "inductorsBtm" 41 | Equation = 1 42 | Material = 1 43 | Body Force = 2 44 | End 45 | 46 | Body 3 47 | Target Bodies(1) = 1 48 | Name = "melt" 49 | Equation = 1 50 | Material = 3 51 | End 52 | 53 | Body 4 54 | Target Bodies(1) = 2 55 | Name = "air" 56 | Equation = 1 57 | Material = 2 58 | End 59 | 60 | Solver 1 61 | Equation = OpenFOAM2Elmer 62 | Procedure = "OpenFOAM2Elmer" "OpenFOAM2ElmerSolver" 63 | 64 | Target Variable 1 = String "OF conductivity" 65 | End 66 | 67 | Solver 2 68 | Equation = "MDynamics" 69 | Variable = AV[AV re:1 AV im:1] 70 | Procedure = "MagnetoDynamics2D" "MagnetoDynamics2DHarmonic" 71 | 72 | Linear System Symmetric = True 73 | 74 | Linear System Solver = "Iterative" 75 | Linear System Preconditioning = ILU0 76 | Linear System Residual Output = 10 77 | Linear System Max Iterations = 500 78 | Linear System Iterative Method = BiCGStabl 79 | Linear System Convergence Tolerance = 1.0e-8 80 | 81 | Steady State Convergence Tolerance = 1e-08 82 | Nonlinear System Max Iterations = 3 83 | Nonlinear System Convergence Tolerance = 1.0e-6 84 | Nonlinear System Relaxation Factor = 1 85 | 86 | Exported Variable 1 = -dg "OF conductivity" 87 | Exported Variable 1 Mask = String "material: mapping" 88 | End 89 | 90 | Solver 3 91 | Equation = "Postprocess" 92 | 93 | Potential Variable = String "AV" 94 | Procedure = "MagnetoDynamics" "MagnetoDynamicsCalcFields" 95 | 96 | Linear System Solver = Iterative 97 | Linear System Iterative Method = CG 98 | Linear System Symmetric = True 99 | Linear System Max Iterations = 100 100 | Linear System Preconditioning = None 101 | Linear System Convergence Tolerance = 1.0e-8 102 | Steady State Convergence Tolerance = 1e-08 103 | 104 | Nonlinear System Max Iterations = 3 105 | Nonlinear System Convergence Tolerance = 1.0e-8 106 | 107 | Calculate Current Density = Logical True 108 | Calculate Magnetic Field Strength = Logical True 109 | Calculate JxB = Logical True 110 | End 111 | 112 | Solver 4 113 | Equation = ElmerToOpenFOAM 114 | Procedure = "Elmer2OpenFOAM" "Elmer2OpenFOAMSolver" 115 | 116 | Target Variable 1 = jxb re 1 117 | Target Variable 2 = jxb re 2 118 | Target Variable 3 = jxb re 3 119 | End 120 | 121 | Equation 1 122 | Name = "MGDyn" 123 | Active Solvers(4) = 1 2 3 4 124 | End 125 | 126 | Material 1 127 | Name = "Copper" 128 | Electric Conductivity = 0 129 | Relative Permittivity = 1 130 | Relative Permeability = 1 131 | End 132 | 133 | Material 2 134 | Name = "Air" 135 | Electric Conductivity = 0 136 | Relative Permittivity = 1 137 | Relative Permeability = 1 138 | End 139 | 140 | Material 3 141 | Name = "Melt" 142 | Electric Conductivity = Equals "OF conductivity" 143 | Relative Permittivity = 1 144 | Relative Permeability = 1 145 | mapping = Logical True 146 | End 147 | 148 | Body Force 1 149 | Name = "Current pos" 150 | Current Density = Real @I@ 151 | End 152 | 153 | Body Force 2 154 | Name = "Current neg" 155 | Current Density = Real -@I@ 156 | End 157 | 158 | Boundary Condition 1 159 | Target Boundaries(1) = 5 160 | AV re {e} = Real 0.0 161 | AV im {e} = Real 0.0 162 | End 163 | 164 | Solver 2 165 | Reference Norm = Real 6.16456343E-04 166 | Reference Norm Tolerance = Real 1E-8 167 | End 168 | -------------------------------------------------------------------------------- /Benchmark/casebase/case.sif: -------------------------------------------------------------------------------- 1 | Header 2 | CHECK KEYWORDS Warn 3 | Mesh DB "." "meshElmer" 4 | Include Path "" 5 | Results Directory "" 6 | End 7 | 8 | $frequency = 9800 9 | $current = 30000000 10 | 11 | Simulation 12 | Max Output Level = 3 13 | Coordinate System = Axi Symmetric 14 | Angular Frequency = $2*pi*9800 15 | 16 | Simulation Type = Steady 17 | Steady State Max Iterations = 10000000 18 | Steady State Min Iterations = 10000000 19 | 20 | Solver Input File = case.sif 21 | Post File = case.vtu 22 | Output Intervals = 2 23 | End 24 | 25 | Constants 26 | Permittivity of Vacuum = 8.8542e-12 27 | Permeability of Vacuum = 1.257e-6 28 | End 29 | 30 | Body 1 31 | Target Bodies(1) = 3 32 | Name = "inductorsTop" 33 | Equation = 1 34 | Material = 1 35 | Body Force = 1 36 | End 37 | 38 | Body 2 39 | Target Bodies(1) = 4 40 | Name = "inductorsBtm" 41 | Equation = 1 42 | Material = 1 43 | Body Force = 2 44 | End 45 | 46 | Body 3 47 | Target Bodies(1) = 1 48 | Name = "melt" 49 | Equation = 1 50 | Material = 3 51 | End 52 | 53 | Body 4 54 | Target Bodies(1) = 2 55 | Name = "air" 56 | Equation = 1 57 | Material = 2 58 | End 59 | 60 | Solver 1 61 | Equation = OpenFOAM2Elmer 62 | Procedure = "OpenFOAM2Elmer" "OpenFOAM2ElmerSolver" 63 | 64 | Target Variable 1 = String "OF conductivity" 65 | End 66 | 67 | Solver 2 68 | Equation = "MDynamics" 69 | Variable = AV[AV re:1 AV im:1] 70 | Procedure = "MagnetoDynamics2D" "MagnetoDynamics2DHarmonic" 71 | 72 | Linear System Symmetric = True 73 | 74 | Linear System Solver = "Iterative" 75 | Linear System Preconditioning = ILU0 76 | Linear System Residual Output = 10 77 | Linear System Max Iterations = 500 78 | Linear System Iterative Method = BiCGStabl 79 | Linear System Convergence Tolerance = 1.0e-8 80 | 81 | Steady State Convergence Tolerance = 1e-08 82 | Nonlinear System Max Iterations = 3 83 | Nonlinear System Convergence Tolerance = 1.0e-6 84 | Nonlinear System Relaxation Factor = 1 85 | 86 | Exported Variable 1 = -dg "OF conductivity" 87 | Exported Variable 1 Mask = String "material: mapping" 88 | End 89 | 90 | Solver 3 91 | Equation = "Postprocess" 92 | 93 | Potential Variable = String "AV" 94 | Procedure = "MagnetoDynamics" "MagnetoDynamicsCalcFields" 95 | 96 | Linear System Solver = Iterative 97 | Linear System Iterative Method = CG 98 | Linear System Symmetric = True 99 | Linear System Max Iterations = 100 100 | Linear System Preconditioning = None 101 | Linear System Convergence Tolerance = 1.0e-8 102 | Steady State Convergence Tolerance = 1e-08 103 | 104 | Nonlinear System Max Iterations = 3 105 | Nonlinear System Convergence Tolerance = 1.0e-8 106 | 107 | Calculate Current Density = Logical True 108 | Calculate Magnetic Field Strength = Logical True 109 | Calculate JxB = Logical True 110 | End 111 | 112 | Solver 4 113 | Equation = ElmerToOpenFOAM 114 | Procedure = "Elmer2OpenFOAM" "Elmer2OpenFOAMSolver" 115 | 116 | Target Variable 1 = jxb re 1 117 | Target Variable 2 = jxb re 2 118 | Target Variable 3 = jxb re 3 119 | End 120 | 121 | Equation 1 122 | Name = "MGDyn" 123 | Active Solvers(4) = 1 2 3 4 124 | End 125 | 126 | Material 1 127 | Name = "Copper" 128 | Electric Conductivity = 0 129 | Relative Permittivity = 1 130 | Relative Permeability = 1 131 | End 132 | 133 | Material 2 134 | Name = "Air" 135 | Electric Conductivity = 0 136 | Relative Permittivity = 1 137 | Relative Permeability = 1 138 | End 139 | 140 | Material 3 141 | Name = "Melt" 142 | Electric Conductivity = Equals "OF conductivity" 143 | Relative Permittivity = 1 144 | Relative Permeability = 1 145 | mapping = Logical True 146 | End 147 | 148 | Body Force 1 149 | Name = "Current pos" 150 | Current Density = Real 30000000 151 | End 152 | 153 | Body Force 2 154 | Name = "Current neg" 155 | Current Density = Real -30000000 156 | End 157 | 158 | Boundary Condition 1 159 | Target Boundaries(1) = 5 160 | AV re {e} = Real 0.0 161 | AV im {e} = Real 0.0 162 | End 163 | 164 | Solver 2 165 | Reference Norm = Real 6.16456343E-04 166 | Reference Norm Tolerance = Real 1E-8 167 | End 168 | -------------------------------------------------------------------------------- /Benchmark/casebase/system/snappyHexMeshDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | o | | 3 | | o o | HELYX-OS | 4 | | o O o | Version: v2.3.1 | 5 | | o o | Web: http://www.engys.com | 6 | | o | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location system; 14 | object snappyHexMeshDict; 15 | } 16 | 17 | castellatedMesh true; 18 | snap true; 19 | addLayers false; 20 | 21 | geometry 22 | { 23 | meltRegion.stl 24 | { 25 | type triSurfaceMesh; 26 | name meltRegion; 27 | 28 | regions 29 | { 30 | left 31 | { 32 | name leftBC; 33 | } 34 | out 35 | { 36 | name outBC; 37 | } 38 | right 39 | { 40 | name rightBC; 41 | } 42 | } 43 | } 44 | } 45 | 46 | castellatedMeshControls 47 | { 48 | features 49 | ( 50 | ); 51 | 52 | refinementSurfaces 53 | { 54 | meltRegion 55 | { 56 | level (0 0); 57 | regions 58 | { 59 | left 60 | { 61 | level (0 0); 62 | } 63 | out 64 | { 65 | level (0 0); 66 | } 67 | right 68 | { 69 | level (0 0); 70 | } 71 | } 72 | } 73 | } 74 | 75 | refinementRegions 76 | { 77 | } 78 | 79 | locationInMesh (0.01 0 0); 80 | maxLocalCells 10000000; 81 | maxGlobalCells 20000000; 82 | minRefinementCells 0; 83 | nCellsBetweenLevels 1; 84 | resolveFeatureAngle 30.0; 85 | allowFreeStandingZoneFaces true; 86 | planarAngle 30.0; 87 | maxLoadUnbalance 0.1; 88 | } 89 | 90 | snapControls 91 | { 92 | nSolveIter 50; 93 | nSmoothPatch 5; 94 | tolerance 1.0; 95 | nRelaxIter 10; 96 | nFeatureSnapIter 20; 97 | implicitFeatureSnap false; 98 | explicitFeatureSnap true; 99 | multiRegionFeatureSnap false; 100 | } 101 | 102 | addLayersControls 103 | { 104 | layers 105 | { 106 | } 107 | 108 | relativeSizes true; 109 | expansionRatio 1.0; 110 | finalLayerThickness 0.3; 111 | minThickness 0.25; 112 | nGrow 0; 113 | featureAngle 130.0; 114 | slipFeatureAngle 30.0; 115 | nRelaxIter 5; 116 | nSmoothSurfaceNormals 1; 117 | nSmoothNormals 3; 118 | nSmoothThickness 10; 119 | maxFaceThicknessRatio 0.5; 120 | maxThicknessToMedialRatio 0.3; 121 | minMedialAxisAngle 90; 122 | nBufferCellsNoExtrude 0; 123 | nLayerIter 50; 124 | nRelaxedIter 20; 125 | } 126 | 127 | meshQualityControls 128 | { 129 | maxNonOrtho 65.0; 130 | maxBoundarySkewness 20; 131 | maxInternalSkewness 4; 132 | maxConcave 80.0; 133 | minFlatness 0.5; 134 | minVol 1.0E-13; 135 | minTetQuality 1.0E-15; 136 | minArea -1.0; 137 | minTwist 0.02; 138 | minDeterminant 0.001; 139 | minFaceWeight 0.05; 140 | minVolRatio 0.01; 141 | minTriangleTwist -1.0; 142 | nSmoothScale 4; 143 | errorReduction 0.75; 144 | } 145 | 146 | debug 0; 147 | mergeTolerance 1.0E-6; 148 | autoBlockMesh true; 149 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/case.sif: -------------------------------------------------------------------------------- 1 | Header 2 | CHECK KEYWORDS Warn 3 | Mesh DB "." "meshElmer" 4 | Include Path "" 5 | Results Directory "" 6 | End 7 | 8 | $frequency = 9800 9 | $current = 30000000 10 | 11 | Simulation 12 | Max Output Level = 3 13 | Coordinate System = Axi Symmetric 14 | Angular Frequency = $2*pi*9800 15 | 16 | Simulation Type = Steady 17 | Steady State Max Iterations = 10000000 18 | Steady State Min Iterations = 10000000 19 | 20 | Solver Input File = case.sif 21 | Post File = case.vtu 22 | Output Intervals = 2 23 | End 24 | 25 | Constants 26 | Permittivity of Vacuum = 8.8542e-12 27 | Permeability of Vacuum = 1.257e-6 28 | End 29 | 30 | Body 1 31 | Target Bodies(1) = 3 32 | Name = "inductorsTop" 33 | Equation = 1 34 | Material = 1 35 | Body Force = 1 36 | End 37 | 38 | Body 2 39 | Target Bodies(1) = 4 40 | Name = "inductorsBtm" 41 | Equation = 1 42 | Material = 1 43 | Body Force = 2 44 | End 45 | 46 | Body 3 47 | Target Bodies(1) = 1 48 | Name = "melt" 49 | Equation = 1 50 | Material = 3 51 | End 52 | 53 | Body 4 54 | Target Bodies(1) = 2 55 | Name = "air" 56 | Equation = 1 57 | Material = 2 58 | End 59 | 60 | Solver 1 61 | Equation = OpenFOAM2Elmer 62 | Procedure = "OpenFOAM2Elmer" "OpenFOAM2ElmerSolver" 63 | 64 | Target Variable 1 = String "OF conductivity" 65 | End 66 | 67 | Solver 2 68 | Equation = "MDynamics" 69 | Variable = AV[AV re:1 AV im:1] 70 | Procedure = "MagnetoDynamics2D" "MagnetoDynamics2DHarmonic" 71 | 72 | Linear System Symmetric = True 73 | 74 | Linear System Solver = "Iterative" 75 | Linear System Preconditioning = ILU0 76 | Linear System Residual Output = 10 77 | Linear System Max Iterations = 500 78 | Linear System Iterative Method = BiCGStabl 79 | Linear System Convergence Tolerance = 1.0e-8 80 | 81 | Steady State Convergence Tolerance = 1e-08 82 | Nonlinear System Max Iterations = 3 83 | Nonlinear System Convergence Tolerance = 1.0e-6 84 | Nonlinear System Relaxation Factor = 1 85 | 86 | Exported Variable 1 = -dg "OF conductivity" 87 | Exported Variable 1 Mask = String "material: mapping" 88 | End 89 | 90 | Solver 3 91 | Equation = "Postprocess" 92 | 93 | Potential Variable = String "AV" 94 | Procedure = "MagnetoDynamics" "MagnetoDynamicsCalcFields" 95 | 96 | Linear System Solver = Iterative 97 | Linear System Iterative Method = CG 98 | Linear System Symmetric = True 99 | Linear System Max Iterations = 100 100 | Linear System Preconditioning = None 101 | Linear System Convergence Tolerance = 1.0e-8 102 | Steady State Convergence Tolerance = 1e-08 103 | 104 | Nonlinear System Max Iterations = 3 105 | Nonlinear System Convergence Tolerance = 1.0e-8 106 | 107 | Calculate Current Density = Logical True 108 | Calculate Magnetic Field Strength = Logical True 109 | Calculate JxB = Logical True 110 | End 111 | 112 | Solver 4 113 | Equation = ElmerToOpenFOAM 114 | Procedure = "Elmer2OpenFOAM" "Elmer2OpenFOAMSolver" 115 | 116 | Target Variable 1 = jxb re 1 117 | Target Variable 2 = jxb re 2 118 | Target Variable 3 = jxb re 3 119 | End 120 | 121 | Equation 1 122 | Name = "MGDyn" 123 | Active Solvers(4) = 1 2 3 4 124 | End 125 | 126 | Material 1 127 | Name = "Copper" 128 | Electric Conductivity = 0 129 | Relative Permittivity = 1 130 | Relative Permeability = 1 131 | End 132 | 133 | Material 2 134 | Name = "Air" 135 | Electric Conductivity = 0 136 | Relative Permittivity = 1 137 | Relative Permeability = 1 138 | End 139 | 140 | Material 3 141 | Name = "Melt" 142 | Electric Conductivity = Equals "OF conductivity" 143 | Relative Permittivity = 1 144 | Relative Permeability = 1 145 | mapping = Logical True 146 | End 147 | 148 | Body Force 1 149 | Name = "Current pos" 150 | Current Density = Real 30000000 151 | End 152 | 153 | Body Force 2 154 | Name = "Current neg" 155 | Current Density = Real -30000000 156 | End 157 | 158 | Boundary Condition 1 159 | Target Boundaries(1) = 5 160 | AV re {e} = Real 0.0 161 | AV im {e} = Real 0.0 162 | End 163 | 164 | Solver 2 165 | Reference Norm = Real 6.16456343E-04 166 | Reference Norm Tolerance = Real 1E-8 167 | End 168 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/system/snappyHexMeshDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | o | | 3 | | o o | HELYX-OS | 4 | | o O o | Version: v2.3.1 | 5 | | o o | Web: http://www.engys.com | 6 | | o | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location system; 14 | object snappyHexMeshDict; 15 | } 16 | 17 | castellatedMesh true; 18 | snap true; 19 | addLayers false; 20 | 21 | geometry 22 | { 23 | meltRegion.stl 24 | { 25 | type triSurfaceMesh; 26 | name meltRegion; 27 | 28 | regions 29 | { 30 | left 31 | { 32 | name leftBC; 33 | } 34 | out 35 | { 36 | name outBC; 37 | } 38 | right 39 | { 40 | name rightBC; 41 | } 42 | } 43 | } 44 | } 45 | 46 | castellatedMeshControls 47 | { 48 | features 49 | ( 50 | ); 51 | 52 | refinementSurfaces 53 | { 54 | meltRegion 55 | { 56 | level (0 0); 57 | regions 58 | { 59 | left 60 | { 61 | level (0 0); 62 | } 63 | out 64 | { 65 | level (0 0); 66 | } 67 | right 68 | { 69 | level (0 0); 70 | } 71 | } 72 | } 73 | } 74 | 75 | refinementRegions 76 | { 77 | } 78 | 79 | locationInMesh (0.01 0 0); 80 | maxLocalCells 10000000; 81 | maxGlobalCells 20000000; 82 | minRefinementCells 0; 83 | nCellsBetweenLevels 1; 84 | resolveFeatureAngle 30.0; 85 | allowFreeStandingZoneFaces true; 86 | planarAngle 30.0; 87 | maxLoadUnbalance 0.1; 88 | } 89 | 90 | snapControls 91 | { 92 | nSolveIter 50; 93 | nSmoothPatch 5; 94 | tolerance 1.0; 95 | nRelaxIter 10; 96 | nFeatureSnapIter 20; 97 | implicitFeatureSnap false; 98 | explicitFeatureSnap true; 99 | multiRegionFeatureSnap false; 100 | } 101 | 102 | addLayersControls 103 | { 104 | layers 105 | { 106 | } 107 | 108 | relativeSizes true; 109 | expansionRatio 1.0; 110 | finalLayerThickness 0.3; 111 | minThickness 0.25; 112 | nGrow 0; 113 | featureAngle 130.0; 114 | slipFeatureAngle 30.0; 115 | nRelaxIter 5; 116 | nSmoothSurfaceNormals 1; 117 | nSmoothNormals 3; 118 | nSmoothThickness 10; 119 | maxFaceThicknessRatio 0.5; 120 | maxThicknessToMedialRatio 0.3; 121 | minMedialAxisAngle 90; 122 | nBufferCellsNoExtrude 0; 123 | nLayerIter 50; 124 | nRelaxedIter 20; 125 | } 126 | 127 | meshQualityControls 128 | { 129 | maxNonOrtho 65.0; 130 | maxBoundarySkewness 20; 131 | maxInternalSkewness 4; 132 | maxConcave 80.0; 133 | minFlatness 0.5; 134 | minVol 1.0E-13; 135 | minTetQuality 1.0E-15; 136 | minArea -1.0; 137 | minTwist 0.02; 138 | minDeterminant 0.001; 139 | minFaceWeight 0.05; 140 | minVolRatio 0.01; 141 | minTriangleTwist -1.0; 142 | nSmoothScale 4; 143 | errorReduction 0.75; 144 | } 145 | 146 | debug 0; 147 | mergeTolerance 1.0E-6; 148 | autoBlockMesh true; 149 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/templatedir/case.template: -------------------------------------------------------------------------------- 1 | Header 2 | CHECK KEYWORDS Warn 3 | Mesh DB "." "meshElmer" 4 | Include Path "" 5 | Results Directory "" 6 | End 7 | 8 | $frequency = @f@ 9 | $current = @I@ 10 | 11 | Simulation 12 | Max Output Level = 3 13 | Coordinate System = Axi Symmetric 14 | Angular Frequency = $2*pi*@f@ 15 | 16 | Simulation Type = Steady 17 | Steady State Max Iterations = 10000000 18 | Steady State Min Iterations = 10000000 19 | 20 | Solver Input File = case.sif 21 | Post File = case.vtu 22 | Output Intervals = 2 23 | End 24 | 25 | Constants 26 | Permittivity of Vacuum = 8.8542e-12 27 | Permeability of Vacuum = 1.257e-6 28 | End 29 | 30 | Body 1 31 | Target Bodies(1) = 3 32 | Name = "inductorsTop" 33 | Equation = 1 34 | Material = 1 35 | Body Force = 1 36 | End 37 | 38 | Body 2 39 | Target Bodies(1) = 4 40 | Name = "inductorsBtm" 41 | Equation = 1 42 | Material = 1 43 | Body Force = 2 44 | End 45 | 46 | Body 3 47 | Target Bodies(1) = 1 48 | Name = "melt" 49 | Equation = 1 50 | Material = 3 51 | End 52 | 53 | Body 4 54 | Target Bodies(1) = 2 55 | Name = "air" 56 | Equation = 1 57 | Material = 2 58 | End 59 | 60 | Solver 1 61 | Equation = OpenFOAM2Elmer 62 | Procedure = "OpenFOAM2Elmer" "OpenFOAM2ElmerSolver" 63 | 64 | Target Variable 1 = String "OF conductivity" 65 | End 66 | 67 | Solver 2 68 | Equation = "MDynamics" 69 | Variable = AV[AV re:1 AV im:1] 70 | Procedure = "MagnetoDynamics2D" "MagnetoDynamics2DHarmonic" 71 | 72 | Linear System Symmetric = True 73 | 74 | Linear System Solver = "Iterative" 75 | Linear System Preconditioning = ILU0 76 | Linear System Residual Output = 10 77 | Linear System Max Iterations = 500 78 | Linear System Iterative Method = BiCGStabl 79 | Linear System Convergence Tolerance = 1.0e-8 80 | 81 | Steady State Convergence Tolerance = 1e-08 82 | Nonlinear System Max Iterations = 3 83 | Nonlinear System Convergence Tolerance = 1.0e-6 84 | Nonlinear System Relaxation Factor = 1 85 | 86 | Exported Variable 1 = -dg "OF conductivity" 87 | Exported Variable 1 Mask = String "material: mapping" 88 | End 89 | 90 | Solver 3 91 | Equation = "Postprocess" 92 | 93 | Potential Variable = String "AV" 94 | Procedure = "MagnetoDynamics" "MagnetoDynamicsCalcFields" 95 | 96 | Linear System Solver = Iterative 97 | Linear System Iterative Method = CG 98 | Linear System Symmetric = True 99 | Linear System Max Iterations = 100 100 | Linear System Preconditioning = None 101 | Linear System Convergence Tolerance = 1.0e-8 102 | Steady State Convergence Tolerance = 1e-08 103 | 104 | Nonlinear System Max Iterations = 3 105 | Nonlinear System Convergence Tolerance = 1.0e-8 106 | 107 | Calculate Current Density = Logical True 108 | Calculate Magnetic Field Strength = Logical True 109 | Calculate JxB = Logical True 110 | End 111 | 112 | Solver 4 113 | Equation = ElmerToOpenFOAM 114 | Procedure = "Elmer2OpenFOAM" "Elmer2OpenFOAMSolver" 115 | 116 | Target Variable 1 = jxb re 1 117 | Target Variable 2 = jxb re 2 118 | Target Variable 3 = jxb re 3 119 | End 120 | 121 | Equation 1 122 | Name = "MGDyn" 123 | Active Solvers(4) = 1 2 3 4 124 | End 125 | 126 | Material 1 127 | Name = "Copper" 128 | Electric Conductivity = 0 129 | Relative Permittivity = 1 130 | Relative Permeability = 1 131 | End 132 | 133 | Material 2 134 | Name = "Air" 135 | Electric Conductivity = 0 136 | Relative Permittivity = 1 137 | Relative Permeability = 1 138 | End 139 | 140 | Material 3 141 | Name = "Melt" 142 | Electric Conductivity = Equals "OF conductivity" 143 | Relative Permittivity = 1 144 | Relative Permeability = 1 145 | mapping = Logical True 146 | End 147 | 148 | Body Force 1 149 | Name = "Current pos" 150 | Current Density = Real @I@ 151 | End 152 | 153 | Body Force 2 154 | Name = "Current neg" 155 | Current Density = Real -@I@ 156 | End 157 | 158 | Boundary Condition 1 159 | Target Boundaries(1) = 5 160 | AV re {e} = Real 0.0 161 | AV im {e} = Real 0.0 162 | End 163 | 164 | Solver 2 165 | Reference Norm = Real 6.16456343E-04 166 | Reference Norm Tolerance = Real 1E-8 167 | End 168 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/case.sif: -------------------------------------------------------------------------------- 1 | Header 2 | CHECK KEYWORDS Warn 3 | Mesh DB "." "meshElmer" 4 | Include Path "" 5 | Results Directory "" 6 | End 7 | 8 | $frequency = 9800 9 | $current = 30000000 10 | 11 | Simulation 12 | Max Output Level = 3 13 | Coordinate System = Axi Symmetric 14 | Angular Frequency = $2*pi*9800 15 | 16 | Simulation Type = Steady 17 | Steady State Max Iterations = 10000000 18 | Steady State Min Iterations = 10000000 19 | 20 | Solver Input File = case.sif 21 | Post File = case.vtu 22 | Output Intervals = 2 23 | End 24 | 25 | Constants 26 | Permittivity of Vacuum = 8.8542e-12 27 | Permeability of Vacuum = 1.257e-6 28 | End 29 | 30 | Body 1 31 | Target Bodies(1) = 3 32 | Name = "inductorsTop" 33 | Equation = 1 34 | Material = 1 35 | Body Force = 1 36 | End 37 | 38 | Body 2 39 | Target Bodies(1) = 4 40 | Name = "inductorsBtm" 41 | Equation = 1 42 | Material = 1 43 | Body Force = 2 44 | End 45 | 46 | Body 3 47 | Target Bodies(1) = 1 48 | Name = "melt" 49 | Equation = 1 50 | Material = 3 51 | End 52 | 53 | Body 4 54 | Target Bodies(1) = 2 55 | Name = "air" 56 | Equation = 1 57 | Material = 2 58 | End 59 | 60 | Solver 1 61 | Equation = OpenFOAM2Elmer 62 | Procedure = "OpenFOAM2Elmer" "OpenFOAM2ElmerSolver" 63 | 64 | Target Variable 1 = String "OF conductivity" 65 | End 66 | 67 | Solver 2 68 | Equation = "MDynamics" 69 | Variable = AV[AV re:1 AV im:1] 70 | Procedure = "MagnetoDynamics2D" "MagnetoDynamics2DHarmonic" 71 | 72 | Linear System Symmetric = True 73 | 74 | Linear System Solver = "Iterative" 75 | Linear System Preconditioning = ILU0 76 | Linear System Residual Output = 10 77 | Linear System Max Iterations = 500 78 | Linear System Iterative Method = BiCGStabl 79 | Linear System Convergence Tolerance = 1.0e-8 80 | 81 | Steady State Convergence Tolerance = 1e-08 82 | Nonlinear System Max Iterations = 3 83 | Nonlinear System Convergence Tolerance = 1.0e-6 84 | Nonlinear System Relaxation Factor = 1 85 | 86 | Exported Variable 1 = -dg "OF conductivity" 87 | Exported Variable 1 Mask = String "material: mapping" 88 | End 89 | 90 | Solver 3 91 | Equation = "Postprocess" 92 | 93 | Potential Variable = String "AV" 94 | Procedure = "MagnetoDynamics" "MagnetoDynamicsCalcFields" 95 | 96 | Linear System Solver = Iterative 97 | Linear System Iterative Method = CG 98 | Linear System Symmetric = True 99 | Linear System Max Iterations = 100 100 | Linear System Preconditioning = None 101 | Linear System Convergence Tolerance = 1.0e-8 102 | Steady State Convergence Tolerance = 1e-08 103 | 104 | Nonlinear System Max Iterations = 3 105 | Nonlinear System Convergence Tolerance = 1.0e-8 106 | 107 | Calculate Current Density = Logical True 108 | Calculate Magnetic Field Strength = Logical True 109 | Calculate JxB = Logical True 110 | End 111 | 112 | Solver 4 113 | Equation = ElmerToOpenFOAM 114 | Procedure = "Elmer2OpenFOAM" "Elmer2OpenFOAMSolver" 115 | 116 | Target Variable 1 = jxb re 1 117 | Target Variable 2 = jxb re 2 118 | Target Variable 3 = jxb re 3 119 | End 120 | 121 | Equation 1 122 | Name = "MGDyn" 123 | Active Solvers(4) = 1 2 3 4 124 | End 125 | 126 | Material 1 127 | Name = "Copper" 128 | Electric Conductivity = 0 129 | Relative Permittivity = 1 130 | Relative Permeability = 1 131 | End 132 | 133 | Material 2 134 | Name = "Air" 135 | Electric Conductivity = 0 136 | Relative Permittivity = 1 137 | Relative Permeability = 1 138 | End 139 | 140 | Material 3 141 | Name = "Melt" 142 | Electric Conductivity = Equals "OF conductivity" 143 | Relative Permittivity = 1 144 | Relative Permeability = 1 145 | mapping = Logical True 146 | End 147 | 148 | Body Force 1 149 | Name = "Current pos" 150 | Current Density = Real 30000000 151 | End 152 | 153 | Body Force 2 154 | Name = "Current neg" 155 | Current Density = Real -30000000 156 | End 157 | 158 | Boundary Condition 1 159 | Target Boundaries(1) = 5 160 | AV re {e} = Real 0.0 161 | AV im {e} = Real 0.0 162 | End 163 | 164 | Solver 2 165 | Reference Norm = Real 6.16456343E-04 166 | Reference Norm Tolerance = Real 1E-8 167 | End 168 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/system/snappyHexMeshDict: -------------------------------------------------------------------------------- 1 | /*--------------------------------*- C++ -*----------------------------------*\ 2 | | o | | 3 | | o o | HELYX-OS | 4 | | o O o | Version: v2.3.1 | 5 | | o o | Web: http://www.engys.com | 6 | | o | | 7 | \*---------------------------------------------------------------------------*/ 8 | FoamFile 9 | { 10 | version 2.0; 11 | format ascii; 12 | class dictionary; 13 | location system; 14 | object snappyHexMeshDict; 15 | } 16 | 17 | castellatedMesh true; 18 | snap true; 19 | addLayers false; 20 | 21 | geometry 22 | { 23 | meltRegion.stl 24 | { 25 | type triSurfaceMesh; 26 | name meltRegion; 27 | 28 | regions 29 | { 30 | left 31 | { 32 | name leftBC; 33 | } 34 | out 35 | { 36 | name outBC; 37 | } 38 | right 39 | { 40 | name rightBC; 41 | } 42 | } 43 | } 44 | } 45 | 46 | castellatedMeshControls 47 | { 48 | features 49 | ( 50 | ); 51 | 52 | refinementSurfaces 53 | { 54 | meltRegion 55 | { 56 | level (0 0); 57 | regions 58 | { 59 | left 60 | { 61 | level (0 0); 62 | } 63 | out 64 | { 65 | level (0 0); 66 | } 67 | right 68 | { 69 | level (0 0); 70 | } 71 | } 72 | } 73 | } 74 | 75 | refinementRegions 76 | { 77 | } 78 | 79 | locationInMesh (0.01 0 0); 80 | maxLocalCells 10000000; 81 | maxGlobalCells 20000000; 82 | minRefinementCells 0; 83 | nCellsBetweenLevels 1; 84 | resolveFeatureAngle 30.0; 85 | allowFreeStandingZoneFaces true; 86 | planarAngle 30.0; 87 | maxLoadUnbalance 0.1; 88 | } 89 | 90 | snapControls 91 | { 92 | nSolveIter 50; 93 | nSmoothPatch 5; 94 | tolerance 1.0; 95 | nRelaxIter 10; 96 | nFeatureSnapIter 20; 97 | implicitFeatureSnap false; 98 | explicitFeatureSnap true; 99 | multiRegionFeatureSnap false; 100 | } 101 | 102 | addLayersControls 103 | { 104 | layers 105 | { 106 | } 107 | 108 | relativeSizes true; 109 | expansionRatio 1.0; 110 | finalLayerThickness 0.3; 111 | minThickness 0.25; 112 | nGrow 0; 113 | featureAngle 130.0; 114 | slipFeatureAngle 30.0; 115 | nRelaxIter 5; 116 | nSmoothSurfaceNormals 1; 117 | nSmoothNormals 3; 118 | nSmoothThickness 10; 119 | maxFaceThicknessRatio 0.5; 120 | maxThicknessToMedialRatio 0.3; 121 | minMedialAxisAngle 90; 122 | nBufferCellsNoExtrude 0; 123 | nLayerIter 50; 124 | nRelaxedIter 20; 125 | } 126 | 127 | meshQualityControls 128 | { 129 | maxNonOrtho 65.0; 130 | maxBoundarySkewness 20; 131 | maxInternalSkewness 4; 132 | maxConcave 80.0; 133 | minFlatness 0.5; 134 | minVol 1.0E-13; 135 | minTetQuality 1.0E-15; 136 | minArea -1.0; 137 | minTwist 0.02; 138 | minDeterminant 0.001; 139 | minFaceWeight 0.05; 140 | minVolRatio 0.01; 141 | minTriangleTwist -1.0; 142 | nSmoothScale 4; 143 | errorReduction 0.75; 144 | } 145 | 146 | debug 0; 147 | mergeTolerance 1.0E-6; 148 | autoBlockMesh true; 149 | -------------------------------------------------------------------------------- /Benchmark/casebase/meshElmer/mesh.boundary: -------------------------------------------------------------------------------- 1 | 1 6 1 0 202 5 4 2 | 2 6 2 0 202 6 5 3 | 3 6 3 0 202 7 6 4 | 4 6 5 0 202 8 7 5 | 5 6 7 0 202 9 8 6 | 6 6 9 0 202 10 9 7 | 7 6 11 0 202 11 10 8 | 8 6 13 0 202 12 11 9 | 9 6 15 0 202 13 12 10 | 10 6 17 0 202 14 13 11 | 11 6 19 0 202 15 14 12 | 12 6 21 0 202 16 15 13 | 13 6 23 0 202 17 16 14 | 14 6 25 0 202 18 17 15 | 15 6 27 0 202 19 18 16 | 16 6 29 0 202 20 19 17 | 17 6 31 0 202 21 20 18 | 18 6 33 0 202 22 21 19 | 19 6 35 0 202 23 22 20 | 20 6 37 0 202 24 23 21 | 21 6 39 0 202 25 24 22 | 22 6 41 0 202 26 25 23 | 23 6 43 0 202 27 26 24 | 24 6 45 0 202 28 27 25 | 25 6 47 0 202 29 28 26 | 26 6 49 0 202 30 29 27 | 27 6 51 0 202 31 30 28 | 28 6 53 0 202 32 31 29 | 29 6 55 0 202 33 32 30 | 30 6 57 0 202 34 33 31 | 31 6 59 0 202 35 34 32 | 32 6 61 0 202 36 35 33 | 33 6 63 0 202 37 36 34 | 34 6 65 0 202 38 37 35 | 35 6 67 0 202 39 38 36 | 36 6 69 0 202 40 39 37 | 37 6 71 0 202 41 40 38 | 38 6 73 0 202 42 41 39 | 39 6 75 0 202 43 42 40 | 40 6 77 0 202 44 43 41 | 41 6 79 0 202 45 44 42 | 42 6 81 0 202 46 45 43 | 43 6 83 0 202 47 46 44 | 44 6 85 0 202 48 47 45 | 45 6 87 0 202 49 48 46 | 46 6 89 0 202 50 49 47 | 47 6 91 0 202 51 50 48 | 48 6 93 0 202 52 51 49 | 49 6 95 0 202 53 52 50 | 50 6 97 0 202 54 53 51 | 51 6 99 0 202 55 54 52 | 52 6 101 0 202 56 55 53 | 53 6 103 0 202 57 56 54 | 54 6 105 0 202 58 57 55 | 55 6 107 0 202 1 58 56 | 56 6 146 0 202 2 59 57 | 57 6 108 0 202 59 60 58 | 58 6 163 0 202 60 61 59 | 59 6 109 0 202 61 62 60 | 60 6 161 0 202 62 63 61 | 61 6 110 0 202 63 64 62 | 62 6 166 0 202 64 65 63 | 63 6 111 0 202 65 66 64 | 64 6 167 0 202 66 67 65 | 65 6 112 0 202 67 68 66 | 66 6 168 0 202 68 69 67 | 67 6 113 0 202 69 70 68 | 68 6 169 0 202 70 71 69 | 69 6 114 0 202 71 72 70 | 70 6 170 0 202 72 73 71 | 71 6 115 0 202 73 74 72 | 72 6 171 0 202 74 75 73 | 73 6 116 0 202 75 76 74 | 74 6 172 0 202 76 77 75 | 75 6 117 0 202 77 78 76 | 76 6 173 0 202 78 79 77 | 77 6 118 0 202 79 80 78 | 78 6 174 0 202 80 81 79 | 79 6 119 0 202 81 82 80 | 80 6 176 0 202 82 83 81 | 81 6 120 0 202 83 84 82 | 82 6 107 0 202 84 1 83 | 83 6 121 0 202 3 85 84 | 84 6 122 0 202 85 86 85 | 85 6 180 0 202 86 87 86 | 86 6 123 0 202 87 88 87 | 87 6 182 0 202 88 89 88 | 88 6 124 0 202 89 90 89 | 89 6 183 0 202 90 91 90 | 90 6 125 0 202 91 92 91 | 91 6 184 0 202 92 93 92 | 92 6 126 0 202 93 94 93 | 93 6 185 0 202 94 95 94 | 94 6 127 0 202 95 96 95 | 95 6 186 0 202 96 97 96 | 96 6 128 0 202 97 98 97 | 97 6 187 0 202 98 99 98 | 98 6 129 0 202 99 100 99 | 99 6 188 0 202 100 101 100 | 100 6 130 0 202 101 102 101 | 101 6 189 0 202 102 103 102 | 102 6 131 0 202 103 104 103 | 103 6 190 0 202 104 105 104 | 104 6 132 0 202 105 106 105 | 105 6 191 0 202 106 107 106 | 106 6 133 0 202 107 108 107 | 107 6 192 0 202 108 109 108 | 108 6 134 0 202 109 110 109 | 109 6 193 0 202 110 111 110 | 110 6 135 0 202 111 112 111 | 111 6 194 0 202 112 113 112 | 112 6 136 0 202 113 114 113 | 113 6 195 0 202 114 115 114 | 114 6 137 0 202 115 116 115 | 115 6 196 0 202 116 117 116 | 116 6 138 0 202 117 118 117 | 117 6 197 0 202 118 119 118 | 118 6 139 0 202 119 120 119 | 119 6 198 0 202 120 121 120 | 120 6 140 0 202 121 122 121 | 121 6 199 0 202 122 123 122 | 122 6 141 0 202 123 124 123 | 123 6 200 0 202 124 125 124 | 124 6 142 0 202 125 126 125 | 125 6 201 0 202 126 127 126 | 126 6 143 0 202 127 128 127 | 127 6 144 0 202 128 2 128 | 128 6 1 0 202 4 129 129 | 129 6 147 0 202 129 130 130 | 130 6 203 0 202 130 131 131 | 131 6 148 0 202 131 132 132 | 132 6 204 0 202 132 133 133 | 133 6 149 0 202 133 134 134 | 134 6 205 0 202 134 135 135 | 135 6 150 0 202 135 136 136 | 136 6 206 0 202 136 137 137 | 137 6 151 0 202 137 138 138 | 138 6 207 0 202 138 139 139 | 139 6 152 0 202 139 140 140 | 140 6 208 0 202 140 141 141 | 141 6 153 0 202 141 142 142 | 142 6 209 0 202 142 143 143 | 143 6 154 0 202 143 144 144 | 144 6 210 0 202 144 145 145 | 145 6 155 0 202 145 146 146 | 146 6 211 0 202 146 147 147 | 147 6 156 0 202 147 148 148 | 148 6 212 0 202 148 149 149 | 149 6 157 0 202 149 150 150 | 150 6 213 0 202 150 151 151 | 151 6 158 0 202 151 152 152 | 152 6 121 0 202 152 3 153 | 153 6 2882 0 202 1524 1509 154 | 154 6 3174 0 202 1525 1524 155 | 155 6 2883 0 202 1526 1525 156 | 156 6 2884 0 202 1527 1526 157 | 157 6 2886 0 202 1528 1527 158 | 158 6 3176 0 202 1529 1528 159 | 159 6 2888 0 202 1530 1529 160 | 160 6 2889 0 202 1531 1530 161 | 161 6 2891 0 202 1532 1531 162 | 162 6 3179 0 202 1533 1532 163 | 163 6 2893 0 202 1534 1533 164 | 164 6 2894 0 202 1535 1534 165 | 165 6 2896 0 202 1536 1535 166 | 166 6 2881 0 202 4 1536 167 | 167 6 2900 0 202 1537 1 168 | 168 6 2902 0 202 1538 1537 169 | 169 6 3182 0 202 1539 1538 170 | 170 6 2904 0 202 1540 1539 171 | 171 6 2905 0 202 1541 1540 172 | 172 6 2907 0 202 1542 1541 173 | 173 6 3184 0 202 1543 1542 174 | 174 6 2909 0 202 1544 1543 175 | 175 6 2910 0 202 1545 1544 176 | 176 6 2912 0 202 1546 1545 177 | 177 6 2914 0 202 1508 1546 178 | 178 5 2882 0 202 1509 1547 179 | 179 5 2915 0 202 1547 1548 180 | 180 5 3187 0 202 1548 1549 181 | 181 5 2916 0 202 1549 1550 182 | 182 5 3188 0 202 1550 1551 183 | 183 5 2917 0 202 1551 1552 184 | 184 5 3189 0 202 1552 1553 185 | 185 5 2918 0 202 1553 1554 186 | 186 5 3190 0 202 1554 1555 187 | 187 5 2919 0 202 1555 1556 188 | 188 5 3191 0 202 1556 1557 189 | 189 5 2920 0 202 1557 1558 190 | 190 5 3192 0 202 1558 1559 191 | 191 5 2921 0 202 1559 1560 192 | 192 5 3193 0 202 1560 1561 193 | 193 5 2922 0 202 1561 1562 194 | 194 5 3194 0 202 1562 1563 195 | 195 5 2923 0 202 1563 1564 196 | 196 5 3195 0 202 1564 1565 197 | 197 5 2924 0 202 1565 1566 198 | 198 5 3196 0 202 1566 1567 199 | 199 5 2925 0 202 1567 1568 200 | 200 5 3198 0 202 1568 1569 201 | 201 5 2914 0 202 1569 1508 202 | 202 6 2926 0 202 1507 1570 203 | 203 6 2927 0 202 1570 1571 204 | 204 6 2928 0 202 1571 1572 205 | 205 6 2929 0 202 1572 1573 206 | 206 6 2930 0 202 1573 1574 207 | 207 6 2931 0 202 1574 1575 208 | 208 6 3099 0 202 1575 1576 209 | 209 6 2932 0 202 1576 1577 210 | 210 6 3207 0 202 1577 1578 211 | 211 6 2933 0 202 1578 1579 212 | 212 6 2934 0 202 1579 1580 213 | 213 6 2935 0 202 1580 1581 214 | 214 6 2936 0 202 1581 1582 215 | 215 6 3090 0 202 1582 1507 216 | 216 6 2937 0 202 1506 1583 217 | 217 6 2938 0 202 1583 1584 218 | 218 6 2939 0 202 1584 1585 219 | 219 6 2940 0 202 1585 1586 220 | 220 6 2941 0 202 1586 1587 221 | 221 6 2942 0 202 1587 1588 222 | 222 6 2943 0 202 1588 1589 223 | 223 6 2944 0 202 1589 1590 224 | 224 6 2945 0 202 1590 1591 225 | 225 6 2946 0 202 1591 1592 226 | 226 6 2947 0 202 1592 1593 227 | 227 6 2948 0 202 1593 1594 228 | 228 6 2949 0 202 1594 1595 229 | 229 6 3092 0 202 1595 1506 230 | 230 6 3285 0 202 1505 1596 231 | 231 6 2950 0 202 1596 1597 232 | 232 6 2951 0 202 1597 1598 233 | 233 6 2952 0 202 1598 1599 234 | 234 6 2953 0 202 1599 1600 235 | 235 6 2954 0 202 1600 1601 236 | 236 6 2955 0 202 1601 1602 237 | 237 6 2956 0 202 1602 1603 238 | 238 6 2957 0 202 1603 1604 239 | 239 6 2958 0 202 1604 1605 240 | 240 6 2959 0 202 1605 1606 241 | 241 6 2960 0 202 1606 1607 242 | 242 6 2961 0 202 1607 1505 243 | 243 6 2962 0 202 1504 1608 244 | 244 6 2963 0 202 1608 1609 245 | 245 6 2964 0 202 1609 1610 246 | 246 6 2965 0 202 1610 1611 247 | 247 6 2966 0 202 1611 1612 248 | 248 6 2967 0 202 1612 1613 249 | 249 6 2968 0 202 1613 1614 250 | 250 6 3534 0 202 1614 1615 251 | 251 6 2969 0 202 1615 1616 252 | 252 6 2970 0 202 1616 1617 253 | 253 6 2971 0 202 1617 1618 254 | 254 6 2972 0 202 1618 1619 255 | 255 6 2973 0 202 1619 1504 256 | 256 6 2974 0 202 1503 1620 257 | 257 6 2975 0 202 1620 1621 258 | 258 6 2976 0 202 1621 1622 259 | 259 6 2977 0 202 1622 1623 260 | 260 6 2978 0 202 1623 1624 261 | 261 6 3248 0 202 1624 1625 262 | 262 6 3091 0 202 1625 1626 263 | 263 6 2979 0 202 1626 1627 264 | 264 6 2980 0 202 1627 1628 265 | 265 6 2981 0 202 1628 1629 266 | 266 6 2982 0 202 1629 1630 267 | 267 6 2983 0 202 1630 1631 268 | 268 6 2984 0 202 1631 1632 269 | 269 6 3093 0 202 1632 1503 270 | 270 6 2985 0 202 1502 1633 271 | 271 6 2986 0 202 1633 1634 272 | 272 6 2987 0 202 1634 1635 273 | 273 6 2988 0 202 1635 1636 274 | 274 6 2989 0 202 1636 1637 275 | 275 6 2990 0 202 1637 1638 276 | 276 6 2991 0 202 1638 1639 277 | 277 6 2992 0 202 1639 1640 278 | 278 6 2993 0 202 1640 1641 279 | 279 6 2994 0 202 1641 1642 280 | 280 6 2995 0 202 1642 1643 281 | 281 6 2996 0 202 1643 1644 282 | 282 6 2997 0 202 1644 1645 283 | 283 6 2998 0 202 1645 1502 284 | 284 6 2999 0 202 1501 1646 285 | 285 6 3000 0 202 1646 1647 286 | 286 6 3001 0 202 1647 1648 287 | 287 6 3002 0 202 1648 1649 288 | 288 6 3003 0 202 1649 1650 289 | 289 6 3004 0 202 1650 1651 290 | 290 6 3273 0 202 1651 1652 291 | 291 6 3360 0 202 1652 1653 292 | 292 6 3005 0 202 1653 1654 293 | 293 6 3006 0 202 1654 1655 294 | 294 6 3007 0 202 1655 1656 295 | 295 6 3008 0 202 1656 1657 296 | 296 6 3009 0 202 1657 1658 297 | 297 6 3010 0 202 1658 1501 298 | 298 6 3283 0 202 1500 1659 299 | 299 6 3011 0 202 1659 1660 300 | 300 6 3012 0 202 1660 1661 301 | 301 6 3013 0 202 1661 1662 302 | 302 6 3014 0 202 1662 1663 303 | 303 6 3015 0 202 1663 1664 304 | 304 6 3016 0 202 1664 1665 305 | 305 6 3017 0 202 1665 1666 306 | 306 6 3018 0 202 1666 1667 307 | 307 6 3019 0 202 1667 1668 308 | 308 6 3020 0 202 1668 1669 309 | 309 6 3359 0 202 1669 1670 310 | 310 6 3274 0 202 1670 1500 311 | 311 6 3368 0 202 1499 1671 312 | 312 6 3021 0 202 1671 1672 313 | 313 6 3022 0 202 1672 1673 314 | 314 6 3023 0 202 1673 1674 315 | 315 6 3024 0 202 1674 1675 316 | 316 6 3025 0 202 1675 1676 317 | 317 6 3026 0 202 1676 1677 318 | 318 6 3027 0 202 1677 1678 319 | 319 6 3028 0 202 1678 1679 320 | 320 6 3029 0 202 1679 1680 321 | 321 6 3030 0 202 1680 1681 322 | 322 6 3031 0 202 1681 1682 323 | 323 6 3032 0 202 1682 1499 324 | 324 6 3033 0 202 1498 1683 325 | 325 6 3034 0 202 1683 1684 326 | 326 6 3035 0 202 1684 1685 327 | 327 6 3382 0 202 1685 1686 328 | 328 6 3036 0 202 1686 1687 329 | 329 6 3037 0 202 1687 1688 330 | 330 6 3038 0 202 1688 1689 331 | 331 6 3384 0 202 1689 1690 332 | 332 6 3039 0 202 1690 1691 333 | 333 6 3040 0 202 1691 1692 334 | 334 6 3041 0 202 1692 1693 335 | 335 6 3042 0 202 1693 1694 336 | 336 6 3043 0 202 1694 1498 337 | 337 6 3460 0 202 1497 1695 338 | 338 6 3044 0 202 1695 1696 339 | 339 6 3045 0 202 1696 1697 340 | 340 6 3046 0 202 1697 1698 341 | 341 6 3047 0 202 1698 1699 342 | 342 6 3048 0 202 1699 1700 343 | 343 6 3049 0 202 1700 1701 344 | 344 6 3394 0 202 1701 1702 345 | 345 6 3395 0 202 1702 1703 346 | 346 6 3050 0 202 1703 1704 347 | 347 6 3051 0 202 1704 1705 348 | 348 6 3052 0 202 1705 1706 349 | 349 6 3053 0 202 1706 1497 350 | 350 6 3401 0 202 1496 1707 351 | 351 6 3054 0 202 1707 1708 352 | 352 6 3055 0 202 1708 1709 353 | 353 6 3056 0 202 1709 1710 354 | 354 6 3057 0 202 1710 1711 355 | 355 6 3058 0 202 1711 1712 356 | 356 6 3059 0 202 1712 1713 357 | 357 6 3406 0 202 1713 1714 358 | 358 6 3060 0 202 1714 1715 359 | 359 6 3061 0 202 1715 1716 360 | 360 6 3062 0 202 1716 1717 361 | 361 6 3063 0 202 1717 1718 362 | 362 6 3411 0 202 1718 1496 363 | 363 6 3064 0 202 1495 1719 364 | 364 6 3065 0 202 1719 1720 365 | 365 6 3066 0 202 1720 1721 366 | 366 6 3067 0 202 1721 1722 367 | 367 6 3068 0 202 1722 1723 368 | 368 6 3069 0 202 1723 1724 369 | 369 6 3070 0 202 1724 1725 370 | 370 6 3286 0 202 1725 1726 371 | 371 6 3071 0 202 1726 1727 372 | 372 6 3072 0 202 1727 1728 373 | 373 6 3073 0 202 1728 1729 374 | 374 6 3074 0 202 1729 1730 375 | 375 6 3075 0 202 1730 1495 376 | 376 6 3076 0 202 1494 1731 377 | 377 6 3077 0 202 1731 1732 378 | 378 6 3078 0 202 1732 1733 379 | 379 6 3079 0 202 1733 1734 380 | 380 6 3080 0 202 1734 1735 381 | 381 6 3081 0 202 1735 1736 382 | 382 6 3082 0 202 1736 1737 383 | 383 6 3083 0 202 1737 1738 384 | 384 6 3084 0 202 1738 1739 385 | 385 6 3085 0 202 1739 1740 386 | 386 6 3086 0 202 1740 1741 387 | 387 6 3087 0 202 1741 1742 388 | 388 6 3088 0 202 1742 1743 389 | 389 6 3089 0 202 1743 1494 390 | 390 6 5181 0 202 1510 1744 391 | 391 6 5188 0 202 1744 1745 392 | 392 6 5182 0 202 1745 1746 393 | 393 6 5183 0 202 1746 1747 394 | 394 6 5190 0 202 1747 1748 395 | 395 6 5184 0 202 1748 1749 396 | 396 6 5179 0 202 1749 1750 397 | 397 6 5185 0 202 1750 1751 398 | 398 6 5180 0 202 1751 1510 399 | 399 6 5220 0 202 1511 1752 400 | 400 6 5232 0 202 1752 1753 401 | 401 6 5221 0 202 1753 1754 402 | 402 6 5230 0 202 1754 1755 403 | 403 6 5222 0 202 1755 1756 404 | 404 6 5223 0 202 1756 1757 405 | 405 6 5224 0 202 1757 1758 406 | 406 6 5225 0 202 1758 1759 407 | 407 6 5226 0 202 1759 1511 408 | 408 6 5257 0 202 1512 1760 409 | 409 6 5265 0 202 1760 1761 410 | 410 6 5258 0 202 1761 1762 411 | 411 6 5264 0 202 1762 1763 412 | 412 6 5259 0 202 1763 1764 413 | 413 6 5260 0 202 1764 1765 414 | 414 6 5261 0 202 1765 1766 415 | 415 6 5268 0 202 1766 1767 416 | 416 6 5262 0 202 1767 1512 417 | 417 6 5292 0 202 1513 1768 418 | 418 6 5304 0 202 1768 1769 419 | 419 6 5293 0 202 1769 1770 420 | 420 6 5294 0 202 1770 1771 421 | 421 6 5295 0 202 1771 1772 422 | 422 6 5296 0 202 1772 1773 423 | 423 6 5297 0 202 1773 1774 424 | 424 6 5298 0 202 1774 1775 425 | 425 6 5299 0 202 1775 1513 426 | 426 6 5328 0 202 1514 1776 427 | 427 6 5334 0 202 1776 1777 428 | 428 6 5329 0 202 1777 1778 429 | 429 6 5337 0 202 1778 1779 430 | 430 6 5330 0 202 1779 1780 431 | 431 6 5331 0 202 1780 1781 432 | 432 6 5332 0 202 1781 1782 433 | 433 6 5340 0 202 1782 1783 434 | 434 6 5333 0 202 1783 1514 435 | 435 6 5364 0 202 1515 1784 436 | 436 6 5373 0 202 1784 1785 437 | 437 6 5365 0 202 1785 1786 438 | 438 6 5358 0 202 1786 1787 439 | 439 6 5366 0 202 1787 1788 440 | 440 6 5367 0 202 1788 1789 441 | 441 6 5368 0 202 1789 1790 442 | 442 6 5371 0 202 1790 1791 443 | 443 6 5369 0 202 1791 1515 444 | 444 6 5400 0 202 1516 1792 445 | 445 6 5401 0 202 1792 1793 446 | 446 6 5402 0 202 1793 1794 447 | 447 6 5403 0 202 1794 1795 448 | 448 6 5404 0 202 1795 1796 449 | 449 6 5405 0 202 1796 1797 450 | 450 6 5406 0 202 1797 1798 451 | 451 6 5407 0 202 1798 1799 452 | 452 6 5408 0 202 1799 1516 453 | 453 6 5435 0 202 1517 1800 454 | 454 6 5436 0 202 1800 1801 455 | 455 6 5444 0 202 1801 1802 456 | 456 6 5437 0 202 1802 1803 457 | 457 6 5430 0 202 1803 1804 458 | 458 6 5438 0 202 1804 1805 459 | 459 6 5439 0 202 1805 1806 460 | 460 6 5440 0 202 1806 1807 461 | 461 6 5441 0 202 1807 1517 462 | 462 6 5470 0 202 1518 1808 463 | 463 6 5471 0 202 1808 1809 464 | 464 6 5478 0 202 1809 1810 465 | 465 6 5472 0 202 1810 1811 466 | 466 6 5482 0 202 1811 1812 467 | 467 6 5473 0 202 1812 1813 468 | 468 6 5477 0 202 1813 1814 469 | 469 6 5474 0 202 1814 1815 470 | 470 6 5475 0 202 1815 1518 471 | 471 6 5505 0 202 1519 1816 472 | 472 6 5506 0 202 1816 1817 473 | 473 6 5514 0 202 1817 1818 474 | 474 6 5507 0 202 1818 1819 475 | 475 6 5508 0 202 1819 1820 476 | 476 6 5509 0 202 1820 1821 477 | 477 6 5502 0 202 1821 1822 478 | 478 6 5510 0 202 1822 1823 479 | 479 6 5511 0 202 1823 1519 480 | 480 6 5539 0 202 1520 1824 481 | 481 6 5540 0 202 1824 1825 482 | 482 6 5541 0 202 1825 1826 483 | 483 6 5542 0 202 1826 1827 484 | 484 6 5535 0 202 1827 1828 485 | 485 6 5543 0 202 1828 1829 486 | 486 6 5537 0 202 1829 1830 487 | 487 6 5544 0 202 1830 1831 488 | 488 6 5545 0 202 1831 1520 489 | 489 6 5575 0 202 1521 1832 490 | 490 6 5576 0 202 1832 1833 491 | 491 6 5577 0 202 1833 1834 492 | 492 6 5578 0 202 1834 1835 493 | 493 6 5584 0 202 1835 1836 494 | 494 6 5579 0 202 1836 1837 495 | 495 6 5585 0 202 1837 1838 496 | 496 6 5580 0 202 1838 1839 497 | 497 6 5581 0 202 1839 1521 498 | 498 6 5610 0 202 1522 1840 499 | 499 6 5611 0 202 1840 1841 500 | 500 6 5621 0 202 1841 1842 501 | 501 6 5612 0 202 1842 1843 502 | 502 6 5613 0 202 1843 1844 503 | 503 6 5196 0 202 1844 1845 504 | 504 6 5193 0 202 1845 1846 505 | 505 6 5160 0 202 1846 1847 506 | 506 6 5158 0 202 1847 1522 507 | 507 6 4649 0 202 1523 1848 508 | 508 6 4607 0 202 1848 1849 509 | 509 6 4590 0 202 1849 1850 510 | 510 6 4585 0 202 1850 1851 511 | 511 6 4525 0 202 1851 1852 512 | 512 6 4515 0 202 1852 1853 513 | 513 6 4246 0 202 1853 1854 514 | 514 6 4095 0 202 1854 1855 515 | 515 6 4028 0 202 1855 1523 516 | -------------------------------------------------------------------------------- /Examples/UQ/casebase/meshElmer/mesh.boundary: -------------------------------------------------------------------------------- 1 | 1 6 1 0 202 5 4 2 | 2 6 2 0 202 6 5 3 | 3 6 3 0 202 7 6 4 | 4 6 5 0 202 8 7 5 | 5 6 7 0 202 9 8 6 | 6 6 9 0 202 10 9 7 | 7 6 11 0 202 11 10 8 | 8 6 13 0 202 12 11 9 | 9 6 15 0 202 13 12 10 | 10 6 17 0 202 14 13 11 | 11 6 19 0 202 15 14 12 | 12 6 21 0 202 16 15 13 | 13 6 23 0 202 17 16 14 | 14 6 25 0 202 18 17 15 | 15 6 27 0 202 19 18 16 | 16 6 29 0 202 20 19 17 | 17 6 31 0 202 21 20 18 | 18 6 33 0 202 22 21 19 | 19 6 35 0 202 23 22 20 | 20 6 37 0 202 24 23 21 | 21 6 39 0 202 25 24 22 | 22 6 41 0 202 26 25 23 | 23 6 43 0 202 27 26 24 | 24 6 45 0 202 28 27 25 | 25 6 47 0 202 29 28 26 | 26 6 49 0 202 30 29 27 | 27 6 51 0 202 31 30 28 | 28 6 53 0 202 32 31 29 | 29 6 55 0 202 33 32 30 | 30 6 57 0 202 34 33 31 | 31 6 59 0 202 35 34 32 | 32 6 61 0 202 36 35 33 | 33 6 63 0 202 37 36 34 | 34 6 65 0 202 38 37 35 | 35 6 67 0 202 39 38 36 | 36 6 69 0 202 40 39 37 | 37 6 71 0 202 41 40 38 | 38 6 73 0 202 42 41 39 | 39 6 75 0 202 43 42 40 | 40 6 77 0 202 44 43 41 | 41 6 79 0 202 45 44 42 | 42 6 81 0 202 46 45 43 | 43 6 83 0 202 47 46 44 | 44 6 85 0 202 48 47 45 | 45 6 87 0 202 49 48 46 | 46 6 89 0 202 50 49 47 | 47 6 91 0 202 51 50 48 | 48 6 93 0 202 52 51 49 | 49 6 95 0 202 53 52 50 | 50 6 97 0 202 54 53 51 | 51 6 99 0 202 55 54 52 | 52 6 101 0 202 56 55 53 | 53 6 103 0 202 57 56 54 | 54 6 105 0 202 58 57 55 | 55 6 107 0 202 1 58 56 | 56 6 146 0 202 2 59 57 | 57 6 108 0 202 59 60 58 | 58 6 163 0 202 60 61 59 | 59 6 109 0 202 61 62 60 | 60 6 161 0 202 62 63 61 | 61 6 110 0 202 63 64 62 | 62 6 166 0 202 64 65 63 | 63 6 111 0 202 65 66 64 | 64 6 167 0 202 66 67 65 | 65 6 112 0 202 67 68 66 | 66 6 168 0 202 68 69 67 | 67 6 113 0 202 69 70 68 | 68 6 169 0 202 70 71 69 | 69 6 114 0 202 71 72 70 | 70 6 170 0 202 72 73 71 | 71 6 115 0 202 73 74 72 | 72 6 171 0 202 74 75 73 | 73 6 116 0 202 75 76 74 | 74 6 172 0 202 76 77 75 | 75 6 117 0 202 77 78 76 | 76 6 173 0 202 78 79 77 | 77 6 118 0 202 79 80 78 | 78 6 174 0 202 80 81 79 | 79 6 119 0 202 81 82 80 | 80 6 176 0 202 82 83 81 | 81 6 120 0 202 83 84 82 | 82 6 107 0 202 84 1 83 | 83 6 121 0 202 3 85 84 | 84 6 122 0 202 85 86 85 | 85 6 180 0 202 86 87 86 | 86 6 123 0 202 87 88 87 | 87 6 182 0 202 88 89 88 | 88 6 124 0 202 89 90 89 | 89 6 183 0 202 90 91 90 | 90 6 125 0 202 91 92 91 | 91 6 184 0 202 92 93 92 | 92 6 126 0 202 93 94 93 | 93 6 185 0 202 94 95 94 | 94 6 127 0 202 95 96 95 | 95 6 186 0 202 96 97 96 | 96 6 128 0 202 97 98 97 | 97 6 187 0 202 98 99 98 | 98 6 129 0 202 99 100 99 | 99 6 188 0 202 100 101 100 | 100 6 130 0 202 101 102 101 | 101 6 189 0 202 102 103 102 | 102 6 131 0 202 103 104 103 | 103 6 190 0 202 104 105 104 | 104 6 132 0 202 105 106 105 | 105 6 191 0 202 106 107 106 | 106 6 133 0 202 107 108 107 | 107 6 192 0 202 108 109 108 | 108 6 134 0 202 109 110 109 | 109 6 193 0 202 110 111 110 | 110 6 135 0 202 111 112 111 | 111 6 194 0 202 112 113 112 | 112 6 136 0 202 113 114 113 | 113 6 195 0 202 114 115 114 | 114 6 137 0 202 115 116 115 | 115 6 196 0 202 116 117 116 | 116 6 138 0 202 117 118 117 | 117 6 197 0 202 118 119 118 | 118 6 139 0 202 119 120 119 | 119 6 198 0 202 120 121 120 | 120 6 140 0 202 121 122 121 | 121 6 199 0 202 122 123 122 | 122 6 141 0 202 123 124 123 | 123 6 200 0 202 124 125 124 | 124 6 142 0 202 125 126 125 | 125 6 201 0 202 126 127 126 | 126 6 143 0 202 127 128 127 | 127 6 144 0 202 128 2 128 | 128 6 1 0 202 4 129 129 | 129 6 147 0 202 129 130 130 | 130 6 203 0 202 130 131 131 | 131 6 148 0 202 131 132 132 | 132 6 204 0 202 132 133 133 | 133 6 149 0 202 133 134 134 | 134 6 205 0 202 134 135 135 | 135 6 150 0 202 135 136 136 | 136 6 206 0 202 136 137 137 | 137 6 151 0 202 137 138 138 | 138 6 207 0 202 138 139 139 | 139 6 152 0 202 139 140 140 | 140 6 208 0 202 140 141 141 | 141 6 153 0 202 141 142 142 | 142 6 209 0 202 142 143 143 | 143 6 154 0 202 143 144 144 | 144 6 210 0 202 144 145 145 | 145 6 155 0 202 145 146 146 | 146 6 211 0 202 146 147 147 | 147 6 156 0 202 147 148 148 | 148 6 212 0 202 148 149 149 | 149 6 157 0 202 149 150 150 | 150 6 213 0 202 150 151 151 | 151 6 158 0 202 151 152 152 | 152 6 121 0 202 152 3 153 | 153 6 2882 0 202 1524 1509 154 | 154 6 3174 0 202 1525 1524 155 | 155 6 2883 0 202 1526 1525 156 | 156 6 2884 0 202 1527 1526 157 | 157 6 2886 0 202 1528 1527 158 | 158 6 3176 0 202 1529 1528 159 | 159 6 2888 0 202 1530 1529 160 | 160 6 2889 0 202 1531 1530 161 | 161 6 2891 0 202 1532 1531 162 | 162 6 3179 0 202 1533 1532 163 | 163 6 2893 0 202 1534 1533 164 | 164 6 2894 0 202 1535 1534 165 | 165 6 2896 0 202 1536 1535 166 | 166 6 2881 0 202 4 1536 167 | 167 6 2900 0 202 1537 1 168 | 168 6 2902 0 202 1538 1537 169 | 169 6 3182 0 202 1539 1538 170 | 170 6 2904 0 202 1540 1539 171 | 171 6 2905 0 202 1541 1540 172 | 172 6 2907 0 202 1542 1541 173 | 173 6 3184 0 202 1543 1542 174 | 174 6 2909 0 202 1544 1543 175 | 175 6 2910 0 202 1545 1544 176 | 176 6 2912 0 202 1546 1545 177 | 177 6 2914 0 202 1508 1546 178 | 178 5 2882 0 202 1509 1547 179 | 179 5 2915 0 202 1547 1548 180 | 180 5 3187 0 202 1548 1549 181 | 181 5 2916 0 202 1549 1550 182 | 182 5 3188 0 202 1550 1551 183 | 183 5 2917 0 202 1551 1552 184 | 184 5 3189 0 202 1552 1553 185 | 185 5 2918 0 202 1553 1554 186 | 186 5 3190 0 202 1554 1555 187 | 187 5 2919 0 202 1555 1556 188 | 188 5 3191 0 202 1556 1557 189 | 189 5 2920 0 202 1557 1558 190 | 190 5 3192 0 202 1558 1559 191 | 191 5 2921 0 202 1559 1560 192 | 192 5 3193 0 202 1560 1561 193 | 193 5 2922 0 202 1561 1562 194 | 194 5 3194 0 202 1562 1563 195 | 195 5 2923 0 202 1563 1564 196 | 196 5 3195 0 202 1564 1565 197 | 197 5 2924 0 202 1565 1566 198 | 198 5 3196 0 202 1566 1567 199 | 199 5 2925 0 202 1567 1568 200 | 200 5 3198 0 202 1568 1569 201 | 201 5 2914 0 202 1569 1508 202 | 202 6 2926 0 202 1507 1570 203 | 203 6 2927 0 202 1570 1571 204 | 204 6 2928 0 202 1571 1572 205 | 205 6 2929 0 202 1572 1573 206 | 206 6 2930 0 202 1573 1574 207 | 207 6 2931 0 202 1574 1575 208 | 208 6 3099 0 202 1575 1576 209 | 209 6 2932 0 202 1576 1577 210 | 210 6 3207 0 202 1577 1578 211 | 211 6 2933 0 202 1578 1579 212 | 212 6 2934 0 202 1579 1580 213 | 213 6 2935 0 202 1580 1581 214 | 214 6 2936 0 202 1581 1582 215 | 215 6 3090 0 202 1582 1507 216 | 216 6 2937 0 202 1506 1583 217 | 217 6 2938 0 202 1583 1584 218 | 218 6 2939 0 202 1584 1585 219 | 219 6 2940 0 202 1585 1586 220 | 220 6 2941 0 202 1586 1587 221 | 221 6 2942 0 202 1587 1588 222 | 222 6 2943 0 202 1588 1589 223 | 223 6 2944 0 202 1589 1590 224 | 224 6 2945 0 202 1590 1591 225 | 225 6 2946 0 202 1591 1592 226 | 226 6 2947 0 202 1592 1593 227 | 227 6 2948 0 202 1593 1594 228 | 228 6 2949 0 202 1594 1595 229 | 229 6 3092 0 202 1595 1506 230 | 230 6 3285 0 202 1505 1596 231 | 231 6 2950 0 202 1596 1597 232 | 232 6 2951 0 202 1597 1598 233 | 233 6 2952 0 202 1598 1599 234 | 234 6 2953 0 202 1599 1600 235 | 235 6 2954 0 202 1600 1601 236 | 236 6 2955 0 202 1601 1602 237 | 237 6 2956 0 202 1602 1603 238 | 238 6 2957 0 202 1603 1604 239 | 239 6 2958 0 202 1604 1605 240 | 240 6 2959 0 202 1605 1606 241 | 241 6 2960 0 202 1606 1607 242 | 242 6 2961 0 202 1607 1505 243 | 243 6 2962 0 202 1504 1608 244 | 244 6 2963 0 202 1608 1609 245 | 245 6 2964 0 202 1609 1610 246 | 246 6 2965 0 202 1610 1611 247 | 247 6 2966 0 202 1611 1612 248 | 248 6 2967 0 202 1612 1613 249 | 249 6 2968 0 202 1613 1614 250 | 250 6 3534 0 202 1614 1615 251 | 251 6 2969 0 202 1615 1616 252 | 252 6 2970 0 202 1616 1617 253 | 253 6 2971 0 202 1617 1618 254 | 254 6 2972 0 202 1618 1619 255 | 255 6 2973 0 202 1619 1504 256 | 256 6 2974 0 202 1503 1620 257 | 257 6 2975 0 202 1620 1621 258 | 258 6 2976 0 202 1621 1622 259 | 259 6 2977 0 202 1622 1623 260 | 260 6 2978 0 202 1623 1624 261 | 261 6 3248 0 202 1624 1625 262 | 262 6 3091 0 202 1625 1626 263 | 263 6 2979 0 202 1626 1627 264 | 264 6 2980 0 202 1627 1628 265 | 265 6 2981 0 202 1628 1629 266 | 266 6 2982 0 202 1629 1630 267 | 267 6 2983 0 202 1630 1631 268 | 268 6 2984 0 202 1631 1632 269 | 269 6 3093 0 202 1632 1503 270 | 270 6 2985 0 202 1502 1633 271 | 271 6 2986 0 202 1633 1634 272 | 272 6 2987 0 202 1634 1635 273 | 273 6 2988 0 202 1635 1636 274 | 274 6 2989 0 202 1636 1637 275 | 275 6 2990 0 202 1637 1638 276 | 276 6 2991 0 202 1638 1639 277 | 277 6 2992 0 202 1639 1640 278 | 278 6 2993 0 202 1640 1641 279 | 279 6 2994 0 202 1641 1642 280 | 280 6 2995 0 202 1642 1643 281 | 281 6 2996 0 202 1643 1644 282 | 282 6 2997 0 202 1644 1645 283 | 283 6 2998 0 202 1645 1502 284 | 284 6 2999 0 202 1501 1646 285 | 285 6 3000 0 202 1646 1647 286 | 286 6 3001 0 202 1647 1648 287 | 287 6 3002 0 202 1648 1649 288 | 288 6 3003 0 202 1649 1650 289 | 289 6 3004 0 202 1650 1651 290 | 290 6 3273 0 202 1651 1652 291 | 291 6 3360 0 202 1652 1653 292 | 292 6 3005 0 202 1653 1654 293 | 293 6 3006 0 202 1654 1655 294 | 294 6 3007 0 202 1655 1656 295 | 295 6 3008 0 202 1656 1657 296 | 296 6 3009 0 202 1657 1658 297 | 297 6 3010 0 202 1658 1501 298 | 298 6 3283 0 202 1500 1659 299 | 299 6 3011 0 202 1659 1660 300 | 300 6 3012 0 202 1660 1661 301 | 301 6 3013 0 202 1661 1662 302 | 302 6 3014 0 202 1662 1663 303 | 303 6 3015 0 202 1663 1664 304 | 304 6 3016 0 202 1664 1665 305 | 305 6 3017 0 202 1665 1666 306 | 306 6 3018 0 202 1666 1667 307 | 307 6 3019 0 202 1667 1668 308 | 308 6 3020 0 202 1668 1669 309 | 309 6 3359 0 202 1669 1670 310 | 310 6 3274 0 202 1670 1500 311 | 311 6 3368 0 202 1499 1671 312 | 312 6 3021 0 202 1671 1672 313 | 313 6 3022 0 202 1672 1673 314 | 314 6 3023 0 202 1673 1674 315 | 315 6 3024 0 202 1674 1675 316 | 316 6 3025 0 202 1675 1676 317 | 317 6 3026 0 202 1676 1677 318 | 318 6 3027 0 202 1677 1678 319 | 319 6 3028 0 202 1678 1679 320 | 320 6 3029 0 202 1679 1680 321 | 321 6 3030 0 202 1680 1681 322 | 322 6 3031 0 202 1681 1682 323 | 323 6 3032 0 202 1682 1499 324 | 324 6 3033 0 202 1498 1683 325 | 325 6 3034 0 202 1683 1684 326 | 326 6 3035 0 202 1684 1685 327 | 327 6 3382 0 202 1685 1686 328 | 328 6 3036 0 202 1686 1687 329 | 329 6 3037 0 202 1687 1688 330 | 330 6 3038 0 202 1688 1689 331 | 331 6 3384 0 202 1689 1690 332 | 332 6 3039 0 202 1690 1691 333 | 333 6 3040 0 202 1691 1692 334 | 334 6 3041 0 202 1692 1693 335 | 335 6 3042 0 202 1693 1694 336 | 336 6 3043 0 202 1694 1498 337 | 337 6 3460 0 202 1497 1695 338 | 338 6 3044 0 202 1695 1696 339 | 339 6 3045 0 202 1696 1697 340 | 340 6 3046 0 202 1697 1698 341 | 341 6 3047 0 202 1698 1699 342 | 342 6 3048 0 202 1699 1700 343 | 343 6 3049 0 202 1700 1701 344 | 344 6 3394 0 202 1701 1702 345 | 345 6 3395 0 202 1702 1703 346 | 346 6 3050 0 202 1703 1704 347 | 347 6 3051 0 202 1704 1705 348 | 348 6 3052 0 202 1705 1706 349 | 349 6 3053 0 202 1706 1497 350 | 350 6 3401 0 202 1496 1707 351 | 351 6 3054 0 202 1707 1708 352 | 352 6 3055 0 202 1708 1709 353 | 353 6 3056 0 202 1709 1710 354 | 354 6 3057 0 202 1710 1711 355 | 355 6 3058 0 202 1711 1712 356 | 356 6 3059 0 202 1712 1713 357 | 357 6 3406 0 202 1713 1714 358 | 358 6 3060 0 202 1714 1715 359 | 359 6 3061 0 202 1715 1716 360 | 360 6 3062 0 202 1716 1717 361 | 361 6 3063 0 202 1717 1718 362 | 362 6 3411 0 202 1718 1496 363 | 363 6 3064 0 202 1495 1719 364 | 364 6 3065 0 202 1719 1720 365 | 365 6 3066 0 202 1720 1721 366 | 366 6 3067 0 202 1721 1722 367 | 367 6 3068 0 202 1722 1723 368 | 368 6 3069 0 202 1723 1724 369 | 369 6 3070 0 202 1724 1725 370 | 370 6 3286 0 202 1725 1726 371 | 371 6 3071 0 202 1726 1727 372 | 372 6 3072 0 202 1727 1728 373 | 373 6 3073 0 202 1728 1729 374 | 374 6 3074 0 202 1729 1730 375 | 375 6 3075 0 202 1730 1495 376 | 376 6 3076 0 202 1494 1731 377 | 377 6 3077 0 202 1731 1732 378 | 378 6 3078 0 202 1732 1733 379 | 379 6 3079 0 202 1733 1734 380 | 380 6 3080 0 202 1734 1735 381 | 381 6 3081 0 202 1735 1736 382 | 382 6 3082 0 202 1736 1737 383 | 383 6 3083 0 202 1737 1738 384 | 384 6 3084 0 202 1738 1739 385 | 385 6 3085 0 202 1739 1740 386 | 386 6 3086 0 202 1740 1741 387 | 387 6 3087 0 202 1741 1742 388 | 388 6 3088 0 202 1742 1743 389 | 389 6 3089 0 202 1743 1494 390 | 390 6 5181 0 202 1510 1744 391 | 391 6 5188 0 202 1744 1745 392 | 392 6 5182 0 202 1745 1746 393 | 393 6 5183 0 202 1746 1747 394 | 394 6 5190 0 202 1747 1748 395 | 395 6 5184 0 202 1748 1749 396 | 396 6 5179 0 202 1749 1750 397 | 397 6 5185 0 202 1750 1751 398 | 398 6 5180 0 202 1751 1510 399 | 399 6 5220 0 202 1511 1752 400 | 400 6 5232 0 202 1752 1753 401 | 401 6 5221 0 202 1753 1754 402 | 402 6 5230 0 202 1754 1755 403 | 403 6 5222 0 202 1755 1756 404 | 404 6 5223 0 202 1756 1757 405 | 405 6 5224 0 202 1757 1758 406 | 406 6 5225 0 202 1758 1759 407 | 407 6 5226 0 202 1759 1511 408 | 408 6 5257 0 202 1512 1760 409 | 409 6 5265 0 202 1760 1761 410 | 410 6 5258 0 202 1761 1762 411 | 411 6 5264 0 202 1762 1763 412 | 412 6 5259 0 202 1763 1764 413 | 413 6 5260 0 202 1764 1765 414 | 414 6 5261 0 202 1765 1766 415 | 415 6 5268 0 202 1766 1767 416 | 416 6 5262 0 202 1767 1512 417 | 417 6 5292 0 202 1513 1768 418 | 418 6 5304 0 202 1768 1769 419 | 419 6 5293 0 202 1769 1770 420 | 420 6 5294 0 202 1770 1771 421 | 421 6 5295 0 202 1771 1772 422 | 422 6 5296 0 202 1772 1773 423 | 423 6 5297 0 202 1773 1774 424 | 424 6 5298 0 202 1774 1775 425 | 425 6 5299 0 202 1775 1513 426 | 426 6 5328 0 202 1514 1776 427 | 427 6 5334 0 202 1776 1777 428 | 428 6 5329 0 202 1777 1778 429 | 429 6 5337 0 202 1778 1779 430 | 430 6 5330 0 202 1779 1780 431 | 431 6 5331 0 202 1780 1781 432 | 432 6 5332 0 202 1781 1782 433 | 433 6 5340 0 202 1782 1783 434 | 434 6 5333 0 202 1783 1514 435 | 435 6 5364 0 202 1515 1784 436 | 436 6 5373 0 202 1784 1785 437 | 437 6 5365 0 202 1785 1786 438 | 438 6 5358 0 202 1786 1787 439 | 439 6 5366 0 202 1787 1788 440 | 440 6 5367 0 202 1788 1789 441 | 441 6 5368 0 202 1789 1790 442 | 442 6 5371 0 202 1790 1791 443 | 443 6 5369 0 202 1791 1515 444 | 444 6 5400 0 202 1516 1792 445 | 445 6 5401 0 202 1792 1793 446 | 446 6 5402 0 202 1793 1794 447 | 447 6 5403 0 202 1794 1795 448 | 448 6 5404 0 202 1795 1796 449 | 449 6 5405 0 202 1796 1797 450 | 450 6 5406 0 202 1797 1798 451 | 451 6 5407 0 202 1798 1799 452 | 452 6 5408 0 202 1799 1516 453 | 453 6 5435 0 202 1517 1800 454 | 454 6 5436 0 202 1800 1801 455 | 455 6 5444 0 202 1801 1802 456 | 456 6 5437 0 202 1802 1803 457 | 457 6 5430 0 202 1803 1804 458 | 458 6 5438 0 202 1804 1805 459 | 459 6 5439 0 202 1805 1806 460 | 460 6 5440 0 202 1806 1807 461 | 461 6 5441 0 202 1807 1517 462 | 462 6 5470 0 202 1518 1808 463 | 463 6 5471 0 202 1808 1809 464 | 464 6 5478 0 202 1809 1810 465 | 465 6 5472 0 202 1810 1811 466 | 466 6 5482 0 202 1811 1812 467 | 467 6 5473 0 202 1812 1813 468 | 468 6 5477 0 202 1813 1814 469 | 469 6 5474 0 202 1814 1815 470 | 470 6 5475 0 202 1815 1518 471 | 471 6 5505 0 202 1519 1816 472 | 472 6 5506 0 202 1816 1817 473 | 473 6 5514 0 202 1817 1818 474 | 474 6 5507 0 202 1818 1819 475 | 475 6 5508 0 202 1819 1820 476 | 476 6 5509 0 202 1820 1821 477 | 477 6 5502 0 202 1821 1822 478 | 478 6 5510 0 202 1822 1823 479 | 479 6 5511 0 202 1823 1519 480 | 480 6 5539 0 202 1520 1824 481 | 481 6 5540 0 202 1824 1825 482 | 482 6 5541 0 202 1825 1826 483 | 483 6 5542 0 202 1826 1827 484 | 484 6 5535 0 202 1827 1828 485 | 485 6 5543 0 202 1828 1829 486 | 486 6 5537 0 202 1829 1830 487 | 487 6 5544 0 202 1830 1831 488 | 488 6 5545 0 202 1831 1520 489 | 489 6 5575 0 202 1521 1832 490 | 490 6 5576 0 202 1832 1833 491 | 491 6 5577 0 202 1833 1834 492 | 492 6 5578 0 202 1834 1835 493 | 493 6 5584 0 202 1835 1836 494 | 494 6 5579 0 202 1836 1837 495 | 495 6 5585 0 202 1837 1838 496 | 496 6 5580 0 202 1838 1839 497 | 497 6 5581 0 202 1839 1521 498 | 498 6 5610 0 202 1522 1840 499 | 499 6 5611 0 202 1840 1841 500 | 500 6 5621 0 202 1841 1842 501 | 501 6 5612 0 202 1842 1843 502 | 502 6 5613 0 202 1843 1844 503 | 503 6 5196 0 202 1844 1845 504 | 504 6 5193 0 202 1845 1846 505 | 505 6 5160 0 202 1846 1847 506 | 506 6 5158 0 202 1847 1522 507 | 507 6 4649 0 202 1523 1848 508 | 508 6 4607 0 202 1848 1849 509 | 509 6 4590 0 202 1849 1850 510 | 510 6 4585 0 202 1850 1851 511 | 511 6 4525 0 202 1851 1852 512 | 512 6 4515 0 202 1852 1853 513 | 513 6 4246 0 202 1853 1854 514 | 514 6 4095 0 202 1854 1855 515 | 515 6 4028 0 202 1855 1523 516 | -------------------------------------------------------------------------------- /Examples/Surrogate Modeling/casebase/meshElmer/mesh.boundary: -------------------------------------------------------------------------------- 1 | 1 6 1 0 202 5 4 2 | 2 6 2 0 202 6 5 3 | 3 6 3 0 202 7 6 4 | 4 6 5 0 202 8 7 5 | 5 6 7 0 202 9 8 6 | 6 6 9 0 202 10 9 7 | 7 6 11 0 202 11 10 8 | 8 6 13 0 202 12 11 9 | 9 6 15 0 202 13 12 10 | 10 6 17 0 202 14 13 11 | 11 6 19 0 202 15 14 12 | 12 6 21 0 202 16 15 13 | 13 6 23 0 202 17 16 14 | 14 6 25 0 202 18 17 15 | 15 6 27 0 202 19 18 16 | 16 6 29 0 202 20 19 17 | 17 6 31 0 202 21 20 18 | 18 6 33 0 202 22 21 19 | 19 6 35 0 202 23 22 20 | 20 6 37 0 202 24 23 21 | 21 6 39 0 202 25 24 22 | 22 6 41 0 202 26 25 23 | 23 6 43 0 202 27 26 24 | 24 6 45 0 202 28 27 25 | 25 6 47 0 202 29 28 26 | 26 6 49 0 202 30 29 27 | 27 6 51 0 202 31 30 28 | 28 6 53 0 202 32 31 29 | 29 6 55 0 202 33 32 30 | 30 6 57 0 202 34 33 31 | 31 6 59 0 202 35 34 32 | 32 6 61 0 202 36 35 33 | 33 6 63 0 202 37 36 34 | 34 6 65 0 202 38 37 35 | 35 6 67 0 202 39 38 36 | 36 6 69 0 202 40 39 37 | 37 6 71 0 202 41 40 38 | 38 6 73 0 202 42 41 39 | 39 6 75 0 202 43 42 40 | 40 6 77 0 202 44 43 41 | 41 6 79 0 202 45 44 42 | 42 6 81 0 202 46 45 43 | 43 6 83 0 202 47 46 44 | 44 6 85 0 202 48 47 45 | 45 6 87 0 202 49 48 46 | 46 6 89 0 202 50 49 47 | 47 6 91 0 202 51 50 48 | 48 6 93 0 202 52 51 49 | 49 6 95 0 202 53 52 50 | 50 6 97 0 202 54 53 51 | 51 6 99 0 202 55 54 52 | 52 6 101 0 202 56 55 53 | 53 6 103 0 202 57 56 54 | 54 6 105 0 202 58 57 55 | 55 6 107 0 202 1 58 56 | 56 6 146 0 202 2 59 57 | 57 6 108 0 202 59 60 58 | 58 6 163 0 202 60 61 59 | 59 6 109 0 202 61 62 60 | 60 6 161 0 202 62 63 61 | 61 6 110 0 202 63 64 62 | 62 6 166 0 202 64 65 63 | 63 6 111 0 202 65 66 64 | 64 6 167 0 202 66 67 65 | 65 6 112 0 202 67 68 66 | 66 6 168 0 202 68 69 67 | 67 6 113 0 202 69 70 68 | 68 6 169 0 202 70 71 69 | 69 6 114 0 202 71 72 70 | 70 6 170 0 202 72 73 71 | 71 6 115 0 202 73 74 72 | 72 6 171 0 202 74 75 73 | 73 6 116 0 202 75 76 74 | 74 6 172 0 202 76 77 75 | 75 6 117 0 202 77 78 76 | 76 6 173 0 202 78 79 77 | 77 6 118 0 202 79 80 78 | 78 6 174 0 202 80 81 79 | 79 6 119 0 202 81 82 80 | 80 6 176 0 202 82 83 81 | 81 6 120 0 202 83 84 82 | 82 6 107 0 202 84 1 83 | 83 6 121 0 202 3 85 84 | 84 6 122 0 202 85 86 85 | 85 6 180 0 202 86 87 86 | 86 6 123 0 202 87 88 87 | 87 6 182 0 202 88 89 88 | 88 6 124 0 202 89 90 89 | 89 6 183 0 202 90 91 90 | 90 6 125 0 202 91 92 91 | 91 6 184 0 202 92 93 92 | 92 6 126 0 202 93 94 93 | 93 6 185 0 202 94 95 94 | 94 6 127 0 202 95 96 95 | 95 6 186 0 202 96 97 96 | 96 6 128 0 202 97 98 97 | 97 6 187 0 202 98 99 98 | 98 6 129 0 202 99 100 99 | 99 6 188 0 202 100 101 100 | 100 6 130 0 202 101 102 101 | 101 6 189 0 202 102 103 102 | 102 6 131 0 202 103 104 103 | 103 6 190 0 202 104 105 104 | 104 6 132 0 202 105 106 105 | 105 6 191 0 202 106 107 106 | 106 6 133 0 202 107 108 107 | 107 6 192 0 202 108 109 108 | 108 6 134 0 202 109 110 109 | 109 6 193 0 202 110 111 110 | 110 6 135 0 202 111 112 111 | 111 6 194 0 202 112 113 112 | 112 6 136 0 202 113 114 113 | 113 6 195 0 202 114 115 114 | 114 6 137 0 202 115 116 115 | 115 6 196 0 202 116 117 116 | 116 6 138 0 202 117 118 117 | 117 6 197 0 202 118 119 118 | 118 6 139 0 202 119 120 119 | 119 6 198 0 202 120 121 120 | 120 6 140 0 202 121 122 121 | 121 6 199 0 202 122 123 122 | 122 6 141 0 202 123 124 123 | 123 6 200 0 202 124 125 124 | 124 6 142 0 202 125 126 125 | 125 6 201 0 202 126 127 126 | 126 6 143 0 202 127 128 127 | 127 6 144 0 202 128 2 128 | 128 6 1 0 202 4 129 129 | 129 6 147 0 202 129 130 130 | 130 6 203 0 202 130 131 131 | 131 6 148 0 202 131 132 132 | 132 6 204 0 202 132 133 133 | 133 6 149 0 202 133 134 134 | 134 6 205 0 202 134 135 135 | 135 6 150 0 202 135 136 136 | 136 6 206 0 202 136 137 137 | 137 6 151 0 202 137 138 138 | 138 6 207 0 202 138 139 139 | 139 6 152 0 202 139 140 140 | 140 6 208 0 202 140 141 141 | 141 6 153 0 202 141 142 142 | 142 6 209 0 202 142 143 143 | 143 6 154 0 202 143 144 144 | 144 6 210 0 202 144 145 145 | 145 6 155 0 202 145 146 146 | 146 6 211 0 202 146 147 147 | 147 6 156 0 202 147 148 148 | 148 6 212 0 202 148 149 149 | 149 6 157 0 202 149 150 150 | 150 6 213 0 202 150 151 151 | 151 6 158 0 202 151 152 152 | 152 6 121 0 202 152 3 153 | 153 6 2882 0 202 1524 1509 154 | 154 6 3174 0 202 1525 1524 155 | 155 6 2883 0 202 1526 1525 156 | 156 6 2884 0 202 1527 1526 157 | 157 6 2886 0 202 1528 1527 158 | 158 6 3176 0 202 1529 1528 159 | 159 6 2888 0 202 1530 1529 160 | 160 6 2889 0 202 1531 1530 161 | 161 6 2891 0 202 1532 1531 162 | 162 6 3179 0 202 1533 1532 163 | 163 6 2893 0 202 1534 1533 164 | 164 6 2894 0 202 1535 1534 165 | 165 6 2896 0 202 1536 1535 166 | 166 6 2881 0 202 4 1536 167 | 167 6 2900 0 202 1537 1 168 | 168 6 2902 0 202 1538 1537 169 | 169 6 3182 0 202 1539 1538 170 | 170 6 2904 0 202 1540 1539 171 | 171 6 2905 0 202 1541 1540 172 | 172 6 2907 0 202 1542 1541 173 | 173 6 3184 0 202 1543 1542 174 | 174 6 2909 0 202 1544 1543 175 | 175 6 2910 0 202 1545 1544 176 | 176 6 2912 0 202 1546 1545 177 | 177 6 2914 0 202 1508 1546 178 | 178 5 2882 0 202 1509 1547 179 | 179 5 2915 0 202 1547 1548 180 | 180 5 3187 0 202 1548 1549 181 | 181 5 2916 0 202 1549 1550 182 | 182 5 3188 0 202 1550 1551 183 | 183 5 2917 0 202 1551 1552 184 | 184 5 3189 0 202 1552 1553 185 | 185 5 2918 0 202 1553 1554 186 | 186 5 3190 0 202 1554 1555 187 | 187 5 2919 0 202 1555 1556 188 | 188 5 3191 0 202 1556 1557 189 | 189 5 2920 0 202 1557 1558 190 | 190 5 3192 0 202 1558 1559 191 | 191 5 2921 0 202 1559 1560 192 | 192 5 3193 0 202 1560 1561 193 | 193 5 2922 0 202 1561 1562 194 | 194 5 3194 0 202 1562 1563 195 | 195 5 2923 0 202 1563 1564 196 | 196 5 3195 0 202 1564 1565 197 | 197 5 2924 0 202 1565 1566 198 | 198 5 3196 0 202 1566 1567 199 | 199 5 2925 0 202 1567 1568 200 | 200 5 3198 0 202 1568 1569 201 | 201 5 2914 0 202 1569 1508 202 | 202 6 2926 0 202 1507 1570 203 | 203 6 2927 0 202 1570 1571 204 | 204 6 2928 0 202 1571 1572 205 | 205 6 2929 0 202 1572 1573 206 | 206 6 2930 0 202 1573 1574 207 | 207 6 2931 0 202 1574 1575 208 | 208 6 3099 0 202 1575 1576 209 | 209 6 2932 0 202 1576 1577 210 | 210 6 3207 0 202 1577 1578 211 | 211 6 2933 0 202 1578 1579 212 | 212 6 2934 0 202 1579 1580 213 | 213 6 2935 0 202 1580 1581 214 | 214 6 2936 0 202 1581 1582 215 | 215 6 3090 0 202 1582 1507 216 | 216 6 2937 0 202 1506 1583 217 | 217 6 2938 0 202 1583 1584 218 | 218 6 2939 0 202 1584 1585 219 | 219 6 2940 0 202 1585 1586 220 | 220 6 2941 0 202 1586 1587 221 | 221 6 2942 0 202 1587 1588 222 | 222 6 2943 0 202 1588 1589 223 | 223 6 2944 0 202 1589 1590 224 | 224 6 2945 0 202 1590 1591 225 | 225 6 2946 0 202 1591 1592 226 | 226 6 2947 0 202 1592 1593 227 | 227 6 2948 0 202 1593 1594 228 | 228 6 2949 0 202 1594 1595 229 | 229 6 3092 0 202 1595 1506 230 | 230 6 3285 0 202 1505 1596 231 | 231 6 2950 0 202 1596 1597 232 | 232 6 2951 0 202 1597 1598 233 | 233 6 2952 0 202 1598 1599 234 | 234 6 2953 0 202 1599 1600 235 | 235 6 2954 0 202 1600 1601 236 | 236 6 2955 0 202 1601 1602 237 | 237 6 2956 0 202 1602 1603 238 | 238 6 2957 0 202 1603 1604 239 | 239 6 2958 0 202 1604 1605 240 | 240 6 2959 0 202 1605 1606 241 | 241 6 2960 0 202 1606 1607 242 | 242 6 2961 0 202 1607 1505 243 | 243 6 2962 0 202 1504 1608 244 | 244 6 2963 0 202 1608 1609 245 | 245 6 2964 0 202 1609 1610 246 | 246 6 2965 0 202 1610 1611 247 | 247 6 2966 0 202 1611 1612 248 | 248 6 2967 0 202 1612 1613 249 | 249 6 2968 0 202 1613 1614 250 | 250 6 3534 0 202 1614 1615 251 | 251 6 2969 0 202 1615 1616 252 | 252 6 2970 0 202 1616 1617 253 | 253 6 2971 0 202 1617 1618 254 | 254 6 2972 0 202 1618 1619 255 | 255 6 2973 0 202 1619 1504 256 | 256 6 2974 0 202 1503 1620 257 | 257 6 2975 0 202 1620 1621 258 | 258 6 2976 0 202 1621 1622 259 | 259 6 2977 0 202 1622 1623 260 | 260 6 2978 0 202 1623 1624 261 | 261 6 3248 0 202 1624 1625 262 | 262 6 3091 0 202 1625 1626 263 | 263 6 2979 0 202 1626 1627 264 | 264 6 2980 0 202 1627 1628 265 | 265 6 2981 0 202 1628 1629 266 | 266 6 2982 0 202 1629 1630 267 | 267 6 2983 0 202 1630 1631 268 | 268 6 2984 0 202 1631 1632 269 | 269 6 3093 0 202 1632 1503 270 | 270 6 2985 0 202 1502 1633 271 | 271 6 2986 0 202 1633 1634 272 | 272 6 2987 0 202 1634 1635 273 | 273 6 2988 0 202 1635 1636 274 | 274 6 2989 0 202 1636 1637 275 | 275 6 2990 0 202 1637 1638 276 | 276 6 2991 0 202 1638 1639 277 | 277 6 2992 0 202 1639 1640 278 | 278 6 2993 0 202 1640 1641 279 | 279 6 2994 0 202 1641 1642 280 | 280 6 2995 0 202 1642 1643 281 | 281 6 2996 0 202 1643 1644 282 | 282 6 2997 0 202 1644 1645 283 | 283 6 2998 0 202 1645 1502 284 | 284 6 2999 0 202 1501 1646 285 | 285 6 3000 0 202 1646 1647 286 | 286 6 3001 0 202 1647 1648 287 | 287 6 3002 0 202 1648 1649 288 | 288 6 3003 0 202 1649 1650 289 | 289 6 3004 0 202 1650 1651 290 | 290 6 3273 0 202 1651 1652 291 | 291 6 3360 0 202 1652 1653 292 | 292 6 3005 0 202 1653 1654 293 | 293 6 3006 0 202 1654 1655 294 | 294 6 3007 0 202 1655 1656 295 | 295 6 3008 0 202 1656 1657 296 | 296 6 3009 0 202 1657 1658 297 | 297 6 3010 0 202 1658 1501 298 | 298 6 3283 0 202 1500 1659 299 | 299 6 3011 0 202 1659 1660 300 | 300 6 3012 0 202 1660 1661 301 | 301 6 3013 0 202 1661 1662 302 | 302 6 3014 0 202 1662 1663 303 | 303 6 3015 0 202 1663 1664 304 | 304 6 3016 0 202 1664 1665 305 | 305 6 3017 0 202 1665 1666 306 | 306 6 3018 0 202 1666 1667 307 | 307 6 3019 0 202 1667 1668 308 | 308 6 3020 0 202 1668 1669 309 | 309 6 3359 0 202 1669 1670 310 | 310 6 3274 0 202 1670 1500 311 | 311 6 3368 0 202 1499 1671 312 | 312 6 3021 0 202 1671 1672 313 | 313 6 3022 0 202 1672 1673 314 | 314 6 3023 0 202 1673 1674 315 | 315 6 3024 0 202 1674 1675 316 | 316 6 3025 0 202 1675 1676 317 | 317 6 3026 0 202 1676 1677 318 | 318 6 3027 0 202 1677 1678 319 | 319 6 3028 0 202 1678 1679 320 | 320 6 3029 0 202 1679 1680 321 | 321 6 3030 0 202 1680 1681 322 | 322 6 3031 0 202 1681 1682 323 | 323 6 3032 0 202 1682 1499 324 | 324 6 3033 0 202 1498 1683 325 | 325 6 3034 0 202 1683 1684 326 | 326 6 3035 0 202 1684 1685 327 | 327 6 3382 0 202 1685 1686 328 | 328 6 3036 0 202 1686 1687 329 | 329 6 3037 0 202 1687 1688 330 | 330 6 3038 0 202 1688 1689 331 | 331 6 3384 0 202 1689 1690 332 | 332 6 3039 0 202 1690 1691 333 | 333 6 3040 0 202 1691 1692 334 | 334 6 3041 0 202 1692 1693 335 | 335 6 3042 0 202 1693 1694 336 | 336 6 3043 0 202 1694 1498 337 | 337 6 3460 0 202 1497 1695 338 | 338 6 3044 0 202 1695 1696 339 | 339 6 3045 0 202 1696 1697 340 | 340 6 3046 0 202 1697 1698 341 | 341 6 3047 0 202 1698 1699 342 | 342 6 3048 0 202 1699 1700 343 | 343 6 3049 0 202 1700 1701 344 | 344 6 3394 0 202 1701 1702 345 | 345 6 3395 0 202 1702 1703 346 | 346 6 3050 0 202 1703 1704 347 | 347 6 3051 0 202 1704 1705 348 | 348 6 3052 0 202 1705 1706 349 | 349 6 3053 0 202 1706 1497 350 | 350 6 3401 0 202 1496 1707 351 | 351 6 3054 0 202 1707 1708 352 | 352 6 3055 0 202 1708 1709 353 | 353 6 3056 0 202 1709 1710 354 | 354 6 3057 0 202 1710 1711 355 | 355 6 3058 0 202 1711 1712 356 | 356 6 3059 0 202 1712 1713 357 | 357 6 3406 0 202 1713 1714 358 | 358 6 3060 0 202 1714 1715 359 | 359 6 3061 0 202 1715 1716 360 | 360 6 3062 0 202 1716 1717 361 | 361 6 3063 0 202 1717 1718 362 | 362 6 3411 0 202 1718 1496 363 | 363 6 3064 0 202 1495 1719 364 | 364 6 3065 0 202 1719 1720 365 | 365 6 3066 0 202 1720 1721 366 | 366 6 3067 0 202 1721 1722 367 | 367 6 3068 0 202 1722 1723 368 | 368 6 3069 0 202 1723 1724 369 | 369 6 3070 0 202 1724 1725 370 | 370 6 3286 0 202 1725 1726 371 | 371 6 3071 0 202 1726 1727 372 | 372 6 3072 0 202 1727 1728 373 | 373 6 3073 0 202 1728 1729 374 | 374 6 3074 0 202 1729 1730 375 | 375 6 3075 0 202 1730 1495 376 | 376 6 3076 0 202 1494 1731 377 | 377 6 3077 0 202 1731 1732 378 | 378 6 3078 0 202 1732 1733 379 | 379 6 3079 0 202 1733 1734 380 | 380 6 3080 0 202 1734 1735 381 | 381 6 3081 0 202 1735 1736 382 | 382 6 3082 0 202 1736 1737 383 | 383 6 3083 0 202 1737 1738 384 | 384 6 3084 0 202 1738 1739 385 | 385 6 3085 0 202 1739 1740 386 | 386 6 3086 0 202 1740 1741 387 | 387 6 3087 0 202 1741 1742 388 | 388 6 3088 0 202 1742 1743 389 | 389 6 3089 0 202 1743 1494 390 | 390 6 5181 0 202 1510 1744 391 | 391 6 5188 0 202 1744 1745 392 | 392 6 5182 0 202 1745 1746 393 | 393 6 5183 0 202 1746 1747 394 | 394 6 5190 0 202 1747 1748 395 | 395 6 5184 0 202 1748 1749 396 | 396 6 5179 0 202 1749 1750 397 | 397 6 5185 0 202 1750 1751 398 | 398 6 5180 0 202 1751 1510 399 | 399 6 5220 0 202 1511 1752 400 | 400 6 5232 0 202 1752 1753 401 | 401 6 5221 0 202 1753 1754 402 | 402 6 5230 0 202 1754 1755 403 | 403 6 5222 0 202 1755 1756 404 | 404 6 5223 0 202 1756 1757 405 | 405 6 5224 0 202 1757 1758 406 | 406 6 5225 0 202 1758 1759 407 | 407 6 5226 0 202 1759 1511 408 | 408 6 5257 0 202 1512 1760 409 | 409 6 5265 0 202 1760 1761 410 | 410 6 5258 0 202 1761 1762 411 | 411 6 5264 0 202 1762 1763 412 | 412 6 5259 0 202 1763 1764 413 | 413 6 5260 0 202 1764 1765 414 | 414 6 5261 0 202 1765 1766 415 | 415 6 5268 0 202 1766 1767 416 | 416 6 5262 0 202 1767 1512 417 | 417 6 5292 0 202 1513 1768 418 | 418 6 5304 0 202 1768 1769 419 | 419 6 5293 0 202 1769 1770 420 | 420 6 5294 0 202 1770 1771 421 | 421 6 5295 0 202 1771 1772 422 | 422 6 5296 0 202 1772 1773 423 | 423 6 5297 0 202 1773 1774 424 | 424 6 5298 0 202 1774 1775 425 | 425 6 5299 0 202 1775 1513 426 | 426 6 5328 0 202 1514 1776 427 | 427 6 5334 0 202 1776 1777 428 | 428 6 5329 0 202 1777 1778 429 | 429 6 5337 0 202 1778 1779 430 | 430 6 5330 0 202 1779 1780 431 | 431 6 5331 0 202 1780 1781 432 | 432 6 5332 0 202 1781 1782 433 | 433 6 5340 0 202 1782 1783 434 | 434 6 5333 0 202 1783 1514 435 | 435 6 5364 0 202 1515 1784 436 | 436 6 5373 0 202 1784 1785 437 | 437 6 5365 0 202 1785 1786 438 | 438 6 5358 0 202 1786 1787 439 | 439 6 5366 0 202 1787 1788 440 | 440 6 5367 0 202 1788 1789 441 | 441 6 5368 0 202 1789 1790 442 | 442 6 5371 0 202 1790 1791 443 | 443 6 5369 0 202 1791 1515 444 | 444 6 5400 0 202 1516 1792 445 | 445 6 5401 0 202 1792 1793 446 | 446 6 5402 0 202 1793 1794 447 | 447 6 5403 0 202 1794 1795 448 | 448 6 5404 0 202 1795 1796 449 | 449 6 5405 0 202 1796 1797 450 | 450 6 5406 0 202 1797 1798 451 | 451 6 5407 0 202 1798 1799 452 | 452 6 5408 0 202 1799 1516 453 | 453 6 5435 0 202 1517 1800 454 | 454 6 5436 0 202 1800 1801 455 | 455 6 5444 0 202 1801 1802 456 | 456 6 5437 0 202 1802 1803 457 | 457 6 5430 0 202 1803 1804 458 | 458 6 5438 0 202 1804 1805 459 | 459 6 5439 0 202 1805 1806 460 | 460 6 5440 0 202 1806 1807 461 | 461 6 5441 0 202 1807 1517 462 | 462 6 5470 0 202 1518 1808 463 | 463 6 5471 0 202 1808 1809 464 | 464 6 5478 0 202 1809 1810 465 | 465 6 5472 0 202 1810 1811 466 | 466 6 5482 0 202 1811 1812 467 | 467 6 5473 0 202 1812 1813 468 | 468 6 5477 0 202 1813 1814 469 | 469 6 5474 0 202 1814 1815 470 | 470 6 5475 0 202 1815 1518 471 | 471 6 5505 0 202 1519 1816 472 | 472 6 5506 0 202 1816 1817 473 | 473 6 5514 0 202 1817 1818 474 | 474 6 5507 0 202 1818 1819 475 | 475 6 5508 0 202 1819 1820 476 | 476 6 5509 0 202 1820 1821 477 | 477 6 5502 0 202 1821 1822 478 | 478 6 5510 0 202 1822 1823 479 | 479 6 5511 0 202 1823 1519 480 | 480 6 5539 0 202 1520 1824 481 | 481 6 5540 0 202 1824 1825 482 | 482 6 5541 0 202 1825 1826 483 | 483 6 5542 0 202 1826 1827 484 | 484 6 5535 0 202 1827 1828 485 | 485 6 5543 0 202 1828 1829 486 | 486 6 5537 0 202 1829 1830 487 | 487 6 5544 0 202 1830 1831 488 | 488 6 5545 0 202 1831 1520 489 | 489 6 5575 0 202 1521 1832 490 | 490 6 5576 0 202 1832 1833 491 | 491 6 5577 0 202 1833 1834 492 | 492 6 5578 0 202 1834 1835 493 | 493 6 5584 0 202 1835 1836 494 | 494 6 5579 0 202 1836 1837 495 | 495 6 5585 0 202 1837 1838 496 | 496 6 5580 0 202 1838 1839 497 | 497 6 5581 0 202 1839 1521 498 | 498 6 5610 0 202 1522 1840 499 | 499 6 5611 0 202 1840 1841 500 | 500 6 5621 0 202 1841 1842 501 | 501 6 5612 0 202 1842 1843 502 | 502 6 5613 0 202 1843 1844 503 | 503 6 5196 0 202 1844 1845 504 | 504 6 5193 0 202 1845 1846 505 | 505 6 5160 0 202 1846 1847 506 | 506 6 5158 0 202 1847 1522 507 | 507 6 4649 0 202 1523 1848 508 | 508 6 4607 0 202 1848 1849 509 | 509 6 4590 0 202 1849 1850 510 | 510 6 4585 0 202 1850 1851 511 | 511 6 4525 0 202 1851 1852 512 | 512 6 4515 0 202 1852 1853 513 | 513 6 4246 0 202 1853 1854 514 | 514 6 4095 0 202 1854 1855 515 | 515 6 4028 0 202 1855 1523 516 | --------------------------------------------------------------------------------