public

README.md

@@ -8,7 +8,11 @@ For more information, please refer to the following: (https://maziarraissi.githu
 
 ## Note
 
-This work is currently under review for publication and the codes and data will be publicly available as soon as the review process concludes.
+The required data (to be copied in the Data directory) and some Matlab scripts (to be copied in the Figures directory) for plotting purposes are provided in the following link:
+
+   - [Data and Figures](https://bit.ly/2NRB65U)
+
+In addition to the Data and Figures directories, the Results folder is currently empty and will be automatically populated after running the corresponding examples provided in the Source and Scripts directories.
 
 ## Citation
 

Scripts/Cylinder2D_flower_case_studies_1.sh

@@ -0,0 +1,35 @@
+#!/bin/bash
+
+FAIL=0
+
+CUDA_VISIBLE_DEVICES=0 python Cylinder2D_flower_systematic.py 201 15000 > Cylinder2D_flower_201_15000_stdout.txt &
+CUDA_VISIBLE_DEVICES=1 python Cylinder2D_flower_systematic.py 101 15000 > Cylinder2D_flower_101_15000_stdout.txt &
+CUDA_VISIBLE_DEVICES=2 python Cylinder2D_flower_systematic.py 51 15000 > Cylinder2D_flower_51_15000_stdout.txt &
+CUDA_VISIBLE_DEVICES=3 python Cylinder2D_flower_systematic.py 26 15000 > Cylinder2D_flower_26_15000_stdout.txt &
+CUDA_VISIBLE_DEVICES=4 python Cylinder2D_flower_systematic.py 201 10000 > Cylinder2D_flower_201_10000_stdout.txt &
+CUDA_VISIBLE_DEVICES=5 python Cylinder2D_flower_systematic.py 101 10000 > Cylinder2D_flower_101_10000_stdout.txt &
+CUDA_VISIBLE_DEVICES=6 python Cylinder2D_flower_systematic.py 51 10000 > Cylinder2D_flower_51_10000_stdout.txt &
+CUDA_VISIBLE_DEVICES=7 python Cylinder2D_flower_systematic.py 26 10000 > Cylinder2D_flower_26_10000_stdout.txt&
+CUDA_VISIBLE_DEVICES=8 python Cylinder2D_flower_systematic.py 201 5000 > Cylinder2D_flower_201_5000_stdout.txt &
+CUDA_VISIBLE_DEVICES=9 python Cylinder2D_flower_systematic.py 101 5000 > Cylinder2D_flower_101_5000_stdout.txt &
+CUDA_VISIBLE_DEVICES=10 python Cylinder2D_flower_systematic.py 51 5000 > Cylinder2D_flower_51_5000_stdout.txt &
+CUDA_VISIBLE_DEVICES=11 python Cylinder2D_flower_systematic.py 26 5000 > Cylinder2D_flower_26_5000_stdout.txt &
+CUDA_VISIBLE_DEVICES=12 python Cylinder2D_flower_systematic.py 201 2500 > Cylinder2D_flower_201_2500_stdout.txt &
+CUDA_VISIBLE_DEVICES=13 python Cylinder2D_flower_systematic.py 101 2500 > Cylinder2D_flower_101_2500_stdout.txt &
+CUDA_VISIBLE_DEVICES=14 python Cylinder2D_flower_systematic.py 51 2500 > Cylinder2D_flower_51_2500_stdout.txt &
+CUDA_VISIBLE_DEVICES=15 python Cylinder2D_flower_systematic.py 26 2500 > Cylinder2D_flower_26_2500_stdout.txt &
+
+for job in `jobs -p`
+do
+echo $job
+    wait $job || let "FAIL+=1"
+done
+
+echo $FAIL
+
+if [ "$FAIL" == "0" ];
+then
+echo "YAY!"
+else
+echo "FAIL! ($FAIL)"
+fi

Scripts/Cylinder2D_flower_case_studies_2.sh

@@ -0,0 +1,35 @@
+#!/bin/bash
+
+FAIL=0
+
+CUDA_VISIBLE_DEVICES=0 python Cylinder2D_flower_systematic.py 26 2500 > Cylinder2D_flower_26_2500_stdout.txt &
+CUDA_VISIBLE_DEVICES=1 python Cylinder2D_flower_systematic.py 13 2500 > Cylinder2D_flower_13_2500_stdout.txt &
+CUDA_VISIBLE_DEVICES=2 python Cylinder2D_flower_systematic.py 7 2500 > Cylinder2D_flower_7_2500_stdout.txt &
+CUDA_VISIBLE_DEVICES=3 python Cylinder2D_flower_systematic.py 3 2500 > Cylinder2D_flower_3_2500_stdout.txt &
+CUDA_VISIBLE_DEVICES=4 python Cylinder2D_flower_systematic.py 26 1500 > Cylinder2D_flower_26_1500_stdout.txt &
+CUDA_VISIBLE_DEVICES=5 python Cylinder2D_flower_systematic.py 13 1500 > Cylinder2D_flower_13_1500_stdout.txt &
+CUDA_VISIBLE_DEVICES=6 python Cylinder2D_flower_systematic.py 7 1500 > Cylinder2D_flower_7_1500_stdout.txt &
+CUDA_VISIBLE_DEVICES=7 python Cylinder2D_flower_systematic.py 3 1500 > Cylinder2D_flower_3_1500_stdout.txt &
+CUDA_VISIBLE_DEVICES=8 python Cylinder2D_flower_systematic.py 26 500 > Cylinder2D_flower_26_500_stdout.txt &
+CUDA_VISIBLE_DEVICES=9 python Cylinder2D_flower_systematic.py 13 500 > Cylinder2D_flower_13_500_stdout.txt &
+CUDA_VISIBLE_DEVICES=10 python Cylinder2D_flower_systematic.py 7 500 > Cylinder2D_flower_7_500_stdout.txt &
+CUDA_VISIBLE_DEVICES=11 python Cylinder2D_flower_systematic.py 3 500 > Cylinder2D_flower_3_500_stdout.txt &
+CUDA_VISIBLE_DEVICES=12 python Cylinder2D_flower_systematic.py 26 250 > Cylinder2D_flower_26_250_stdout.txt &
+CUDA_VISIBLE_DEVICES=13 python Cylinder2D_flower_systematic.py 13 250 > Cylinder2D_flower_13_250_stdout.txt &
+CUDA_VISIBLE_DEVICES=14 python Cylinder2D_flower_systematic.py 7 250 > Cylinder2D_flower_7_250_stdout.txt &
+CUDA_VISIBLE_DEVICES=15 python Cylinder2D_flower_systematic.py 3 250 > Cylinder2D_flower_3_250_stdout.txt &
+
+for job in `jobs -p`
+do
+echo $job
+    wait $job || let "FAIL+=1"
+done
+
+echo $FAIL
+
+if [ "$FAIL" == "0" ];
+then
+echo "YAY!"
+else
+echo "FAIL! ($FAIL)"
+fi

Scripts/Cylinder2D_flower_case_studies_noise_1.sh

@@ -0,0 +1,35 @@
+#!/bin/bash
+
+FAIL=0
+
+CUDA_VISIBLE_DEVICES=0 python Cylinder2D_flower_systematic_noise.py 0.01 > Cylinder2D_flower_noise_01_stdout.txt &
+CUDA_VISIBLE_DEVICES=1 python Cylinder2D_flower_systematic_noise.py 0.02 > Cylinder2D_flower_noise_02_stdout.txt &
+CUDA_VISIBLE_DEVICES=2 python Cylinder2D_flower_systematic_noise.py 0.04 > Cylinder2D_flower_noise_04_stdout.txt &
+CUDA_VISIBLE_DEVICES=3 python Cylinder2D_flower_systematic_noise.py 0.06 > Cylinder2D_flower_noise_06_stdout.txt &
+CUDA_VISIBLE_DEVICES=4 python Cylinder2D_flower_systematic_noise.py 0.08 > Cylinder2D_flower_noise_08_stdout.txt &
+CUDA_VISIBLE_DEVICES=5 python Cylinder2D_flower_systematic_noise.py 0.10 > Cylinder2D_flower_noise_10_stdout.txt &
+CUDA_VISIBLE_DEVICES=6 python Cylinder2D_flower_systematic_noise.py 0.12 > Cylinder2D_flower_noise_12_stdout.txt &
+CUDA_VISIBLE_DEVICES=7 python Cylinder2D_flower_systematic_noise.py 0.14 > Cylinder2D_flower_noise_14_stdout.txt &
+CUDA_VISIBLE_DEVICES=8 python Cylinder2D_flower_systematic_noise.py 0.16 > Cylinder2D_flower_noise_16_stdout.txt &
+CUDA_VISIBLE_DEVICES=9 python Cylinder2D_flower_systematic_noise.py 0.18 > Cylinder2D_flower_noise_18_stdout.txt &
+CUDA_VISIBLE_DEVICES=10 python Cylinder2D_flower_systematic_noise.py 0.20 > Cylinder2D_flower_noise_20_stdout.txt &
+CUDA_VISIBLE_DEVICES=11 python Cylinder2D_flower_systematic_noise.py 0.22 > Cylinder2D_flower_noise_22_stdout.txt &
+CUDA_VISIBLE_DEVICES=12 python Cylinder2D_flower_systematic_noise.py 0.24 > Cylinder2D_flower_noise_24_stdout.txt &
+CUDA_VISIBLE_DEVICES=13 python Cylinder2D_flower_systematic_noise.py 0.26 > Cylinder2D_flower_noise_26_stdout.txt &
+CUDA_VISIBLE_DEVICES=14 python Cylinder2D_flower_systematic_noise.py 0.28 > Cylinder2D_flower_noise_28_stdout.txt &
+CUDA_VISIBLE_DEVICES=15 python Cylinder2D_flower_systematic_noise.py 0.30 > Cylinder2D_flower_noise_30_stdout.txt &
+
+for job in `jobs -p`
+do
+echo $job
+    wait $job || let "FAIL+=1"
+done
+
+echo $FAIL
+
+if [ "$FAIL" == "0" ];
+then
+echo "YAY!"
+else
+echo "FAIL! ($FAIL)"
+fi

Scripts/Cylinder2D_flower_case_studies_noise_2.sh

@@ -0,0 +1,35 @@
+#!/bin/bash
+
+FAIL=0
+
+CUDA_VISIBLE_DEVICES=0 python Cylinder2D_flower_systematic_noise.py 0.32 > Cylinder2D_flower_noise_32_stdout.txt &
+CUDA_VISIBLE_DEVICES=1 python Cylinder2D_flower_systematic_noise.py 0.34 > Cylinder2D_flower_noise_34_stdout.txt &
+CUDA_VISIBLE_DEVICES=2 python Cylinder2D_flower_systematic_noise.py 0.36 > Cylinder2D_flower_noise_36_stdout.txt &
+CUDA_VISIBLE_DEVICES=3 python Cylinder2D_flower_systematic_noise.py 0.38 > Cylinder2D_flower_noise_38_stdout.txt &
+CUDA_VISIBLE_DEVICES=4 python Cylinder2D_flower_systematic_noise.py 0.40 > Cylinder2D_flower_noise_40_stdout.txt &
+CUDA_VISIBLE_DEVICES=5 python Cylinder2D_flower_systematic_noise.py 0.42 > Cylinder2D_flower_noise_42_stdout.txt &
+CUDA_VISIBLE_DEVICES=6 python Cylinder2D_flower_systematic_noise.py 0.44 > Cylinder2D_flower_noise_44_stdout.txt &
+CUDA_VISIBLE_DEVICES=7 python Cylinder2D_flower_systematic_noise.py 0.46 > Cylinder2D_flower_noise_46_stdout.txt &
+CUDA_VISIBLE_DEVICES=8 python Cylinder2D_flower_systematic_noise.py 0.48 > Cylinder2D_flower_noise_48_stdout.txt &
+CUDA_VISIBLE_DEVICES=9 python Cylinder2D_flower_systematic_noise.py 0.50 > Cylinder2D_flower_noise_50_stdout.txt &
+CUDA_VISIBLE_DEVICES=10 python Cylinder2D_flower_systematic_noise.py 0.52 > Cylinder2D_flower_noise_52_stdout.txt &
+CUDA_VISIBLE_DEVICES=11 python Cylinder2D_flower_systematic_noise.py 0.54 > Cylinder2D_flower_noise_54_stdout.txt &
+CUDA_VISIBLE_DEVICES=12 python Cylinder2D_flower_systematic_noise.py 0.56 > Cylinder2D_flower_noise_56_stdout.txt &
+CUDA_VISIBLE_DEVICES=13 python Cylinder2D_flower_systematic_noise.py 0.58 > Cylinder2D_flower_noise_58_stdout.txt &
+CUDA_VISIBLE_DEVICES=14 python Cylinder2D_flower_systematic_noise.py 0.60 > Cylinder2D_flower_noise_60_stdout.txt &
+CUDA_VISIBLE_DEVICES=15 python Cylinder2D_flower_systematic_noise.py 0.62 > Cylinder2D_flower_noise_62_stdout.txt &
+
+for job in `jobs -p`
+do
+echo $job
+    wait $job || let "FAIL+=1"
+done
+
+echo $FAIL
+
+if [ "$FAIL" == "0" ];
+then
+echo "YAY!"
+else
+echo "FAIL! ($FAIL)"
+fi

Scripts/Cylinder2D_flower_case_studies_noise_3.sh

@@ -0,0 +1,35 @@
+#!/bin/bash
+
+FAIL=0
+
+CUDA_VISIBLE_DEVICES=0 python Cylinder2D_flower_systematic_noise.py 0.64 > Cylinder2D_flower_noise_64_stdout.txt &
+CUDA_VISIBLE_DEVICES=1 python Cylinder2D_flower_systematic_noise.py 0.66 > Cylinder2D_flower_noise_66_stdout.txt &
+CUDA_VISIBLE_DEVICES=2 python Cylinder2D_flower_systematic_noise.py 0.68 > Cylinder2D_flower_noise_68_stdout.txt &
+CUDA_VISIBLE_DEVICES=3 python Cylinder2D_flower_systematic_noise.py 0.70 > Cylinder2D_flower_noise_70_stdout.txt &
+CUDA_VISIBLE_DEVICES=4 python Cylinder2D_flower_systematic_noise.py 0.72 > Cylinder2D_flower_noise_72_stdout.txt &
+CUDA_VISIBLE_DEVICES=5 python Cylinder2D_flower_systematic_noise.py 0.74 > Cylinder2D_flower_noise_74_stdout.txt &
+CUDA_VISIBLE_DEVICES=6 python Cylinder2D_flower_systematic_noise.py 0.76 > Cylinder2D_flower_noise_76_stdout.txt &
+CUDA_VISIBLE_DEVICES=7 python Cylinder2D_flower_systematic_noise.py 0.78 > Cylinder2D_flower_noise_78_stdout.txt &
+CUDA_VISIBLE_DEVICES=8 python Cylinder2D_flower_systematic_noise.py 0.80 > Cylinder2D_flower_noise_80_stdout.txt &
+CUDA_VISIBLE_DEVICES=9 python Cylinder2D_flower_systematic_noise.py 0.82 > Cylinder2D_flower_noise_82_stdout.txt &
+CUDA_VISIBLE_DEVICES=10 python Cylinder2D_flower_systematic_noise.py 0.84 > Cylinder2D_flower_noise_84_stdout.txt &
+CUDA_VISIBLE_DEVICES=11 python Cylinder2D_flower_systematic_noise.py 0.86 > Cylinder2D_flower_noise_86_stdout.txt &
+CUDA_VISIBLE_DEVICES=12 python Cylinder2D_flower_systematic_noise.py 0.88 > Cylinder2D_flower_noise_88_stdout.txt &
+CUDA_VISIBLE_DEVICES=13 python Cylinder2D_flower_systematic_noise.py 0.90 > Cylinder2D_flower_noise_90_stdout.txt &
+CUDA_VISIBLE_DEVICES=14 python Cylinder2D_flower_systematic_noise.py 0.92 > Cylinder2D_flower_noise_92_stdout.txt &
+CUDA_VISIBLE_DEVICES=15 python Cylinder2D_flower_systematic_noise.py 0.94 > Cylinder2D_flower_noise_94_stdout.txt &
+
+for job in `jobs -p`
+do
+echo $job
+    wait $job || let "FAIL+=1"
+done
+
+echo $FAIL
+
+if [ "$FAIL" == "0" ];
+then
+echo "YAY!"
+else
+echo "FAIL! ($FAIL)"
+fi

Scripts/Cylinder2D_flower_case_studies_noise_4.sh

@@ -0,0 +1,35 @@
+#!/bin/bash
+
+FAIL=0
+
+CUDA_VISIBLE_DEVICES=0 python Cylinder2D_flower_systematic_noise.py 0.98 > Cylinder2D_flower_noise_98_stdout.txt &
+CUDA_VISIBLE_DEVICES=1 python Cylinder2D_flower_systematic_noise.py 1.00 > Cylinder2D_flower_noise_100_stdout.txt &
+CUDA_VISIBLE_DEVICES=2 python Cylinder2D_flower_systematic_noise.py 1.02 > Cylinder2D_flower_noise_102_stdout.txt &
+CUDA_VISIBLE_DEVICES=3 python Cylinder2D_flower_systematic_noise.py 1.04 > Cylinder2D_flower_noise_104_stdout.txt &
+CUDA_VISIBLE_DEVICES=4 python Cylinder2D_flower_systematic_noise.py 1.06 > Cylinder2D_flower_noise_106_stdout.txt &
+CUDA_VISIBLE_DEVICES=5 python Cylinder2D_flower_systematic_noise.py 1.08 > Cylinder2D_flower_noise_108_stdout.txt &
+CUDA_VISIBLE_DEVICES=6 python Cylinder2D_flower_systematic_noise.py 1.10 > Cylinder2D_flower_noise_110_stdout.txt &
+CUDA_VISIBLE_DEVICES=7 python Cylinder2D_flower_systematic_noise.py 1.12 > Cylinder2D_flower_noise_112_stdout.txt &
+CUDA_VISIBLE_DEVICES=8 python Cylinder2D_flower_systematic_noise.py 1.14 > Cylinder2D_flower_noise_114_stdout.txt &
+CUDA_VISIBLE_DEVICES=9 python Cylinder2D_flower_systematic_noise.py 1.16 > Cylinder2D_flower_noise_116_stdout.txt &
+CUDA_VISIBLE_DEVICES=10 python Cylinder2D_flower_systematic_noise.py 1.18 > Cylinder2D_flower_noise_118_stdout.txt &
+CUDA_VISIBLE_DEVICES=11 python Cylinder2D_flower_systematic_noise.py 1.20 > Cylinder2D_flower_noise_120_stdout.txt &
+CUDA_VISIBLE_DEVICES=12 python Cylinder2D_flower_systematic_noise.py 1.22 > Cylinder2D_flower_noise_122_stdout.txt &
+CUDA_VISIBLE_DEVICES=13 python Cylinder2D_flower_systematic_noise.py 1.24 > Cylinder2D_flower_noise_124_stdout.txt &
+CUDA_VISIBLE_DEVICES=14 python Cylinder2D_flower_systematic_noise.py 1.26 > Cylinder2D_flower_noise_126_stdout.txt &
+CUDA_VISIBLE_DEVICES=15 python Cylinder2D_flower_systematic_noise.py 1.28 > Cylinder2D_flower_noise_128_stdout.txt &
+
+for job in `jobs -p`
+do
+echo $job
+    wait $job || let "FAIL+=1"
+done
+
+echo $FAIL
+
+if [ "$FAIL" == "0" ];
+then
+echo "YAY!"
+else
+echo "FAIL! ($FAIL)"
+fi

Scripts/Cylinder2D_flower_case_studies_noise_5.sh

@@ -0,0 +1,35 @@
+#!/bin/bash
+
+FAIL=0
+
+CUDA_VISIBLE_DEVICES=0 python Cylinder2D_flower_systematic_noise.py 1.30 > Cylinder2D_flower_noise_130_stdout.txt &
+CUDA_VISIBLE_DEVICES=1 python Cylinder2D_flower_systematic_noise.py 1.32 > Cylinder2D_flower_noise_132_stdout.txt &
+CUDA_VISIBLE_DEVICES=2 python Cylinder2D_flower_systematic_noise.py 1.34 > Cylinder2D_flower_noise_134_stdout.txt &
+CUDA_VISIBLE_DEVICES=3 python Cylinder2D_flower_systematic_noise.py 1.36 > Cylinder2D_flower_noise_136_stdout.txt &
+CUDA_VISIBLE_DEVICES=4 python Cylinder2D_flower_systematic_noise.py 1.38 > Cylinder2D_flower_noise_138_stdout.txt &
+CUDA_VISIBLE_DEVICES=5 python Cylinder2D_flower_systematic_noise.py 1.40 > Cylinder2D_flower_noise_140_stdout.txt &
+CUDA_VISIBLE_DEVICES=6 python Cylinder2D_flower_systematic_noise.py 1.42 > Cylinder2D_flower_noise_142_stdout.txt &
+CUDA_VISIBLE_DEVICES=7 python Cylinder2D_flower_systematic_noise.py 1.44 > Cylinder2D_flower_noise_144_stdout.txt &
+CUDA_VISIBLE_DEVICES=8 python Cylinder2D_flower_systematic_noise.py 1.46 > Cylinder2D_flower_noise_146_stdout.txt &
+CUDA_VISIBLE_DEVICES=9 python Cylinder2D_flower_systematic_noise.py 1.48 > Cylinder2D_flower_noise_148_stdout.txt &
+CUDA_VISIBLE_DEVICES=10 python Cylinder2D_flower_systematic_noise.py 1.50 > Cylinder2D_flower_noise_150_stdout.txt &
+CUDA_VISIBLE_DEVICES=11 python Cylinder2D_flower_systematic_noise.py 1.52 > Cylinder2D_flower_noise_152_stdout.txt &
+CUDA_VISIBLE_DEVICES=12 python Cylinder2D_flower_systematic_noise.py 1.54 > Cylinder2D_flower_noise_154_stdout.txt &
+CUDA_VISIBLE_DEVICES=13 python Cylinder2D_flower_systematic_noise.py 1.56 > Cylinder2D_flower_noise_156_stdout.txt &
+CUDA_VISIBLE_DEVICES=14 python Cylinder2D_flower_systematic_noise.py 1.58 > Cylinder2D_flower_noise_158_stdout.txt &
+CUDA_VISIBLE_DEVICES=15 python Cylinder2D_flower_systematic_noise.py 1.60 > Cylinder2D_flower_noise_160_stdout.txt &
+
+for job in `jobs -p`
+do
+echo $job
+    wait $job || let "FAIL+=1"
+done
+
+echo $FAIL
+
+if [ "$FAIL" == "0" ];
+then
+echo "YAY!"
+else
+echo "FAIL! ($FAIL)"
+fi

Scripts/kitchen_sink.sh

@@ -0,0 +1,30 @@
+#!/bin/bash
+
+FAIL=0
+
+CUDA_VISIBLE_DEVICES=0 python Stenosis2D.py > Stenosis2D_stdout.txt &
+CUDA_VISIBLE_DEVICES=1 python Cylinder2D.py > Cylinder2D_stdout.txt &
+CUDA_VISIBLE_DEVICES=2 python Aneurysm3D.py > Aneurysm3D_stdout.txt &
+CUDA_VISIBLE_DEVICES=3 python Cylinder2D_Pec_Re.py > Cylinder2D_Pec_Re_stdout.txt &
+CUDA_VISIBLE_DEVICES=4 python Stenosis2D_Pec_Re.py > Stenosis2D_Pec_Re_stdout.txt &
+CUDA_VISIBLE_DEVICES=5 python Davinci.py > Davinci_stdout.txt &
+CUDA_VISIBLE_DEVICES=6 python Cylinder3D.py > Cylinder3D_stdout.txt &
+CUDA_VISIBLE_DEVICES=7 python Cylinder2D_flower.py > Cylinder2D_flower_stdout.txt &
+CUDA_VISIBLE_DEVICES=8 python Cylinder2D_No_Slip.py > Cylinder2D_No_Slip_stdout.txt &
+CUDA_VISIBLE_DEVICES=9 python Cylinder2D_flower_convergence_plot.py > Cylinder2D_flower_convergence_plot_stdout.txt &
+CUDA_VISIBLE_DEVICES=10 python Aneurysm3D_Wall_Stresses.py > Aneurysm3D_Wall_Stresses_stdout.txt &
+
+for job in `jobs -p`
+do
+echo $job
+    wait $job || let "FAIL+=1"
+done
+
+echo $FAIL
+
+if [ "$FAIL" == "0" ];
+then
+echo "YAY!"
+else
+echo "FAIL! ($FAIL)"
+fi