Skip to content

Commit 02f5388

Browse files
committed
Proposal for a test case for cosine gusts in z-direction on a 3D mesh.
Test duration: 0.83 min Not sure where the mesh and the restart files should go, so this will break the testing until I commit the missing files.
1 parent 2259d67 commit 02f5388

2 files changed

Lines changed: 181 additions & 0 deletions

File tree

Lines changed: 172 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,172 @@
1+
% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
2+
%
3+
% Solver type (EULER, NAVIER_STOKES, RANS,
4+
% INC_EULER, INC_NAVIER_STOKES, INC_RANS,
5+
% NEMO_EULER, NEMO_NAVIER_STOKES,
6+
% FEM_EULER, FEM_NAVIER_STOKES, FEM_RANS, FEM_LES,
7+
% HEAT_EQUATION_FVM, ELASTICITY)
8+
SOLVER= EULER
9+
%
10+
% Mathematical problem (DIRECT, CONTINUOUS_ADJOINT)
11+
MATH_PROBLEM= DIRECT
12+
%
13+
% Restart solution (NO, YES)
14+
RESTART_SOL= YES
15+
%
16+
% Iteration number to begin unsteady restarts
17+
RESTART_ITER= 72
18+
%
19+
% ------------------------------- SOLVER CONTROL ------------------------------%
20+
%
21+
% Maximum number of inner iterations
22+
INNER_ITER= 30
23+
%
24+
% Min value of the residual (log10 of the residual)
25+
CONV_RESIDUAL_MINVAL= -6
26+
%
27+
% Start convergence criteria at iteration number
28+
CONV_STARTITER= 0
29+
%
30+
% ------------------------- UNSTEADY SIMULATION -------------------------------%
31+
%
32+
TIME_DOMAIN=YES
33+
TIME_MARCHING= DUAL_TIME_STEPPING-2ND_ORDER
34+
TIME_STEP= 0.001
35+
TIME_ITER= 80
36+
%
37+
% --------------------------- GUST SIMULATION ---------------------------------%
38+
%
39+
% The gust simulation requires the GRID_MOVEMENT flag to be set to YES.
40+
% and the GRID_MOVEMENT_KIND to be set to GUST or any of the other options.
41+
% Apply a wind gust (NO, YES)
42+
GRID_MOVEMENT=GUST
43+
WIND_GUST= YES
44+
GUST_TYPE= ONE_M_COSINE
45+
GUST_DIR= Z_DIR
46+
GUST_WAVELENGTH= 5.0
47+
GUST_PERIODS= 1.0
48+
% Gust amplitude corresponds to ~2.0 deg AoA
49+
GUST_AMPL= 2.37
50+
%
51+
GUST_BEGIN_TIME= 0.0
52+
GUST_BEGIN_LOC=-10.0
53+
% -------------------- COMPRESSIBLE FREE-STREAM DEFINITION --------------------%
54+
%
55+
% Mach number (non-dimensional, based on the free-stream values)
56+
MACH_NUMBER= 0.2
57+
%
58+
% Angle of attack (degrees, only for compressible flows)
59+
AOA= 0.0
60+
%
61+
% Side-slip angle (degrees, only for compressible flows)
62+
SIDESLIP_ANGLE= 0.0
63+
%
64+
% Free-stream option to choose between density and temperature (default) for
65+
% initializing the solution (TEMPERATURE_FS, DENSITY_FS)
66+
FREESTREAM_OPTION= DENSITY_FS
67+
%
68+
% Free-stream pressure (101325.0 N/m^2, 2116.216 psf by default)
69+
FREESTREAM_PRESSURE= 101325.0
70+
%
71+
% Free-stream density (1.2886 Kg/m^3, 0.0025 slug/ft^3 by default)
72+
FREESTREAM_DENSITY= 1.225
73+
%
74+
% Free-stream temperature (288.15 K, 518.67 R by default)
75+
FREESTREAM_TEMPERATURE= 288.15
76+
%
77+
% ---------------------- REFERENCE VALUE DEFINITION ---------------------------%
78+
%
79+
% Reference origin for moment computation (m or in)
80+
REF_ORIGIN_MOMENT_X = 0.25
81+
REF_ORIGIN_MOMENT_Y = 0.00
82+
REF_ORIGIN_MOMENT_Z = 0.00
83+
%
84+
% Reference length for moment non-dimensional coefficients (m or in)
85+
REF_LENGTH= 1.0
86+
%
87+
% Reference area for non-dimensional force coefficients (0 implies automatic
88+
% calculation) (m^2 or in^2)
89+
REF_AREA= 3.0
90+
%
91+
% Aircraft semi-span (0 implies automatic calculation) (m or in)
92+
SEMI_SPAN= 1.5
93+
%
94+
% Flow non-dimensionalization (DIMENSIONAL, FREESTREAM_PRESS_EQ_ONE,
95+
% FREESTREAM_VEL_EQ_MACH, FREESTREAM_VEL_EQ_ONE)
96+
REF_DIMENSIONALIZATION= DIMENSIONAL
97+
98+
% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
99+
%
100+
% Euler wall boundary marker(s) (NONE = no marker)
101+
% Implementation identical to MARKER_SYM.
102+
MARKER_EULER= ( wing, tip)
103+
%
104+
% Far-field boundary marker(s) (NONE = no marker)
105+
MARKER_FAR= ( far_away )
106+
107+
% ------------------------ SURFACES IDENTIFICATION ----------------------------%
108+
%
109+
% Marker(s) of the surface in the surface flow solution file
110+
MARKER_PLOTTING = ( wing, tip )
111+
%
112+
% Marker(s) of the surface where the non-dimensional coefficients are evaluated.
113+
MARKER_MONITORING = ( wing, tip )
114+
115+
% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
116+
%
117+
% CFL number (initial value for the adaptive CFL number)
118+
CFL_NUMBER= 1000.0
119+
%
120+
% Adaptive CFL number (NO, YES)
121+
CFL_ADAPT= YES
122+
%
123+
% -------------------------- MULTIGRID PARAMETERS -----------------------------%
124+
%
125+
% Multi-grid levels (0 = no multi-grid)
126+
MGLEVEL= 3
127+
%
128+
% Multi-grid cycle (V_CYCLE, W_CYCLE, FULLMG_CYCLE)
129+
MGCYCLE= W_CYCLE
130+
131+
% -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
132+
%
133+
% Convective numerical method (JST, JST_KE, JST_MAT, LAX-FRIEDRICH, CUSP, ROE, AUSM,
134+
% AUSMPLUSUP, AUSMPLUSUP2, AUSMPWPLUS, HLLC, TURKEL_PREC,
135+
% SW, MSW, FDS, SLAU, SLAU2, L2ROE, LMROE)
136+
CONV_NUM_METHOD_FLOW= JST
137+
%
138+
% Max number of iterations of the linear solver for the implicit formulation
139+
LINEAR_SOLVER_ITER= 20
140+
%
141+
% ------------------------- SCREEN/HISTORY VOLUME OUTPUT --------------------------%
142+
%
143+
% Screen output fields (use 'SU2_CFD -d <config_file>' to view list of available fields)
144+
SCREEN_OUTPUT= (TIME_ITER, INNER_ITER, RMS_DENSITY, LIFT, MOMENT_X, MOMENT_Y, MOMENT_Z)
145+
%
146+
% History output groups (use 'SU2_CFD -d <config_file>' to view list of available fields)
147+
HISTORY_OUTPUT= (ITER, RMS_RES, AERO_COEFF, CAUCHY, WALL_TIME)
148+
%
149+
% ------------------------- INPUT/OUTPUT FILE INFORMATION --------------------------%
150+
%
151+
% Mesh input file
152+
MESH_FILENAME= mesh_rectangular_wing.su2
153+
%
154+
% Mesh input file format (SU2, CGNS)
155+
MESH_FORMAT= SU2
156+
%
157+
% Restart flow input file
158+
SOLUTION_FILENAME= restart_gust.dat
159+
RESTART_FILENAME= restart_gust.dat
160+
SURFACE_FILENAME= surface_gust
161+
VOLUME_FILENAME= volume_gust
162+
CONV_FILENAME= history_gust
163+
%
164+
% Output tabular file format (TECPLOT, CSV)
165+
TABULAR_FORMAT= CSV
166+
%
167+
% Files to output
168+
% Possible formats : (TECPLOT, TECPLOT_BINARY, SURFACE_TECPLOT,
169+
% SURFACE_TECPLOT_BINARY, CSV, SURFACE_CSV, PARAVIEW, PARAVIEW_BINARY, SURFACE_PARAVIEW,
170+
% SURFACE_PARAVIEW_BINARY, MESH, RESTART_BINARY, RESTART_ASCII, CGNS, SURFACE_CGNS, STL)
171+
OUTPUT_FILES= (RESTART, RESTART_ASCII, TECPLOT, SURFACE_TECPLOT)
172+

TestCases/hybrid_regression.py

Lines changed: 9 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -478,6 +478,15 @@ def main():
478478
sine_gust.test_vals = [-1.977520, 3.481804, -0.012402, -0.007454]
479479
sine_gust.unsteady = True
480480
test_list.append(sine_gust)
481+
482+
# Cosine gust in z-direction
483+
cosine_gust = TestCase('cosine_gust_zdir')
484+
cosine_gust.cfg_dir = "gust"
485+
cosine_gust.cfg_file = "cosine_gust_zdir.cfg"
486+
cosine_gust.test_iter = 79
487+
cosine_gust.test_vals = [-2.418813, 0.004650, -0.001878, -0.000637, -0.000271]
488+
cosine_gust.unsteady = True
489+
test_list.append(cosine_gust)
481490

482491
# Aeroelastic
483492
aeroelastic = TestCase('aeroelastic')

0 commit comments

Comments
 (0)