CFDEMperiodicTest/parCFDDEMrun.sh at main · acatwithacomputer/CFDEMperiodicTest · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
#!/bin/bash

#===================================================================#
# allrun script for testcase as part of test routine
# run settlingTest CFD part
# Christoph Goniva - Feb. 2011
#===================================================================#

#- source CFDEM env vars
. ~/.bashrc

#- include functions
source $CFDEM_SRC_DIR/lagrangian/cfdemParticle/etc/functions.sh

#--------------------------------------------------------------------------------#
#- define variables
casePath="$(dirname "$(readlink -f ${BASH_SOURCE[0]})")"
logpath=$casePath
headerText="cfdemSolverIB_periodicTest"
logfileName="log_$headerText"
solverName="cfdemSolverIB"
nrProcs="4"
machineFileName="none"   # yourMachinefileName | none
debugMode="off"          # on | off| strict
testHarnessPath="$CFDEM_TEST_HARNESS_PATH"
runOctave="false"
postproc="true"
#--------------------------------------------------------------------------------#

#- call function to run a parallel CFD-DEM case
parCFDDEMrun $logpath $logfileName $casePath $headerText $solverName $nrProcs $machineFileName $debugMode

if [ $runOctave == "true" ]
  then

    cd $casePath/CFD/octave
    octave --no-gui postproc.m
    evince pos_y_two_part_rec_glow.eps
    evince vel_y_two_part_rec_glow.eps
    #display pos_y_two_part_rec_glow.png &
    #display vel_y_two_part_rec_glow.png &
fi

if [ $postproc == "true" ]
  then
    #- get VTK data from liggghts dump file
    cd $casePath/DEM/post
    python2 -i $CFDEM_LPP_DIR/lpp.py  dump.liggghts_run

    #- get VTK data from CFD sim
    cd $casePath/CFD
    for ((i = 0 ; i < $nrProcs ; i++));
    do
        cp -r ./processor${i}/constant/polyMesh/ ./processor${i}/polyMesh
    done
    reconstructParMesh -mergeTol 1e-06
    reconstructPar -noLagrangian
    foamToVTK

    #- start paraview
    #paraview

    #- keep terminal open (if started in new terminal)
    echo "...press enter to clean up case"
    echo "press Ctr+C to keep data"
    read
fi

#- copy log file to test harness
#cp ../../$logfileName $testHarnessPath

#- clean up case
#keepDEMrestart="false"
#cleanCFDEMcase $casePath/CFD $keepDEMrestart