Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

hit formatting of "complete geometry" input files #135

Open
wants to merge 1 commit into
base: devel
Choose a base branch
from
Open

hit formatting of "complete geometry" input files #135

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
27 changes: 10 additions & 17 deletions examples/sps/multiapp/complete_geometry/electrothermal/engineering_scale_electrothermal.i
View file
Open in desktop
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]
sapitts marked this conversation as resolved.
Show resolved Hide resolved
[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