Skip to content

Commit

Permalink
hit formatting of diemodel input files Refs idaholab#128
Browse files Browse the repository at this point in the history
  • Loading branch information
Armando Pabon committed Jun 24, 2024
1 parent c3a6a23 commit d0a7328
Show file tree
Hide file tree
Showing 7 changed files with 107 additions and 94 deletions.
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
## Units in the input file: m-Pa-s-K

initial_temperature=300 #roughly 600C where the pyrometer kicks in
initial_temperature = 300 #roughly 600C where the pyrometer kicks in

[GlobalParams]
order = SECOND
Expand Down Expand Up @@ -46,7 +46,8 @@ initial_temperature=300 #roughly 600C where the pyrometer kicks in
[sigma_aeh]
initial_condition = 2.0e-10 #in units eV/((nV)^2-s-nm)
[]
[microapp_potential] #converted to microapp electronVolts units
[microapp_potential]
#converted to microapp electronVolts units
[]
[E_x]
order = FIRST
Expand Down Expand Up @@ -143,7 +144,7 @@ initial_temperature=300 #roughly 600C where the pyrometer kicks in
variable = heat_transfer_radiation
boundary = 'outer_die_wall'
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 300.0' #roughly room temperature, which is probably too cold
expression = '-boltzmann*epsilon*(temperature^4-temperature_farfield^4)'
[]
Expand Down Expand Up @@ -386,7 +387,6 @@ initial_temperature=300 #roughly 600C where the pyrometer kicks in
[]
[]


[Materials]
[graphite_density]
type = ADGenericConstantMaterial
Expand Down Expand Up @@ -457,16 +457,16 @@ initial_temperature=300 #roughly 600C where the pyrometer kicks in
[]
[electrical_conductivity]
type = ADParsedMaterial
# coupled_variables = 'sigma_aeh'
# expression = 'sigma_aeh*1.602e8' #converts to units of J/(V^2-m-s)
# coupled_variables = 'sigma_aeh'
# expression = '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 @@ -537,13 +537,12 @@ initial_temperature=300 #roughly 600C where the pyrometer kicks in
[]
[]


[Outputs]
csv = true
exodus = true
perf_graph = true
[ckpt]
type =Checkpoint
type = Checkpoint
time_step_interval = 1
num_files = 2
[]
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -59,7 +59,8 @@ initial_temperature = 300 #roughly 600C where the pyrometer kicks in
# order = FIRST
# block = 'upper_plunger lower_plunger die_wall'
# []
[microapp_potential] #converted to microapp electronVolts units
[microapp_potential]
#converted to microapp electronVolts units
[]
[E_x]
order = FIRST
Expand All @@ -84,23 +85,25 @@ initial_temperature = 300 #roughly 600C where the pyrometer kicks in
[]

[Physics]
[SolidMechanics/QuasiStatic]
[graphite]
strain = FINITE
add_variables = true
use_automatic_differentiation = true
generate_output = 'strain_xx strain_xy strain_yy strain_zz stress_xx stress_xy stress_yy stress_zz'
extra_vector_tags = 'ref'
eigenstrain_names = 'graphite_thermal_expansion'
block = 'upper_plunger lower_plunger die_wall'
[]
[yttria]
strain = FINITE
add_variables = true
use_automatic_differentiation = true
generate_output = 'strain_xx strain_xy strain_yy strain_zz stress_xx stress_xy stress_yy stress_zz'
extra_vector_tags = 'ref'
block = 'powder_compact'
[SolidMechanics]
[QuasiStatic]
[graphite]
strain = FINITE
add_variables = true
use_automatic_differentiation = true
generate_output = 'strain_xx strain_xy strain_yy strain_zz stress_xx stress_xy stress_yy stress_zz'
extra_vector_tags = 'ref'
eigenstrain_names = 'graphite_thermal_expansion'
block = 'upper_plunger lower_plunger die_wall'
[]
[yttria]
strain = FINITE
add_variables = true
use_automatic_differentiation = true
generate_output = 'strain_xx strain_xy strain_yy strain_zz stress_xx stress_xy stress_yy stress_zz'
extra_vector_tags = 'ref'
block = 'powder_compact'
[]
[]
[]
[]
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -59,7 +59,8 @@ initial_temperature = 300 #roughly 600C where the pyrometer kicks in
# order = FIRST
# block = 'upper_plunger lower_plunger die_wall'
# []
[microapp_potential] #converted to microapp electronVolts units
[microapp_potential]
#converted to microapp electronVolts units
[]
[E_x]
order = FIRST
Expand All @@ -84,24 +85,26 @@ initial_temperature = 300 #roughly 600C where the pyrometer kicks in
[]

[Physics]
[SolidMechanics/QuasiStatic]
[graphite]
strain = FINITE
add_variables = true
use_automatic_differentiation = true
generate_output = 'strain_xx strain_xy strain_yy strain_zz stress_xx stress_xy stress_yy stress_zz'
extra_vector_tags = 'ref'
eigenstrain_names = 'graphite_thermal_expansion'
block = 'upper_plunger lower_plunger die_wall'
[]
[yttria]
strain = FINITE
add_variables = true
use_automatic_differentiation = true
generate_output = 'strain_xx strain_xy strain_yy strain_zz stress_xx stress_xy stress_yy stress_zz'
extra_vector_tags = 'ref'
eigenstrain_names = 'yttria_thermal_expansion'
block = 'powder_compact'
[SolidMechanics]
[QuasiStatic]
[graphite]
strain = FINITE
add_variables = true
use_automatic_differentiation = true
generate_output = 'strain_xx strain_xy strain_yy strain_zz stress_xx stress_xy stress_yy stress_zz'
extra_vector_tags = 'ref'
eigenstrain_names = 'graphite_thermal_expansion'
block = 'upper_plunger lower_plunger die_wall'
[]
[yttria]
strain = FINITE
add_variables = true
use_automatic_differentiation = true
generate_output = 'strain_xx strain_xy strain_yy strain_zz stress_xx stress_xy stress_yy stress_zz'
extra_vector_tags = 'ref'
eigenstrain_names = 'yttria_thermal_expansion'
block = 'powder_compact'
[]
[]
[]
[]
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -49,24 +49,26 @@ initial_temperature = 300 #roughly 600C where the pyrometer kicks in
[]

[Physics]
[SolidMechanics/QuasiStatic]
[graphite]
strain = FINITE
add_variables = true
use_automatic_differentiation = true
generate_output = 'strain_xx strain_xy strain_yy strain_zz stress_xx stress_xy stress_yy stress_zz'
extra_vector_tags = 'ref'
eigenstrain_names = 'graphite_thermal_expansion'
block = 'upper_plunger lower_plunger die_wall'
[]
[yttria]
strain = FINITE
add_variables = true
use_automatic_differentiation = true
generate_output = 'plastic_strain_yy strain_xx strain_xy strain_yy strain_zz stress_xx stress_xy stress_yy stress_zz'
extra_vector_tags = 'ref'
eigenstrain_names = 'yttria_thermal_expansion'
block = 'powder_compact'
[SolidMechanics]
[QuasiStatic]
[graphite]
strain = FINITE
add_variables = true
use_automatic_differentiation = true
generate_output = 'strain_xx strain_xy strain_yy strain_zz stress_xx stress_xy stress_yy stress_zz'
extra_vector_tags = 'ref'
eigenstrain_names = 'graphite_thermal_expansion'
block = 'upper_plunger lower_plunger die_wall'
[]
[yttria]
strain = FINITE
add_variables = true
use_automatic_differentiation = true
generate_output = 'plastic_strain_yy strain_xx strain_xy strain_yy strain_zz stress_xx stress_xy stress_yy stress_zz'
extra_vector_tags = 'ref'
eigenstrain_names = 'yttria_thermal_expansion'
block = 'powder_compact'
[]
[]
[]
[]
Expand Down Expand Up @@ -223,7 +225,6 @@ initial_temperature = 300 #roughly 600C where the pyrometer kicks in
# relative_tolerance = 1e-20
# absolute_tolerance = 1e-8
# max_inelastic_increment = 0.000001

[]
[yttria_thermal_expansion]
type = ADComputeThermalExpansionEigenstrain
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -100,7 +100,8 @@ initial_temperature = 300 #roughly 600C where the pyrometer kicks in
# order = FIRST
# block = 'upper_plunger'
# []
[microapp_potential] #converted to microapp electronVolts units
[microapp_potential]
#converted to microapp electronVolts units
[]
[E_x]
order = FIRST
Expand All @@ -121,15 +122,17 @@ initial_temperature = 300 #roughly 600C where the pyrometer kicks in
[]

[Physics]
[SolidMechanics/QuasiStatic]
[graphite]
strain = FINITE
add_variables = true
use_automatic_differentiation = true
generate_output = 'strain_xx strain_xy strain_yy strain_zz stress_xx stress_xy stress_yy stress_zz'
extra_vector_tags = 'ref'
# eigenstrain_names = 'graphite_thermal_expansion'
block = 'upper_plunger'
[SolidMechanics]
[QuasiStatic]
[graphite]
strain = FINITE
add_variables = true
use_automatic_differentiation = true
generate_output = 'strain_xx strain_xy strain_yy strain_zz stress_xx stress_xy stress_yy stress_zz'
extra_vector_tags = 'ref'
# eigenstrain_names = 'graphite_thermal_expansion'
block = 'upper_plunger'
[]
[]
[]
[]
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -80,14 +80,16 @@ initial_temperature = 600 #roughly 600C where the pyrometer kicks in
[]

[Physics]
[SolidMechanics/QuasiStatic]
[graphite]
strain = FINITE
add_variables = true
use_automatic_differentiation = true
generate_output = 'plastic_strain_yy strain_xx strain_xy strain_yy strain_zz stress_xx stress_xy stress_yy stress_zz'
extra_vector_tags = 'ref'
eigenstrain_names = 'thermal_expansion'
[SolidMechanics]
[QuasiStatic]
[graphite]
strain = FINITE
add_variables = true
use_automatic_differentiation = true
generate_output = 'plastic_strain_yy strain_xx strain_xy strain_yy strain_zz stress_xx stress_xy stress_yy stress_zz'
extra_vector_tags = 'ref'
eigenstrain_names = 'thermal_expansion'
[]
[]
[]
[]
Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -63,7 +63,8 @@ initial_temperature = 300 #roughly 600C where the pyrometer kicks in
[specific_heat_capacity_va]
initial_condition = 842.2 # at 1500K #568.73 at 1000K #447.281 # at 293K
[]
[density_va] # initial condition set in ICs block below
[density_va]
# initial condition set in ICs block below
[]
[heat_transfer_radiation]
[]
Expand All @@ -82,16 +83,17 @@ initial_temperature = 300 #roughly 600C where the pyrometer kicks in
[]

[Physics]
[SolidMechanics/QuasiStatic]
[graphite]
strain = FINITE
incremental = true
add_variables = true
use_automatic_differentiation = true
generate_output = 'strain_xx strain_xy strain_yy strain_zz stress_xx stress_xy stress_yy '
'stress_zz plastic_strain_yy'
extra_vector_tags = 'ref'
eigenstrain_names = 'thermal_expansion'
[SolidMechanics]
[QuasiStatic]
[graphite]
strain = FINITE
incremental = true
add_variables = true
use_automatic_differentiation = true
generate_output = 'strain_xx strain_xy strain_yy strain_zz stress_xx stress_xy stress_yy stress_zz plastic_strain_yy'
extra_vector_tags = 'ref'
eigenstrain_names = 'thermal_expansion'
[]
[]
[]
[]
Expand Down

0 comments on commit d0a7328

Please sign in to comment.