TRACE Model Development

WORKING GROUP MEETING SPRING 2016TASK 2 BREAKOUT SESSION

ANN ARBOR, MI

Emory Brown

• Task Description

• Timeline

• Pump

• Simple Models

• Revised Geometry and Test Train Model

• Current Challenges and Questions

Outline

June 6, 2016

1

Task 2.2 Overview

June 6, 2016

2

Flow of information

June 6, 2016

3

Commercial PWR

TREAT Reactor

TWERL Loop

Neutronic boundary conditions

Commercial PWR

Transient CHF data

TRTL Loop

Task 2.2.11: Modeling of benchmark test with TRACE

• OSU will develop a TRACE model for one of the benchmark tests performed using the U.S. NRC code TRACE. Modeling of the benchmark test will be done blindly, based on the design package put together as a part of task 2.2.3. The data will not be made available until the modeling and results have been completed.

June 6, 2016

4

Task 2.2.11 specific objectives

• Start with TWERL model.

• Known geometries.

• TRTL model will be nearly identical

• Compare model with RELAP5-3D

• Finalize TRTL model

• Model benchmark simulation and compare to experimental benchmark data (2.2.13)

June 6, 2016

5

TRACE Model Development Timeline

June 6, 2016

6

NovOct Dec FebJan Mar Apr

Gain access to TRACE, SNAP, and AptPlot

Simple models to familiarize with code.

Gain access to TWERL Solidworks model

Gradually increase complexity and identify issues and gaps in knowledge

Found Theory and Model Guidelines manuals

Focus on test section and compare to RELAP-5

June 6, 2016

7

1.5 10 · .

Pump Impeller Inertia *

118 ·

.

Pump Motor Inertia *

* Wylie et al.

Pump Information

· · ·

Flow Rate Head Efficiency Power Normalized Curve

GPM m3/s ft m m2/s2 % KW(in) KW(heat) q/ω h/ω2

0.24 1.48E‐05 59.64 18.18 178.32 0.1% 2.649 2.646 2.599E‐05 9.112E‐06

2.88 1.81E‐04 59.27 18.07 177.23 2.0% 1.610 1.578 3.179E‐04 9.056E‐06

5.88 3.71E‐04 58.78 17.92 175.75 6.0% 1.089 1.024 6.506E‐04 8.981E‐06

9.18 5.79E‐04 58.04 17.69 173.53 12.0% 0.839 0.738 1.015E‐03 8.867E‐06

13.68 8.63E‐04 56.46 17.21 168.82 18.0% 0.811 0.665 1.513E‐03 8.627E‐06

17.61 1.11E‐03 54.69 16.67 163.53 22.0% 0.827 0.645 1.947E‐03 8.357E‐06

20.86 1.32E‐03 52.73 16.07 157.66 24.0% 0.866 0.658 2.307E‐03 8.057E‐06

RATED 25.03 1.58E‐03 50.14 15.28 149.93 26.0% 0.912 0.675 2.768E‐03 7.661E‐0629.24 1.84E‐03 47.81 14.57 142.96 27.5% 0.961 0.697 3.233E‐03 7.305E‐06

31.65 2.00E‐03 45.99 14.02 137.50 27.8% 0.990 0.714 3.500E‐03 7.026E‐06

33.72 2.13E‐03 43.96 13.40 131.44 27.8% 1.008 0.728 3.729E‐03 6.717E‐06

37.10 2.34E‐03 41.11 12.53 122.93 27.5% 1.048 0.760 4.103E‐03 6.282E‐06

40.36 2.55E‐03 38.64 11.78 115.53 26.7% 1.104 0.809 4.463E‐03 5.904E‐06

42.10 2.66E‐03 37.13 11.32 111.02 26.1% 1.132 0.836 4.655E‐03 5.673E‐06

44.03 2.78E‐03 35.17 10.72 105.17 25.5% 1.148 0.855 4.868E‐03 5.374E‐06

46.78 2.95E‐03 32.12 9.79 96.03 24.5% 1.159 0.875 5.172E‐03 4.907E‐06

50.07 3.16E‐03 27.87 8.49 83.33 23.5% 1.122 0.858 5.536E‐03 4.258E‐06

gc 9.81m/s2 QR 1.58E‐03 m3/sS.G. 1 [unitless] HR 149.93 m2/s2

RPM3450RPM 5.790 in‐lbf ΩR 361.3 rad/sec361 rad/sec 0.482 ft‐lbf HPR, liquid 0.316945556 hp

Impeller Size 3.5 in 0.654 kg‐m2/s2 τR 0.654173293 Pa‐m3

IP 1.04E‐03 kg‐m2

IM 5.98E‐04 kg‐m2

Effective M.O.I 1.6350E‐03 kg‐m2

Pump Information

June 6, 2016

8

Pump Information

June 6, 2016

9

0%

5%

10%

15%

20%

25%

30%

0 10 20 30 40 50 60

Head

[ft]

Flow Rate [GPM]

Pump Curve ‐ 3.5" Impeller @ 67 Hz

Head

Efficiency

Simple loops…

• Learning the ropes of TRACE, SNAP, and APTPLOT. • Creating incrementally complex heated water loops to verify

solutions to analytical thermal hydraulic problems.

June 6, 2016

10

t[s]=m kg

Q kJs

CpkJ

kg K∆T K

4.619782.31173 4.18 12

99.83

Creating a TRACE model of TWERL• Extract physical dimensions from

Solidworks model. • Pump information provided from

vendor pump curve. • Translate into TRACE model• Run initial case similar to simplified

model. Q=(m/t)(Cp)(∆T)

June 6, 2016

11

4/21: Model revamp with lessons learned

June 6, 2016

12

AA

0.076200 m

0.340030 m

0.021819 m

0.027051 m

0.027864 m

0.024689 m

0.027864 m 0.026848 m

0.057150 m 0.005715 m

SECTION A-ASCALE 1 : 2

AA

0.136624 m

0.223510 m

0.079055 m

C

A

B

SECTION A-ASCALE 1 : 2

0.003200 m

DETAIL CSCALE 2 : 1

Top of 3rd Section

0.003200 m

DETAIL ASCALE 2 : 1

Bottom of 1st Section

0.003200 m DETAIL B

SCALE 2 : 1

Bottom of fueled section.Top is the top of the last pellet

0.052502 m

0.039051 m

0.026645 m

0.025425 m

0.024549 m

0.116442 m

Test Section Geometry

June 6, 2016

13

Comp. Cell Length (m) dZ (m) In Edge Flow Area (m2) HD (m) Out Edge Flow Area 

(m2) HD (m) Changing Diameter?

Annular? Dinner (m) Douter

(m) VolumeVolume 

average flow area

110

1 0.0762 0.0762 1 0.00057472 0.027051 2 0.00057472 0.027051 0 4.37937E‐05 0.0005747202 0.34003 0.34003 2 0.000609785 0.027864 3 0.000609785 0.027864 0 0.000207345 0.0006097853 0.021819 0.021819 3 0.000478737 0.024689 4 0.000478737 0.024689 0 1.04456E‐05 0.0004787374 0.026848 0.026848 4 0.000609785 0.027864 5 0.000609785 0.027864 0 1.63715E‐05 0.0006097855 0.05715 0.05715 5 0.000248287 0.01778 6 0.000216617 0.005715 1 X 0.00635 0.01778 0.000247683 0.004333917

120 1 0.329828 0.329828 1 0.000127286 0.003188 2 0.000127286 0.003188 0 X 0.009526 0.0159 4.19824E‐05 0.000127286

1211 0.136624 0.136624 1 0.000127674 0.0032 2 0.000127674 0.0032 0 X 0.00950 0.01590 1.74434E‐05 0.0001276742 0.22351 0.22351 2 0.000127674 0.0032 3 0.000127674 0.0032 0 X 0.00950 0.01590 2.85365E‐05 0.0001276743 0.079055 0.079055 3 0.000127674 0.0032 4 0.000127674 0.0032 0 X 0.00950 0.01590 1.00933E‐05 0.000127674

122

1 0.441325 0.441325 1 0.000137688 0.003393 2 0.000137688 0.003393 0 X 0.009524 0.01631 6.07651E‐05 0.000137688

2 0.022225 0.022225 2 0.000071601 0.001981 3 0.000546175 0.008784 1 X * See attached figure 6.20899E‐06 0.000279370

3 0.069056 0.069056 3 0.000387948 0.022225 4 0.000387948 0.022225 0 2.67901E‐05 0.000387948123 1 0.572294 0.572294 1 0.000507705 0.025425 2 0.000507705 0.025425 0 0.000290557 0.000507705

130

1 0.0257285 0.0257285 1 0.000507705 0.025425 2 0.002164919 0.052502 1 3.19121E‐05 0.0012403412 0.026645 0.026645 2 0.002164919 0.052502 3 0.002164919 0.052502 0 5.76843E‐05 0.0021649193 0.1031195 0.1031195 3 0.002164919 0.052502 4 0.002164919 0.052502 0 0.000223245 0.0021649191 0.024549 1 0.000557598 0.026645 2 0.000557598 0.026645 0 1.36885E‐05 0.000557598

131 1 0.50165 0.50165 1 0.002123227 0.051994 2 0.002123227 0.051994 0.001065117 0.002123227

140

1 0.023749 1 0.000557598 0.026645 2 0.000557598 0.026645 0 1.32424E‐05 0.0005575982 0.0598473 2 0.000557598 0.026645 3 0.000557598 0.026645 0 3.33707E‐05 0.0005575983 0.060541 0.060541 3 0.000557598 0.026645 4 0.000557598 0.026645 0 3.37576E‐05 0.0005575984 0.020231 0.020231 4 0.000478737 0.024689 5 0.000478737 0.024689 0 9.68533E‐06 0.0004787375 0.103308 0.103308 5 0.000609785 0.027864 6 0.000609785 0.027864 0 6.29957E‐05 0.0006097856 0.021819 0.021819 6 0.000478737 0.024689 7 0.000478737 0.024689 0 1.04456E‐05 0.000478737

1501 1.692453 1.692453 1 0.000609785 0.027864 2 0.000609785 0.027864 0 0.001032033 0.0006097852 0.021819 0.021819 2 0.000478737 0.024689 3 0.000478737 0.024689 0 1.04456E‐05 0.0004787373 0.103124 0.103124 3 0.000609785 0.027864 4 0.000609785 0.027864 0 6.28835E‐05 0.000609785

160

1 0.021819 0.021819 1 0.000478729 0.0246888 2 0.000478729 0.0246888 0 1.04454E‐05 0.0004787292 0.144805 0.144805 2 0.000609785 0.027864 3 0.000609785 0.027864 0 8.82999E‐05 0.0006097853 0.042141 3 0.000609785 0.027864 4 0.000609785 0.027864 0 2.5697E‐05 0.0006097854 0.0762 0.0762 4 0.000609785 0.027864 5 0.000609785 0.027864 0 4.64656E‐05 0.0006097855 0.042141 5 0.000609785 0.027864 6 0.000609785 0.027864 0 2.5697E‐05 0.0006097856 0.0762 0.0762 6 0.000609785 0.027864 7 0.000609785 0.027864 0 4.64656E‐05 0.000609785

180 1 0.914385 0.914385 1 0.008107 0.101598 2 0.008107 0.101598 0 0.00741292 0.008107000

Test Section Parameters

June 6, 2016

14

Fill/Break Component Length (m) Volume (m3) Pressure Liquid Temp 

(K) Vapor Temp (K) TypeLiquid Velocity (m/s)

Vapor Velocity (m/s)

1 Fill 0.05715 1.84E‐06 1.5513E+07 550 550 Constant Velocity 0.9 0.911 Break 0.081867502 2.63E‐06 1.5513E+07 550 550 No tables 0.9 0.9

Heat Structure

Component Axial Cell Inner Surface BC Outer Surface BC

Radial Node 1 (K)

Radial Node 2 (K)

Radial Node 3 (K)

Radial Node 4 (K)

Material Region

Inner Length (m)

Outer Length (m) Material

21

1 130:1 Flux: 0.0 550 550 550 550 1 3.200E‐03 3.797E‐03 Zircaloy2 130:2 Flux: 0.0 550 550 550 550 2 3.797E‐03 3.848E‐03 Gap Gases3 130:3 Flux: 0.0 550 550 550 550 3 3.848E‐03 7.950E‐03 Mixed Oxide4 130:4 Flux: 0.0 550 550 550 5505 130:5 Flux: 0.0 550 550 550 550 Gas Gap HTC 6300 W/m2/K6 130:6 Flux: 0.0 550 550 550 5507 130:7 Flux: 0.0 550 550 550 5508 130:8 Flux: 0.0 550 550 550 5509 130:9 Flux: 0.0 550 550 550 550

10 130:10 Flux: 0.0 550 550 550 550

Power Con

nection

Component Axial Cell Relative Power Shape Radial Node Radial Locations (m)

Radial Power Shape

31

1 1 1 3.20E‐03 02 1 2 3.797E‐03 03 1 3 3.848E‐03 14 1 4 7.950E‐03 15 16 1 Power Option Constant Power7 1 Initial Power 5000 W8 19 1

10 1

Test Train Entrance Nozzle (110)

June 6, 2016

15

True Geometry Calculated Geometry As-Modeled Geometry

Exit Nozzle Detail (122)

June 6, 2016

True Geometry Calculated Geometry As-Modeled Geometry

Test Train Results

June 6, 2016

17

• Steady State Calculations

• Compare to RELAP-5 simulations to determine if physics is being modeled correctly

• Change parameters to verify phenomena is being modeled as expected

Test Train Results

June 6, 2016

18

0%

2%

4%

6%

8%

10%

12%

14%

570

575

580

585

590

595

600

0 5 10 15 20 25 30 35 40

Temperature

[K]

Cell

Temperature Void Frac Subcooled Boiling P = 15.5 MPaTsat = 620 K

Q = 5 kWMass flow = 0.0416 kg/s

Test Train Results

June 6, 2016

19

Subcooled Boiling P = 15.5 MPaTsat = 620 K

Q = 10 kWMass flow = 0.0208 kg/s

0%

10%

20%

30%

40%

50%

60%

570

575

580

585

590

595

600

605

610

615

620

625

0 5 10 15 20 25 30 35 40

Temperature

[K]

Cell

Temperature Void Frac

Model Specific Questions

June 6, 2016

20

Pipes Pump• Inner vs outer. Does it

matter?

• Appropriate power profile(axial and radial)

Heat Structure• Effective Moment of

Inertia verification

• 2-phase curves required?

• Dealing with Tee junction

• (Pipe with side junction)

• Liquid-gas interface in component 131

• Dealing with abrupt geometry changes. When to smooth out

High Level Challenges

June 6, 2016

21

TRACE SNAP• Is there a manual?

• How to do axial plots? (no spatial variable to plot against)

• MATLAB Interfacing worth it? (lots of preprocessing just to organize data)

APTPlot• Model Notes?

• Improving workflow between SNAP and TRACE

• Error message explanations and causes.

• Understanding numerical instabilities.

• Access to resources beyond that of manuals.

Thank you for your time.

Questions?

June 6, 2016

22