From 60716a4d1d501e22ae626aa9c7a9af962513692c Mon Sep 17 00:00:00 2001 From: L-Nafaryus Date: Sat, 20 Feb 2021 14:59:59 +0500 Subject: [PATCH] Mod: try to use potentialFoam --- run.sh | 14 +++++++++++ src/genmesh.py | 2 +- src/prefoam.py | 13 ++++++++++ src/simple-cubic/0/U | 2 +- src/simple-cubic/0/p | 4 ++-- src/simple-cubic/constant/transportProperties | 2 +- .../constant/turbulenceProperties | 11 +-------- src/simple-cubic/mesh.py | 2 +- src/simple-cubic/system/controlDict | 2 +- src/simple-cubic/system/fvSchemes | 4 ++++ src/simple-cubic/system/fvSolution | 24 +++++++++++++++---- 11 files changed, 58 insertions(+), 22 deletions(-) create mode 100755 run.sh create mode 100644 src/prefoam.py diff --git a/run.sh b/run.sh new file mode 100755 index 0000000..d31994d --- /dev/null +++ b/run.sh @@ -0,0 +1,14 @@ +#!/usr/bin/env bash + +path=build/simple-cubic/0.1 + +# python src/genmesh.py +# python src/prefoam.py + +ideasUnvFoam -case $path mesh.unv +checkMesh -case $path +foamDictionary -case $path constant/polyMesh/boundary -entry entry0.wall.type -set wall +#potentialFoam -case $path +decomposePar -case $path +mpirun -np 4 --oversubscribe simpleFoam -parallel -case $path > $path/simpleFoam.log + diff --git a/src/genmesh.py b/src/genmesh.py index f67abcd..2314b8d 100644 --- a/src/genmesh.py +++ b/src/genmesh.py @@ -10,7 +10,7 @@ if not os.path.exists(build): ### -alpha = [0.1, 0.15, 0.2] +alpha = [0.1] #, 0.15, 0.2] simpleCubic = os.path.join(src, "simple-cubic/main.py") # Body-centered cubic diff --git a/src/prefoam.py b/src/prefoam.py new file mode 100644 index 0000000..7c27f05 --- /dev/null +++ b/src/prefoam.py @@ -0,0 +1,13 @@ +import os, shutil + +src = os.getcwd() +build = os.path.join(src, "../build") + +if not os.path.exists(build): + os.makedirs(build) + +foamCase = [ "0", "constant", "system" ] + +for d in foamCase: + shutil.copytree("{}/simple-cubic/{}".format(src, d), + "{}/simple-cubic/0.1/{}".format(build, d)) diff --git a/src/simple-cubic/0/U b/src/simple-cubic/0/U index 1279edb..cde1cf0 100644 --- a/src/simple-cubic/0/U +++ b/src/simple-cubic/0/U @@ -23,7 +23,7 @@ boundaryField inlet { type pressureInletVelocity; - valut uniform (0 0 0); + value uniform (0 0 0); } outlet diff --git a/src/simple-cubic/0/p b/src/simple-cubic/0/p index b7d670d..ab26b2b 100644 --- a/src/simple-cubic/0/p +++ b/src/simple-cubic/0/p @@ -16,14 +16,14 @@ FoamFile dimensions [0 2 -2 0 0 0 0]; -internaleField uniform 0; +internalField uniform 0; boundaryField { inlet { type fixedValue; - value uniform 1; + value uniform 0.001; } outlet diff --git a/src/simple-cubic/constant/transportProperties b/src/simple-cubic/constant/transportProperties index b514dfd..ef5fb7b 100644 --- a/src/simple-cubic/constant/transportProperties +++ b/src/simple-cubic/constant/transportProperties @@ -17,6 +17,6 @@ FoamFile transportModel Newtonian; -nu 1e-05; +nu 1e-06; // ************************************************************************* // diff --git a/src/simple-cubic/constant/turbulenceProperties b/src/simple-cubic/constant/turbulenceProperties index b6bbfab..9ec63ee 100644 --- a/src/simple-cubic/constant/turbulenceProperties +++ b/src/simple-cubic/constant/turbulenceProperties @@ -15,18 +15,9 @@ FoamFile } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // -simulationType RAS; +simulationType laminar; -RAS -{ - // Tested with kEpsilon, realizableKE, kOmega, kOmegaSST, - // ShihQuadraticKE, LienCubicKE. - RASModel kEpsilon; - turbulence on; - - printCoeffs on; -} // ************************************************************************* // diff --git a/src/simple-cubic/mesh.py b/src/simple-cubic/mesh.py index 31b31fc..afa666e 100644 --- a/src/simple-cubic/mesh.py +++ b/src/simple-cubic/mesh.py @@ -17,7 +17,7 @@ def create(geomObj, bc): param.SetCheckChartBoundary( 0 ) param.SetMinSize( 0.01 ) param.SetMaxSize( 0.1 ) - param.SetFineness( 3 ) + param.SetFineness( 4 ) #param.SetGrowthRate( 0.1 ) #param.SetNbSegPerEdge( 5 ) #param.SetNbSegPerRadius( 10 ) diff --git a/src/simple-cubic/system/controlDict b/src/simple-cubic/system/controlDict index f7623b2..ad712c8 100644 --- a/src/simple-cubic/system/controlDict +++ b/src/simple-cubic/system/controlDict @@ -23,7 +23,7 @@ startTime 0; stopAt endTime; -endTime 500; +endTime 5000; deltaT 1; diff --git a/src/simple-cubic/system/fvSchemes b/src/simple-cubic/system/fvSchemes index 6d94655..d02dcab 100644 --- a/src/simple-cubic/system/fvSchemes +++ b/src/simple-cubic/system/fvSchemes @@ -23,17 +23,21 @@ ddtSchemes gradSchemes { default Gauss linear; +// grad(Phi) Gauss linear; } divSchemes { default none; div(phi,U) bounded Gauss linearUpwind grad(U); + div((nuEff*dev2(T(grad(U))))) Gauss linear; + div(nonlinearStress) Gauss linear; } laplacianSchemes { default Gauss linear corrected; +// laplacian(1, Phi) Gauss linear corrected; } interpolationSchemes diff --git a/src/simple-cubic/system/fvSolution b/src/simple-cubic/system/fvSolution index 891704c..c8b25d1 100644 --- a/src/simple-cubic/system/fvSolution +++ b/src/simple-cubic/system/fvSolution @@ -33,16 +33,30 @@ solvers tolerance 1e-08; relTol 0.1; } + + Phi + { + solver smoothSolver; + smoother GaussSeidel; + nSweeps 2; + tolerance 1e-08; + relTol 0.1; + } +} + +potentialFlow +{ + nNonOrthogonalCorrectors 3; } SIMPLE { - nNonOrthogonalCorrectors 0; + nNonOrthogonalCorrectors 10; residualControl { - p 1e-5; - U 1e-5; + p 1e-3; + U 1e-3; } } @@ -50,11 +64,11 @@ relaxationFactors { fields { - p 0.3; + p 0.05; } equations { - U 0.7; + U 0.2; } }