Skip to content

Commit

Permalink
Add tests for error catching and calculated conductances.
Browse files Browse the repository at this point in the history
Refs #13
  • Loading branch information
cticenhour committed Jun 30, 2022
1 parent 0a3b70d commit 1ce3839
Show file tree
Hide file tree
Showing 3 changed files with 47 additions and 2 deletions.
Binary file not shown.
18 changes: 18 additions & 0 deletions test/tests/interfacekernels/thermal_conductance/tests
Original file line number Diff line number Diff line change
Expand Up @@ -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'
[../]
[]
Original file line number Diff line number Diff line change
@@ -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
[]
Expand Down Expand Up @@ -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
Expand All @@ -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
Expand Down Expand Up @@ -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]
Expand Down

0 comments on commit 1ce3839

Please sign in to comment.