diff --git a/openfoam2012/continuemany.sh b/openfoam2012/continuemany.sh index fd5f796359be2fe1d70b610b2c63b06dac1c2b3a..e011b218e9fac04b55dd4dc28f7d30346462d7f5 100755 --- a/openfoam2012/continuemany.sh +++ b/openfoam2012/continuemany.sh @@ -9,7 +9,7 @@ echo "continue until ${endTime}" top="$(pwd)" for n in 1 2 do - for logu in -2 -1 0 1 + for logu in -1 do cd "${top}" new="freeCylinderFlow-N${n}-u${logu}" @@ -17,8 +17,8 @@ do then echo "continue ${new}" >> runmany.log cd "${new}" - sed -i 's#^end_time .*;$#end_time '"${end_Time}"';#' parameters - ./continue.sh $n + sed -i 's#^end_Time .*;$#end_Time '"${endTime}"';#' parameters + ./continue.sh else echo "could not continue ${new} because nonexistent" fi diff --git a/openfoam2012/freeCylinderFlow/continue.sh b/openfoam2012/freeCylinderFlow/continue.sh index df6e6dc9310daee0cae2d0575b8e3d6fbaf035bb..506164e80b67171a994fd649febb7cc3bd3fc9c7 100755 --- a/openfoam2012/freeCylinderFlow/continue.sh +++ b/openfoam2012/freeCylinderFlow/continue.sh @@ -1,9 +1,9 @@ #! /usr/bin/env sh set -e # f = force = do not error on nonexisting files -rm -rf 'postProcessing/patchAverage(name=cylinder,mag(grad(T)))' -rm -rf '0/grad(T)' -rm -rf '0/mag(grad(T))' +# rm -rf 'postProcessing/patchAverage(name=cylinder,mag(grad(T)))' +# rm -rf '0/grad(T)' +# rm -rf '0/mag(grad(T))' decomposePar -force foamJob -wait -p buoyantBoussinesqPimpleFoam diff --git a/openfoam2012/freeCylinderFlow/parameters b/openfoam2012/freeCylinderFlow/parameters index c3ca90ab391830e5e07f8300668405adb986ebc0..2f80399248d2d095ba4486dd4d4ef662df934496 100644 --- a/openfoam2012/freeCylinderFlow/parameters +++ b/openfoam2012/freeCylinderFlow/parameters @@ -1,8 +1,8 @@ -N 1; +N 4; inVelocity 5e-2; Tinf 300; Tw 400; Pr_t 1.0; thickness 1; radius 1; -end_Time 5; +end_Time 20; diff --git a/openfoam2012/freeCylinderFlow/system/blockMeshDict b/openfoam2012/freeCylinderFlow/system/blockMeshDict index 799a36298545ea9a80097f4ca9579c4245081b81..f66b17b4641505ff4ce96d05fdf91d07f0ec70de 100644 --- a/openfoam2012/freeCylinderFlow/system/blockMeshDict +++ b/openfoam2012/freeCylinderFlow/system/blockMeshDict @@ -54,8 +54,8 @@ ghi 18; // grading horizontally in front of square gho 30; // grading horizontall behind the square gv 16; // grading vertically below or above square // gvs, ghs not done -gr #eval{pow(1.1,$cr)}; // grading in radial direction in square -/* gr 3; */ +/* gr #eval{min(pow(1.1,$cr),45)}; // grading in radial direction in square */ +gr 45; vertices ( diff --git a/openfoam2012/freeCylinderFlow/system/controlDict b/openfoam2012/freeCylinderFlow/system/controlDict index f9af0036566a5913e48ee9ec1afff25179843373..7797a9ea33e5212b94e2fc41eacabcb1f530033c 100644 --- a/openfoam2012/freeCylinderFlow/system/controlDict +++ b/openfoam2012/freeCylinderFlow/system/controlDict @@ -28,15 +28,15 @@ stopAt endTime; /* endTime 1000; */ endTime $end_Time; // hoping that steady state or regular fluctuation is reached -deltaT 0.01; +deltaT 0.1; adjustTimeStep yes; maxCo 1; // according to some forum, this does not destroy stability but accuracy -maxDeltaT 0.1; // otherwise for low Reynold numbers, time step became > 1 s +maxDeltaT 0.2; // otherwise for low Reynold numbers, time step became > 1 s writeControl runTime; -writeInterval 0.1; +writeInterval 1; purgeWrite 0; diff --git a/openfoam2012/runmany.sh b/openfoam2012/runmany.sh index 90ec09cdd0e8941b0a96837925bd0e0fd007a783..e3ea6203b670ad15d682fd9ec6896dce16182ac5 100755 --- a/openfoam2012/runmany.sh +++ b/openfoam2012/runmany.sh @@ -1,13 +1,14 @@ #! /usr/bin/env sh set -e +rm -f freeCylinderFlow/scripts.log echo "" >> runmany.archive touch runmany.log cat runmany.log >> runmany.archive 2> /dev/null echo "" > runmany.log top="$(pwd)" -for n in 1 2 +for n in 1 #2 4 do - for logu in -2 -1 0 1 + for logu in 0 #-1 # 1 -2 do cd "${top}" new="freeCylinderFlow-N${n}-u${logu}"