diff --git a/test/tests/interfacekernels/thermal_conductance/gold/thermal_interface_calculated_out.e b/test/tests/interfacekernels/thermal_conductance/gold/thermal_interface_calculated_out.e new file mode 100644 index 00000000..038c41a2 Binary files /dev/null and b/test/tests/interfacekernels/thermal_conductance/gold/thermal_interface_calculated_out.e differ diff --git a/test/tests/interfacekernels/thermal_conductance/tests b/test/tests/interfacekernels/thermal_conductance/tests index e02395aa..3c2616fd 100644 --- a/test/tests/interfacekernels/thermal_conductance/tests +++ b/test/tests/interfacekernels/thermal_conductance/tests @@ -40,6 +40,24 @@ ratio_tol = 3e-7 #ratio = 2.80685e-7 cli_args = 'Mesh/file=thermal_interface_jacobian_mesh.e Executioner/end_time=0.1' prereq = 'low_low low_high high_low high_high' + requirement = 'The system shall support the calculation of heat transfer across an interface due to thermal conductivity and electrostatic joule heating with the correct jacobian.' issues = '#7 #10' [../] + [./conductance_error] + type = RunException + input = 'thermal_interface.i' + cli_args = 'InterfaceKernels/thermal_contact_conductance/mean_hardness=1.0' + expect_err = "In thermal_contact_conductance, both user-supplied thermal/electrical conductances" + prereq = 'low_low low_high high_low high_high' + requirement = 'The system shall support error-generation when the user supplies both conductances and mean hardness values (to calculate a conductance) to calculate thermal conductivity across an interface.' + issues = '#7' + [../] + [./conductance_calculated] + type = Exodiff + input = 'thermal_interface.i' + exodiff = 'thermal_interface_calculated_out.e' + cli_args = 'Executioner/end_time=0.2 Outputs/file_base=thermal_interface_calculated_out InterfaceKernels/active="thermal_contact_conductance_calculated electric_contact_conductance" Materials/active="heat_conductor_graphite rho_graphite sigma_graphite heat_conductor_stainless_steel rho_stainless_steel sigma_stainless_steel mean_hardness"' + requirement = 'The system shall support the calculation of thermal and electrical conductances when determining the heat transfer across an interface due to thermal conductivity and electrostatic joule heating.' + issues = '#7' + [../] [] diff --git a/test/tests/interfacekernels/thermal_conductance/thermal_interface.i b/test/tests/interfacekernels/thermal_conductance/thermal_interface.i index 62c197ad..0c2fc5cd 100644 --- a/test/tests/interfacekernels/thermal_conductance/thermal_interface.i +++ b/test/tests/interfacekernels/thermal_conductance/thermal_interface.i @@ -1,3 +1,4 @@ +# Input file block to generate this mesh can be found in thermal_interface_mesh.i [Mesh] file = thermal_interface_regular_mesh.e [] @@ -110,14 +111,18 @@ [] [InterfaceKernels] + active = 'thermal_contact_conductance electric_contact_conductance' + [./thermal_contact_conductance] type = ThermalContactCondition variable = temperature_stainless_steel neighbor_var = temperature_graphite master_potential = potential_stainless_steel neighbor_potential = potential_graphite - master_conductivity = thermal_conductivity - neighbor_conductivity = thermal_conductivity + master_thermal_conductivity = thermal_conductivity + neighbor_thermal_conductivity = thermal_conductivity + master_electrical_conductivity = electrical_conductivity + neighbor_electrical_conductivity = electrical_conductivity user_electrical_contact_conductance = 2.5e5 # as described in Cincotti et al (DOI: 10.1002/aic.11102) user_thermal_contact_conductance = 7 # also from Cincotti et al boundary = ssg_interface @@ -129,9 +134,26 @@ boundary = ssg_interface user_electrical_contact_conductance = 2.5e5 # as described in Cincotti et al (DOI: 10.1002/aic.11102) [../] + + [./thermal_contact_conductance_calculated] + type = ThermalContactCondition + variable = temperature_stainless_steel + neighbor_var = temperature_graphite + master_potential = potential_stainless_steel + neighbor_potential = potential_graphite + master_thermal_conductivity = thermal_conductivity + neighbor_thermal_conductivity = thermal_conductivity + master_electrical_conductivity = electrical_conductivity + neighbor_electrical_conductivity = electrical_conductivity + mean_hardness = graphite_stainless_mean_hardness + mechanical_pressure = 8.52842e10 # resulting in electrical contact conductance = ~1.4715e5, thermal contact conductance = ~3.44689e7 + boundary = ssg_interface + [../] [] [Materials] + active = 'heat_conductor_graphite rho_graphite sigma_graphite heat_conductor_stainless_steel rho_stainless_steel sigma_stainless_steel' + #graphite [./heat_conductor_graphite] type = HeatConductionMaterial @@ -175,6 +197,11 @@ temperature_coefficient = 0 # makes conductivity constant block = stainless_steel [../] + + # harmonic mean of graphite and stainless steel hardness + [./mean_hardness] + type = GraphiteStainlessMeanHardness + [../] [] [Preconditioning]