Thermohydraulic Analysis of the ITER Magnet Systems

Appl. Math. Dept. of STC “SINTEZ”, Efremov Research Institute, St. Petersburg, RF
Thermohydraulic Analysis of the
ITER Magnet Systems
• Toroidal Fields Coils
• Winding & Case
• Cryoplant interface
• Poloidal Field Coils
• Winding
• Cryoplant interface
• Central Solenoid
• Winding
• Cryoplant interface

Numerical Simulation of the TF
Magnet Cooling System: Model
• Hydraulic Layout
(Two closed cooling circuits with the cryolines, feeders, pumps, valves
and heat exchangers for the TF Winding and the TF Case separately )
• Winding cooling circuit
• 1-dimensional dimensional finite element model for the transient
compressed flows of SHe in the conductor channels,
cooling pipes and cryolines
• 7 pairs of channels for the individual modeling of 7
pancakes (half of the TF winding)
• Opposite He flow direction in a double pancake
• Case cooling circuit
• 25+25 (for half of the coil) TF case cooling pipes for two
separate refrigerating circuits controlled by the valves
• Coil structure
Appl. Math. Dept. of STC “SINTEZ”, Efremov Research Institute, St. Petersburg, RF
• 2-dimensional dimensional finite element model for the transient
thermal problem in the TF Coil cross-sections
cross sections
• 32 variable TF Coil cross sections (case & winding) in
poloidal direction
• Space and time distribution of the heat loads along and
across the TF coil
• Space and time distribution of AC losses and nuclear
heating in the pancakes
VINCENTA v.4.1 code

Temperature [K]
5.2
5.1
5
4.9
4.8
4.7
4.6
4.5
turn_1 turn_2
Appl. Math. Dept. of STC “SINTEZ”, Efremov Research Institute, St. Petersburg, RF
Numerical Simulation of the TF
Magnet Cooling System: Results
The general numerical model is treated by VINCENTA v.4.1 code
Winding cooling circuit
4.4
0 50 100 150 200 250 300 350
6.5 s
100 s
130 s
230 s
330 s
430 s
530 s
630 s
725 s
825 s
925 s
1245 s
1590 s
1690 s
1790 s
1800 (0) s
Channel Length [m]
Normal
operation
without
controllability
Temperature evolution of the TF conductor for
pancake #5 during plasma pulse (“odd” pancake,
C63) (no controllability).
Temperature [K]
6.2
6
5.8
5.6
5.4
5.2
5
4.8
4.6
Case cooling circuit
4.4
0 2 4 6 8 10 12 14 16
6.5 s
100 s
130 s
Channel Length [m]
Normal
230 s
330 s
430 s
530 s operation
630 s
725 s
825 s
925 s
1245 s
1590 s
1690 s
1790 s
1800 (0) s
without
controllability
Temperature distribution along the "long"
cooling pass #4 (“front” wall) during 1800 s
plasma pulse (no controllability)
8
7
6
5
Coil structure
Normal
operation
without
controllability
Temperature diagram for section # 9 (middle plane) at 530 s (end of
plasma burning).

Temperature [K]
5.3
5.2
5.1
5
4.9
4.8
4.7
4.6
4.5
4.4
Appl. Math. Dept. of STC “SINTEZ”, Efremov Research Institute, St. Petersburg, RF
Numerical Simulation of the TF
Magnet Cooling System: Results
The general numerical model is treated by VINCENTA v.4.1 code
Winding cooling circuit
turn_1 turn_2
4.3
0 50 100 150 200 250 300 350
6.5 s
100 s
130 s
230 s
330 s
430 s
530 s
630 s
725 s
825 s
925 s
1245 s
1590 s
1690 s
1790 s
1800 (0) s
Channel Length [m]
Normal
operation
with
controllability
Temperature evolution of the TF conductor for
pancake #5 during plasma pulse (“odd” pancake,
C63).
Temperature [K]
6.6
6.4
6.2
6
5.8
5.6
5.4
5.2
5
4.8
4.6
4.4
Case cooling circuit
4.2
0 2 4 6 8 10 12 14 16
6.5 s
100 s
130 s
230 s
330 s
430 s
530 s
630 s
725 s
825 s
925 s
1245 s
1590 s
1690 s
1790 s
1800 (0) s
Channel Length [m]
Normal
operation
with
controllability
Temperature distribution along the "long"
cooling pass #4 (“front” wall) during 1800 s
plasma pulse.
8
7
6
5
Coil structure
Normal
operation
with
controllability
Temperature diagram for section # 9 (middle plane) at 530 s (end of
plasma burning).

Temperature [K]
5.8
5.6
5.4
5.2
5
4.8
4.6
4.4
Appl. Math. Dept. of STC “SINTEZ”, Efremov Research Institute, St. Petersburg, RF
Numerical Simulation of the TF
Magnet Cooling System: Results
The general numerical model is treated by VINCENTA v.4.1 code
Winding cooling circuit
turn_1 turn_2
4.2
0 50 100 150 200 250 300 350
Channel Length [m]
6.5 s
100 s
130 s
230 s
330 s
430 s
530 s
630 s
725 s
825 s
925 s
1245 s
1590 s
1690 s
1790 s
1800 (0) s
Temperature evolution of the TF conductor for
pancake #5 (“odd” pancake, C63) during 9 th
plasma pulse ended by plasma disruption.
Temperature [K]
9.5
9
8.5
8
7.5
7
6.5
6
5.5
5
4.5
Case cooling circuit
4
0 2 4 6 8 10 12 14 16
Channel Length [m]
6.5 s
100 s
130 s
230 s
330 s
430 s
530 s
630 s
725 s
825 s
925 s
1245 s
1590 s
1690 s
1790 s
1800 (0) s
Temperature distribution along the "long"
cooling pass #4 (“front” wall) during 9 th plasma
pulse ended by plasma disruption.
14
12
10
8
6
Coil structure
Operation with plasma disruption
Temperature diagram for section # 9 (middle plane) at 590 s (one
minute pass after 9 th plasma pulse ended by plasma disruption).

Temperature (K)
6
5.8
5.6
5.4
5.2
5
4.8
4.6
4.4
Appl. Math. Dept. of STC “SINTEZ”, Efremov Research Institute, St. Petersburg, RF
Numerical Simulation of the TF
Magnet Cooling System: Results
The general numerical model is treated by VINCENTA v.4.1 code
16 246
4.2
0 50 100 150 200 250
6.5 s
100 s
130 s
230 s
330 s
430 s
530 s
540 s
550 s
560 s
570 s
580 s
590 s
650 s
950 s
1800 (0) s
x (m)
Winding cooling circuit
Pressure (MPa)
0.8
0.75
0.7
0.65
0.6
0.55
0.5
16 246
0.45
0 50 100 150 200 250 300
6.5 s
100 s
130 s
230 s
330 s
430 s
530 s
630 s
725 s
825 s
925 s
1245 s
1590 s
1690 s
1790 s
1800 (0) s
SHe temperature and Pressure evolution in the feeder C73, cryoline C74, heat exchanger C75,
cryoline C76 and feeder C77 during 9 th plasma pulse ended by plasma disruption.
x (m)
Temperature (K)
4
0 20 40 60 80 100 120 140 160
Operation with plasma disruption
8
7.5
7
6.5
6
5.5
5
4.5
6.5 s
100 s
130 s
230 s
330 s
430 s
530 s
540 s
550 s
560 s
570 s
580 s
590 s
650 s
950 s
1800 (0) s
21 131
x (m)
Case cooling circuit
0.5
0 20 40 60 80 100 120 140 160
SHe temperature and pressure evolution in the feeder C53, cryoline C54, heat exchanger
C55, cryoline C56 and feeder C57 during 9 th plasma pulse ended by plasma disruption.
Pressure [MPa]
1.7
1.6
1.5
1.4
1.3
1.2
1.1
1
0.9
0.8
0.7
0.6
21 131
6.5 s
100 s
130 s
230 s
330 s
430 s
530 s
630 s
725 s
825 s
925 s
1245 s
1590 s
1690 s
1790 s
1800 (0) s
x (m)

Appl. Math. Dept. of STC “SINTEZ”, Efremov Research Institute, St. Petersburg, RF
Numerical Simulation of the TF
Magnet Cooling System: Results
The general numerical model is treated by VINCENTA v.4.1 code
Normal operation under Cryoplant Controllability (TF Case circuit circuit
adjusted only)
HEAT LOAD [kW] Normal
PRESSURE [MPa]
20
15
10
5
Winding cooling circuit
0
0 1
Q_cryoplant
Q_winding & pump
Q_average
2 3 4
PULSE
5 6 7 8
0.7
0.65
0.6
0.55
0.5
0.45
0.4
0.35
Evolution of the heat released in the cooling circuit and
absorbed by heat exchanger.
0.3
0 1 2 3 4 5 6 7 8
After the pump
Before the pump
PULSE
Evolution of the He pressure before and after the
pump.
HEAT LOAD [kW]
PRESSURE [MPa]
30
25
20
15
10
5
Case cooling circuit
0
0 1
Q_cryoplant
Q_case & pump
Q_average
2 3 4
PULSE
5 6 7 8
Evolution of the heat released in the cooling circuit and
absorbed by heat exchanger.
1
0.9
0.8
0.7
0.6
0.5
0.4
0.3
0 1 2 3 4 5 6 7 8
After the pump
Before the pump
Evolution of the He pressure before and after the
pump.
PULSE
HEAT LOAD [kW]
Total: Winding & Coil structure
30
25
20
15
10
5
0
0 1
Q_cryoplant (total)
Q_cryoplant (average)
2 3 4
PULSE
5 6 7 8
MASS FLOW RATE (g/s)
Evolution of the total power absorbed by both heat
exchangers.
4000
3500
3000
2500
2000
1500
1000
500
0
0 1 2 3 4 5 6 7 8
PULSE
pump
control valve
bypass valve
Evolution of the He mass flow rate through the pump,
control and bypass valves (TF Case cooling circuit).

Appl. Math. Dept. of STC “SINTEZ”, Efremov Research Institute, St. Petersburg, RF
Numerical Simulation of the PF
Magnet Cooling System: Model
Thermohydraulic Model
• 1-dimensional dimensional finite element model for the
transient compressed flows of SHe in the
conductor channels, cooling pipes and cryolines
• 17 pairs of channels for the individual
modeling of 6 PF coils
• individual distribution of the heat loads in
space and time among the conductors including
AC losses, Nuclear heating, joint losses etc.
• cryoplant interface: centrifugal pump, heat
exchanger, controll & bypass valves
PF conductors
PF1&6 PF2,3&4 PF5
coolant normal / backup inlet 4.7K /4.4K inlet 4.7K inlet 4.7K
Type of Strand NbTi NbTi NbTi
Operating Current (kA) normal / backup 45 / 52 45 / 52 45 / 52
Nominal Peak Field (T) normal / backup 6.0 / 6.4 4.0 5.0
Operating Temperature (K)
normal / backup 5.0 / 4.7 5.0 5.0
Equiv. Disch. Time Constant (s) hot spot 18 18 18
Tcs (K) normal / backup 6.5 / 6.27 6.65 / 6.51 6.60 / 6.51
Iop/Ic normal / backup 0.127 / 0.144 0.365 / 0.422 0.264 / 0.305
Cable diameter (mm) 38.2 34.5 35.4
Central spiral od x id (mm) 12 x 10 12 x 10 12 x 10
Conductor OD (mm) 53.8 x 53.8 52.3 x 53.2 51.9 x 51.9
Jacket steel steel steel
sc strand diam (mm) 0.73 0.73 0.72
sc strand cu : non-cu 1.6 6.9 4.4
3x4x4x5x6 ((3x3x4+1) ((3x3x4+1)
cabling pattern
x4+1)x6 x5+1)x6
sc strand Nr 1440 864 1080
Cu core 2/3/4 stage (mm) 0/0/0 0 / 1.8 / 3.5 0 /1.2 /2.7
ocal Void Fraction (%) in strand bundle 34.5 34.2 34.3
Helium in Annulus (mm2 ) 351.1 277.3 294.8
Helium in strand bundle (mm2) 334.5 261.2 278.6
Total Helium Area (mm2 ) 429.6 355.8 373.3
SC strand total perimeter (m)
3.302
1.981
2.443
[twisted]
[3.390]
[2.034]
[2.508]
A-ncu (mm2 )
229.3
45.2
80.5
[twisted]
[241.3]
[47.6]
[84.8]
A-cu sc str (mm2 )
366.8
312.4
354.3
[twisted]
[386.1]
[328.8]
[373.0]
A-cu extra (mm2 )
0
118.4
67.9
[twisted]
[0]
[120.9]
[69.5]
SC strand weight/m of conductor (kg/m) 4.885 2.931 3.564
? P/L (Pa/m) at 5K@5bar 8 g/s 66.9 69.1 71.4
J cable space (A/mm2) normal / backup 39.26 / 45.37 48.14 / 55.62 45.72 / 52.83
Conductor Cost (IUA/m) 1.26 0.97 1.05
Channel geometry and hydraulic parameters
Channel Cross section
#
area, mm 2
Hydraulic Length,
Comments
ID, mm m
C1-C2 351.1 0.42 196 Cable Spaces of PF1 cables
C3-C5 277.3 0.54 280 Cable Spaces of PF2 cables
C6-C8 277.3 0.54 441 Cable Spaces of PF3 cables
C9-C11 277.3 0.54 144 Cable Spaces of PF4 cables
C12-C14 294.8 0.47 144 Cable Spaces of PF5 cables
C15-C17 351.1 0.42 144 Cable Spaces of PF6 cables
C18- C34 78.5 10 - Central channel of PF1-PF6 cables
C41, C43 5030 80 80 Return & Supply Cryolines
C42 20000 20 20 Heat Exchanger
C35-40, C44-49 1257 40 16 Return & Supply pipes
VINCENTA v.4.code

MASS FLOW RATE (g/s)
2000
1500
1000
500
0
Appl. Math. Dept. of STC “SINTEZ”, Efremov Research Institute, St. Petersburg, RF
Numerical Simulation of the PF
Magnet Cooling System: Results
500
0 0.5 1 1.5 2 2.5 3
PULSE
3.5 4 4.5 5 5.5 6
pump
control valve
bypass valve
Evolution of the SHe mass flow rate through the pump,
the control and by-pass valves
Evolution of the SHe pressure along the PF cooling
loop for the repetitive pulsing mode.
Generalized view of the temperature evolution for
all simulated PF cables.
VINCENTA v.4.1 code

Appl. Math. Dept. of STC “SINTEZ”, Efremov Research Institute, St. Petersburg, RF
Numerical Simulation of the PF
Magnet Cooling System: Results
Temperature (K)
Mass flow rate (g/s)
5.2
5
4.8
4.6
0 50 100 150 200
10
9
8
7
Inlet
Outlet
Time (min)
PF3: “top”
6
0 50 100 150 200
Inlet
Outlet
Time (min)
PF3: “top”
Temperature (K)
Mass flow rate (g/s)
4.85
4.8
4.75
4.7
4.65
4.6
0 50 100 150 200
10
9
8
7
Inlet
Outlet
Time (min)
PF3: “regular”
6
0 50 100 150 200
Inlet
Outlet
Time (min)
PF3: “regular”
4.6
0 50 100 150 200
Evolution of the SHe temperature and mass flow rate at the inlet/outlet of the PF3 conductors conductors
Temperature (K)
Mass flow rate (g/s)
5.1
5
4.9
4.8
4.7
10
9
8
7
Inlet
Outlet
Time (min)
PF3: “bottom”
6
0 50 100 150 200
Inlet
Outlet
Time (min)
PF3: “bottom”
VINCENTA v.4.1 code

Numerical Simulation of the CSMC
& CS Insert Cooling system: Model
Thermohydraulics:
Thermohydraulics
Appl. Math. Dept. of STC “SINTEZ”, Efremov Research Institute, St. Petersburg, RF
• 1-dimensional dimensional finite element model for the
transient compressed flows of SHe in the
conductor channels, cooling pipes and cryolines
• 18 pairs of conductors for modeling of 10+8
solenoid layers
• distribution of AC losses in space and time of
solenoid conductors
• two separate cooling loop for Inner & Outer
Modules controlled by valves
• cryoplant interface: centrifugal pump, heat
exchanger, control & bypass valves
Coil structure:
Joint
• 2-dimensional dimensional finite
element model for the
transient thermal problem in
the CSMC cross-sections
cross sections
• 5 cross sections in
circumference Insert
Supporting
Structure
Inner Module
Outer Module
Buffer Spacer
1m
VINCENTA v.4.1

Mass flow rate [g/s]
Temperature [K]
16
14
12
10
8
6
Appl. Math. Dept. of STC “SINTEZ”, Efremov Research Institute, St. Petersburg, RF
Numerical Simulation of the CSMC
& CS Insert Cooling system : Results
4
0 200 400 600 800 1000 1200 1400 1600 1800 2000
X [m]
0 s (30 min)
32.5 s
37.5 s
1 min
2 min
6 min
10 min
14 min
18 min
22 min
26 min
Generalized view of He mass flow rate evolution along the
conductors of layers number 1 through 18. Option 1.
6
5.8
5.6
5.4
5.2
5
4.8
4.6
4.4
4.2
0 200 400 600 800 1000 1200 1400 1600 1800 2000
X [m]
0 s (30 min)
32.5 s
37.5 s
1 min
2 min
6 min
10 min
14 min
18 min
22 min
26 min
Generalized view of conductor temperature evolution for
layers number 1 through 18. Option 1.
Temperature diagram for CSMC Winding & Buffer zone at 6th minute of pulse during repetitive pulsing mode.
•Option 1: mass flow rate is 370 g/s, effective time constant nτ is 50 ms, and repetitive time of pulsing trep is 30
minutes.
•Option 2: mass flow rate is 185 g/s and another parameters as Option 1.
•Option 3: effective time constant nτ is 100 ms and another parameters as Option 1.
•Option 4: repetitive time of pulsing trep is 20 minutes and another parameters as Option 1.
VINCENTA v.4.1

Pressure (MPa)
1
0.95
0.9
0.85
0.8
0.75
0.7
0.65
0.6
0.55
0.5
Appl. Math. Dept. of STC “SINTEZ”, Efremov Research Institute, St. Petersburg, RF
Validation of the CSMC Model:
Calculating & Experimental Data
Calculation
0 2000 4000 6000 8000 1 10 4
0.45
Time (sec)
Coil Inlet
Coil Outlet
Experiment with the series of current pulses 30 kA, 40 kA, 43 kA and 46 kA
Evolution of pressure at inlet and outlet of
the CSMC for series of current pulses.
Pressure (MPa)
1.2 10 4
1
0.95
0.9
0.85
0.8
0.75
0.7
0.65
0.6
0.55
0.5
Experiment
0 2000 4000 6000 8000 1 10 4
0.45
Time (sec)
CSV_PT_PI_CB40X
CSV_PT_PI_CB40E
1.2 10 4
Temperature (K)
5.8
5.6
5.4
5.2
5
4.8
4.6
4.4
Calculation
0 2000 4000 6000 8000 1 10 4
Time (sec)
Coil Inlet
Coil Outlet
Temperature (K)
1.2 10 4
5.8
5.6
5.4
5.2
4.8
4.6
4.4
0 2000 4000 6000 8000 1 10 4
Time (sec)
Evolution of temperatures at inlet and outlet of
the CSMC for series of current pulses.
VINCENTA v.4.1
5
Experiment
CSV_TC_TB_CB40X
CSV_TC_TB_CB40E
1.2 10

Temperature (K)
6
5.8
5.6
5.4
5.2
5
4.8
4.6
Appl. Math. Dept. of STC “SINTEZ”, Efremov Research Institute, St. Petersburg, RF
Validation of the CSMC Model:
Calculating & Experimental Data
Calculation
0 2000 4000 6000 8000 1 10 4
4.4
Time (sec)
Inner Module Inlet
Inner Module Outlet
Experiment with the series of current pulses 30 kA, 40 kA, 43 kA and 46 kA
1.2 10 4
Experiment
Evolution of temperatures at inlet and
outlet of the inner CSMC module.
Temperature (K)
6
5.8
5.6
5.4
5.2
5
4.8
4.6
0 2000 4000 6000 8000 1 10 4
4.4
Time (sec)
CSV_TC_TB_CS1E
CSV_TC_TB_CS1X
Temperature (K)
1.2 10 4
5.3
5.2
5.1
5
4.9
4.8
4.7
4.6
4.5
Calculation
0 2000 4000 6000 8000 1 10 4
4.4
Time [sec]
Outer Module Inlet
Outer Module Outlet
Temperature (K)
1.2 10 4
5.3
5.2
5.1
4.9
4.8
4.7
4.6
4.5
Experiment
0 2000 4000 6000 8000 1 10 4
4.4
Time (sec)
Evolution of temperatures at inlet and outlet of
the outer CSMC module.
VINCENTA v.4.1
5
CSV_TC_TB_CS2E
CSV_TC_TB_CS2X
1.2 10

Temperature [K]
7
6.5
6
5.5
5
4.5
Appl. Math. Dept. of STC “SINTEZ”, Efremov Research Institute, St. Petersburg, RF
Validation of the CSMC Model:
Calculating & Experimental Data
Calculation
0 500 1000 1500 2000 2500 3000
1A
2A
3A
4A
5A
6A
7A
8A
9A
10A
Experiment with the series of current pulses 30 kA, 40 kA, 43 kA and 46 kA
Time (s)
Experiment
0 500 1000 1500 2000 2500 3000
Evolution of temperatures at outlet of 10 (A)
conductors of the inner CSMC module for 46 kA
current pulse.
Temperature (K)
7
6.5
6
5.5
5
4.5
Time (sec)
MCI_TC_01AO
MCI_TC_02AO
MCI_TC_03AO
MCI_TC_04AO
MCI_TC_05AO
MCI_TC_06AO
MCI_TC_07AO
MCI_TC_08AO
MCI_TC_09AO
MCI_TC_10AO
Temperature [K]
6.2
6
5.8
5.6
5.4
5.2
5
4.8
4.6
Calculation 5.8
Experiment
4.4
0 500 1000 1500 2000 2500 3000
11A
12A
13A
14A
15A
16A
17A
18A
Time (s)
Temperature (K)
4.4
0 500 1000 1500 2000 2500 3000
Evolution of the She temperatures at outlet of 8
(A) conductors of the outer CSMC module
for 46 kA current pulse.
6.2
6
5.6
5.4
5.2
5
4.8
4.6
Time (sec)
MCO_TC_11AO
MCO_TC_12AO
MCO_TC_13AO
MCO_TC_14AO
MCO_TC_15AO
MCO_TC_16AO
MCO_TC_17AO
MCO_TC_18AO
VINCENTA v.4.1

T, K 70
65
60
55
50
45
40
35
30
25
20
15
10
5
Numerical Simulation of the
CICC Conductor Performance
Quench Analysis
1.0s
Appl. Math. Dept. of STC “SINTEZ”, Efremov Research Institute, St. Petersburg, RF
Time=1s
2s
3s
4s
5s
6s
7s
8s
9s
10s
15s
20s
24s
67 68 69 70 71 72 73 74 75 76 77 78 79 80
Length, m
Evolution of the strand temperature along
the cable during the Quench (Current (40-
0KA) and field (13.4-0T) Dump time
constant: 15s)
24s
.
.
∂Tk
∂ ⎛ ∂Tk
⎞ 1 ∂ ⎛ ∂Tk
⎞
Ck
( T ) = qV
k ( x,
r,
t)
+ ⎜ κk
( T ) ⎟ + ⎜ κk
( T ) ⋅ r ⎟
∂t
∂x
⎝ ∂x
⎠ r ∂r
⎝ ∂r
⎠
Cable model treated by
∑ ki
i iVi
k
t x Ai
ρ
Γ
∂ρ ∂ρ
+ =
∂ ∂
V
∑ ki
iVi
fi
iVi
Vi
k
( Pi
iVi
)
t x
Dh
A
i
i
ρ
Γ
∂ρ ∂
2 − 2 ρ
+ + ρ =
+
∂ ∂
conv
ρH
∑ Qmi
+ ∑ Γki
⎛
2
Vi
P ⎞ ⎛
2
∂
i ∂
V ⎞ i m
k
ρ ⎜
i H
⎟
i
ρiV
⎜
i H ⎟
i =
t ⎜
+ −
⎟
+
i x ⎜
+
∂
∂
⎟
⎝ 2 ρ ⎠ ⎝ 2 ⎠ Ai
1 1 2 2 ∂Tm 2 ∂ ⎛
1 1 2 2
( A C + A C ) = − cos θ ⎜ ( A k + A k )
m m
m m
T = T
∂Tk
κ k + hi
⋅
∂n
k
He ( T − T ) = 0
k
∂t
m
∂Tk
κ k = ϕ(
x,
t)
∂n
i
∂x⎝
+ Q
Joule
m
m m
+
∑
i
Q
m m
conv
im
∂Tm
∂x
⎞
⎟ +
⎠
+ Q + Q
∑
n
cond
nm
Ext
m
Stability Analysis
CICC stability diagram for the 1 st cabling stage

Appl. Math. Dept. of STC “SINTEZ”, Efremov Research Institute, St. Petersburg, RF
VINCENTA - Numerical Tools
for Thermohydraulic Analysis of
Cryogenic Systems
• OVERVIEW
The code VINCENTA is designed to enhance the simulation of transient thermohydraulic processes in cryogenic environment of superconducting magnet systems. An advanced,
powerful calculation algorithm enables combined 1D, 2D and quasi-3D thermal calculations for a full range of operating modes, including normal operation, quench, ramping,
cooldown and emergency. A rich set of easy-to-link mathematical models provide a maximum realistic approximations for system geometry, material properties, coolant behavior and
accessory arrangement (piping, valves, pumps, etc).
• OPERATIONAL FEATURES
• comprehensive full-scale mathematical simulation of both the whole system and its components;
• forecasting simulation to help solving constructive problems from many points of view;
• realistic modeling for a variety of cryogenic systems and operation conditions;
• estimation of temperature distribution, coolant parameters, heat exchange, and nonlinear effects;
• high adaptability for each particular application.
• APPLICATIONS
A multipurpose versatility, extensive calculation capability and high performance of the VINCENTA code allows a wide range of most demanding applications, including:
• Fusion and particle accelerator magnets
• MRI-magnets
• Experimental and diagnostic devices for scientific research
• Superconducting generators
• Superconducting cables and joints

Appl. Math. Dept. of STC “SINTEZ”, Efremov Research Institute, St. Petersburg, RF
COND - Numerical Tools for
Thermohydraulic Analysis of SC
Magnet Systems
• OVERVIEW
The code COND is intended to simulate the long-term transient thermohydraulic processes such as cool-down, warm-up etc. in large SC Magnet Systems. An advanced, powerful
calculation algorithm enables 3D thermal calculations for a full range of operating modes, providing a maximum realistic approximations for system geometry, material properties.
• OPERATIONAL FEATURES
• comprehensive full-scale mathematical simulation of both the whole system and its components;
• forecasting simulation to help solving constructive problems from many points of view;
• realistic modeling for a variety of operation conditions;
• estimation of temperature distribution, coolant parameters, heat exchange, and nonlinear effects;
• high adaptability for each particular application.
• APPLICATIONS
A multipurpose versatility, extensive calculation capability and high performance of the COND
code allows a wide range of most demanding applications, including:
• Fusion and particle accelerator magnets
• MRI-magnets
• Experimental devices for scientific research
• Superconducting generators

GLORY - 3D Transient
Thermohydraulic Analysis of SC
Magnet Systems
• OVERVIEW
The code GLORY is a versatile, integrated software for extensive analysis of thermohydraulic processes in superconducting magnet systems. An effective
combination of problem-oriented packages enables numerical simulation of both stationary and transient thermohydraulic processes. The code provides a 3D
dynamic analysis, with estimates of thermal field, heat load, heat transfer, and thermal inertia for cooled magnet structures (cables, coils, supports, etc). Flexible,
maximum realistic models allow a comprehensive forecasting simulation to choose the best design and materials for each particular application.
• HIGHLIGHTS
Appl. Math. Dept. of STC “SINTEZ”, Efremov Research Institute, St. Petersburg, RF
• high-volume data processing;
• full-scale real-time mathematical simulation of both the entire system and separate components;
• enhanced modeling capability;
• easy to adapt for exact needs.
The code provides an analysis of transient behavior of SC magnet systems for a
wide range of geometry's, materials, cooling strategies and operational
conditions. A 1,000÷10,000,000 node finite-element mesh may be used for
calculation depending on the accuracy required.

Appl. Math. Dept. of STC “SINTEZ”, Efremov Research Institute, St. Petersburg, RF
KOMPOT-T. KOMPOT T. 3-D 3 D Thermostatic
Field Computation
• OVERVIEW
The KOMPOT code is the well-proven integrated software designed for numerical simulation and analysis of 3D thermostatic field. An efficient calculation
algorithm enables a versatile thermal field analysis using medium-scale computers. The numerical simulation provides a desired accuracy with the allowance for
complex system geometry and non-linear effects..
• HIGHLIGHTS
• an efficient numerical simulation algorithm capable of precise thermostatic field analysis;
• pre- and post processing of input/output data;
• prolonged intensive wide-range applications.
PERFORMANCE
The numerical simulation algorithm is based on the classical Poisson equation, finite-element
method, and symmetric successive overrelaxation method combined with a polynomial
acceleration of a convergence rate. An effective integration with another special program
packages using the same finite-element mesh makes it possible to obtain the heat releases
distribution from the steady and eddy currents, induced by electrical and transient magnetic
fields. Special calculating procedures allow to perform a combined thermohydraulic analysis
together with the codes VINCENTA and COND, that allows take into account the heat
exchange with a wide range of coolants (liquid and supercritical helium, water, etc) flowing in
cooling channels; take into account the boil-of of coolants on the given surfaces
.
REQUIREMENTS
•1000 ÷ 10,000,000 node mesh is possible for field
analysis;
•mid-runtime full-scale precise calculations are provided on
i586 PC or better