Create MVP for FFD-box writing

Author: TobiKattmann

design/Species_Transport/README.md

@@ -0,0 +1 @@
+The necessary files for this tutorial are over at /incompressible_flow/Inc_Species_Transport

incompressible_flow/Inc_Species_Transport/1__FFD-box-writing/species3_primitiveVenturi.cfg

@@ -140,12 +140,6 @@ WRT_PERFORMANCE= YES
 %
 % -------------------- FREE-FORM DEFORMATION PARAMETERS -----------------------%
 %
-% Tolerance of the Free-Form Deformation point inversion
-FFD_TOLERANCE= 1E-10
-%
-% Maximum number of iterations in the Free-Form Deformation point inversion
-FFD_ITERATIONS= 500
-%
 % FFD box definition: 3D case (FFD_BoxTag, X1, Y1, Z1, X2, Y2, Z2, X3, Y3, Z3, X4, Y4, Z4,
 %                              X5, Y5, Z5, X6, Y6, Z6, X7, Y7, Z7, X8, Y8, Z8)
 %                     2D case (FFD_BoxTag, X1, Y1, 0.0, X2, Y2, 0.0, X3, Y3, 0.0, X4, Y4, 0.0,
@@ -161,71 +155,18 @@ FFD_DEFINITION= (BOX, \
 0.0, 0.0, 0.0, \
 0.0, 0.0, 0.0 )
 %
-% FFD box degree: 3D case (x_degree, y_degree, z_degree)
-%                 2D case (x_degree, y_degree, 0)
+% FFD box degree: 3D case (i_degree, j_degree, k_degree)
+%                 2D case (i_degree, j_degree, 0)
 FFD_DEGREE= (6, 1, 0)
 %
-% Surface grid continuity at the intersection with the faces of the FFD boxes.
-% To keep a particular level of surface continuity, SU2 automatically freezes the right
-% number of control point planes (NO_DERIVATIVE, 1ST_DERIVATIVE, 2ND_DERIVATIVE, USER_INPUT)
-FFD_CONTINUITY= NO_DERIVATIVE
-%
-% Definition of the FFD planes to be frozen in the FFD (x,y,z).
-% Value from 0 FFD degree in that direction. Pick a value larger than degree if you don't want to fix any plane.
-%FFD_FIX_I= (0,2,3)
-%FFD_FIX_J= (0,2,3)
-%FFD_FIX_K= (0,2,3)
-
-% ----------------------- DESIGN VARIABLE PARAMETERS --------------------------%
-%
-% Kind of deformation (NO_DEFORMATION, SCALE_GRID, TRANSLATE_GRID, ROTATE_GRID,
-%                      FFD_SETTING, FFD_NACELLE,
-%                      FFD_CONTROL_POINT, FFD_CAMBER, FFD_THICKNESS, FFD_TWIST
-%                      FFD_CONTROL_POINT_2D, FFD_CAMBER_2D, FFD_THICKNESS_2D, 
-%                      FFD_TWIST_2D, HICKS_HENNE, SURFACE_BUMP, SURFACE_FILE)
 DV_KIND= FFD_SETTING
 %
 % Marker of the surface in which we are going apply the shape deformation
 % NOTE: for deformation the outlet should be a MARKER_SYM to hinder the mesh being ripped apart.
 DV_MARKER= ( wall )
 %
 % Parameters of the shape deformation
-% - NO_DEFORMATION ( 1.0 )
 % - FFD_SETTING ( 1.0 )
 % - FFD_CONTROL_POINT ( FFD_BoxTag, i_Ind, j_Ind, k_Ind, x_Disp, y_Disp, z_Disp )
 DV_PARAM= ( 1.0 )
-%
-% Value of the shape deformation
-DV_VALUE= 0.01
 
-% ------------------------ GRID DEFORMATION PARAMETERS ------------------------%
-%
-% Linear solver or smoother for implicit formulations (FGMRES, RESTARTED_FGMRES, BCGSTAB)
-DEFORM_LINEAR_SOLVER= FGMRES
-%
-% Preconditioner of the Krylov linear solver (ILU, LU_SGS, JACOBI)
-DEFORM_LINEAR_SOLVER_PREC= ILU
-%
-% Number of smoothing iterations for mesh deformation
-DEFORM_LINEAR_SOLVER_ITER= 1000
-%
-% Number of nonlinear deformation iterations (surface deformation increments)
-DEFORM_NONLINEAR_ITER= 1
-%
-% Minimum residual criteria for the linear solver convergence of grid deformation
-DEFORM_LINEAR_SOLVER_ERROR= 1E-14
-%
-% Print the residuals during mesh deformation to the console (YES, NO)
-DEFORM_CONSOLE_OUTPUT= YES
-%
-% Deformation coefficient (linear elasticity limits from -1.0 to 0.5, a larger
-% value is also possible)
-DEFORM_COEFF = 1E6
-%
-% Type of element stiffness imposed for FEA mesh deformation (INVERSE_VOLUME,
-%                                           WALL_DISTANCE, CONSTANT_STIFFNESS)
-DEFORM_STIFFNESS_TYPE= WALL_DISTANCE
-%
-% Deform the grid only close to the surface. It is possible to specify how much
-% of the volumetric grid is going to be deformed in meters or inches (1E6 by default)
-DEFORM_LIMIT = 1E6