corrected sign convention in tut4

sandialabs · Dec 12, 2023 · 1413321 · 1413321
1 parent 4c4f60b
commit 1413321
Showing 1 changed file with 4 additions and 4 deletions.
diff --git a/examples/tutorial_4_Pioneer.ipynb b/examples/tutorial_4_Pioneer.ipynb
@@ -772,11 +772,11 @@
     "P_radiated = ((1/2)*(Rad_res * np.abs(Vel)**2 ).squeeze().sum('omega').item()) # after Falnes Eq. 6.4\n",
     "P_excited= (1/4)*(Fex*np.conjugate(Vel) + np.conjugate(Fex)*Vel ).squeeze().sum('omega').item() # after Falnes Eq. 6.3\n",
     "\n",
-    "power_flows = {'Optimal Excitation' : 2* (np.real(P_max_absorbable)),    #positive because the only inflow\n",
-    "            'Radiated': -1*(np.real(P_radiated)), #negative because \"out\"flow\n",
-    "            'Actual Excitation': -1*(np.real(P_excited)), #negative because \"out\"flow\n",
+    "power_flows = {'Optimal Excitation' : -2* (np.real(P_max_absorbable)),    \n",
+    "            'Radiated': -1*(np.real(P_radiated)), \n",
+    "            'Actual Excitation': -1*(np.real(P_excited)), \n",
     "}\n",
-    "power_flows['Unused Potential'] =  -1*power_flows['Optimal Excitation'] - power_flows['Actual Excitation']\n",
+    "power_flows['Unused Potential'] =  power_flows['Optimal Excitation'] - power_flows['Actual Excitation']\n",
     "power_flows['Absorbed'] =  power_flows['Actual Excitation'] - power_flows['Radiated']  # after Falnes Eq. 6.2 \n",
     "#absorbed by WEC-PTO (difference between actual excitation power and actual radiated power needs to be absorbed by PTO)\n",
     "\n"