-
Notifications
You must be signed in to change notification settings - Fork 81
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
- Loading branch information
Showing
491 changed files
with
2,822,337 additions
and
4,190 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,8 @@ | ||
--- | ||
layout: default | ||
--- | ||
|
||
<div class="container text-center"> | ||
<p>The page you are looking for cannot be found.</p> | ||
<h2>404</h2> | ||
</div> |
File renamed without changes
File renamed without changes
File renamed without changes
File renamed without changes
File renamed without changes
File renamed without changes
File renamed without changes
367 changes: 367 additions & 0 deletions
367
Inviscid_2D_Unconstrained_NACA0012/inv_NACA0012_basic.cfg
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Original file line number | Diff line number | Diff line change |
---|---|---|
@@ -0,0 +1,367 @@ | ||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% | ||
% % | ||
% SU2 configuration file % | ||
% Case description: Transonic inviscid optimization of a NACA0012 airfoil % | ||
% Author: Francisco Palacios % | ||
% Institution: Stanford University % | ||
% Date: 2013.09.29 % | ||
% File Version 5.0.0 "Raven" % | ||
% % | ||
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% | ||
|
||
% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------% | ||
% | ||
% Physical governing equations (EULER, NAVIER_STOKES, | ||
% WAVE_EQUATION, HEAT_EQUATION, FEM_ELASTICITY, | ||
% POISSON_EQUATION) | ||
PHYSICAL_PROBLEM= EULER | ||
% | ||
% Mathematical problem (DIRECT, CONTINUOUS_ADJOINT) | ||
MATH_PROBLEM= DIRECT | ||
% | ||
% Restart solution (NO, YES) | ||
RESTART_SOL= YES | ||
|
||
% -------------------- COMPRESSIBLE FREE-STREAM DEFINITION --------------------% | ||
% | ||
% Mach number (non-dimensional, based on the free-stream values) | ||
MACH_NUMBER= 0.8 | ||
% | ||
% Angle of attack (degrees) | ||
AOA= 1.25 | ||
% | ||
% Free-stream pressure (101325.0 N/m^2 by default, only Euler flows) | ||
FREESTREAM_PRESSURE= 101325.0 | ||
% | ||
% Free-stream temperature (288.15 K by default) | ||
FREESTREAM_TEMPERATURE= 288.15 | ||
|
||
% ---------------------- REFERENCE VALUE DEFINITION ---------------------------% | ||
% | ||
% Reference origin for moment computation | ||
REF_ORIGIN_MOMENT_X = 0.25 | ||
REF_ORIGIN_MOMENT_Y = 0.00 | ||
REF_ORIGIN_MOMENT_Z = 0.00 | ||
% | ||
% Reference length for pitching, rolling, and yawing non-dimensional moment | ||
REF_LENGTH= 1.0 | ||
% | ||
% Reference area for force coefficients (0 implies automatic calculation) | ||
REF_AREA= 1.0 | ||
% | ||
% Flow non-dimensionalization (DIMENSIONAL, FREESTREAM_PRESS_EQ_ONE, | ||
% FREESTREAM_VEL_EQ_MACH, FREESTREAM_VEL_EQ_ONE) | ||
REF_DIMENSIONALIZATION= FREESTREAM_PRESS_EQ_ONE | ||
|
||
% ----------------------- BOUNDARY CONDITION DEFINITION -----------------------% | ||
% | ||
% Marker of the Euler boundary (0 = no marker) | ||
MARKER_EULER= ( airfoil ) | ||
% | ||
% Marker of the far field (0 = no marker) | ||
MARKER_FAR= ( farfield ) | ||
|
||
% ------------------------ SURFACES IDENTIFICATION ----------------------------% | ||
% | ||
% Marker of the surface which is going to be plotted or designed | ||
MARKER_PLOTTING= ( airfoil ) | ||
% | ||
% Marker of the surface where the functional (Cd, Cl, etc.) will be evaluated | ||
MARKER_MONITORING= ( airfoil ) | ||
|
||
% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------% | ||
% | ||
% Numerical method for spatial gradients (GREEN_GAUSS, WEIGHTED_LEAST_SQUARES) | ||
NUM_METHOD_GRAD= GREEN_GAUSS | ||
% | ||
% Courant-Friedrichs-Lewy condition of the finest grid | ||
CFL_NUMBER= 10.0 | ||
% | ||
% 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 ) | ||
% | ||
% Runge-Kutta alpha coefficients | ||
RK_ALPHA_COEFF= ( 0.66667, 0.66667, 1.000000 ) | ||
% | ||
% Number of total iterations | ||
EXT_ITER= 1000 | ||
|
||
% ------------------------ LINEAR SOLVER DEFINITION ---------------------------% | ||
% | ||
% Linear solver for the implicit (or discrete adjoint) formulation (LU_SGS, | ||
% SYM_GAUSS_SEIDEL, BCGSTAB, GMRES) | ||
LINEAR_SOLVER= FGMRES | ||
% | ||
% Preconditioner of the Krylov linear solver (NONE, JACOBI, LINELET, LUSGS) | ||
LINEAR_SOLVER_PREC= LU_SGS | ||
% | ||
% Min error of the linear solver for the implicit formulation | ||
LINEAR_SOLVER_ERROR= 1E-4 | ||
% | ||
% Max number of iterations of the linear solver for the implicit formulation | ||
LINEAR_SOLVER_ITER= 2 | ||
|
||
% -------------------------- MULTIGRID PARAMETERS -----------------------------% | ||
% | ||
% Multi-Grid Levels (0 = no multi-grid) | ||
MGLEVEL= 2 | ||
% | ||
% Multi-grid cycle (V_CYCLE, W_CYCLE, FULLMG_CYCLE) | ||
MGCYCLE= V_CYCLE | ||
% | ||
% Multi-Grid PreSmoothing Level | ||
MG_PRE_SMOOTH= ( 1, 2, 3, 3 ) | ||
% | ||
% Multi-Grid PostSmoothing Level | ||
MG_POST_SMOOTH= ( 0, 0, 0, 0 ) | ||
% | ||
% Jacobi implicit smoothing of the correction | ||
MG_CORRECTION_SMOOTH= ( 0, 0, 0, 0 ) | ||
% | ||
% Damping factor for the residual restriction | ||
MG_DAMP_RESTRICTION= 1.0 | ||
% | ||
% Damping factor for the correction prolongation | ||
MG_DAMP_PROLONGATION= 1.0 | ||
|
||
% --------------------- FLOW NUMERICAL METHOD DEFINITION ----------------------% | ||
% Convective numerical method (JST, LAX-FRIEDRICH, ROE-1ST_ORDER, | ||
% ROE-2ND_ORDER) | ||
CONV_NUM_METHOD_FLOW= JST | ||
% | ||
% Slope limiter (VENKATAKRISHNAN) | ||
SLOPE_LIMITER_FLOW= VENKATAKRISHNAN | ||
% | ||
% 2nd and 4th order artificial dissipation coefficients | ||
JST_SENSOR_COEFF= ( 0.5, 0.02 ) | ||
% | ||
% Time discretization (RUNGE-KUTTA_EXPLICIT, EULER_IMPLICIT, EULER_EXPLICIT) | ||
TIME_DISCRE_FLOW= EULER_IMPLICIT | ||
|
||
% ---------------- ADJOINT-FLOW NUMERICAL METHOD DEFINITION -------------------% | ||
% Adjoint problem boundary condition (DRAG, LIFT, SIDEFORCE, MOMENT_X, | ||
% MOMENT_Y, MOMENT_Z, EFFICIENCY, | ||
% EQUIVALENT_AREA, NEARFIELD_PRESSURE, | ||
% FORCE_X, FORCE_Y, FORCE_Z, THRUST, | ||
% TORQUE) | ||
OBJECTIVE_FUNCTION= DRAG | ||
% | ||
% Convective numerical method (JST, LAX-FRIEDRICH, ROE-1ST_ORDER, | ||
% ROE-2ND_ORDER) | ||
CONV_NUM_METHOD_ADJFLOW= JST | ||
% | ||
% Slope limiter (VENKATAKRISHNAN, SHARP_EDGES) | ||
SLOPE_LIMITER_ADJFLOW= VENKATAKRISHNAN | ||
% | ||
% 2nd, and 4th order artificial dissipation coefficients | ||
ADJ_JST_SENSOR_COEFF= ( 0.0, 0.02 ) | ||
% | ||
% Time discretization (RUNGE-KUTTA_EXPLICIT, EULER_IMPLICIT) | ||
TIME_DISCRE_ADJFLOW= EULER_IMPLICIT | ||
% | ||
% Reduction factor of the CFL coefficient in the adjoint problem | ||
CFL_REDUCTION_ADJFLOW= 0.8 | ||
% | ||
% Limit value for the adjoint variable | ||
LIMIT_ADJFLOW= 1E6 | ||
|
||
% ----------------------- GEOMETRY EVALUATION PARAMETERS ----------------------% | ||
% | ||
% Marker(s) of the surface where geometrical based function will be evaluated | ||
GEO_MARKER= ( airfoil ) | ||
% | ||
% Description of the geometry to be analyzed (AIRFOIL, WING, FUSELAGE) | ||
GEO_DESCRIPTION= AIRFOIL | ||
% | ||
% Geometrical evaluation mode (FUNCTION, GRADIENT) | ||
GEO_MODE= FUNCTION | ||
|
||
% ----------------------- DESIGN VARIABLE PARAMETERS --------------------------% | ||
% | ||
% Kind of deformation (FFD_SETTING, HICKS_HENNE, HICKS_HENNE_NORMAL, PARABOLIC, | ||
% HICKS_HENNE_SHOCK, NACA_4DIGITS, DISPLACEMENT, ROTATION, | ||
% FFD_CONTROL_POINT, FFD_DIHEDRAL_ANGLE, FFD_TWIST_ANGLE, | ||
% FFD_ROTATION) | ||
DV_KIND= HICKS_HENNE | ||
% | ||
% Marker of the surface in which we are going apply the shape deformation | ||
DV_MARKER= ( airfoil ) | ||
% | ||
% Parameters of the shape deformation | ||
% - HICKS_HENNE_FAMILY ( Lower(0)/Upper(1) side, x_Loc ) | ||
% - NACA_4DIGITS ( 1st digit, 2nd digit, 3rd and 4th digit ) | ||
% - PARABOLIC ( 1st digit, 2nd and 3rd digit ) | ||
% - DISPLACEMENT ( x_Disp, y_Disp, z_Disp ) | ||
% - ROTATION ( x_Orig, y_Orig, z_Orig, x_End, y_End, z_End ) | ||
DV_PARAM= ( 1, 0.5 ) | ||
% | ||
% Value of the shape deformation deformation | ||
DV_VALUE= 1.0 | ||
|
||
% ------------------------ GRID DEFORMATION PARAMETERS ------------------------% | ||
% | ||
% Number of smoothing iterations for FEA mesh deformation | ||
DEFORM_LINEAR_ITER= 500 | ||
% | ||
% Number of nonlinear deformation iterations (surface deformation increments) | ||
DEFORM_NONLINEAR_ITER= 1 | ||
% | ||
% Print the residuals during mesh deformation to the console (YES, NO) | ||
DEFORM_CONSOLE_OUTPUT= YES | ||
% | ||
% Factor to multiply smallest cell volume for deform tolerance (0.001 default) | ||
DEFORM_TOL_FACTOR = 0.001 | ||
% | ||
% Type of element stiffness imposed for FEA mesh deformation (INVERSE_VOLUME, | ||
% WALL_DISTANCE, CONSTANT_STIFFNESS) | ||
DEFORM_STIFFNESS_TYPE= INVERSE_VOLUME | ||
|
||
% --------------------------- CONVERGENCE PARAMETERS --------------------------% | ||
% Convergence criteria (CAUCHY, RESIDUAL) | ||
% | ||
CONV_CRITERIA= RESIDUAL | ||
% | ||
% Residual reduction (order of magnitude with respect to the initial value) | ||
RESIDUAL_REDUCTION= 6 | ||
% | ||
% Min value of the residual (log10 of the residual) | ||
RESIDUAL_MINVAL= -13 | ||
% | ||
% Start Cauchy criteria at iteration number | ||
STARTCONV_ITER= 10 | ||
% | ||
% Number of elements to apply the criteria | ||
CAUCHY_ELEMS= 100 | ||
% | ||
% Epsilon to control the series convergence | ||
CAUCHY_EPS= 1E-6 | ||
% | ||
% Direct function to apply the convergence criteria (LIFT, DRAG, NEARFIELD_PRESS) | ||
CAUCHY_FUNC_FLOW= DRAG | ||
% | ||
% Adjoint function to apply the convergence criteria (SENS_GEOMETRY, SENS_MACH) | ||
CAUCHY_FUNC_ADJFLOW= SENS_GEOMETRY | ||
|
||
% ------------------------- INPUT/OUTPUT INFORMATION --------------------------% | ||
% | ||
% Mesh input file | ||
MESH_FILENAME= mesh_NACA0012_inv.su2 | ||
% | ||
% Mesh input file format (SU2, CGNS, NETCDF_ASCII) | ||
MESH_FORMAT= SU2 | ||
% | ||
% Mesh output file | ||
MESH_OUT_FILENAME= mesh_out.su2 | ||
% | ||
% Restart flow input file | ||
SOLUTION_FLOW_FILENAME= solution_flow.dat | ||
% | ||
% Restart adjoint input file | ||
SOLUTION_ADJ_FILENAME= solution_adj.dat | ||
% | ||
% Output file format (PARAVIEW, TECPLOT) | ||
OUTPUT_FORMAT= TECPLOT | ||
% | ||
% Output file convergence history (w/o extension) | ||
CONV_FILENAME= history | ||
% | ||
% Output file restart flow | ||
RESTART_FLOW_FILENAME= restart_flow.dat | ||
% | ||
% Output file restart adjoint | ||
RESTART_ADJ_FILENAME= restart_adj.dat | ||
% | ||
% Output file flow (w/o extension) variables | ||
VOLUME_FLOW_FILENAME= flow | ||
% | ||
% Output file adjoint (w/o extension) variables | ||
VOLUME_ADJ_FILENAME= adjoint | ||
% | ||
% Output Objective function gradient (using continuous adjoint) | ||
GRAD_OBJFUNC_FILENAME= of_grad.dat | ||
% | ||
% Output file surface flow coefficient (w/o extension) | ||
SURFACE_FLOW_FILENAME= surface_flow | ||
% | ||
% Output file surface adjoint coefficient (w/o extension) | ||
SURFACE_ADJ_FILENAME= surface_adjoint | ||
% | ||
% Writing solution file frequency | ||
WRT_SOL_FREQ= 250 | ||
% | ||
% Writing solution file frequency for physical time steps (dual time) | ||
WRT_SOL_FREQ_DUALTIME= 1 | ||
% | ||
% Writing convergence history frequency | ||
WRT_CON_FREQ= 1 | ||
% | ||
% Writing convergence history frequency (dual time, only written to screen) | ||
WRT_CON_FREQ_DUALTIME= 10 | ||
% | ||
% Output rind layers in the solution files | ||
WRT_HALO= NO | ||
|
||
% --------------------- OPTIMAL SHAPE DESIGN DEFINITION -----------------------% | ||
% Available flow based objective functions or constraint functions | ||
% DRAG, LIFT, SIDEFORCE, EFFICIENCY, | ||
% FORCE_X, FORCE_Y, FORCE_Z, | ||
% MOMENT_X, MOMENT_Y, MOMENT_Z, | ||
% THRUST, TORQUE, FIGURE_OF_MERIT, | ||
% EQUIVALENT_AREA, NEARFIELD_PRESSURE, | ||
% | ||
% Available geometrical based objective functions or constraint functions | ||
% AIRFOIL_AREA, AIRFOIL_THICKNESS, AIRFOIL_CHORD, AIRFOIL_TOC, AIRFOIL_AOA, | ||
% WING_VOLUME, WING_MIN_THICKNESS, WING_MAX_THICKNESS, WING_MAX_CHORD, WING_MIN_TOC, WING_MAX_TWIST, WING_MAX_CURVATURE, WING_MAX_DIHEDRAL | ||
% STATION#_WIDTH, STATION#_AREA, STATION#_THICKNESS, STATION#_CHORD, STATION#_TOC, | ||
% STATION#_TWIST (where # is the index of the station defined in GEO_LOCATION_STATIONS) | ||
% | ||
% Available design variables | ||
% HICKS_HENNE ( 1, Scale | Mark. List | Lower(0)/Upper(1) side, x_Loc ) | ||
% SPHERICAL ( 3, Scale | Mark. List | ControlPoint_Index, Theta_Disp, R_Disp ) | ||
% NACA_4DIGITS ( 4, Scale | Mark. List | 1st digit, 2nd digit, 3rd and 4th digit ) | ||
% DISPLACEMENT ( 5, Scale | Mark. List | x_Disp, y_Disp, z_Disp ) | ||
% ROTATION ( 6, Scale | Mark. List | x_Axis, y_Axis, z_Axis, x_Turn, y_Turn, z_Turn ) | ||
% FFD_CONTROL_POINT ( 7, Scale | Mark. List | FFD_BoxTag, i_Ind, j_Ind, k_Ind, x_Mov, y_Mov, z_Mov ) | ||
% FFD_TWIST ( 9, Scale | Mark. List | FFD_BoxTag, x_Orig, y_Orig, z_Orig, x_End, y_End, z_End ) | ||
% FFD_ROTATION ( 10, Scale | Mark. List | FFD_BoxTag, x_Orig, y_Orig, z_Orig, x_End, y_End, z_End ) | ||
% FFD_CAMBER ( 11, Scale | Mark. List | FFD_BoxTag, i_Ind, j_Ind ) | ||
% FFD_THICKNESS ( 12, Scale | Mark. List | FFD_BoxTag, i_Ind, j_Ind ) | ||
% FFD_VOLUME ( 13, Scale | Mark. List | FFD_BoxTag, i_Ind, j_Ind ) | ||
% FOURIER ( 14, Scale | Mark. List | Lower(0)/Upper(1) side, index, cos(0)/sin(1) ) | ||
% | ||
% Optimization objective function with scaling factor, separated by semicolons. | ||
% To include quadratic penalty function: use OPT_CONSTRAINT option syntax within the OPT_OBJECTIVE list. | ||
% ex= Objective * Scale | ||
OPT_OBJECTIVE= DRAG | ||
% | ||
% Optimization constraint functions with pushing factors (affects its value, not the gradient | ||
% in the python scripts), separated by semicolons | ||
% ex= (Objective = Value ) * Scale, use '>','<','=' | ||
OPT_CONSTRAINT= NONE | ||
% | ||
% Factor to reduce the norm of the gradient (affects the objective function and gradient in the python scripts) | ||
% In general, a norm of the gradient ~1E-6 is desired. | ||
OPT_GRADIENT_FACTOR= 1E-6 | ||
% | ||
% Factor to relax or accelerate the optimizer convergence (affects the line search in SU2_DEF) | ||
% In general, surface deformations of 0.01'' or 0.0001m are desirable | ||
OPT_RELAX_FACTOR= 1E3 | ||
% | ||
% Maximum number of optimizer iterations | ||
OPT_ITERATIONS= 100 | ||
% | ||
% Requested accuracy | ||
OPT_ACCURACY= 1E-10 | ||
% | ||
% Upper bound for each design variable | ||
OPT_BOUND_UPPER= 0.1 | ||
% | ||
% Lower bound for each design variable | ||
OPT_BOUND_LOWER= -0.1 | ||
% | ||
% Optimization design variables, separated by semicolons | ||
DEFINITION_DV= ( 1, 1.0 | airfoil | 0, 0.05 ); ( 1, 1.0 | airfoil | 0, 0.10 ); ( 1, 1.0 | airfoil | 0, 0.15 ); ( 1, 1.0 | airfoil | 0, 0.20 ); ( 1, 1.0 | airfoil | 0, 0.25 ); ( 1, 1.0 | airfoil | 0, 0.30 ); ( 1, 1.0 | airfoil | 0, 0.35 ); ( 1, 1.0 | airfoil | 0, 0.40 ); ( 1, 1.0 | airfoil | 0, 0.45 ); ( 1, 1.0 | airfoil | 0, 0.50 ); ( 1, 1.0 | airfoil | 0, 0.55 ); ( 1, 1.0 | airfoil | 0, 0.60 ); ( 1, 1.0 | airfoil | 0, 0.65 ); ( 1, 1.0 | airfoil | 0, 0.70 ); ( 1, 1.0 | airfoil | 0, 0.75 ); ( 1, 1.0 | airfoil | 0, 0.80 ); ( 1, 1.0 | airfoil | 0, 0.85 ); ( 1, 1.0 | airfoil | 0, 0.90 ); ( 1, 1.0 | airfoil | 0, 0.95 ); ( 1, 1.0 | airfoil | 1, 0.05 ); ( 1, 1.0 | airfoil | 1, 0.10 ); ( 1, 1.0 | airfoil | 1, 0.15 ); ( 1, 1.0 | airfoil | 1, 0.20 ); ( 1, 1.0 | airfoil | 1, 0.25 ); ( 1, 1.0 | airfoil | 1, 0.30 ); ( 1, 1.0 | airfoil | 1, 0.35 ); ( 1, 1.0 | airfoil | 1, 0.40 ); ( 1, 1.0 | airfoil | 1, 0.45 ); ( 1, 1.0 | airfoil | 1, 0.50 ); ( 1, 1.0 | airfoil | 1, 0.55 ); ( 1, 1.0 | airfoil | 1, 0.60 ); ( 1, 1.0 | airfoil | 1, 0.65 ); ( 1, 1.0 | airfoil | 1, 0.70 ); ( 1, 1.0 | airfoil | 1, 0.75 ); ( 1, 1.0 | airfoil | 1, 0.80 ); ( 1, 1.0 | airfoil | 1, 0.85 ); ( 1, 1.0 | airfoil | 1, 0.90 ); ( 1, 1.0 | airfoil | 1, 0.95 ) |
Oops, something went wrong.