13 const Real
SIGMA = 5.670374419E-8;
21 InputParameters params = validParams<InterfaceMaterial>();
22 params.addClassDescription(
"This material constructs the necessary coefficients and properties "
23 "for SideSetHeatTransferKernel.");
24 params.addParam<FunctionName>(
"conductivity", 0.0,
"Heat conductivity in W/m/K.");
25 params.addParam<FunctionName>(
"conductivity_temperature_function",
26 "Heat conductivity in W/m/K as a function of temperature.");
29 "Coupled Temperature of gap, used for computing temperature dependent conductivity.");
30 params.addParam<FunctionName>(
"gap_length", 1.0,
"Total width of gap in m.");
31 params.addParam<FunctionName>(
"Tbulk", 300,
"Bulk temperature of gap in K.");
32 params.addParam<FunctionName>(
33 "h_master", 0.0,
"Convective heat transfer coefficient (master face) in W/m^2/K.");
34 params.addParam<FunctionName>(
35 "h_neighbor", 0.0,
"Convective heat transfer coefficient (neighbor face) in W/m^2/K.");
36 params.addParam<FunctionName>(
"emissivity_master", 0.0,
"Master face emissivity.");
37 params.addParam<FunctionName>(
"emissivity_neighbor", 0.0,
"Neighbor face emissivity.");
38 params.addParam<FunctionName>(
"reflectivity_master",
39 "Master face reflectivity, uses (1-emissivity) if not provided.");
40 params.addParam<FunctionName>(
"reflectivity_neighbor",
41 "Neighbor face reflectivity, uses (1-emissivity) if not provided.");
46 : InterfaceMaterial(parameters),
47 _kgap(isParamValid(
"conductivity_temperature_function")
48 ? getFunction(
"conductivity_temperature_function")
49 : getFunction(
"conductivity")),
50 _Tk(isCoupled(
"gap_temperature") ? &coupledValue(
"gap_temperature") : nullptr),
51 _dgap(getFunction(
"gap_length")),
52 _Tb(getFunction(
"Tbulk")),
53 _hp(getFunction(
"h_master")),
54 _hm(getFunction(
"h_neighbor")),
55 _eps_p(getFunction(
"emissivity_master")),
56 _eps_m(getFunction(
"emissivity_neighbor")),
57 _rho_p(isParamValid(
"reflectivity_master") ? &getFunction(
"reflectivity_master") : nullptr),
58 _rho_m(isParamValid(
"reflectivity_neighbor") ? &getFunction(
"reflectivity_neighbor") : nullptr),
59 _cond(declareProperty<Real>(
"gap_conductance")),
60 _Tbulk(declareProperty<Real>(
"gap_Tbulk")),
61 _h_master(declareProperty<Real>(
"gap_h_master")),
62 _h_neighbor(declareProperty<Real>(
"gap_h_neighbor")),
63 _emmissivity_eff_master(declareProperty<Real>(
"gap_emissivity_eff_master")),
64 _emmissivity_eff_neighbor(declareProperty<Real>(
"gap_emissivity_eff_neighbor")),
67 if ((parameters.isParamSetByUser(
"conductivity") ||
68 isParamValid(
"conductivity_temperature_function")) &&
69 !parameters.isParamSetByUser(
"gap_length"))
70 paramError(
"gap_length",
"gap_length not set, but conduction term requested.");
71 if (parameters.isParamSetByUser(
"gap_length") &&
72 !(parameters.isParamSetByUser(
"conductivity") ||
73 isParamValid(
"conductivity_temperature_function")))
74 paramError(
"conductivity",
"conductivity not set, but conduction term requested.");
75 if (isParamValid(
"conductivity_temperature_function") &&
76 parameters.isParamSetByUser(
"conductivity"))
77 paramError(
"conductivity",
78 "Cannot specify both conductivity and conductivity_temperature_function.");
79 if (isParamValid(
"conductivity_temperature_function") && !
_Tk)
80 paramError(
"gap_temperature",
81 "Variable specification for temp needed if specifying a temperature dependent "
84 if (parameters.isParamSetByUser(
"h_master") && !parameters.isParamSetByUser(
"h_neighbor"))
85 paramError(
"h_neighbor",
"h_neighbor not set, but convection term requested.");
86 if (parameters.isParamSetByUser(
"h_neighbor") && !parameters.isParamSetByUser(
"h_master"))
87 paramError(
"h_master",
"h_master not set, but convection term requested.");
89 if (parameters.isParamSetByUser(
"emissivity_master") &&
90 !parameters.isParamSetByUser(
"emissivity_neighbor"))
91 paramError(
"emissivity_neighbor",
"emissivity_neighbor not set, but radiation term requested");
92 if (parameters.isParamSetByUser(
"emissivity_neighbor") &&
93 !parameters.isParamSetByUser(
"emissivity_master"))
94 paramError(
"emissivity_master",
"emissivity_master not set, but radiation term requested");
101 Real tqp = (
_Tk ? (*_Tk)[_qp] : _t);
102 _cond[_qp] =
_kgap.value(tqp, _q_point[_qp]) /
_dgap.value(_t, _q_point[_qp]);
105 _Tbulk[_qp] =
_Tb.value(_t, _q_point[_qp]);
110 Real rhop = (
_rho_p ? (*_rho_p).value(_t, _q_point[_qp]) : 1.0 -
_eps_p.value(_t, _q_point[_qp]));
111 Real rhom = (
_rho_m ? (*_rho_m).value(_t, _q_point[_qp]) : 1.0 -
_eps_m.value(_t, _q_point[_qp]));