Adding grid vel term to Euler wall/symmetry BC.

Author: economon

SU2_CFD/src/solvers/CEulerSolver.cpp

@@ -6874,6 +6874,12 @@ void CEulerSolver::BC_Sym_Plane(CGeometry      *geometry,
             velocity in normal direction: v_r = v - 2 (v*n)n ---*/
       ProjVelocity_i = nodes->GetProjVel(iPoint,UnitNormal);
 
+      /*--- Adjustment to v.n due to grid movement. ---*/
+      if (dynamic_grid) {
+          for (iDim = 0; iDim < nDim; iDim++)
+            ProjVelocity_i -= geometry->nodes->GetGridVel(iPoint)[iDim]*UnitNormal[iDim];
+      }
+  
       for (iDim = 0; iDim < nDim; iDim++)
         V_reflected[iDim+1] = nodes->GetVelocity(iPoint,iDim) - 2.0 * ProjVelocity_i*UnitNormal[iDim];
 

TestCases/rotating/caradonna_tung/rot_caradonna_tung.cfg

@@ -4,7 +4,7 @@
 % Case description: Caradonna & Tung rotor in hover (still air)                %
 % Author: Thomas D. Economon                                                   %
 % Institution: Stanford University                                             %
-% Date: 2013.09.29                                                             %
+% Date: 2020.05.24                                                             %
 % File Version 7.0.4 "Blackbird"                                                %
 %                                                                              %
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
@@ -91,39 +91,39 @@ MARKER_MONITORING= ( blade_1, blade_2 )
 NUM_METHOD_GRAD= GREEN_GAUSS
 %
 % Courant-Friedrichs-Lewy condition of the finest grid
-CFL_NUMBER= 6.0
+CFL_NUMBER= 1e3
 %
 % Adaptive CFL number (NO, YES)
 CFL_ADAPT= NO
 %
 % Parameters of the adaptive CFL number (factor down, factor up, CFL min value,
 %                                        CFL max value )
-CFL_ADAPT_PARAM= ( 1.5, 0.5, 1.0, 100.0 )
+CFL_ADAPT_PARAM= ( 0.1, 2.0, 15, 1e6 )
 %
 % Runge-Kutta alpha coefficients
 RK_ALPHA_COEFF= ( 0.66667, 0.66667, 1.000000 )
 %
 % Number of total iterations
-EXT_ITER= 99999
+ITER= 99999
 
 % ------------------------ LINEAR SOLVER DEFINITION ---------------------------%
 %
 % Linear solver for implicit formulations (BCGSTAB, FGMRES)
 LINEAR_SOLVER= FGMRES
 %
 % Preconditioner of the Krylov linear solver (JACOBI, LINELET, LU_SGS)
-LINEAR_SOLVER_PREC= LU_SGS
+LINEAR_SOLVER_PREC= ILU
 %
 % Minimum error of the linear solver for implicit formulations
-LINEAR_SOLVER_ERROR= 1E-4
+LINEAR_SOLVER_ERROR= 1E-10
 %
 % Max number of iterations of the linear solver for the implicit formulation
-LINEAR_SOLVER_ITER= 4
+LINEAR_SOLVER_ITER= 10
 
 % -------------------------- MULTIGRID PARAMETERS -----------------------------%
 %
 % Multi-Grid Levels (0 = no multi-grid)
-MGLEVEL= 2
+MGLEVEL= 3
 %
 % Multi-grid cycle (V_CYCLE, W_CYCLE, FULLMG_CYCLE)
 MGCYCLE= W_CYCLE
@@ -168,7 +168,6 @@ TIME_DISCRE_FLOW= EULER_IMPLICIT
 %
 CONV_CRITERIA= RESIDUAL
 %
-%
 % Min value of the residual (log10 of the residual)
 CONV_RESIDUAL_MINVAL= -12
 %
@@ -231,3 +230,5 @@ WRT_SOL_FREQ= 100
 %
 % Writing convergence history frequency
 WRT_CON_FREQ= 1
+
+SCREEN_OUTPUT=(INNER_ITER, WALL_TIME, RMS_DENSITY, LIFT, DRAG, AVG_CFL)