boundaryWallFunctionsProfile

#!/bin/sh
cd ${0%/*} || exit 1 # 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
# Test if gnuplot exists on the system
command -v gnuplot >/dev/null 2>&1 || {
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

実行年月：06/19/2020