su2code · bigfooted · Dec 22, 2025 · Dec 27, 2025 · Jan 7, 2026 · Jan 11, 2026
diff --git a/multiphysics/TFC_python/adiabatic/inlet.dat b/multiphysics/TFC_python/adiabatic/inlet.dat
@@ -0,0 +1,35 @@
+NMARK= 1
+MARKER_TAG= inlet
+NROW=30
+NCOL=9
+# COORD-X               COORD-Y                 TEMPERATURE             VELOCITY                NORMAL-X                NORMAL-Y                TKE                     DISSIPATION             SPECIES_0               
+-0.1 0.0 673 52.113575 1 0 0.022036247 2249.2866 0
+-0.1 0.00072056672 673 52.113575 1 0 0.022036247 2249.2866 0
+-0.1 0.0014123108 673 52.130169 1 0 0.021853773 2227.1838 0
+-0.1 0.0020763851 673 52.155247 1 0 0.021966202 2223.8691 0
+-0.1 0.0027138963 673 52.185356 1 0 0.025456807 2238.5793 0
+-0.1 0.003325907 673 52.181305 1 0 0.11307145 2565.6211 0
+-0.1 0.0039134375 673 51.904926 1 0 0.91227305 3924.4668 0
+-0.1 0.0044774669 673 51.180588 1 0 3.0400403 5234.1846 0
+-0.1 0.0050189346 673 50.250782 1 0 5.6927605 5843.0571 0
+-0.1 0.0055387439 673 49.246353 1 0 8.354352 6278.2769 0
+-0.1 0.006037761 673 48.157909 1 0 11.022179 6693.3311 0
+-0.1 0.0065168175 673 46.984482 1 0 13.741647 7145.3198 0
+-0.1 0.0069767111 673 45.774849 1 0 16.35457 7609.7793 0
+-0.1 0.0074182097 673 44.565525 1 0 18.708464 8072.2192 0
+-0.1 0.0078420481 673 43.370502 1 0 20.779472 8543.3164 0
+-0.1 0.0082489327 673 42.1898 1 0 22.594801 9048.4883 0
+-0.1 0.0086395424 673 41.016678 1 0 24.199833 9620.7598 0
+-0.1 0.0090145273 673 39.840576 1 0 25.640738 10297.504 0
+-0.1 0.0093745124 673 38.646294 1 0 26.959324 11122.591 0
+-0.1 0.0097200992 673 37.411186 1 0 28.190733 12151.757 0
+-0.1 0.010051861 673 36.11153 1 0 29.362619 13455.807 0
+-0.1 0.010370354 673 34.725197 1 0 30.4884 15132.082 0
+-0.1 0.010676106 673 33.222515 1 0 31.5849 17335.744 0
+-0.1 0.010969629 673 31.561712 1 0 32.676102 20327.223 0
+-0.1 0.01125141 673 29.680511 1 0 33.802517 24576.211 0
+-0.1 0.011521921 673 27.479441 1 0 35.046871 31017.684 0
+-0.1 0.011781611 673 24.78553 1 0 36.619263 41792.531 0
+-0.1 0.012030913 673 21.265036 1 0 39.229969 62961.266 0
+-0.1 0.012270243 673 15.721132 1 0 44.032578 119438.84 0
+-0.1 0.0125 673 0 1 0 1e-10 211634.48 0
diff --git a/multiphysics/TFC_python/adiabatic/psi.cfg b/multiphysics/TFC_python/adiabatic/psi.cfg
@@ -0,0 +1,166 @@
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+%                                                                              %
+% SU2 configuration file                                                       %
+% Case description: Turbulent premixed high pressure combustion chamber.       %
+% Author: N. Beishuizen                                                        %
+% Institution: Bosch Thermotechniek B.V.                                       %
+% Date: 2025/01/01                                                             %
+% File Version 8.0 "Harrier"                                                   %
+%                                                                              %
+%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+% ------------- DIRECT, ADJOINT, AND LINEARIZED PROBLEM DEFINITION ------------%
+%
+SOLVER= INC_RANS
+KIND_TURB_MODEL= SST
+SST_OPTIONS= V1994m
+
+RESTART_SOL= YES
+
+% temperature and density 
+FREESTREAM_TEMPERATURE = 673
+FREESTREAM_DENSITY = 2.55
+% ---------------- INCOMPRESSIBLE FLOW CONDITION DEFINITION -------------------%
+%
+INC_DENSITY_MODEL= VARIABLE
+INC_DENSITY_INIT= 2.55
+%
+INC_VELOCITY_INIT= (40.00, 0.0, 0.0 )
+%
+INC_ENERGY_EQUATION= NO
+INC_TEMPERATURE_INIT= 673.0
+%
+INC_NONDIM= DIMENSIONAL
+%
+% -------------------- FLUID PROPERTIES ------------------------------------- %
+%
+FLUID_MODEL= INC_IDEAL_GAS
+%
+CONDUCTIVITY_MODEL= CONSTANT_PRANDTL
+%THERMAL_CONDUCTIVITY_CONSTANT= 0.0357
+%
+PRANDTL_LAM= 0.72
+TURBULENT_CONDUCTIVITY_MODEL= CONSTANT_PRANDTL_TURB
+%TURBULENT_CONDUCTIVITY_MODEL= NONE 
+PRANDTL_TURB= 0.90
+%
+%VISCOSITY_MODEL= SUTHERLAND
+VISCOSITY_MODEL= POLYNOMIAL_VISCOSITY 
+% polynomial viscosity for methane-air (phi=0.50) 
+MU_POLYCOEFFS= (1.097e-5, 3.437e-8, 3.332e-12, 0, 0)
+
+MU_CONSTANT=  1.716E-5
+MU_REF = 1.716e-5
+MU_T_REF= 273.15
+SUTHERLAND_CONSTANT = 110.4
+
+SPECIFIC_HEAT_CP = 1350
+%
+% -------------------- BOUNDARY CONDITION DEFINITION --------------------------%
+%
+MARKER_HEATFLUX= ( wall_top, 0.0, wall_side, 0.0, wall_pipe, 0.0, wall_out,0.0, wall_side_out,0.0 )
+
+% note, case is axisymmetric
+MARKER_SYM= ( symmetry )
+AXISYMMETRIC= YES
+%
+SPECIFIED_INLET_PROFILE= YES
+INLET_MATCHING_TOLERANCE=1e-4
+INLET_FILENAME= inlet.dat
+INC_INLET_TYPE=  VELOCITY_INLET
+INC_INLET_DAMPING= 0.1
+MARKER_INLET= ( inlet, 673, 40.0, 1.0,  0.0, 0.0)
+MARKER_INLET_TURBULENT = (inlet, 0.10, 15)
+MARKER_INLET_SPECIES= (inlet, 0.0)
+%
+INC_OUTLET_TYPE= PRESSURE_OUTLET
+INC_OUTLET_DAMPING= 0.1
+MARKER_OUTLET= ( outlet, 0.0 )
+%
+% ------------- COMMON PARAMETERS DEFINING THE NUMERICAL METHOD ---------------%
+%
+NUM_METHOD_GRAD= WEIGHTED_LEAST_SQUARES
+%
+CFL_NUMBER= 10.0
+CFL_REDUCTION_SPECIES= 1.0
+CFL_REDUCTION_TURB= 1.0
+CFL_ADAPT= YES
+CFL_ADAPT_PARAM= ( 0.95, 1.01, 1.0, 250, 1.0e-4, 0)
+
+%
+ITER= 1
+%
+% ------------------------ LINEAR SOLVER DEFINITION ---------------------------%
+%
+LINEAR_SOLVER= FGMRES
+LINEAR_SOLVER_PREC= ILU
+LINEAR_SOLVER_ERROR= 1E-6
+LINEAR_SOLVER_ITER= 10
+
+% -------------------- FLOW NUMERICAL METHOD DEFINITION -----------------------%
+%
+CONV_NUM_METHOD_FLOW= FDS
+MUSCL_FLOW= NO
+SLOPE_LIMITER_FLOW = NONE
+TIME_DISCRE_FLOW= EULER_IMPLICIT
+%
+% -------------------- SCALAR TRANSPORT ---------------------------------------%
+%
+KIND_SCALAR_MODEL= SPECIES_TRANSPORT
+%DIFFUSIVITY_MODEL= CONSTANT_DIFFUSIVITY
+DIFFUSIVITY_MODEL= CONSTANT_SCHMIDT
+%DIFFUSIVITY_MODEL= CONSTANT_LEWIS
+SCHMIDT_NUMBER_LAMINAR= 1.0
+CONSTANT_LEWIS_NUMBER= 1.0 1.0
+DIFFUSIVITY_CONSTANT= 5.0e-5
+
+% according to the paper
+SCHMIDT_NUMBER_TURBULENT= 0.7
+
+%
+CONV_NUM_METHOD_SPECIES= BOUNDED_SCALAR
+MUSCL_SPECIES= NO
+SLOPE_LIMITER_SPECIES = NONE
+%
+TIME_DISCRE_SPECIES= EULER_IMPLICIT
+%
+SPECIES_INIT= 0.0
+SPECIES_CLIPPING= YES
+SPECIES_CLIPPING_MIN= 0.0
+SPECIES_CLIPPING_MAX= 1.0
+%
+% -------------------- TURBULENT TRANSPORT ---------------------------------------%
+%
+CONV_NUM_METHOD_TURB= BOUNDED_SCALAR
+MUSCL_TURB= NO
+%
+% --------------------------- CONVERGENCE PARAMETERS --------------------------%
+%
+CONV_FIELD= RMS_PRESSURE, RMS_VELOCITY-X, RMS_VELOCITY-Y, RMS_TKE, RMS_SPECIES
+CONV_RESIDUAL_MINVAL= -12
+CONV_STARTITER= 10
+%
+% ------------------------- INPUT/OUTPUT INFORMATION --------------------------%
+%
+MESH_FILENAME= psi.su2
+%
+SCREEN_OUTPUT= INNER_ITER WALL_TIME \
+               RMS_PRESSURE RMS_VELOCITY-X RMS_VELOCITY-Y RMS_TKE RMS_DISSIPATION RMS_SPECIES_0
+
+SCREEN_WRT_FREQ_INNER= 1
+%
+HISTORY_OUTPUT= ITER RMS_RES LINSOL SPECIES_COEFF SPECIES_COEFF_SURF
+CONV_FILENAME= history
+MARKER_ANALYZE= gas_inlet, air_axial_inlet, outlet
+MARKER_ANALYZE_AVERAGE= AREA
+%
+OUTPUT_FILES= RESTART, PARAVIEW_MULTIBLOCK
+VOLUME_OUTPUT= RESIDUAL, PRIMITIVE, SPECIES_UDS_0
+OUTPUT_WRT_FREQ= 100
+%
+READ_BINARY_RESTART= YES
+RESTART_FILENAME= restart
+SOLUTION_FILENAME= solution
+%
+WRT_PERFORMANCE= YES
+PYTHON_CUSTOM_SOURCE= YES