Skip to content

Commit 1ac7e15

Browse files
authored
Merge pull request #9 from AlbertoGuardone/patch-1
NICFD tutorial
2 parents 9038ab4 + ab0042d commit 1ac7e15

File tree

10 files changed

+10668
-0
lines changed

10 files changed

+10668
-0
lines changed

NICFD_nozzle/NICFD_nozzle.cfg

Lines changed: 260 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,260 @@
1+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2+
% %
3+
% SU2 configuration file %
4+
% Case description: Non-ideal compressible fluid flow in a converging- %
5+
% diverging supersonic nozzle for siloxane fluid MDM %
6+
% Author: Alberto Guardone %
7+
% Institution: Politecnico di Milano %
8+
% Date: 2019.05.03 %
9+
% File Version 6.2.0 "Falcon" %
10+
% %
11+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
12+
13+
% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
14+
%
15+
% Physical governing equations (EULER, NAVIER_STOKES,
16+
% FEM_EULER, FEM_NAVIER_STOKES, FEM_RANS, FEM_LES,
17+
% WAVE_EQUATION, HEAT_EQUATION, FEM_ELASTICITY,
18+
% POISSON_EQUATION)
19+
PHYSICAL_PROBLEM= NAVIER_STOKES
20+
%
21+
% Specify turbulence model (NONE, SA, SA_NEG, SST, SA_E, SA_COMP, SA_E_COMP)
22+
KIND_TURB_MODEL= SST
23+
%
24+
% Mathematical problem (DIRECT, CONTINUOUS_ADJOINT, DISCRETE_ADJOINT)
25+
MATH_PROBLEM= DIRECT
26+
%
27+
% Regime type (COMPRESSIBLE, INCOMPRESSIBLE)
28+
REGIME_TYPE= COMPRESSIBLE
29+
%
30+
% Restart solution (NO, YES)
31+
RESTART_SOL= NO
32+
%
33+
% System of measurements (SI, US)
34+
% International system of units (SI): ( meters, kilograms, Kelvins,
35+
% Newtons = kg m/s^2, Pascals = N/m^2,
36+
% Density = kg/m^3, Speed = m/s,
37+
% Equiv. Area = m^2 )
38+
% United States customary units (US): ( inches, slug, Rankines, lbf = slug ft/s^2,
39+
% psf = lbf/ft^2, Density = slug/ft^3,
40+
% Speed = ft/s, Equiv. Area = ft^2 )
41+
SYSTEM_MEASUREMENTS= SI
42+
%
43+
% -------------------- COMPRESSIBLE FREE-STREAM DEFINITION --------------------%
44+
%
45+
% Mach number (non-dimensional, based on the free-stream values)
46+
MACH_NUMBER= 1E-9
47+
%
48+
% Angle of attack (degrees, only for compressible flows)
49+
AOA= 0.0
50+
%
51+
% Side-slip angle (degrees, only for compressible flows)
52+
SIDESLIP_ANGLE= 0.0
53+
%
54+
% Init option to choose between Reynolds (default) or thermodynamics quantities
55+
% for initializing the solution (REYNOLDS, TD_CONDITIONS)
56+
INIT_OPTION= TD_CONDITIONS
57+
%
58+
% Free-stream option to choose between density and temperature (default) for
59+
% initializing the solution (TEMPERATURE_FS, DENSITY_FS)
60+
FREESTREAM_OPTION= TEMPERATURE_FS
61+
%
62+
% Free-stream pressure (101325.0 N/m^2, 2116.216 psf by default)
63+
FREESTREAM_PRESSURE= 904388
64+
%
65+
% Free-stream temperature (288.15 K, 518.67 R by default)
66+
FREESTREAM_TEMPERATURE= 542.13
67+
%
68+
% Compressible flow non-dimensionalization (DIMENSIONAL, FREESTREAM_PRESS_EQ_ONE,
69+
% FREESTREAM_VEL_EQ_MACH, FREESTREAM_VEL_EQ_ONE)
70+
REF_DIMENSIONALIZATION= DIMENSIONAL
71+
72+
% ---- IDEAL GAS, POLYTROPIC, VAN DER WAALS AND PENG ROBINSON CONSTANTS -------%
73+
%
74+
% Fluid model (STANDARD_AIR, IDEAL_GAS, VW_GAS, PR_GAS,
75+
% CONSTANT_DENSITY, INC_IDEAL_GAS, INC_IDEAL_GAS_POLY)
76+
FLUID_MODEL= PR_GAS
77+
%
78+
% Ratio of specific heats (1.4 default and the value is hardcoded
79+
% for the model STANDARD_AIR, compressible only)
80+
GAMMA_VALUE= 1.01767
81+
%
82+
% Specific gas constant (287.058 J/kg*K default and this value is hardcoded
83+
% for the model STANDARD_AIR, compressible only)
84+
GAS_CONSTANT= 35.17
85+
%
86+
% Critical Temperature (131.00 K by default)
87+
CRITICAL_TEMPERATURE= 565.3609
88+
%
89+
% Critical Pressure (3588550.0 N/m^2 by default)
90+
CRITICAL_PRESSURE= 1437500
91+
%
92+
% Acentric factor (0.035 (air))
93+
ACENTRIC_FACTOR= 0.524
94+
95+
% --------------------------- VISCOSITY MODEL ---------------------------------%
96+
%
97+
% Viscosity model (SUTHERLAND, CONSTANT_VISCOSITY, POLYNOMIAL_VISCOSITY).
98+
VISCOSITY_MODEL= CONSTANT_VISCOSITY
99+
%
100+
% Molecular Viscosity that would be constant (1.716E-5 by default)
101+
MU_CONSTANT= 1.21409E-05
102+
103+
% --------------------------- THERMAL CONDUCTIVITY MODEL ----------------------%
104+
%
105+
% Laminar Conductivity model (CONSTANT_CONDUCTIVITY, CONSTANT_PRANDTL,
106+
% POLYNOMIAL_CONDUCTIVITY).
107+
CONDUCTIVITY_MODEL= CONSTANT_CONDUCTIVITY
108+
%
109+
% Molecular Thermal Conductivity that would be constant (0.0257 by default)
110+
KT_CONSTANT= 0.030542828
111+
112+
% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
113+
%
114+
% Navier-Stokes (no-slip), constant heat flux wall marker(s) (NONE = no marker)
115+
% Format: ( marker name, constant heat flux (J/m^2), ... )
116+
MARKER_HEATFLUX= ( WALL, 0.0 )
117+
%
118+
% Symmetry boundary marker(s) (NONE = no marker)
119+
MARKER_SYM= ( SYMMETRY )
120+
%
121+
% Riemann boundary marker(s) (NONE = no marker)
122+
% Format: (marker, data kind flag, list of data)
123+
MARKER_RIEMANN= ( INFLOW, TOTAL_CONDITIONS_PT, 904388, 542.13, 1.0, 0.0, 0.0, OUTFLOW, STATIC_PRESSURE, 200000.0, 0.0, 0.0, 0.0, 0.0 )
124+
125+
% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
126+
%
127+
% Numerical method for spatial gradients (GREEN_GAUSS, WEIGHTED_LEAST_SQUARES)
128+
NUM_METHOD_GRAD= GREEN_GAUSS
129+
%
130+
% CFL number (initial value for the adaptive CFL number)
131+
CFL_NUMBER= 10.0
132+
%
133+
% Adaptive CFL number (NO, YES)
134+
CFL_ADAPT= YES
135+
%
136+
% Parameters of the adaptive CFL number (factor down, factor up, CFL min value,
137+
% CFL max value )
138+
CFL_ADAPT_PARAM= ( 1.5, 0.5, 10.0, 1000.0 )
139+
%
140+
% Maximum Delta Time in local time stepping simulations
141+
MAX_DELTA_TIME= 1E6
142+
143+
% ----------- SLOPE LIMITER AND DISSIPATION SENSOR DEFINITION -----------------%
144+
%
145+
% Monotonic Upwind Scheme for Conservation Laws (TVD) in the flow equations.
146+
% Required for 2nd order upwind schemes (NO, YES)
147+
MUSCL_FLOW= YES
148+
%
149+
% Slope limiter (NONE, VENKATAKRISHNAN, VENKATAKRISHNAN_WANG,
150+
% BARTH_JESPERSEN, VAN_ALBADA_EDGE)
151+
SLOPE_LIMITER_FLOW= NONE
152+
%
153+
% Monotonic Upwind Scheme for Conservation Laws (TVD) in the turbulence equations.
154+
% Required for 2nd order upwind schemes (NO, YES)
155+
MUSCL_TURB= NO
156+
157+
% ------------------------ LINEAR SOLVER DEFINITION ---------------------------%
158+
%
159+
% Linear solver or smoother for implicit formulations (BCGSTAB, FGMRES, SMOOTHER_JACOBI,
160+
% SMOOTHER_ILU, SMOOTHER_LUSGS,
161+
% SMOOTHER_LINELET)
162+
LINEAR_SOLVER= FGMRES
163+
%
164+
% Preconditioner of the Krylov linear solver (ILU, LU_SGS, LINELET, JACOBI)
165+
LINEAR_SOLVER_PREC= ILU
166+
%
167+
% Linael solver ILU preconditioner fill-in level (0 by default)
168+
LINEAR_SOLVER_ILU_FILL_IN= 0
169+
%
170+
% Minimum error of the linear solver for implicit formulations
171+
LINEAR_SOLVER_ERROR= 1E-6
172+
%
173+
% Max number of iterations of the linear solver for the implicit formulation
174+
LINEAR_SOLVER_ITER= 5
175+
176+
% -------------------------- MULTIGRID PARAMETERS -----------------------------%
177+
%
178+
% Multi-grid levels (0 = no multi-grid)
179+
MGLEVEL= 0
180+
181+
% -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
182+
%
183+
% Convective numerical method (JST, LAX-FRIEDRICH, CUSP, ROE, AUSM, AUSMPLUSUP, AUSMPLUSUP2, HLLC,
184+
% TURKEL_PREC, MSW, FDS)
185+
CONV_NUM_METHOD_FLOW= ROE
186+
%
187+
% Entropy fix coefficient (0.0 implies no entropy fixing, 1.0 implies scalar
188+
% artificial dissipation)
189+
ENTROPY_FIX_COEFF= 0.1
190+
%
191+
% Time discretization (RUNGE-KUTTA_EXPLICIT, EULER_IMPLICIT, EULER_EXPLICIT)
192+
TIME_DISCRE_FLOW= EULER_IMPLICIT
193+
%
194+
% Relaxation coefficient
195+
RELAXATION_FACTOR_FLOW= 0.95
196+
197+
% -------------------- TURBULENT NUMERICAL METHOD DEFINITION ------------------%
198+
%
199+
% Convective numerical method (SCALAR_UPWIND)
200+
CONV_NUM_METHOD_TURB= SCALAR_UPWIND
201+
%
202+
% Time discretization (EULER_IMPLICIT)
203+
TIME_DISCRE_TURB= EULER_IMPLICIT
204+
%
205+
% Reduction factor of the CFL coefficient in the turbulence problem
206+
CFL_REDUCTION_TURB= 1.0
207+
%
208+
% Relaxation coefficient
209+
RELAXATION_FACTOR_TURB= 0.95
210+
211+
% --------------------------- CONVERGENCE PARAMETERS --------------------------%
212+
%
213+
% Number of total iterations
214+
EXT_ITER= 1000
215+
%
216+
% Convergence criteria (CAUCHY, RESIDUAL)
217+
%
218+
CONV_CRITERIA= RESIDUAL
219+
%
220+
% Residual reduction (order of magnitude with respect to the initial value)
221+
RESIDUAL_REDUCTION= 6
222+
%
223+
% Min value of the residual (log10 of the residual)
224+
RESIDUAL_MINVAL= -24
225+
%
226+
% Start convergence criteria at iteration number
227+
STARTCONV_ITER= 10
228+
229+
% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
230+
%
231+
% Mesh input file
232+
MESH_FILENAME= NICFD_nozzle.su2
233+
%
234+
% Mesh input file format (SU2, CGNS)
235+
MESH_FORMAT= SU2
236+
%
237+
% Mesh output file
238+
MESH_OUT_FILENAME= mesh_out.su2
239+
%
240+
% Restart flow input file
241+
SOLUTION_FLOW_FILENAME= solution_flow.dat
242+
%
243+
% Output file format (TECPLOT, TECPLOT_BINARY, PARAVIEW, PARAVIEW_BINARY,
244+
% FIELDVIEW, FIELDVIEW_BINARY)
245+
OUTPUT_FORMAT= PARAVIEW
246+
%
247+
% Output file convergence history (w/o extension)
248+
CONV_FILENAME= history
249+
%
250+
% Output file restart flow
251+
RESTART_FLOW_FILENAME= restart_flow.dat
252+
%
253+
% Output file flow (w/o extension) variables
254+
VOLUME_FLOW_FILENAME= flow
255+
%
256+
% Output file surface flow coefficient (w/o extension)
257+
SURFACE_FLOW_FILENAME= surface_flow
258+
%
259+
% Writing solution file frequency
260+
WRT_SOL_FREQ= 1000

0 commit comments

Comments
 (0)