Skip to content

Commit

Permalink
hit formatting of complete geometry input files
Browse files Browse the repository at this point in the history 
Refs #128
  • Loading branch information
Armando Pabon committed Jul 3, 2024
1 parent 137e7c9 commit a197872
Show file tree
Hide file tree
Showing 2 changed files with 164 additions and 172 deletions.
27 changes: 10 additions & 17 deletions examples/sps/multiapp/complete_geometry/electrothermal/engineering_scale_electrothermal.i
View file
Original file line number Diff line number Diff line change
Expand Up @@ -36,7 +36,7 @@
# Reference for graphite, stainless steel: Cincotti et al, DOI 10.1002/aic.11102
# Assorted references for yttria, listed as comments in input file

initial_temperature=293 #roughly 600C where the pyrometer kicks in
initial_temperature = 293 #roughly 600C where the pyrometer kicks in
#initial_porosity=0.36 #Maximum random jammed packing, Donev et al (2004) Science Magazine

[Mesh]
Expand Down Expand Up @@ -88,7 +88,8 @@ initial_temperature=293 #roughly 600C where the pyrometer kicks in
initial_condition = 2.0e-10 #in units eV/((nV)^2-s-nm)
block = 'powder_compact'
[]
[microapp_potential] #converted to microapp electronVolts units
[microapp_potential]
#converted to microapp electronVolts units
block = 'powder_compact'
[]
[E_x]
Expand All @@ -107,7 +108,6 @@ initial_temperature=293 #roughly 600C where the pyrometer kicks in
block = 'powder_compact'
[]


# [T_infinity]
# initial_condition = ${initial_temperature}
# []
Expand Down Expand Up @@ -384,7 +384,7 @@ initial_temperature=293 #roughly 600C where the pyrometer kicks in
variable = heat_transfer_radiation
boundary = 'outer_radiative_spacers outer_die_wall radiative_upper_plunger radiative_lower_plunger'
coupled_variables = 'temperature'
constant_names = 'boltzmann epsilon temperature_farfield' #published emissivity for graphite is 0.85
constant_names = 'boltzmann epsilon temperature_farfield' #published emissivity for graphite is 0.85
constant_expressions = '5.67e-8 0.85 293.0' #roughly room temperature, which is probably too cold
expression = '-boltzmann*epsilon*(temperature^4-temperature_farfield^4)'
[]
Expand All @@ -393,7 +393,7 @@ initial_temperature=293 #roughly 600C where the pyrometer kicks in
variable = heat_transfer_radiation
boundary = 'outer_radiative_stainless_steel'
coupled_variables = 'temperature_stainless_steel'
constant_names = 'boltzmann epsilon temperature_farfield' #published emissivity for graphite is 0.85
constant_names = 'boltzmann epsilon temperature_farfield' #published emissivity for graphite is 0.85
constant_expressions = '5.67e-8 0.4 293.0' #roughly room temperature, which is probably too cold
expression = '-boltzmann*epsilon*(temperature_stainless_steel^4-temperature_farfield^4)'
[]
Expand Down Expand Up @@ -620,10 +620,8 @@ initial_temperature=293 #roughly 600C where the pyrometer kicks in
gap_conductivity_function_variable = temperature
normal_smoothing_distance = 0.1
[]
[]

##Thermal Contact between gapped graphite die components
[ThermalContact]
[upper_plunger_spacer_gap_thermal]
type = GapHeatTransfer
primary = spacer_facing_upper_plunger
Expand Down Expand Up @@ -678,7 +676,7 @@ initial_temperature=293 #roughly 600C where the pyrometer kicks in
[]
[upper_plunger_diewall_gap_thermal]
type = GapHeatTransfer
primary = inner_die_wall ### paired temperature doesn't show on inner die wall, but temperature profile looks reasonable
primary = inner_die_wall ### paired temperature doesn't show on inner die wall, but temperature profile looks reasonable
secondary = die_wall_facing_upper_plunger
variable = temperature
quadrature = true
Expand All @@ -691,7 +689,7 @@ initial_temperature=293 #roughly 600C where the pyrometer kicks in
[]
[lower_plunger_diewall_gap_thermal]
type = GapHeatTransfer
primary = inner_die_wall ### paired temperature doesn't show on inner die wall, but temperature profile looks reasonable
primary = inner_die_wall ### paired temperature doesn't show on inner die wall, but temperature profile looks reasonable
secondary = die_wall_facing_lower_plunger
variable = temperature
quadrature = true
Expand All @@ -702,10 +700,8 @@ initial_temperature=293 #roughly 600C where the pyrometer kicks in
gap_conductivity_function_variable = temperature
normal_smoothing_distance = 0.1
[]
[]

## Thermal Contact between touching components of powder and die
[ThermalContact]
[upper_plunger_powder_thermal]
type = GapHeatTransfer
primary = bottom_upper_plunger
Expand Down Expand Up @@ -760,7 +756,6 @@ initial_temperature=293 #roughly 600C where the pyrometer kicks in
[]
[]


[Materials]
## graphite blocks
[graphite_thermal]
Expand Down Expand Up @@ -836,16 +831,16 @@ initial_temperature=293 #roughly 600C where the pyrometer kicks in
[]
[electrical_conductivity]
type = ADParsedMaterial
# coupled_variables = 'yttria_sigma_aeh'
# expression = 'yttria_sigma_aeh*1.602e8' #converts to units of J/(V^2-m-s)
# coupled_variables = 'yttria_sigma_aeh'
# expression = 'yttria_sigma_aeh*1.602e8' #converts to units of J/(V^2-m-s)
property_name = 'electrical_conductivity'
output_properties = electrical_conductivity
outputs = 'exodus csv'
block = powder_compact
# type = ADDerivativeParsedMaterial
# property_name = electrical_conductivity
coupled_variables = 'temperature'
constant_names = 'Q_elec kB prefactor_solid initial_porosity'
constant_names = 'Q_elec kB prefactor_solid initial_porosity'
constant_expressions = '1.61 8.617343e-5 1.25e-4 0.38'
expression = '(1-initial_porosity) * prefactor_solid * exp(-Q_elec/kB/temperature) * 1.602e8' # in eV/(nV^2 s nm) per chat with Larry, last term converts to units of J/(V^2-m-s)
[]
Expand Down Expand Up @@ -895,7 +890,6 @@ initial_temperature=293 #roughly 600C where the pyrometer kicks in
[]
[]


[Postprocessors]
[temperature_pp]
type = AverageNodalVariableValue
Expand Down Expand Up @@ -939,7 +933,6 @@ initial_temperature=293 #roughly 600C where the pyrometer kicks in
[]
[]


[Outputs]
csv = true
exodus = true
Expand Down
Loading

0 comments on commit a197872

Please sign in to comment.