boundaryWallFunctionsProfile

#!/bin/sh
cd “${0%/*}” || exit # Run from this directory
. ${WM_PROJECT_DIR:?}/bin/tools/RunFunctions # Tutorial run functions
#——————————————————————————

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 -s $e $(getApplication)

# extract y+, U+
# note: both must be added to foamLog.db
runApplication -s $e foamLog log.boundaryFoam.$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

if notTest “$@”
then
# create validation plot

# Require gnuplot
command -v gnuplot >/dev/null || {
echo “gnuplot not found – skipping graph creation” 1>&2
exit 1
}

graphName=”OF_vs_ANALYTICAL.eps”

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

実行年月：08/01/2020