|
[Sponsors] |
December 2, 2015, 13:03 |
SU2 Validation for Hypersonic Flows
|
#1 |
New Member
Mike
Join Date: Oct 2015
Posts: 11
Rep Power: 11 |
I'm working towards my masters thesis looking at a hypersonic (M=5) ramjet inlet and I'm trying to validate the SU2 code against some wind tunnel data we have taken. The ultimate goal is to run an optimisation study and assess if the code can be used to replace some of the commercial codes we have access to such as Fluent and Cobalt. I'd be very happy to share my findings later on.
Unfortunately, after using the code for a couple weeks we've found it very difficult to achieve convergence - it frequently diverges and exits even with a simplified axisymmetric model after only ~20-50 iterations. The grid works well with Cobalt and the results match very well with the experiment. Therefore it must be either a code/cfg file problem. Could any developers or people with experience using the code please take a look at the settings I'm using and suggest any potential improvements. Cheers, Mike. Code:
% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------% % % Physical governing equations (EULER, NAVIER_STOKES, % TNE2_EULER, TNE2_NAVIER_STOKES, % WAVE_EQUATION, HEAT_EQUATION, LINEAR_ELASTICITY, % POISSON_EQUATION) PHYSICAL_PROBLEM= NAVIER_STOKES % % Specify turbulent model (NONE, SA, SA_NEG, SST) KIND_TURB_MODEL= SST % % Mathematical problem (DIRECT, ADJOINT, LINEARIZED) MATH_PROBLEM= DIRECT % % Restart solution (NO, YES) RESTART_SOL= NO % -------------------- COMPRESSIBLE FREE-STREAM DEFINITION --------------------% % % Mach number (non-dimensional, based on the free-stream values) MACH_NUMBER= 5.0 % % Angle of attack (degrees) AoA= 0.0 % % Free-stream temperature (288.15 K by default) FREESTREAM_TEMPERATURE= 223.25 % % Reynolds number (non-dimensional, based on the free-stream values) REYNOLDS_NUMBER= 42.4858E6 % % Reynolds length (1 m by default) REYNOLDS_LENGTH= 1.0 % ---- IDEAL GAS, POLYTROPIC, VAN DER WAALS AND PENG ROBINSON CONSTANTS -------% % % Different gas model (STANDARD_AIR, IDEAL_GAS, VW_GAS, PR_GAS) FLUID_MODEL= STANDARD_AIR % % Ratio of specific heats (1.4 default and the value is hardcoded % for the model STANDARD_AIR) GAMMA_VALUE= 1.4 % % Specific gas constant (287.058 J/kg*K default and this value is hardcoded % for the model STANDARD_AIR) GAS_CONSTANT= 287.058 % % Critical Temperature (131.00 K by default) CRITICAL_TEMPERATURE= 131.00 % % Critical Pressure (3588550.0 N/m^2 by default) CRITICAL_PRESSURE= 3588550.0 % % Critical Density (263.0 Kg/m3 by default) CRITICAL_DENSITY= 263.0 % % Acentric factor (0.035 (air)) ACENTRIC_FACTOR= 0.035 % --------------------------- VISCOSITY MODEL ---------------------------------% % % Viscosity model (SUTHERLAND, CONSTANT_VISCOSITY). VISCOSITY_MODEL= SUTHERLAND % % Molecular Viscosity that would be constant (1.716E-5 by default) MU_CONSTANT= 1.716E-5 % % Sutherland Viscosity Ref (1.716E-5 default value for AIR SI) MU_REF= 1.716E-5 % % Sutherland Temperature Ref (273.15 K default value for AIR SI) MU_T_REF= 273.15 % % Sutherland constant (110.4 default value for AIR SI) SUTHERLAND_CONSTANT= 110.4 % ---------------------- 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_MOMENT= 1.0 % % Reference area for force coefficients (0 implies automatic calculation) REF_AREA= 1.0 % % Reference element length for computing the slope limiter epsilon REF_ELEM_LENGTH= 0.000007 % -------------------- BOUNDARY CONDITION DEFINITION --------------------------% % % Navier-Stokes wall boundary marker(s) (NONE = no marker) MARKER_HEATFLUX= ( WALL_UPPER, 0.0, LOWER_WALL 0.0 ) % % Supersonic inlet boundary marker(s) (NONE = no marker) % Total Conditions: (inlet marker, temperature, static pressure, velocity_x, % velocity_y, velocity_z, ... ), i.e. all variables specified. MARKER_SUPERSONIC_INLET= ( INLET, 223.25, 26500.0, 1497.513356, 0.0, 0.0 ) % MARKER_SYM= ( AXIS ) % AXISYMMETRIC= YES MARKER_SUPERSONIC_OUTLET= ( OUTLET, FARFIELD ) % % Marker(s) of the surface to be plotted or designed MARKER_PLOTTING= ( LOWER_WALL ) % % Marker(s) of the surface where the functional (Cd, Cl, etc.) will be evaluated MARKER_MONITORING= ( LOWER_WALL ) % ------------- 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 % % 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.25, 2, 100 ) % % Number of total iterations EXT_ITER= 2000 % % ------------------------ LINEAR SOLVER DEFINITION ---------------------------% % % Linear solver for the implicit (or discrete adjoint) formulation (BCGSTAB, FGMRES) LINEAR_SOLVER= FGMRES % % Min error of the linear solver for the implicit formulation LINEAR_SOLVER_ERROR= 1E-6 % % Max number of iterations of the linear solver for the implicit formulation LINEAR_SOLVER_ITER= 10 % -------------------------- MULTIGRID PARAMETERS -----------------------------% % % Multi-Grid Levels (0 = no multi-grid) MGLEVEL= 0 % % Multi-grid cycle (V_CYCLE, W_CYCLE, FULLMG_CYCLE) MGCYCLE= W_CYCLE % % Multi-grid pre-smoothing level MG_PRE_SMOOTH= ( 1, 2, 2, 2 ) % % Multi-grid post-smoothing level MG_POST_SMOOTH= ( 2, 2, 2, 2 ) % % Jacobi implicit smoothing of the correction MG_CORRECTION_SMOOTH= ( 0, 0, 0, 0 ) % % Damping factor for the residual restriction MG_DAMP_RESTRICTION= 0.7 % % Damping factor for the correction prolongation MG_DAMP_PROLONGATION= 0.7 % -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------% % % Convective numerical method (JST, LAX-FRIEDRICH, CUSP, ROE, AUSM, HLLC, % TURKEL_PREC, MSW) CONV_NUM_METHOD_FLOW= AUSM % % Spatial numerical order integration (1ST_ORDER, 2ND_ORDER, 2ND_ORDER_LIMITER) % SPATIAL_ORDER_FLOW= 2ND_ORDER % % Slope limiter (VENKATAKRISHNAN) SLOPE_LIMITER_FLOW= VENKATAKRISHNAN % % Coefficient for the limiter (smooth regions) LIMITER_COEFF= 100 % % 1st, 2nd and 4th order artificial dissipation coefficients AD_COEFF_FLOW= ( 0.15, 0.5, 0.02 ) % % Time discretization (RUNGE-KUTTA_EXPLICIT, EULER_IMPLICIT, EULER_EXPLICIT) TIME_DISCRE_FLOW= EULER_IMPLICIT % -------------------- TURBULENT NUMERICAL METHOD DEFINITION ------------------% % % Convective numerical method (SCALAR_UPWIND) CONV_NUM_METHOD_TURB= SCALAR_UPWIND % % Spatial numerical order integration (1ST_ORDER, 2ND_ORDER, 2ND_ORDER_LIMITER) % SPATIAL_ORDER_TURB= 1ST_ORDER % % Slope limiter (VENKATAKRISHNAN, BARTH) SLOPE_LIMITER_TURB= VENKATAKRISHNAN % % Time discretization (EULER_IMPLICIT) TIME_DISCRE_TURB= EULER_IMPLICIT % --------------------------- CONVERGENCE PARAMETERS --------------------------% % % Convergence criteria (CAUCHY, RESIDUAL) % CONV_CRITERIA= RESIDUAL % % Residual reduction (order of magnitude with respect to the initial value) RESIDUAL_REDUCTION= 4 % % Min value of the residual (log10 of the residual) RESIDUAL_MINVAL= -10 % % Start convergence 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 % % Function to apply the criteria (LIFT, DRAG, NEARFIELD_PRESS, SENS_GEOMETRY, % SENS_MACH, DELTA_LIFT, DELTA_DRAG) CAUCHY_FUNC_FLOW= DRAG |
|
December 7, 2015, 17:53 |
|
#2 |
Senior Member
Heather Kline
Join Date: Jun 2013
Posts: 309
Rep Power: 14 |
Thanks for your question.
The first thing you should do is lower the CFL number (it is currently 10, which is high unless you have a restarted solution or a very stable problem - 1.0 is more common). If you want to have higher CFL numbers as the solution continues, set CFL_ADAPT =YES. |
|
April 8, 2016, 00:03 |
|
#3 |
Member
Paul Zhang
Join Date: Feb 2011
Posts: 44
Rep Power: 15 |
Hello,
I am modeling hypersonic flows using SU2. I was wondering if you've solved the issue. Paul |
|
Tags |
convergence, hypersonic, inlet |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
interFoam - validation for bubble/droplet flows in microfluidics | duongquaphim | OpenFOAM Verification & Validation | 54 | February 14, 2023 00:46 |
Reacting Flows in SU2? | rking | SU2 | 3 | May 16, 2016 11:55 |
CFX problem in ubuntu (linux) | Vigneshramaero | CFX | 0 | July 13, 2012 11:22 |
CFX-Pre problem, pls help!!! | cth_yao | CFX | 0 | February 17, 2012 01:52 |
Validation tests for 3-D flows | Alexey | Main CFD Forum | 1 | January 6, 2000 00:48 |