...
 
Commits (3)
......@@ -27,8 +27,8 @@
[Functions]
[./time_stepper]
type = PiecewiseLinear
x = '-100 -45 -44 -1 0 2 3 10 11 1e5'
y = ' 0.1 0.1 1 1 1 0.2 0.5 0.5 2 2'
x = '-100 -45 -44 -1 0 2 2.2 10 11 100 101 1e5'
y = ' 0.1 0.1 1 1 0.1 0.1 0.2 0.2 1 1 2 2'
[../]
[./pump_coastdown]
......@@ -37,17 +37,13 @@
4.00E+00 4.40E+00 4.80E+00 5.20E+00 5.60E+00 6.00E+00 6.40E+00 6.80E+00 7.20E+00 7.60E+00
8.000E+00 1.000E+01 2.000E+01 3.000E+01 4.000E+01 5.000E+01 6.000E+01 7.000E+01 8.000E+01 9.000E+01
1.000E+02 1.100E+02 1.200E+02 1.300E+02 1.400E+02 1.500E+02 1.600E+02 1.700E+02 1.800E+02 1.900E+02
2.000E+02 2.100E+02 2.200E+02 2.300E+02 2.400E+02 2.500E+02 2.600E+02 2.700E+02 2.800E+02 2.900E+02
3.000E+02 3.100E+02 3.200E+02 3.300E+02 3.400E+02 3.500E+02 3.600E+02 3.700E+02 3.800E+02 3.900E+02
4.000E+02 4.100E+02 4.200E+02 1.00E+05'
2.000E+02 2.100E+02 2.200E+02 2.300E+02 2.400E+02 2.500E+02 2.600E+02 2.700E+02 2.800E+02 1.00E+05'
y ='1.000E+00 1.000E+00 9.246E-01 8.573E-01 7.972E-01 7.432E-01 6.944E-01 6.504E-01 6.104E-01 5.739E-01 5.407E-01
5.102E-01 4.823E-01 4.565E-01 4.328E-01 4.109E-01 3.906E-01 3.718E-01 3.543E-01 3.380E-01 3.228E-01
3.086E-01 2.500E-01 1.111E-01 6.250E-02 4.000E-02 2.778E-02 2.041E-02 1.563E-02 1.235E-02 1.000E-02
8.264E-03 6.944E-03 5.917E-03 5.102E-03 4.444E-03 3.906E-03 3.460E-03 3.086E-03 2.770E-03 2.500E-03
2.268E-03 2.066E-03 1.890E-03 1.736E-03 1.600E-03 1.479E-03 1.372E-03 1.276E-03 1.189E-03 1.111E-03
1.041E-03 9.766E-04 9.083E-04 8.651E-04 8.163E-04 7.716E-04 7.305E-04 6.925E-04 6.575E-04 6.250E-04
5.949E-04 5.669E-04 0 0'
scale_factor = 1e3
2.268E-03 2.066E-03 1.890E-03 1.736E-03 1.600E-03 1.479E-03 1.372E-03 1.276E-03 0 0'
scale_factor = 1e4
[../]
[]
......@@ -61,11 +57,11 @@
[]
[Components]
[./pipe1] # core inlet pipe
[./pipe1]
type = PBOneDFluidComponent
eos = eos
position = '0 -1 0'
orientation = '0 1 0'
position = '0 1 0'
orientation = '0 -1 0'
A = 0.2376
Dh = 0.55
length = 1
......@@ -76,7 +72,7 @@
type = PBOneDFluidComponent
eos = eos
position = '0 0 3'
orientation = '0 -1 0'
orientation = '0 1 0'
A = 0.1590
Dh = 0.45
......@@ -123,17 +119,19 @@
[../]
[./Branch4]
type = PBSingleJunction
type = PBBranch
inputs = 'pipe3(out)'
outputs = 'IHX(primary_in)'
K = '0.0 0.0'
Area = 0.50265
eos = eos
[../]
[./Branch5]
type = PBBranch
inputs = 'IHX(primary_out)'
outputs = 'pipe4(in)'
K = '0. 0.' # Form loss coefficient at pump inlet and outlet
Area = 0.2376 # Reference pump flow area
K = '0.0 0.0'
Area = 0.50265
eos = eos
[../]
......@@ -145,8 +143,8 @@
K = '0. 0.' # Form loss coefficient at pump inlet and outlet
Area = 0.2376 # Reference pump flow area
initial_P = 1.5e5 # Initial pressure
Head = 1e3 # Pump head, Pa
# Head_fn = pump_coastdown
# Head = 1e4 # Pump head, Pa
Head_fn = pump_coastdown
[../]
......@@ -192,6 +190,11 @@
variable = velocity
input = pipe1(out)
[../]
[./IHX_Tin]
type = ComponentBoundaryVariableValue
variable = temperature
input = IHX(primary_in)
[../]
[./IHX_Tout]
type = ComponentBoundaryVariableValue
variable = temperature
......@@ -231,7 +234,7 @@
[../]
start_time = -50
num_steps = 1000
end_time = 0
end_time = 500
nl_rel_tol = 1e-6
nl_abs_tol = 1e-6
......@@ -252,7 +255,7 @@
use_displaced = true
execute_on = 'initial timestep_end'
sequence = false
interval = 10
interval = 10
[../]
[./console]
type = Console
......
......@@ -27,8 +27,8 @@
[Functions]
[./time_stepper]
type = PiecewiseLinear
x = '-100 -45 -44 -1 0 2 3 10 11 1e5'
y = ' 0.1 0.1 1 1 1 0.2 0.5 0.5 2 2'
x = '-100 -45 -44 -1 0 2 2.2 10 11 100 101 1e5'
y = ' 0.1 0.1 1 1 0.1 0.1 0.2 0.2 1 1 2 2'
[../]
[./pump_coastdown]
......@@ -37,17 +37,13 @@
4.00E+00 4.40E+00 4.80E+00 5.20E+00 5.60E+00 6.00E+00 6.40E+00 6.80E+00 7.20E+00 7.60E+00
8.000E+00 1.000E+01 2.000E+01 3.000E+01 4.000E+01 5.000E+01 6.000E+01 7.000E+01 8.000E+01 9.000E+01
1.000E+02 1.100E+02 1.200E+02 1.300E+02 1.400E+02 1.500E+02 1.600E+02 1.700E+02 1.800E+02 1.900E+02
2.000E+02 2.100E+02 2.200E+02 2.300E+02 2.400E+02 2.500E+02 2.600E+02 2.700E+02 2.800E+02 2.900E+02
3.000E+02 3.100E+02 3.200E+02 3.300E+02 3.400E+02 3.500E+02 3.600E+02 3.700E+02 3.800E+02 3.900E+02
4.000E+02 4.100E+02 4.200E+02 1.00E+05'
2.000E+02 2.100E+02 2.200E+02 2.300E+02 2.400E+02 2.500E+02 2.600E+02 2.700E+02 2.800E+02 1.00E+05'
y ='1.000E+00 1.000E+00 9.246E-01 8.573E-01 7.972E-01 7.432E-01 6.944E-01 6.504E-01 6.104E-01 5.739E-01 5.407E-01
5.102E-01 4.823E-01 4.565E-01 4.328E-01 4.109E-01 3.906E-01 3.718E-01 3.543E-01 3.380E-01 3.228E-01
3.086E-01 2.500E-01 1.111E-01 6.250E-02 4.000E-02 2.778E-02 2.041E-02 1.563E-02 1.235E-02 1.000E-02
8.264E-03 6.944E-03 5.917E-03 5.102E-03 4.444E-03 3.906E-03 3.460E-03 3.086E-03 2.770E-03 2.500E-03
2.268E-03 2.066E-03 1.890E-03 1.736E-03 1.600E-03 1.479E-03 1.372E-03 1.276E-03 1.189E-03 1.111E-03
1.041E-03 9.766E-04 9.083E-04 8.651E-04 8.163E-04 7.716E-04 7.305E-04 6.925E-04 6.575E-04 6.250E-04
5.949E-04 5.669E-04 0 0'
scale_factor = 1e3
2.268E-03 2.066E-03 1.890E-03 1.736E-03 1.600E-03 1.479E-03 1.372E-03 1.276E-03 0 0'
scale_factor = 1e4
[../]
[]
......@@ -61,11 +57,11 @@
[]
[Components]
[./pipe1] # core inlet pipe
[./pipe1]
type = PBOneDFluidComponent
eos = eos
position = '0 -1 0'
orientation = '0 1 0'
position = '0 1 0'
orientation = '0 -1 0'
A = 0.2376
Dh = 0.55
length = 1
......@@ -76,7 +72,7 @@
type = PBOneDFluidComponent
eos = eos
position = '0 0 3'
orientation = '0 -1 0'
orientation = '0 1 0'
A = 0.1590
Dh = 0.45
......@@ -123,17 +119,19 @@
[../]
[./Branch4]
type = PBSingleJunction
type = PBBranch
inputs = 'pipe3(out)'
outputs = 'IHX(primary_in)'
K = '0.0 0.0'
Area = 0.50265
eos = eos
[../]
[./Branch5]
type = PBBranch
inputs = 'IHX(primary_out)'
outputs = 'pipe4(in)'
K = '0. 0.' # Form loss coefficient at pump inlet and outlet
Area = 0.2376 # Reference pump flow area
K = '0.0 0.0'
Area = 0.50265
eos = eos
[../]
......@@ -145,8 +143,8 @@
K = '0. 0.' # Form loss coefficient at pump inlet and outlet
Area = 0.2376 # Reference pump flow area
initial_P = 1.5e5 # Initial pressure
Head = 1e3 # Pump head, Pa
# Head_fn = pump_coastdown
# Head = 1e4 # Pump head, Pa
Head_fn = pump_coastdown
[../]
......@@ -192,6 +190,11 @@
variable = velocity
input = pipe1(out)
[../]
[./IHX_Tin]
type = ComponentBoundaryVariableValue
variable = temperature
input = IHX(primary_in)
[../]
[./IHX_Tout]
type = ComponentBoundaryVariableValue
variable = temperature
......@@ -226,7 +229,7 @@
[../]
start_time = -50
num_steps = 1000
end_time = 0
end_time = 500
nl_rel_tol = 1e-6
nl_abs_tol = 1e-6
......@@ -247,7 +250,7 @@
use_displaced = true
execute_on = 'initial timestep_end'
sequence = false
interval = 10
interval = 10
[../]
[./console]
type = Console
......