-
Notifications
You must be signed in to change notification settings - Fork 0
/
RadMom1DInput_Third_Order.cc
129 lines (105 loc) · 3.88 KB
/
RadMom1DInput_Third_Order.cc
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
/*!\file RadMom1DInput.cc
\brief Subroutines for 1D Radiative Moment Input Classes. */
/* Include RADMOM header files */
#include "RadMom1DInput.h" /* Include 1D RadMom input parameter header file. */
// static variables
template<>
const int RadMom1D_Input_Parameters<RadMom1D_cState_Third_Order,
RadMom1D_pState_Third_Order>::i_Residual_Variable;
template<>
const int RadMom1D_Input_Parameters<RadMom1D_cState_Third_Order,
RadMom1D_pState_Third_Order>::Number_of_Residual_Norms;
/******************************************************//**
* Routine: Set_Default_Input_Parameters
*
* Assigns default values to the input parameters.
*
********************************************************/
template<>
void RadMom1D_Input_Parameters<RadMom1D_cState_Third_Order,RadMom1D_pState_Third_Order>::Set_Default_Input_Parameters(void) {
strcpy(Input_File_Name, "RadMom1D.in");
strcpy(Time_Integration_Type, "Explicit_Euler");
i_Time_Integration = TIME_STEPPING_EXPLICIT_EULER;
Time_Accurate = 0;
Local_Time_Stepping = SCALAR_LOCAL_TIME_STEPPING;
CFL_Number = 0.3;
Time_Max = ZERO;
// Reconstruction type:
strcpy(Reconstruction_Type, "Least_Squares");
i_Reconstruction = RECONSTRUCTION_LEAST_SQUARES;
// Limiter type:
strcpy(Limiter_Type, "Barth_Jespersen");
i_Limiter = LIMITER_BARTH_JESPERSEN;
// Flux function:
strcpy(Flux_Function_Type, "HLLE");
i_Flux_Function = FLUX_FUNCTION_HLLE;
// Moment function:
strcpy(Moment_Closure_Type, "P3");
i_Moment_Closure = MOMENT_CLOSURE_P3;
// Initial conditions:
strcpy(ICs_Type, "Uniform");
strcpy(ICs_Intensity_or_Temperature, "Temperature");
i_ICs = IC_UNIFORM;
i_ICs_Flags = RADMOM_ICS_TEMPERATURE;
Intensity = ZERO;
Temperature = THOUSAND; //[K]
// Pressure = PRESSURE_STDATM;//[Pa]
xco = 0.01;
xh2o = 0.2;
xco2 = 0.1;
xo2 = 0.0;
fsoot = 0.0;
Reference_Temp = Temperature;
Case = 0;
// Gas parameters
i_AbsorptionModel = MEDIUM1D_ABSORB_GRAY;
strcpy(AbsorptionModel, "Gray");
i_ScatteringFunc = RADIATION_SCATTER_ISO;
strcpy(ScatteringFunc, "Isotropic");
AbsorptionCoef = ONE;
ScatteringCoef = ZERO;
OptThickness = ONE;
Albedo = ZERO;
// boundary conditions
EastWallTemp = ZERO;
WestWallTemp = ZERO;
EastWallEmiss = ZERO;
WestWallEmiss = ZERO;
// Grid parameters:
Box_Width = ONE;
X_Min = ZERO;
X_Max = ONE;
Number_of_Cells_Idir = 100;
Number_of_Ghost_Cells = 2;
Number_of_Blocks_Idir = 1;
Number_of_Blocks = 1;
Number_of_Processors = 1;
// Boundary conditions:
strcpy(Boundary_Conditions_Specified,"OFF");
BCs_Specified = OFF;
BC_East = BC_GRAY_WALL;
BC_West = BC_GRAY_WALL;
strcpy(Boundary_Conditions_Enforcement,"WEAK");
BCs_Enforcement = WEAK_BCS;
strcpy(BC_East_Type,"None");
strcpy(BC_West_Type,"None");
BC_East = BC_NONE;
BC_West = BC_NONE;
// Default output file names and parameters:
strcpy(Output_File_Name, "outputfile.dat");
strcpy(Grid_File_Name, "gridfile.grid");
strcpy(Restart_File_Name, "restartfile.soln");
strcpy(Output_Format_Type, "Tecplot");
i_Output_Format = IO_TECPLOT;
Restart_Solution_Save_Frequency = 10000;
Output_Progress_Frequency = 50;
// Input_file parameters:
strcpy(Next_Control_Parameter, " ");
Line_Number = 0;
// Freezing_Limiter
Freeze_Limiter = 0;
Freeze_Limiter_Residual_Level = 1e-4;
Min_Residual_Level = TOLER;
// Setup conserved and medium state
SetupInputState();
}