Merge pull request #1656 from HL-ZHI/master

Fix missing phase angle in calculation of grid velocity for rigid pitching motion

Author: pcarruscag

Common/src/grid_movement/CVolumetricMovement.cpp

@@ -2066,9 +2066,9 @@ void CVolumetricMovement::Rigid_Pitching(CGeometry *geometry, CConfig *config, u
 
   /*--- Angular velocity at the new time ---*/
 
-  alphaDot[0] = -Omega[0]*Ampl[0]*cos(Omega[0]*time_new);
-  alphaDot[1] = -Omega[1]*Ampl[1]*cos(Omega[1]*time_new);
-  alphaDot[2] = -Omega[2]*Ampl[2]*cos(Omega[2]*time_new);
+  alphaDot[0] = -Omega[0]*Ampl[0]*cos(Omega[0]*time_new + Phase[0]);
+  alphaDot[1] = -Omega[1]*Ampl[1]*cos(Omega[1]*time_new + Phase[1]);
+  alphaDot[2] = -Omega[2]*Ampl[2]*cos(Omega[2]*time_new + Phase[2]);
 
   if (rank == MASTER_NODE && iter == 0) {
       cout << " Pitching frequency: (" << Omega[0] << ", " << Omega[1];