Change unsteady CHT Testcase options.

Author: TobiKattmann

multiphysics/unsteady_cht/cht_2d_3cylinders.cfg

@@ -1,61 +1,50 @@
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %                                                                               %
 % SU2 configuration file                                                        %
-% Case description: 2D cylinder array with CHT couplings                        %
+% Case description: 2D transient cylinder array with CHT couplings              %
 % Author: O. Burghardt, T. Economon                                             %
 % Institution: Chair for Scientific Computing, TU Kaiserslautern                %
-% Date: August 8, 2019                                                          %
-% File Version 6.0.1 "Falcon"                                                   %
+% Date: March 18, 2021                                                          %
+% File Version 7.1.1 "Blackbird"                                                %
 %                                                                               %
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
 %
-% Physical governing equations (EULER, NAVIER_STOKES,
-%                               WAVE_EQUATION, HEAT_EQUATION, FEM_ELASTICITY,
-%                               POISSON_EQUATION)
 SOLVER= MULTIPHYSICS
 %
-% Mathematical problem (DIRECT, CONTINUOUS_ADJOINT, DISCRETE_ADJOINT)
-MATH_PROBLEM= DIRECT
-%
-%
 CONFIG_LIST = (flow_cylinder.cfg, solid_cylinder1.cfg, solid_cylinder2.cfg, solid_cylinder3.cfg)
 %
-%
 MARKER_ZONE_INTERFACE= (cylinder_outer1, cylinder_inner1, cylinder_outer2, cylinder_inner2, cylinder_outer3, cylinder_inner3)
-%
-%
 MARKER_CHT_INTERFACE= (cylinder_outer1, cylinder_inner1, cylinder_outer2, cylinder_inner2, cylinder_outer3, cylinder_inner3)
 %
-%
 TIME_DOMAIN = YES
-%
-%
 TIME_MARCHING= DUAL_TIME_STEPPING-2ND_ORDER
-%
-%
 TIME_STEP= 0.05
-%
-%
 MAX_TIME= 100.0
-%
-%
 TIME_ITER= 2000
 %
-% Number of total iterations
+% Number of total outer iterations
 OUTER_ITER = 100
 %
-% Courant-Friedrichs-Lewy condition of the finest grid
+% Courant-Friedrichs-Lewy condition
 CFL_NUMBER= 100.0
 %
-% Mesh input file
-MESH_FILENAME= mesh_cht_3cyl.su2
+% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
 %
-% Mesh input file format (SU2, CGNS, NETCDF_ASCII)
-MESH_FORMAT= SU2
+MESH_FILENAME= mesh_cht_3cyl.su2
 %
+SCREEN_OUTPUT= ( TIME_ITER, OUTR_ITER, \
+ BGS_PRESSURE[0], BGS_TEMPERATURE[0], BGS_TEMPERATURE[1], BGS_TEMPERATURE[2], BGS_TEMPERATURE[3], \
+ DRAG[0], SURFACE_STATIC_TEMPERATURE[0], AVG_TEMPERATURE[1], TOTAL_HEATFLUX[1] )
+SCREEN_WRT_FREQ_OUTER= 25
 %
 HISTORY_OUTPUT= (ITER, BGS_RES[0], BGS_RES[1], BGS_RES[2], BGS_RES[3], HEAT[0], AERO_COEFF[0])
 %
-%
 OUTPUT_FILES=(RESTART, PARAVIEW_MULTIBLOCK)
+%
+RESTART_SOL= NO
+READ_BINARY_RESTART= YES
+SOLUTION_FILENAME= solution
+RESTART_FILENAME= solution
+%
+VOLUME_FILENAME= flow
+

multiphysics/unsteady_cht/flow_cylinder.cfg

@@ -1,10 +1,10 @@
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 %                                                                              %
 % SU2 configuration file                                                       %
-% Case description: Steady incompressible laminar flow around heated cylinders %
+% Case description: Transient incomp laminar flow around heated cylinders      %
 %                                                                              %
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
-
+%
 % ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
 %
 % Physical governing equations (EULER, NAVIER_STOKES,
@@ -15,40 +15,11 @@ SOLVER= INC_NAVIER_STOKES
 % If Navier-Stokes, kind of turbulent model (NONE, SA)
 KIND_TURB_MODEL= NONE
 %
-% Restart solution (NO, YES)
-RESTART_SOL= NO
-%
-% Objective function in gradient evaluation   (DRAG, LIFT, SIDEFORCE, MOMENT_X,
-%                                             MOMENT_Y, MOMENT_Z, EFFICIENCY,
-%                                             EQUIVALENT_AREA, NEARFIELD_PRESSURE,
-%                                             FORCE_X, FORCE_Y, FORCE_Z, THRUST,
-%                                             TORQUE, TOTAL_HEATFLUX,
-%                                             MAXIMUM_HEATFLUX, INVERSE_DESIGN_PRESSURE,
-%                                             INVERSE_DESIGN_HEATFLUX, SURFACE_TOTAL_PRESSURE,
-%                                             SURFACE_MASSFLOW, SURFACE_STATIC_PRESSURE, SURFACE_MACH)
-% For a weighted sum of objectives: separate by commas, add OBJECTIVE_WEIGHT and MARKER_MONITORING in matching order.
-OBJECTIVE_FUNCTION= TOTAL_HEATFLUX
-%
-% List of weighting values when using more than one OBJECTIVE_FUNCTION. Separate by commas and match with MARKER_MONITORING.
-OBJECTIVE_WEIGHT = 1.0
-%
-% Read binary restart files (YES, NO)
-READ_BINARY_RESTART = YES
-%
-% Data written to history file
-HISTORY_OUTPUT=(ITER, RMS_RES, HEAT)
-
 % -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
 %
 % Farfield boundary marker(s) (NONE = no marker)
 MARKER_FAR= ( farfield )
 %
-% Marker(s) of the surface to be plotted or designed
-MARKER_PLOTTING= (cylinder_outer1, cylinder_outer2, cylinder_outer3)
-%
-% Marker(s) of the surface where the functional (Cd, Cl, etc.) will be evaluated
-MARKER_MONITORING= (cylinder_outer1, cylinder_outer2, cylinder_outer3)
-
 % ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
 %
 % Density model within the incompressible flow solver.
@@ -74,7 +45,7 @@ INC_TEMPERATURE_INIT= 288.15
 % INITIAL_VALUES (default), REFERENCE_VALUES, or DIMENSIONAL.
 % INC_*_REF values are ignored unless REFERENCE_VALUES is chosen.
 INC_NONDIM= DIMENSIONAL
-
+%
 % ---- IDEAL GAS, POLYTROPIC, VAN DER WAALS AND PENG ROBINSON CONSTANTS -------%
 %
 % Fluid model (STANDARD_AIR, IDEAL_GAS, VW_GAS, PR_GAS,
@@ -88,7 +59,7 @@ SPECIFIC_HEAT_CP= 1004.703
 % Molecular weight for an incompressible ideal gas (28.96 g/mol (air) default)
 % Incompressible fluids with energy eqn. only (CONSTANT_DENSITY, INC_IDEAL_GAS).
 MOLECULAR_WEIGHT= 28.96
-
+%
 % --------------------------- VISCOSITY MODEL ---------------------------------%
 %
 % Viscosity model (SUTHERLAND, CONSTANT_VISCOSITY).
@@ -119,12 +90,12 @@ PRANDTL_LAM= 0.72
 %
 % Turbulent Prandtl number (0.9 (air), only for CONSTANT_PRANDTL)
 PRANDTL_TURB= 0.90
-
+%
 % ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
 %
 % Numerical method for spatial gradients (GREEN_GAUSS, WEIGHTED_LEAST_SQUARES)
 NUM_METHOD_GRAD= GREEN_GAUSS
-
+%
 % ------------------------ LINEAR SOLVER DEFINITION ---------------------------%
 %
 % Linear solver or smoother for implicit formulations (BCGSTAB, FGMRES, SMOOTHER_JACOBI,
@@ -142,8 +113,8 @@ LINEAR_SOLVER_ILU_FILL_IN= 0
 LINEAR_SOLVER_ERROR= 1E-15
 %
 % Max number of iterations of the linear solver for the implicit formulation
-LINEAR_SOLVER_ITER= 5
-
+LINEAR_SOLVER_ITER= 20
+%
 % -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
 %
 % Convective numerical method (JST, LAX-FRIEDRICH, CUSP, ROE, AUSM, HLLC,
@@ -160,7 +131,7 @@ SLOPE_LIMITER_FLOW= NONE
 %
 % Time discretization (RUNGE-KUTTA_EXPLICIT, EULER_IMPLICIT, EULER_EXPLICIT)
 TIME_DISCRE_FLOW= EULER_IMPLICIT
-
+%
 % --------------------------- CONVERGENCE PARAMETERS --------------------------%
 %
 % Min value of the residual (log10 of the residual)
@@ -174,64 +145,12 @@ CONV_CAUCHY_ELEMS= 100
 %
 % Epsilon to control the series convergence
 CONV_CAUCHY_EPS= 1E-6
-
-% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
-%
-% Restart flow input file
-SOLUTION_FILENAME= solution_flow.dat
-%
-% Restart adjoint input file
-SOLUTION_ADJ_FILENAME= solution_adj.dat
-%
-% Output file format (TECPLOT, TECPLOT_BINARY, PARAVIEW,
-%                     FIELDVIEW, FIELDVIEW_BINARY)
-TABULAR_FORMAT= TECPLOT
-%
-% Output file convergence history (w/o extension)
-CONV_FILENAME= history
-%
-% Output file with the forces breakdown
-BREAKDOWN_FILENAME= forces_breakdown.dat
-%
-% Output file restart flow
-RESTART_FILENAME= restart_flow.dat
-%
-% Output file restart adjoint
-RESTART_ADJ_FILENAME= restart_adj.dat
-%
-% Output file flow (w/o extension) variables
-VOLUME_FILENAME= flow
-%
-% Output file adjoint (w/o extension) variables
-VOLUME_ADJ_FILENAME= adjoint
 %
-% Output Objective function
-VALUE_OBJFUNC_FILENAME= of_eval.dat
-%
-% Output objective function gradient (using continuous adjoint)
-GRAD_OBJFUNC_FILENAME= of_grad.dat
-%
-% Output file surface flow coefficient (w/o extension)
-SURFACE_FILENAME= surface_flow
-%
-% Output file surface adjoint coefficient (w/o extension)
-SURFACE_ADJ_FILENAME= surface_adjoint
-%
-% ------------------------ GRID DEFORMATION PARAMETERS ------------------------%
-%
-% Linear solver or smoother for implicit formulations (FGMRES, RESTARTED_FGMRES, BCGSTAB)
-DEFORM_LINEAR_SOLVER= FGMRES
-%
-% Number of smoothing iterations for mesh deformation
-DEFORM_LINEAR_SOLVER_ITER= 200
-%
-% Number of nonlinear deformation iterations (surface deformation increments)
-DEFORM_NONLINEAR_ITER= 1
+% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
 %
-% Print the residuals during mesh deformation to the console (YES, NO)
-DEFORM_CONSOLE_OUTPUT= YES
+MARKER_PLOTTING= (cylinder_outer1, cylinder_outer2, cylinder_outer3)
 %
-% Type of element stiffness imposed for FEA mesh deformation (INVERSE_VOLUME,
-%                                          WALL_DISTANCE, CONSTANT_STIFFNESS)
-DEFORM_STIFFNESS_TYPE= INVERSE_VOLUME
+MARKER_MONITORING= (cylinder_outer1, cylinder_outer2, cylinder_outer3)
 %
+MARKER_ANALYZE= (cylinder_outer1, cylinder_outer2, cylinder_outer3)
+