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

Replace flux transfer in HPMR #285

Open
wants to merge 4 commits into
base: devel
Choose a base branch
from
Open

Replace flux transfer in HPMR #285

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
c4ef587
Manually editing the flux for transfer is brittle, and UO transfer is…
aprilnovak Sep 20, 2023
bed0f85
Use 3 cm mesh size.
aprilnovak Sep 20, 2023
ceb765b
Regold Sockeye test
aprilnovak Sep 20, 2023
70431b4
Actually need more than 1 time step, because the retain-solution feat…
aprilnovak Oct 5, 2023
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
2 changes: 1 addition & 1 deletion microreactors/mrad/heat_pipe_failure/HPMR_sockeye_tr.i
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -65,7 +65,7 @@ length_cond = 90.0e-2

# Mesh density
# The dimensions are nicely divisible by 3 cm mesh.
nelem_base_evap = 50
nelem_base_evap = 60
nelem_base_adia = 10
nelem_base_cond = 30
mesh_density = 3
Expand Down
21 changes: 2 additions & 19 deletions microreactors/mrad/heat_pipe_failure/HPMR_thermo_tr.i
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -5,14 +5,6 @@
## Thermal Only Physics ##
################################################################################

R_clad_o = 0.0105 # heat pipe outer radius
R_hp_hole = 0.0107 # heat pipe + gap
num_sides = 12 # number of sides of heat pipe as a result of mesh polygonization
alpha = '${fparse 2 * pi / num_sides}'
perimeter_correction = '${fparse 0.5 * alpha / sin(0.5 * alpha)}' # polygonization correction factor for perimeter
area_correction = '${fparse sqrt(alpha / sin(alpha))}' # polygonization correction factor for area
corr_factor = '${fparse R_hp_hole / R_clad_o * area_correction / perimeter_correction}'

[Problem]
restart_file_base = '../steady/HPMR_dfem_griffin_ss_out_bison0_cp/LATEST'
force_restart = true
Expand Down Expand Up @@ -83,9 +75,6 @@ corr_factor = '${fparse R_hp_hole / R_clad_o * area_correction / perimeter_corre
[flux_uo] #auxvariable to hold heat pipe surface flux from UserObject
block = 'reflector_quad monolith'
[]
[flux_uo_corr] #auxvariable to hold corrected flux_uo
block = 'reflector_quad monolith'
[]
[hp_temp_aux]
block = 'reflector_quad monolith'
[]
Expand Down Expand Up @@ -133,12 +122,6 @@ corr_factor = '${fparse R_hp_hole / R_clad_o * area_correction / perimeter_corre
variable = flux_uo
user_object = flux_uo
[]
[flux_uo_corr]
type = NormalizationAux
variable = flux_uo_corr
source_variable = flux_uo
normal_factor = '${fparse corr_factor}'
[]
[]

[BCs]
Expand Down Expand Up @@ -269,9 +252,9 @@ corr_factor = '${fparse R_hp_hole / R_clad_o * area_correction / perimeter_corre
execute_on = 'initial timestep_begin'
[]
[to_sockeye_flux]
type = MultiAppGeneralFieldNearestLocationTransfer
type = MultiAppGeneralFieldUserObjectTransfer
to_multi_app = sockeye
source_variable = flux_uo_corr
source_user_object = flux_uo
variable = master_flux
execute_on = 'initial timestep_begin'
[]
Expand Down
2 changes: 1 addition & 1 deletion microreactors/mrad/legacy/steady/MP_FC_ss_sockeye.i
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -26,7 +26,7 @@ length_cond = 90.0e-2

# Mesh density
# The dimensions are nicely divisible by 3 cm mesh.
nelem_base_evap = 50
nelem_base_evap = 60
nelem_base_adia = 10
nelem_base_cond = 30
mesh_density = 3
Expand Down
2 changes: 1 addition & 1 deletion microreactors/mrad/load_following/HPMR_sockeye_tr.i
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -65,7 +65,7 @@ length_cond = 90.0e-2

# Mesh density
# The dimensions are nicely divisible by 3 cm mesh.
nelem_base_evap = 50
nelem_base_evap = 60
nelem_base_adia = 10
nelem_base_cond = 30
mesh_density = 3
Expand Down
21 changes: 2 additions & 19 deletions microreactors/mrad/load_following/HPMR_thermo_tr.i
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -5,14 +5,6 @@
## Thermal Only Physics ##
################################################################################

R_clad_o = 0.0105 # heat pipe outer radius
R_hp_hole = 0.0107 # heat pipe + gap
num_sides = 12 # number of sides of heat pipe as a result of mesh polygonization
alpha = '${fparse 2 * pi / num_sides}'
perimeter_correction = '${fparse 0.5 * alpha / sin(0.5 * alpha)}' # polygonization correction factor for perimeter
area_correction = '${fparse sqrt(alpha / sin(alpha))}' # polygonization correction factor for area
corr_factor = '${fparse R_hp_hole / R_clad_o * area_correction / perimeter_correction}'

[Problem]
restart_file_base = '../steady/HPMR_dfem_griffin_ss_out_bison0_cp/LATEST'
force_restart = true
Expand Down Expand Up @@ -83,9 +75,6 @@ corr_factor = '${fparse R_hp_hole / R_clad_o * area_correction / perimeter_corre
[flux_uo] #auxvariable to hold heat pipe surface flux from UserObject
block = 'reflector_quad monolith'
[]
[flux_uo_corr] #auxvariable to hold corrected flux_uo
block = 'reflector_quad monolith'
[]
[hp_temp_aux]
block = 'reflector_quad monolith'
[]
Expand Down Expand Up @@ -133,12 +122,6 @@ corr_factor = '${fparse R_hp_hole / R_clad_o * area_correction / perimeter_corre
variable = flux_uo
user_object = flux_uo
[]
[flux_uo_corr]
type = NormalizationAux
variable = flux_uo_corr
source_variable = flux_uo
normal_factor = '${fparse corr_factor}'
[]
[]

[BCs]
Expand Down Expand Up @@ -269,9 +252,9 @@ corr_factor = '${fparse R_hp_hole / R_clad_o * area_correction / perimeter_corre
execute_on = 'initial timestep_begin'
[]
[to_sockeye_flux]
type = MultiAppGeneralFieldNearestLocationTransfer
type = MultiAppGeneralFieldUserObjectTransfer
to_multi_app = sockeye
source_variable = flux_uo_corr
source_user_object = flux_uo
variable = master_flux
execute_on = 'initial timestep_begin'
[]
Expand Down
5 changes: 2 additions & 3 deletions microreactors/mrad/steady/HPMR_sockeye_ss.i
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -64,7 +64,7 @@ length_cond = 90.0e-2

# Mesh density
# The dimensions are nicely divisible by 3 cm mesh.
nelem_base_evap = 50
nelem_base_evap = 60
nelem_base_adia = 10
nelem_base_cond = 30
mesh_density = 3
Expand Down Expand Up @@ -427,8 +427,7 @@ q_evap = '${fparse Q_hp / S_evap}'

start_time = -5e4 # negative start time so we can start running from t = 0
end_time = 0
dtmin = 1
dt = 1000
dt = 5e3
[]

[Outputs]
Expand Down
24 changes: 3 additions & 21 deletions microreactors/mrad/steady/HPMR_thermo_ss.i
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -5,14 +5,6 @@
## Thermal Only Physics ##
################################################################################

R_clad_o = 0.0105 # heat pipe outer radius
R_hp_hole = 0.0107 # heat pipe + gap
num_sides = 12 # number of sides of heat pipe as a result of mesh polygonization
alpha = '${fparse 2 * pi / num_sides}'
perimeter_correction = '${fparse 0.5 * alpha / sin(0.5 * alpha)}' # polygonization correction factor for perimeter
area_correction = '${fparse sqrt(alpha / sin(alpha))}' # polygonization correction factor for area
corr_factor = '${fparse R_hp_hole / R_clad_o * area_correction / perimeter_correction}'

[GlobalParams]
flux_conversion_factor = 1
[]
Expand Down Expand Up @@ -82,9 +74,6 @@ corr_factor = '${fparse R_hp_hole / R_clad_o * area_correction / perimeter_corre
[flux_uo] #auxvariable to hold heat pipe surface flux from UserObject
block = 'reflector_quad monolith'
[]
[flux_uo_corr] #auxvariable to hold corrected flux_uo
block = 'reflector_quad monolith'
[]
[hp_temp_aux]
block = 'reflector_quad monolith'
[]
Expand Down Expand Up @@ -132,12 +121,6 @@ corr_factor = '${fparse R_hp_hole / R_clad_o * area_correction / perimeter_corre
variable = flux_uo
user_object = flux_uo
[]
[flux_uo_corr]
type = NormalizationAux
variable = flux_uo_corr
source_variable = flux_uo
normal_factor = '${fparse corr_factor}'
[]
[]

[BCs]
Expand Down Expand Up @@ -271,10 +254,10 @@ corr_factor = '${fparse R_hp_hole / R_clad_o * area_correction / perimeter_corre
[]
[to_sockeye_flux]
type = MultiAppGeneralFieldUserObjectTransfer
variable = master_flux
to_multi_app = sockeye
execute_on = 'timestep_begin'
source_user_object = flux_uo
variable = master_flux
execute_on = 'timestep_begin'
[]
[]

Expand Down Expand Up @@ -302,8 +285,7 @@ corr_factor = '${fparse R_hp_hole / R_clad_o * area_correction / perimeter_corre

start_time = -5e4 # negative start time so we can start running from t = 0
end_time = 0
dtmin = 1
num_steps = 1
dt = 5e3
[]

[Postprocessors]
Expand Down
Loading