#!/bin/sh
cd ${0%/*} || exit 1 # Run from this directory# Source tutorial run functions
. $WM_PROJECT_DIR/bin/tools/RunFunctionsrunApplication 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/transportPropertiesrunApplication -s $e $(getApplication)
# extract y+, U+
# note: both must be added to foamLog.db
runApplication -s $e foamLog log.boundaryFoam.$eif [ -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
firm -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
figraphName=”OF_vs_ANALYTICAL.eps”
echo “Creating graph of u+ vs y+ to $graphName”
gnuplot
実行年月:07/01/2021