QTF Example (#86)

* Modifies the MOST example to include the QTFs

* Fixes a typo in Mesh.cal

* runs mostIO before starting the tests

* comments the turbSim test as it a part of mostIO

* fixes a bug

* Adds the QTF results and modifies gitignore

* limits the test to the constant case only

* Small change so tests will rerun

* Update wecSimInputFile.m

* Update wecSimInputFile.m

disables the QTF calculation in the example case

* Delete Generalized_Body_Modes/gbmWEC-Sim.mat

---------

Co-authored-by: jtgrasb <87095491+jtgrasb@users.noreply.github.com>

Author: MShabara

.gitignore

@@ -95,6 +95,4 @@ WECCCOMP/**/SS*.fig
 Passive_Yaw/**/IrrYaw.mat
 Passive_Yaw/**/figYawIrr.fig
 yaw*.mat
-
-
-
+!MOST/hydroData/VolturnUS15MW_nemoh/Results/QTF/*.dat

MOST/hydroData/VolturnUS15MW_nemoh/Mesh.cal

@@ -0,0 +1,10 @@
+Mesh/Volturn_nemoh
+0
+ 0. 0. 
+ 0.000000 0.000000 0
+980
+ 2
+ 0. 
+ 1.
+1000.
+9.81

MOST/hydroData/VolturnUS15MW_nemoh/Nemoh.cal

@@ -31,7 +31,18 @@
 0				! Show pressure
 0	0 180		! Kochin function 		! Number of directions of calculation (0 for no calculations), Min and Max (degrees)
 0	50 400 400	! Free surface elevation 	! Number of points in x direction (0 for no calcutions) and y direction and dimensions of domain in x and y direction
-0					! Response Amplitude Operator (RAO), 0 no calculation, 1 calculated
+1					! Response Amplitude Operator (RAO), 0 no calculation, 1 calculated
 1					! output freq type, 1,2,3=[rad/s,Hz,s]
 --- QTF----------------------------------------------------------------------------------------------------------------------------
-0         				! QTF flag, 1 is calculated
+1         			! QTF flag, 1 is calculated 
+65	0.0628 4.082 	! Number of radial frequencies, Min, and Max values for the QTF computation
+1	  			  	! 1 Unidirection, Bidirection 0 
+2         			! Contrib, 1 DUOK, 2 DUOK+HASBO, 3 Full QTF (DUOK+HASBO+HASFS+ASYMP)
+NA 					! Name of free surface meshfile (Only for Contrib 3), type 'NA' if not applicable 
+0 	0	0			! Free surface QTF parameters: Re Nre NBessel (for Contrib 3)
+0 					! 1 Includes Hydrostatic terms of the quadratic first order motion, -[K]xi2_tilde
+1					! For QTFposProc, output freq type, 1,2,3=[rad/s,Hz,s]
+1         			! For QTFposProc, 1 includes DUOK in total QTFs, 0 otherwise
+1         			! For QTFposProc, 1 includes HASBO in total QTFs, 0 otherwise 
+0         			! For QTFposProc, 1 includes HASFS+ASYMP in total QTFs, 0 otherwise
+-----------------------------------------------------------------------------------------------------------------------------------

MOST/hydroData/VolturnUS15MW_nemoh/Results/QTF/OUT_QTFM_N_0.dat

@@ -0,0 +1,4 @@
+2				! Gauss quadrature (GQ) surface integration, N^2 GQ Nodes, specify N(1,4)
+0.001			! eps_zmin for determine minimum z of flow and source points of panel, zmin=eps_zmin*body_diameter
+1 				! 0 GAUSS ELIM.; 1 LU DECOMP.: 2 GMRES	!Linear system solver
+10 1e-5 1000  	! Restart parameter, Relative Tolerance, max iter -> additional input for GMRES

MOST/hydroData/VolturnUS15MW_nemoh/Results/QTF/OUT_QTFP_N_0.dat

@@ -44,11 +44,18 @@ function runBEMIO(testCase)
         end
         function runTurbSim(testCase)
             % Only re-run turbSim before testing. The *.mat file is too
-            % large to commit to the repo. 
+            % large to commit to the repo.
             cd(fullfile(testCase.mostDataDir,'turbSim'));
             run_turbsim();
             cd(testCase.testDir)
         end
+        function runMOSTIO(testCase)
+            % Test the pre-processing here so that it does not interfere
+            % with the precise comparison of results.
+            cd(testCase.mostDataDir);
+            mostIO
+            cd(testCase.testDir)
+        end
         function runConstantTest(testCase)
             cd(fullfile(testCase.testDir,'constant'))
             runLoadConstant;
@@ -199,13 +206,5 @@ function turbulent_windSpeed(testCase)
             fprintf(['Wind speed, Diff = '     ...
                      num2str(max(abs(org-new))) '\n']);
         end
-
-        function runMOSTIO(testCase)
-            % Test the pre-processing here so that it does not interfere
-            % with the precise comparison of results.
-            cd(testCase.mostDataDir);
-            mostIO
-            cd(testCase.testDir)
-        end
     end
 end

MOST/hydroData/VolturnUS15MW_nemoh/input.txt

@@ -39,6 +39,7 @@
                          0.00E+00	8.92E+06	0.00E+00	1.68E+10	0.00E+00	0.00E+00
                         -8.92E+06	0.00E+00	0.00E+00	0.00E+00	1.68E+10	0.00E+00
                          0.00E+00	0.00E+00	0.00E+00	0.00E+00	0.00E+00	4.80E+10];
+body(1).QTFs = 0;
 
 %%  Mooring class
 mooring(1) = mooringClass('mooring1');                                                                     % Initialize mooringClass

MOST/hydroData/VolturnUS15MW_nemoh/input_solver.txt

@@ -39,6 +39,7 @@
                          0.00E+00	8.92E+06	0.00E+00	1.68E+10	0.00E+00	0.00E+00
                         -8.92E+06	0.00E+00	0.00E+00	0.00E+00	1.68E+10	0.00E+00
                          0.00E+00	0.00E+00	0.00E+00	0.00E+00	0.00E+00	4.80E+10];
+body(1).QTFs = 0;
 
 %%  Mooring class
 mooring(1) = mooringClass('mooring1');                                                                     % Initialize mooringClass