From 99ecc75a50aaf2faef6d1ac712de3cea977023b3 Mon Sep 17 00:00:00 2001 From: Samuel Ayala Date: Fri, 24 May 2024 17:19:47 -0400 Subject: [PATCH] add very fine mesh --- data/theodorsen.py | 9 +- mesh/infinite_rectangular_56x5.x | 260 +++++++++++++++++++++++++++++++ tests/test_uvlm_theodorsen.cpp | 17 +- 3 files changed, 271 insertions(+), 15 deletions(-) create mode 100644 mesh/infinite_rectangular_56x5.x diff --git a/data/theodorsen.py b/data/theodorsen.py index 090e236..4f247b7 100644 --- a/data/theodorsen.py +++ b/data/theodorsen.py @@ -68,11 +68,10 @@ def uvlm_data(filename): ) files = [ - "cl_data" - # "cl_data_01", - # "cl_data_03", - # "cl_data_05", - # "cl_data_07", + #"cl_data" + "cl_data_025", + "cl_data_050", + "cl_data_075", ] for file in files: diff --git a/mesh/infinite_rectangular_56x5.x b/mesh/infinite_rectangular_56x5.x new file mode 100644 index 0000000..5a9e0b2 --- /dev/null +++ b/mesh/infinite_rectangular_56x5.x @@ -0,0 +1,260 @@ + 1 + 57 6 1 + 0.0000000000000000 0.017857142857142905 0.035714285714285809 0.053571428571428603 + 0.071428571428571508 0.089285714285714302 0.10714285714285721 0.12500000000000000 + 0.14285714285714290 0.16071428571428581 0.17857142857142860 0.19642857142857151 + 0.21428571428571430 0.23214285714285721 0.25000000000000000 0.26785714285714290 + 0.28571428571428581 0.30357142857142860 0.32142857142857151 0.33928571428571430 + 0.35714285714285721 0.37500000000000000 0.39285714285714290 0.41071428571428570 + 0.42857142857142860 0.44642857142857151 0.46428571428571430 0.48214285714285721 + 0.50000000000000000 0.51785714285714290 0.53571428571428581 0.55357142857142860 + 0.57142857142857140 0.58928571428571430 0.60714285714285721 0.62500000000000000 + 0.64285714285714290 0.66071428571428581 0.67857142857142860 0.69642857142857140 + 0.71428571428571430 0.73214285714285721 0.75000000000000000 0.76785714285714290 + 0.78571428571428570 0.80357142857142860 0.82142857142857140 0.83928571428571430 + 0.85714285714285721 0.87500000000000000 0.89285714285714290 0.91071428571428570 + 0.92857142857142860 0.94642857142857140 0.96428571428571430 0.98214285714285710 + 1.0000000000000000 0.0000000000000000 0.017857142857142884 0.035714285714285768 + 0.053571428571428638 0.071428571428571438 0.089285714285714302 0.10714285714285718 + 0.12499999999999997 0.14285714285714288 0.16071428571428575 0.17857142857142860 + 0.19642857142857145 0.21428571428571433 0.23214285714285715 0.25000000000000000 + 0.26785714285714290 0.28571428571428575 0.30357142857142866 0.32142857142857151 + 0.33928571428571441 0.35714285714285732 0.37500000000000011 0.39285714285714290 + 0.41071428571428570 0.42857142857142871 0.44642857142857140 0.46428571428571430 + 0.48214285714285721 0.50000000000000000 0.51785714285714290 0.53571428571428581 + 0.55357142857142860 0.57142857142857140 0.58928571428571430 0.60714285714285721 + 0.62500000000000000 0.64285714285714268 0.66071428571428581 0.67857142857142860 + 0.69642857142857129 0.71428571428571441 0.73214285714285721 0.75000000000000000 + 0.76785714285714290 0.78571428571428570 0.80357142857142860 0.82142857142857129 + 0.83928571428571441 0.85714285714285698 0.87499999999999978 0.89285714285714290 + 0.91071428571428570 0.92857142857142860 0.94642857142857140 0.96428571428571430 + 0.98214285714285710 1.0000000000000000 0.0000000000000000 0.017857142857142860 + 0.035714285714285719 0.053571428571428693 0.071428571428571369 0.089285714285714371 + 0.10714285714285721 0.12500000000000000 0.14285714285714285 0.16071428571428573 + 0.17857142857142874 0.19642857142857137 0.21428571428571436 0.23214285714285721 + 0.25000000000000000 0.26785714285714285 0.28571428571428570 0.30357142857142871 + 0.32142857142857140 0.33928571428571430 0.35714285714285732 0.37500000000000011 + 0.39285714285714279 0.41071428571428570 0.42857142857142871 0.44642857142857145 + 0.46428571428571436 0.48214285714285721 0.50000000000000000 0.51785714285714290 + 0.53571428571428570 0.55357142857142849 0.57142857142857129 0.58928571428571419 + 0.60714285714285721 0.62500000000000000 0.64285714285714290 0.66071428571428570 + 0.67857142857142849 0.69642857142857151 0.71428571428571441 0.73214285714285743 + 0.75000000000000000 0.76785714285714268 0.78571428571428548 0.80357142857142827 + 0.82142857142857151 0.83928571428571430 0.85714285714285698 0.87500000000000000 + 0.89285714285714290 0.91071428571428570 0.92857142857142849 0.94642857142857129 + 0.96428571428571441 0.98214285714285721 1.0000000000000000 0.0000000000000000 + 0.017857142857142839 0.035714285714285678 0.053571428571428742 0.071428571428571314 + 0.089285714285714343 0.10714285714285721 0.12500000000000000 0.14285714285714285 + 0.16071428571428567 0.17857142857142869 0.19642857142857131 0.21428571428571436 + 0.23214285714285721 0.25000000000000000 0.26785714285714285 0.28571428571428570 + 0.30357142857142871 0.32142857142857129 0.33928571428571436 0.35714285714285721 + 0.37500000000000000 0.39285714285714285 0.41071428571428564 0.42857142857142871 + 0.44642857142857129 0.46428571428571452 0.48214285714285721 0.50000000000000000 + 0.51785714285714279 0.53571428571428570 0.55357142857142849 0.57142857142857129 + 0.58928571428571441 0.60714285714285721 0.62500000000000000 0.64285714285714279 + 0.66071428571428570 0.67857142857142849 0.69642857142857129 0.71428571428571441 + 0.73214285714285721 0.75000000000000000 0.76785714285714279 0.78571428571428559 + 0.80357142857142871 0.82142857142857106 0.83928571428571441 0.85714285714285721 + 0.87500000000000000 0.89285714285714302 0.91071428571428581 0.92857142857142849 + 0.94642857142857129 0.96428571428571441 0.98214285714285721 1.0000000000000000 + 0.0000000000000000 0.017857142857142815 0.035714285714285629 0.053571428571428784 + 0.071428571428571244 0.089285714285714385 0.10714285714285721 0.12500000000000000 + 0.14285714285714282 0.16071428571428564 0.17857142857142877 0.19642857142857123 + 0.21428571428571438 0.23214285714285721 0.25000000000000000 0.26785714285714279 + 0.28571428571428564 0.30357142857142877 0.32142857142857123 0.33928571428571441 + 0.35714285714285721 0.37500000000000000 0.39285714285714279 0.41071428571428559 + 0.42857142857142877 0.44642857142857123 0.46428571428571441 0.48214285714285721 + 0.50000000000000000 0.51785714285714279 0.53571428571428559 0.55357142857142838 + 0.57142857142857117 0.58928571428571441 0.60714285714285721 0.62500000000000000 + 0.64285714285714279 0.66071428571428559 0.67857142857142838 0.69642857142857117 + 0.71428571428571441 0.73214285714285721 0.75000000000000000 0.76785714285714279 + 0.78571428571428559 0.80357142857142838 0.82142857142857117 0.83928571428571441 + 0.85714285714285721 0.87500000000000000 0.89285714285714279 0.91071428571428559 + 0.92857142857142838 0.94642857142857117 0.96428571428571441 0.98214285714285721 + 1.0000000000000000 0.0000000000000000 0.017857142857142794 0.035714285714285587 + 0.053571428571428825 0.071428571428571175 0.089285714285714413 0.10714285714285721 + 0.12500000000000000 0.14285714285714279 0.16071428571428559 0.17857142857142883 + 0.19642857142857117 0.21428571428571441 0.23214285714285721 0.25000000000000000 + 0.26785714285714279 0.28571428571428559 0.30357142857142883 0.32142857142857117 + 0.33928571428571441 0.35714285714285721 0.37500000000000000 0.39285714285714279 + 0.41071428571428559 0.42857142857142883 0.44642857142857117 0.46428571428571441 + 0.48214285714285721 0.50000000000000000 0.51785714285714279 0.53571428571428559 + 0.55357142857142838 0.57142857142857117 0.58928571428571441 0.60714285714285721 + 0.62500000000000000 0.64285714285714279 0.66071428571428559 0.67857142857142838 + 0.69642857142857117 0.71428571428571441 0.73214285714285721 0.75000000000000000 + 0.76785714285714279 0.78571428571428559 0.80357142857142838 0.82142857142857117 + 0.83928571428571441 0.85714285714285721 0.87500000000000000 0.89285714285714279 + 0.91071428571428559 0.92857142857142838 0.94642857142857117 0.96428571428571441 + 0.98214285714285721 1.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 8000.0000000000000 8000.0000000000000 8000.0000000000009 + 8000.0000000000000 8000.0000000000000 8000.0000000000018 8000.0000000000018 + 8000.0000000000018 8000.0000000000018 8000.0000000000036 8000.0000000000018 + 8000.0000000000018 8000.0000000000018 8000.0000000000018 8000.0000000000009 + 8000.0000000000009 8000.0000000000009 8000.0000000000045 8000.0000000000009 + 8000.0000000000009 8000.0000000000009 8000.0000000000027 8000.0000000000045 + 8000.0000000000045 8000.0000000000036 8000.0000000000036 8000.0000000000036 + 8000.0000000000036 8000.0000000000036 8000.0000000000036 8000.0000000000036 + 8000.0000000000036 8000.0000000000036 8000.0000000000027 8000.0000000000045 + 8000.0000000000045 8000.0000000000045 8000.0000000000064 8000.0000000000064 + 8000.0000000000064 8000.0000000000064 8000.0000000000064 8000.0000000000064 + 8000.0000000000045 8000.0000000000064 8000.0000000000045 8000.0000000000045 + 8000.0000000000045 8000.0000000000064 8000.0000000000064 8000.0000000000064 + 8000.0000000000082 8000.0000000000073 8000.0000000000073 8000.0000000000073 + 8000.0000000000073 8000.0000000000073 16000.000000000000 16000.000000000000 + 16000.000000000000 16000.000000000000 16000.000000000000 16000.000000000000 + 16000.000000000000 16000.000000000000 16000.000000000000 16000.000000000000 + 16000.000000000000 16000.000000000000 16000.000000000000 16000.000000000000 + 16000.000000000000 16000.000000000002 16000.000000000002 16000.000000000002 + 15999.999999999998 15999.999999999995 15999.999999999998 15999.999999999998 + 16000.000000000002 16000.000000000002 16000.000000000002 15999.999999999995 + 15999.999999999995 15999.999999999998 15999.999999999998 15999.999999999998 + 15999.999999999995 15999.999999999995 15999.999999999995 16000.000000000002 + 15999.999999999998 15999.999999999998 15999.999999999998 15999.999999999995 + 15999.999999999995 15999.999999999995 15999.999999999995 15999.999999999995 + 15999.999999999996 15999.999999999996 15999.999999999996 15999.999999999996 + 15999.999999999996 15999.999999999996 15999.999999999996 15999.999999999996 + 15999.999999999996 15999.999999999996 15999.999999999996 15999.999999999996 + 15999.999999999996 15999.999999999996 15999.999999999996 24000.000000000000 + 24000.000000000000 24000.000000000000 23999.999999999996 24000.000000000000 + 24000.000000000004 24000.000000000000 24000.000000000000 24000.000000000000 + 24000.000000000000 23999.999999999996 24000.000000000000 24000.000000000000 + 24000.000000000000 24000.000000000000 24000.000000000000 24000.000000000000 + 24000.000000000000 24000.000000000000 24000.000000000000 24000.000000000000 + 24000.000000000000 24000.000000000000 24000.000000000000 24000.000000000000 + 24000.000000000000 24000.000000000000 24000.000000000000 23999.999999999996 + 23999.999999999996 23999.999999999996 23999.999999999996 24000.000000000004 + 23999.999999999996 23999.999999999996 23999.999999999996 23999.999999999996 + 23999.999999999996 24000.000000000004 24000.000000000004 24000.000000000004 + 24000.000000000004 24000.000000000004 24000.000000000004 24000.000000000004 + 24000.000000000004 24000.000000000004 24000.000000000004 24000.000000000004 + 24000.000000000004 24000.000000000004 24000.000000000004 24000.000000000004 + 24000.000000000004 24000.000000000004 24000.000000000004 24000.000000000004 + 32000.000000000000 32000.000000000000 31999.999999999993 32000.000000000000 + 32000.000000000004 32000.000000000004 32000.000000000000 32000.000000000000 + 32000.000000000000 31999.999999999993 32000.000000000000 32000.000000000000 + 32000.000000000000 32000.000000000000 32000.000000000000 32000.000000000000 + 32000.000000000000 32000.000000000000 32000.000000000000 32000.000000000000 + 32000.000000000000 32000.000000000000 32000.000000000000 32000.000000000000 + 32000.000000000000 32000.000000000000 32000.000000000000 32000.000000000000 + 32000.000000000000 32000.000000000000 32000.000000000000 32000.000000000000 + 32000.000000000000 32000.000000000000 32000.000000000000 32000.000000000000 + 32000.000000000000 32000.000000000000 32000.000000000000 32000.000000000000 + 32000.000000000000 32000.000000000000 32000.000000000000 32000.000000000000 + 32000.000000000000 32000.000000000000 32000.000000000000 32000.000000000000 + 32000.000000000000 32000.000000000000 32000.000000000000 32000.000000000000 + 32000.000000000000 32000.000000000000 32000.000000000000 32000.000000000000 + 32000.000000000000 40000.000000000000 40000.000000000000 40000.000000000000 + 40000.000000000000 40000.000000000000 40000.000000000000 40000.000000000000 + 40000.000000000000 40000.000000000000 40000.000000000000 40000.000000000000 + 40000.000000000000 40000.000000000000 40000.000000000000 40000.000000000000 + 40000.000000000000 40000.000000000000 40000.000000000000 40000.000000000000 + 40000.000000000000 40000.000000000000 40000.000000000000 40000.000000000000 + 40000.000000000000 40000.000000000000 40000.000000000000 40000.000000000000 + 40000.000000000000 40000.000000000000 40000.000000000000 40000.000000000000 + 40000.000000000000 40000.000000000000 40000.000000000000 40000.000000000000 + 40000.000000000000 40000.000000000000 40000.000000000000 40000.000000000000 + 40000.000000000000 40000.000000000000 40000.000000000000 40000.000000000000 + 40000.000000000000 40000.000000000000 40000.000000000000 40000.000000000000 + 40000.000000000000 40000.000000000000 40000.000000000000 40000.000000000000 + 40000.000000000000 40000.000000000000 40000.000000000000 40000.000000000000 + 40000.000000000000 40000.000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 0.0000000000000000 0.0000000000000000 + 0.0000000000000000 0.0000000000000000 diff --git a/tests/test_uvlm_theodorsen.cpp b/tests/test_uvlm_theodorsen.cpp index 9f0943f..e1be41a 100644 --- a/tests/test_uvlm_theodorsen.cpp +++ b/tests/test_uvlm_theodorsen.cpp @@ -15,7 +15,7 @@ #include "vlm_utils.hpp" #include "vlm_executor.hpp" -#define DEBUG_DISPLACEMENT_DATA +// #define DEBUG_DISPLACEMENT_DATA using namespace vlm; using namespace linalg::ostream_overloads; @@ -107,8 +107,8 @@ void print_buffer(const T* start, u64 size) { int main() { const tiny::ScopedTimer timer("UVLM TOTAL"); - const u64 ni = 40; - const u64 nj = 10; + const u64 ni = 56; + const u64 nj = 5; // vlm::Executor::instance(1); //const std::vector meshes = {"../../../../mesh/rectangular_5x10.x"}; const std::vector meshes = {"../../../../mesh/infinite_rectangular_" + std::to_string(ni) + "x" + std::to_string(nj) + ".x"}; @@ -123,7 +123,7 @@ int main() { const f32 cycles = 3.0f; const f32 u_inf = 1.0f; // freestream velocity const f32 amplitude = 0.1f; // amplitude of the wing motion - const f32 k = 0.7; // reduced frequency + const f32 k = 1.0; // reduced frequency const f32 omega = k * 2.0f * u_inf / (2*b); const f32 t_final = cycles * 2.0f * PI_f / omega; // 4 periods //const f32 t_final = 5.0f; @@ -230,12 +230,12 @@ int main() { } } + std::ofstream cl_data("cl_data.txt"); + cl_data << k << "\n"; + #ifdef DEBUG_DISPLACEMENT_DATA std::ofstream wing_data("wing_data.txt"); std::ofstream wake_data("wake_data.txt"); - std::ofstream cl_data("cl_data.txt"); - - cl_data << k << "\n"; wing_data << mesh->nc << " " << mesh->ns << "\n"; wing_data << vec_t.size() - 1 << "\n\n"; @@ -264,7 +264,6 @@ int main() { const u64 wake_start = (mesh->nc + mesh->nw - i) * (mesh->ns + 1); const u64 wake_end = mesh->nb_vertices_total(); - // std::cout << "Buffer size: " << mesh->v.x.size() << " | " << wake_start << " | " << wake_end << std::endl; dump_buffer(wake_data, mesh->v.x.data() + wake_start, mesh->v.x.data() + wake_end); dump_buffer(wake_data, mesh->v.y.data() + wake_start, mesh->v.y.data() + wake_end); @@ -297,9 +296,7 @@ int main() { const f32 cl_unsteady = backend->compute_coefficient_unsteady_cl(freestream, velocities, dt, mesh->s_ref, 0, mesh->ns); std::printf("t: %f, CL: %f\n", t, cl_unsteady); - #ifdef DEBUG_DISPLACEMENT_DATA cl_data << t << " " << mesh->v.z[0] << " " << cl_unsteady << " " << std::sin(omega * t) << "\n"; - #endif } //backend->wake_rollup(dt); backend->shed_gamma(); // shed before moving & incrementing currentnw