Skip to content

Commit 98c7a4c

Browse files
bigfootedbigfooted
committed
add testcase
1 parent a1a0d7b commit 98c7a4c

File tree

5 files changed

+447
-1
lines changed

5 files changed

+447
-1
lines changed

‎Common/src/CConfig.cpp‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -513,7 +513,7 @@ void CConfig::addWallSpeciesType(const string& name, unsigned short & nMarker_Wa
513513
assert(option_map.find(name) == option_map.end());
514514
all_options.insert(pair<string, bool>(name, true));
515515

516-
//COptionBase* val = new COptionStringValuesList<WALL_SPECIES_TYPE*>(name, nMarker_Wall_Species, Marker_Wall_Species,
516+
//COptionBase* val = new COptionStringValuesList<<static_cast>(unsigned short)WALL_SPECIES_TYPE*>(name, nMarker_Wall_Species, Marker_Wall_Species,
517517
COptionBase* val = new COptionStringValuesList<unsigned short*>(name, nMarker_Wall_Species, Marker_Wall_Species,
518518
wall_species_type, nSpecies_per_Wall);
519519
option_map.insert(pair<string, COptionBase *>(name, val));

‎SU2_CFD/include/solvers/CSpeciesSolver.hpp‎

Lines changed: 12 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -118,6 +118,18 @@ class CSpeciesSolver : public CScalarSolver<CSpeciesVariable> {
118118
void BC_Isothermal_Wall(CGeometry* geometry, CSolver** solver_container, CNumerics* conv_numerics,
119119
CNumerics* visc_numerics, CConfig* config, unsigned short val_marker) override;
120120

121+
/*!
122+
* \brief Impose the Navier-Stokes wall boundary condition.
123+
* \param[in] geometry - Geometrical definition of the problem.
124+
* \param[in] solver_container - Container vector with all the solutions.
125+
* \param[in] conv_numerics - Description of the numerical method.
126+
* \param[in] visc_numerics - Description of the numerical method.
127+
* \param[in] config - Definition of the particular problem.
128+
* \param[in] val_marker - Surface marker where the boundary condition is applied.
129+
*/
130+
void BC_HeatFlux_Wall(CGeometry* geometry, CSolver** solver_container, CNumerics* conv_numerics,
131+
CNumerics* visc_numerics, CConfig* config, unsigned short val_marker) override;
132+
121133
/*!
122134
* \brief Generic implementation of the isothermal wall also covering CHT cases,
123135
* for which the wall temperature is given by GetConjugateHeatVariable.

‎SU2_CFD/src/solvers/CSpeciesSolver.cpp‎

Lines changed: 6 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -451,6 +451,12 @@ void CSpeciesSolver::BC_Isothermal_Wall(CGeometry* geometry, CSolver** solver_co
451451
BC_Isothermal_Wall_Generic(geometry, solver_container, conv_numerics, visc_numerics, config, val_marker);
452452
}
453453

454+
void CSpeciesSolver::BC_HeatFlux_Wall(CGeometry* geometry, CSolver** solver_container,
455+
CNumerics* conv_numerics, CNumerics* visc_numerics, CConfig* config,
456+
unsigned short val_marker) {
457+
BC_Isothermal_Wall_Generic(geometry, solver_container, conv_numerics, visc_numerics, config, val_marker);
458+
}
459+
454460
void CSpeciesSolver::BC_Isothermal_Wall_Generic(CGeometry* geometry, CSolver** solver_container,
455461
CNumerics* conv_numerics, CNumerics* visc_numerics,
456462
CConfig* config, unsigned short val_marker, bool cht_mode) {

‎TestCases/py_wrapper/psi.cfg‎

Lines changed: 172 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,172 @@
1+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
2+
% %
3+
% SU2 configuration file %
4+
% Case description: Turbulent premixed high pressure combustion chamber. %
5+
% Author: N. Beishuizen %
6+
% Institution: Bosch Thermotechniek B.V. %
7+
% Date: 2025/01/01 %
8+
% File Version 8.0 "Harrier" %
9+
% %
10+
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
11+
12+
% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
13+
%
14+
SOLVER= INC_RANS
15+
KIND_TURB_MODEL= SST
16+
SST_OPTIONS= V1994m
17+
18+
RESTART_SOL= YES
19+
20+
MGLEVEL= 0
21+
%
22+
23+
% temperature and density
24+
FREESTREAM_TEMPERATURE = 673
25+
FREESTREAM_DENSITY = 2.55
26+
% ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
27+
%
28+
%INC_DENSITY_MODEL= CONSTANT
29+
INC_DENSITY_MODEL= VARIABLE
30+
INC_DENSITY_INIT= 2.55
31+
%
32+
INC_VELOCITY_INIT= (40.00, 0.0, 0.0 )
33+
%
34+
INC_ENERGY_EQUATION= YES
35+
INC_TEMPERATURE_INIT= 673.0
36+
%
37+
INC_NONDIM= DIMENSIONAL
38+
%
39+
% -------------------- FLUID PROPERTIES ------------------------------------- %
40+
%
41+
FLUID_MODEL= INC_IDEAL_GAS
42+
%
43+
%CONDUCTIVITY_MODEL= CONSTANT_CONDUCTIVITY
44+
CONDUCTIVITY_MODEL= CONSTANT_PRANDTL
45+
THERMAL_CONDUCTIVITY_CONSTANT= 0.0357
46+
%
47+
PRANDTL_LAM= 0.72
48+
TURBULENT_CONDUCTIVITY_MODEL= NONE
49+
PRANDTL_TURB= 0.90
50+
%
51+
VISCOSITY_MODEL= SUTHERLAND
52+
MU_CONSTANT= 1.716E-5
53+
MU_REF = 1.716e-5
54+
MU_T_REF= 273.15
55+
SUTHERLAND_CONSTANT = 110.4
56+
57+
SPECIFIC_HEAT_CP = 1150
58+
%
59+
% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
60+
%
61+
MARKER_HEATFLUX= ( wall_side, 0.0, wall_pipe, 0.0, wall_out, 0.0 )
62+
MARKER_ISOTHERMAL= ( wall_top, 400.0 )
63+
64+
% note, case is axisymmetric
65+
MARKER_SYM= ( symmetry )
66+
AXISYMMETRIC= YES
67+
%
68+
SPECIFIED_INLET_PROFILE= NO
69+
INLET_MATCHING_TOLERANCE=1e-4
70+
INLET_FILENAME= inlet.dat
71+
%INLET_INTERPOLATION_FUNCTION= LINEAR_1D
72+
INC_INLET_TYPE= VELOCITY_INLET
73+
INC_INLET_DAMPING= 0.01
74+
MARKER_INLET= ( inlet, 673, 40.0, 1.0, 0.0, 0.0)
75+
MARKER_INLET_TURBULENT = (inlet, 0.10, 15)
76+
77+
78+
MARKER_INLET_SPECIES= (inlet, 0.0, 400.0)
79+
%WALL_SPECIES_TYPE = wall_side, FLUX, VALUE, wall_pipe, FLUX, VALUE,wall_out, FLUX, VALUE,wall_top, FLUX, VALUE
80+
WALL_SPECIES_TYPE = wall_side, 0 , 1 , wall_pipe, 0 , 1 ,wall_out 0 , 1 , wall_top 0 , 1
81+
MARKER_WALL_SPECIES= wall_side, 0.0, 300.0, wall_pipe, 0.0, 400.0, wall_out, 0.0, 500.0, wall_top, 0.0, 300.0
82+
83+
%
84+
INC_OUTLET_TYPE= PRESSURE_OUTLET
85+
INC_OUTLET_DAMPING= 0.01
86+
MARKER_OUTLET= ( outlet, 0.0 )
87+
%
88+
% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
89+
%
90+
NUM_METHOD_GRAD= WEIGHTED_LEAST_SQUARES
91+
%
92+
CFL_NUMBER= 10.0
93+
CFL_REDUCTION_SPECIES= 1.0
94+
CFL_REDUCTION_TURB= 1.0
95+
%
96+
ITER= 100
97+
%
98+
% ------------------------ LINEAR SOLVER DEFINITION ---------------------------%
99+
%
100+
LINEAR_SOLVER= FGMRES
101+
LINEAR_SOLVER_PREC= ILU
102+
LINEAR_SOLVER_ERROR= 1E-12
103+
LINEAR_SOLVER_ITER= 25
104+
105+
% -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
106+
%
107+
CONV_NUM_METHOD_FLOW= FDS
108+
MUSCL_FLOW= NO
109+
SLOPE_LIMITER_FLOW = NONE
110+
TIME_DISCRE_FLOW= EULER_IMPLICIT
111+
%
112+
% -------------------- SCALAR TRANSPORT ---------------------------------------%
113+
%
114+
KIND_SCALAR_MODEL= SPECIES_TRANSPORT
115+
DIFFUSIVITY_MODEL= CONSTANT_DIFFUSIVITY
116+
SCHMIDT_NUMBER_LAMINAR= 1.0
117+
DIFFUSIVITY_CONSTANT= 7.56e-5
118+
119+
% according to the paper
120+
SCHMIDT_NUMBER_TURBULENT= 0.7
121+
122+
%
123+
CONV_NUM_METHOD_SPECIES= BOUNDED_SCALAR
124+
MUSCL_SPECIES= NO
125+
SLOPE_LIMITER_SPECIES = NONE
126+
%
127+
TIME_DISCRE_SPECIES= EULER_IMPLICIT
128+
%
129+
SPECIES_INIT= 0.0 0.0
130+
%SPECIES_CLIPPING= YES
131+
%SPECIES_CLIPPING_MIN= 0.0
132+
%SPECIES_CLIPPING_MAX= 1.0
133+
%
134+
% -------------------- TURBULENT TRANSPORT ---------------------------------------%
135+
%
136+
CONV_NUM_METHOD_TURB= BOUNDED_SCALAR
137+
MUSCL_TURB= NO
138+
139+
%
140+
% --------------------------- CONVERGENCE PARAMETERS --------------------------%
141+
%
142+
CONV_FIELD= RMS_PRESSURE, RMS_VELOCITY-X, RMS_VELOCITY-Y, RMS_TKE, RMS_SPECIES
143+
CONV_RESIDUAL_MINVAL= -18
144+
CONV_STARTITER= 10
145+
%
146+
% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
147+
%
148+
%MESH_FILENAME= psi_coarse.su2
149+
MESH_FILENAME= psi.su2
150+
%MESH_FILENAME= psi_fine.su2
151+
%
152+
SCREEN_OUTPUT= INNER_ITER WALL_TIME \
153+
RMS_TEMPERATURE RMS_PRESSURE RMS_VELOCITY-X RMS_VELOCITY-Y RMS_TKE RMS_DISSIPATION RMS_SPECIES_0 \
154+
LINSOL_ITER LINSOL_RESIDUAL \
155+
LINSOL_ITER_TURB LINSOL_RESIDUAL_TURB \
156+
LINSOL_ITER_SPECIES LINSOL_RESIDUAL_SPECIES
157+
SCREEN_WRT_FREQ_INNER= 1
158+
%
159+
HISTORY_OUTPUT= ITER RMS_RES LINSOL SPECIES_COEFF SPECIES_COEFF_SURF
160+
CONV_FILENAME= history
161+
MARKER_ANALYZE= gas_inlet, air_axial_inlet, outlet
162+
MARKER_ANALYZE_AVERAGE= AREA
163+
%
164+
OUTPUT_FILES= RESTART_ASCII, PARAVIEW_MULTIBLOCK
165+
VOLUME_OUTPUT= RESIDUAL, PRIMITIVE
166+
OUTPUT_WRT_FREQ= 1000
167+
%
168+
READ_BINARY_RESTART= NO
169+
RESTART_FILENAME= restart
170+
SOLUTION_FILENAME= solution
171+
%
172+
WRT_PERFORMANCE= YES

0 commit comments

Comments
 (0)