[themrdjj]

I'm trying to run a 3D simulation for the first time (compressible Euler, steady state). The example/tutorial cases work fine. The mesh also works fine, no errors there. However, SU2 crashes right after "Begin Solver". It does so without displaying any error messages and I can't find any log files either.

What's the approach to debugging such an issue? It's hard to reduce the problem to something simpler since the mesh is quite complex.

Also, there is one thing that I was unsure about while setting up the config file: I need my freestream flow to be along the y-axism not the x-axis, so I set the sideslip angle to 90 deg. Is this a good approach?

[jcownbey]

If SU2 is crashing before the first iteration, then either something is wrong in the config file or you have a meshing issue. If you need to the free stream to be in the y-direction, it may be easier to reorient your mesh so the flow is in the x-direction. SU2 should have an option to change the freestream vector (e.g. flow from [1 0 0] to [0 1 0]).

[themrdjj]

The mesh itself is fine I think, also when running SU2_MSH there are no errors. Reorienting the mesh would be a lot of work. Also, SU2 displays the following, which makes me think that it's fine.

-- Initial and free-stream conditions:
+------------------------------------------------------------------------------+
|                  Name|    Dim. value|    Ref. value|      Unit|Non-dim. value|
+------------------------------------------------------------------------------+
|       Static Pressure|       23842.3|             1|        Pa|       23842.3|
|               Density|       0.37959|             1|    kg/m^3|       0.37959|
|           Temperature|       218.808|             1|         K|       218.808|
|          Total Energy|        181103|             1|   m^2/s^2|        181103|
|            Velocity-X|  1.34363e-014|             1|       m/s|  1.34363e-014|
|            Velocity-Y|       219.438|             1|       m/s|       219.438|
|            Velocity-Z|             0|             1|       m/s|             0|
|    Velocity Magnitude|       219.438|             1|       m/s|       219.438|
+------------------------------------------------------------------------------+
|           Mach Number|             -|             -|         -|          0.74|
+------------------------------------------------------------------------------+

Is there really no way to see where and why exactly it stops running?

[themrdjj]

I have now tried simulating an empty channel using the same .cfg file (of course adapted to the new mesh). It still crashes.

Perhaps someone could tell me what could be wrong here? Here, the mesh is aligned with the x-axis.

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%                                                                              %
% SU2 configuration file                                                       %
% Case description: Empty channel                                              %
% Author:                                                                      %
% Institution:                                                                 %
% Date: __________                                                             %
% File Version 7.0.0 "Blackbird"                                               %
%                                                                              %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
%
% Solver type (EULER, NAVIER_STOKES, RANS,
%                               INC_EULER, INC_NAVIER_STOKES, INC_RANS
%                               FEM_EULER, FEM_NAVIER_STOKES, FEM_RANS, FEM_LES,
%                               HEAT_EQUATION_FVM, ELASTICITY)
SOLVER= EULER
%
% Mathematical problem (DIRECT, CONTINUOUS_ADJOINT, DISCRETE_ADJOINT)
MATH_PROBLEM= DIRECT
%
% Restart solution (NO, YES)
RESTART_SOL= NO
%
% System of measurements (SI, US)
SYSTEM_MEASUREMENTS= SI


% -------------------- COMPRESSIBLE FREE-STREAM DEFINITION --------------------%
%
% Mach number (non-dimensional, based on the free-stream values)
MACH_NUMBER= 0.74
%
% Angle of attack (degrees, only for compressible flows)
AOA= 0
%
% Side-slip angle (degrees, only for compressible flows)
SIDESLIP_ANGLE= 0
%
% Free-stream pressure (101325.0 N/m^2, 2116.216 psf by default)
%FREESTREAM_PRESSURE= 23842.3
%
% Free-stream temperature (288.15 K, 518.67 R by default)
%FREESTREAM_TEMPERATURE= 218.808

% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
%
% Euler wall boundary marker(s) (NONE = no marker)
% Implementation identical to MARKER_SYM.
MARKER_EULER= ( walls )
%
%
% Far-field boundary marker(s) (NONE = no marker)
%MARKER_FAR= ( inlet, outlet, euler_wall, side_wall, top_bottom_wall )
%
% Symmetry boundary marker(s) (NONE = no marker)
% Implementation identical to MARKER_EULER.
%MARKER_SYM= ( symmetry )
%
% Inlet boundary type (TOTAL_CONDITIONS, MASS_FLOW)
INLET_TYPE= TOTAL_CONDITIONS
%
% Inlet boundary marker(s) with the following formats (NONE = no marker)
% Total Conditions: (inlet marker, total temp, total pressure, flow_direction_x,
%           flow_direction_y, flow_direction_z, ... ) where flow_direction is
%           a unit vector.
% Mass Flow: (inlet marker, density, velocity magnitude, flow_direction_x,
%           flow_direction_y, flow_direction_z, ... ) where flow_direction is
%           a unit vector.
MARKER_INLET= ( inlet, 242.77, 34300.5 , 1.0, 0, 0)
%
% Outlet boundary marker(s) (NONE = no marker)
% Compressible: ( outlet marker, back pressure (static thermodynamic), ... )
MARKER_OUTLET= ( outlet, 23842.3)

% ------------------------- INPUT/OUTPUT FILE INFORMATION --------------------------%
%
% Mesh input file
MESH_FILENAME= channel_empty_x.su2
%
% Mesh input file format (SU2, CGNS)
MESH_FORMAT= SU2
%
% Output tabular file format (TECPLOT, CSV)
TABULAR_FORMAT= TECPLOT
%
% Files to output
OUTPUT_FILES= (RESTART, PARAVIEW, SURFACE_PARAVIEW)

This is the output:

-------------------------------------------------------------------------
|    ___ _   _ ___                                                      |
|   / __| | | |_  )   Release 7.0.0  "Blackbird"                        |
|   \__ \ |_| |/ /                                                      |
|   |___/\___//___|   Suite (Computational Fluid Dynamics Code)         |
|                                                                       |
-------------------------------------------------------------------------
| SU2 Project Website: https://su2code.github.io                        |
|                                                                       |
| The SU2 Project is maintained by the SU2 Foundation                   |
| (http://su2foundation.org)                                            |
-------------------------------------------------------------------------
| Copyright 2012-2019, SU2 Contributors                                 |
|                                                                       |
| SU2 is free software; you can redistribute it and/or                  |
| modify it under the terms of the GNU Lesser General Public            |
| License as published by the Free Software Foundation; either          |
| version 2.1 of the License, or (at your option) any later version.    |
|                                                                       |
| SU2 is distributed in the hope that it will be useful,                |
| but WITHOUT ANY WARRANTY; without even the implied warranty of        |
| MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU      |
| Lesser General Public License for more details.                       |
|                                                                       |
| You should have received a copy of the GNU Lesser General Public      |
| License along with SU2. If not, see <http://www.gnu.org/licenses/>.   |
-------------------------------------------------------------------------

Parsing config file for zone 0

----------------- Physical Case Definition ( Zone 0 ) -------------------
Compressible Euler equations.
Mach number: 0.74.
Angle of attack (AoA): 0 deg, and angle of sideslip (AoS): 0 deg.
No restart solution, use the values at infinity (freestream).
Dimensional simulation.
The reference area is 1 m^2.
The semi-span will be computed using the max y(3D) value.
The reference length is 1 m.
Input mesh file name: channel_empty_x.su2

--------------- Space Numerical Integration ( Zone 0 ) ------------------
Gradient for upwind reconstruction: inverse-distance weighted Least-Squares.
Gradient for viscous and source terms: inverse-distance weighted Least-Squares.

--------------- Time Numerical Integration  ( Zone 0 ) ------------------
Local time stepping (steady state simulation).
Euler implicit method for the flow equations.
FGMRES is used for solving the linear system.
Using a ILU(0) preconditioning.
Convergence criteria of the linear solver: 1e-006.
Max number of linear iterations: 10.
No CFL adaptation.
Courant-Friedrichs-Lewy number:     1.25

------------------ Convergence Criteria  ( Zone 0 ) ---------------------
Maximum number of solver subiterations: 1000.
Begin convergence monitoring at iteration 5.
Residual minimum value: 1e-14.
Cauchy series min. value: 1e-010.
Number of Cauchy elements: 100.

-------------------- Output Information ( Zone 0 ) ----------------------
Writing solution files every 250 iterations.
Writing the convergence history file every 1 inner iterations.
Writing the screen convergence history every 1 inner iterations.
The tabular file format is Tecplot (.dat).
Convergence history file name: history.
Forces breakdown file name: forces_breakdown.dat.
Surface file name: surface.
Volume file name: vol_solution.
Restart file name: restart.dat.

------------- Config File Boundary Information ( Zone 0 ) ---------------
+-----------------------------------------------------------------------+
|                        Marker Type|                        Marker Name|
+-----------------------------------------------------------------------+
|                         Euler wall|                              walls|
+-----------------------------------------------------------------------+
|                     Inlet boundary|                              inlet|
+-----------------------------------------------------------------------+
|                    Outlet boundary|                             outlet|
+-----------------------------------------------------------------------+

-------------------- Output Preprocessing ( Zone 0 ) --------------------
Euler/Navier-Stokes/RANS output structure.
Screen output fields: INNER_ITER, RMS_DENSITY, RMS_MOMENTUM-X, RMS_MOMENTUM-Y, RMS_ENERGY
History output group(s): ITER, RMS_RES
Convergence field(s): RMS_DENSITY
Volume output fields: COORDINATES, SOLUTION, PRIMITIVE

------------------- Geometry Preprocessing ( Zone 0 ) -------------------
Three dimensional problem.
3248 grid points.
12732 volume elements.
3 surface markers.
3634 boundary elements in index 0 (Marker = walls).
214 boundary elements in index 1 (Marker = inlet).
228 boundary elements in index 2 (Marker = outlet).
12732 tetrahedra.
Setting point connectivity.
Renumbering points (Reverse Cuthill McKee Ordering).
Recomputing point connectivity.
Setting element connectivity.
Checking the numerical grid orientation.
There has been a re-orientation of the TETRAHEDRON volume elements.
There has been a re-orientation of the TRIANGLE surface elements.
Identifying edges and vertices.
Computing centers of gravity.
Setting the control volume structure.
Volume of the computational grid: 500.
Searching for the closest normal neighbors to the surfaces.
Storing a mapping from global to local point index.
Compute the surface curvature.
Max K: 27.2517. Mean K: 1.96196. Standard deviation K: 4.87474.
Checking for periodicity.
Computing mesh quality statistics for the dual control volumes.
+--------------------------------------------------------------+
|           Mesh Quality Metric|        Minimum|        Maximum|
+--------------------------------------------------------------+
|    Orthogonality Angle (deg.)|        34.5549|        86.0477|
|     CV Face Area Aspect Ratio|        1.11965|        20.6257|
|           CV Sub-Volume Ratio|        1.00418|        25.1827|
+--------------------------------------------------------------+
Finding max control volume width.
Semi-span length = -1e+010 m.
Wetted area = 0 m^2.
Area projection in the x-plane = 0 m^2, y-plane = 0 m^2, z-plane = 0 m^2.
Max. coordinate in the x-direction = -1e+010 m, y-direction = -1e+010 m, z-direction = -1e+010 m.
Min. coordinate in the x-direction = 1e+010 m, y-direction = 1e+010 m, z-direction = 1e+010 m.
Checking if Euler & Symmetry markers are straight/plane:
Boundary marker walls is NOT a single plane.

-------------------- Solver Preprocessing ( Zone 0 ) --------------------
Inviscid flow: Computing density based on free-stream
temperature and pressure using the ideal gas law.
Force coefficients computed using free-stream values.

-- Models:
+------------------------------------------------------------------------------+
|          Viscosity Model|        Conductivity Model|              Fluid Model|
+------------------------------------------------------------------------------+
|                        -|                         -|             STANDARD_AIR|
+------------------------------------------------------------------------------+
-- Fluid properties:
+------------------------------------------------------------------------------+
|                  Name|    Dim. value|    Ref. value|      Unit|Non-dim. value|
+------------------------------------------------------------------------------+
|          Gas Constant|       287.058|             1|  N.m/kg.K|       287.058|
|      Spec. Heat Ratio|             -|             -|         -|           1.4|
+------------------------------------------------------------------------------+
-- Initial and free-stream conditions:
+------------------------------------------------------------------------------+
|                  Name|    Dim. value|    Ref. value|      Unit|Non-dim. value|
+------------------------------------------------------------------------------+
|       Static Pressure|        101325|             1|        Pa|        101325|
|               Density|       1.22498|             1|    kg/m^3|       1.22498|
|           Temperature|        288.15|             1|         K|        288.15|
|          Total Energy|        238496|             1|   m^2/s^2|        238496|
|            Velocity-X|        251.82|             1|       m/s|        251.82|
|            Velocity-Y|             0|             1|       m/s|             0|
|            Velocity-Z|             0|             1|       m/s|             0|
|    Velocity Magnitude|        251.82|             1|       m/s|        251.82|
+------------------------------------------------------------------------------+
|           Mach Number|             -|             -|         -|          0.74|
+------------------------------------------------------------------------------+
Initialize Jacobian structure (Euler). MG level: 0.

------------------- Numerics Preprocessing ( Zone 0 ) -------------------
No centered scheme.

----------------- Integration Preprocessing ( Zone 0 ) ------------------

------------------- Iteration Preprocessing ( Zone 0 ) ------------------
Euler/Navier-Stokes/RANS fluid iteration.

------------------------------ Begin Solver -----------------------------

Simulation Run using the Single-zone Driver

[jcownbey]

your boundary conditions are incorrect. You have inlet and outlet option set twice. Do you have a picture of your mesh for your boundary conditions?

[talbring]

You are missing the definition of the numerical schemes. The option CONV_NUM_METHOD_FLOW must be present. We could think about throwing an error in that case ...

[themrdjj]

Originally Posted by jcownbey

your boundary conditions are incorrect. You have inlet and outlet option set twice. Do you have a picture of your mesh for your boundary conditions?

I only have one MARKER_INLET and one MARKER_OUTLET, right? I can't get a picture of the mesh right now but the domain is just a simple cuboid.

[themrdjj]

Originally Posted by talbring

You are missing the definition of the numerical schemes. The option CONV_NUM_METHOD_FLOW must be present. We could think about throwing an error in that case ...

Ahhh that explains, it's working now! I definitely think there should be an error message. Especially because nowhere on the website could I find an indication of which inputs have a default value and which don't. I will open a feature request on Github.