boundaryWallFunctionsProfile

#!/bin/sh
cd ${0%/*} || exit 1 # Run from this directory

# Source tutorial run functions
. $WM_PROJECT_DIR/bin/tools/RunFunctions

runApplication blockMesh

exponents=”2 3 4 5 6 7 8″

echo “# yPlus vs uPlus” > yPlus_vs_uPlus

for e in $exponents
do
echo ” Setting nu to 1e-$e”

sed “s/XXX/$e/g” constant/transportProperties.template \r
> constant/transportProperties

runApplication `getApplication`

mv log.boundaryFoam log.boundaryFoam_$e

# extract y+, U+
# note: both must be added to foamLog.db
runApplication foamLog log.boundaryFoam_$e

mv log.foamLog log.foamLog_$e

if [ -e logs/yPlus_0 ]
then
yPlus=`awk < logs/yPlus_0 'END{print $2}'` uPlus=`awk < logs/uPlus_0 'END{print $2}'` echo "$yPlus $uPlus" >> yPlus_vs_uPlus
fi

rm -rf logs
done

# create validation plot
# test if gnuplot exists on the system
if ! which gnuplot > /dev/null 2>&1
then
echo “gnuplot not found – skipping graph creation” >&2
exit 1
fi

graphName=”OF_vs_ANAYTICAL.eps”

echo “Creating graph of u+ vs y+ to $graphName”
gnuplot

実行年月:12/20/2015