Channel Under Transversely Uniform and Non-uniform Heating Omar S. - - PowerPoint PPT Presentation

β–Ά
channel under transversely uniform and
SMART_READER_LITE
LIVE PREVIEW

Channel Under Transversely Uniform and Non-uniform Heating Omar S. - - PowerPoint PPT Presentation

18 th IGORR Conference, Sydney, Australia, Dec 3 rd -7 th , 2017. Onset of Flow Instability in a Rectangular Channel Under Transversely Uniform and Non-uniform Heating Omar S. Al-Yahia, Taewoo Kim, Daeseong Jo * *Corresponding author:


slide-1
SLIDE 1

18th IGORR Conference, Sydney, Australia, Dec 3rd -7th, 2017.

School of Mechanical Engineering, Kyungpook National University

Omar S. Al-Yahia, Taewoo Kim, Daeseong Jo*

80 Daehak-ro, Buk-gu, Daegu, 41566, Republic of Korea

Onset of Flow Instability in a Rectangular Channel Under Transversely Uniform and Non-uniform Heating

*Corresponding author: djo@knu.ac.kr

slide-2
SLIDE 2

18th IGORR Conference, Sydney, Australia, Dec 3rd -7th, 2017.

CONTENTS

Introduction Experimental setup Data reduction 1 Results and discussion Conclusion 2 3 4 5

slide-3
SLIDE 3

18th IGORR Conference, Sydney, Australia, Dec 3rd -7th, 2017.

  • 1. INTRODUCTION (1/2)

The minimum point on pressure drop-mass flux curve is referred to the OFI incipience

Axial distance Wall Temperature Wall Temperature Heat flux ONB ONB At fixed position Constant heat flux

1- Visualization 2- Wall temperature measurement ONB observation OFI observation

slide-4
SLIDE 4

18th IGORR Conference, Sydney, Australia, Dec 3rd -7th, 2017.

  • 1. INTRODUCTION (2/2)

Transverse heat flux distribution in the plate type fuel research reactor

  • ONB is local phenomena depending on the local heat flux and wall temperature.
  • OFI depends on the total thermal power deposited in the flow channel
  • In the plate type fuel research reactors, the power distribution is non-uniform along the axial direction

as well as the transverse direction Study Objective:

  • Investigate the effect of transverse

power distribution on the ONB and OFI incipience.

  • Compare the thermal hydraulic

behavior of ONB and OFI between uniform and non-uniform heat flux distribution.

*Jo, D., Seo, C.G., 2015. Effects of transverse power distribution on thermal hydraulic analysis. Progress in Nuclear Energy 81, 16-21.

slide-5
SLIDE 5

18th IGORR Conference, Sydney, Australia, Dec 3rd -7th, 2017.

Cross sectional view Front view with TC locations

Non-Uniform test section Uniform test section

Front view with TC locations

  • Two SUS316L cartridge

heaters

  • Aluminum block
  • 10 single thermocouples (TCs)

8 TCs distributed axially 2 TCs distributed transversely

  • Two SUS316L cartridge

heaters

  • Copper block
  • Air gap
  • Aluminum block
  • 6 single thermocouples (TCs)
  • 7 double thermocouples (TCs)

Cross sectional view

  • 2. EXPERIMENT SETUP (1/4)
slide-6
SLIDE 6

18th IGORR Conference, Sydney, Australia, Dec 3rd -7th, 2017.

  • 2. EXPERIMENT SETUP (2/4)

Schematic diagram for the experimental facility

Flowmeter Supply Pump Sink V3 V4 Main Pump

TCs

Ξ”P

T P T P

Preheater V1 V2 Condensing Tank Water Reservoir Heat Exchanger High Speed Camera Test Section Check Valve

T D

V5 Drain Line Pump Control Panel Personal Computer Data Acquisition System

D T P

Pressure Measurement Temperature Measurement Density Measurement Main Pipe Supply Pipe Power Cable Signal Cable

Power Control Panel Power Transformer

slide-7
SLIDE 7

18th IGORR Conference, Sydney, Australia, Dec 3rd -7th, 2017.

  • 2. EXPERIMENT SETUP (3/4)
slide-8
SLIDE 8

18th IGORR Conference, Sydney, Australia, Dec 3rd -7th, 2017.

  • 2. EXPERIMENT SETUP (4/4)

Degassing Adjust the mass flow rate Adjust the inlet temperature Apply the heat step wisely Record the data

Experimental procedure Test conditions

Parameter Value Flow rate [kg/s] 0.030-0.130 Heat flux [kW/m2] 100-800 Power distribution Uniform/Non-uniform Inlet temperature [C] 35-65 Pressure atm~ Hydraulic diameter [m] 0.004504

(1) Adjust the mass flow rate (2) Adjust the power

(1) Increase the power step wisely

(2) Decrease the mass flowrate step wisely

slide-9
SLIDE 9

18th IGORR Conference, Sydney, Australia, Dec 3rd -7th, 2017.

  • 3. Data reduction

𝑼𝒙 = 𝑼𝑼𝑫 βˆ’ 𝒓"π’Žπ’‘π’…π’– 𝒍

𝒓"π’Žπ’‘π’… = 𝒓"π’ƒπ’˜π’‰ = 𝑹𝒇 π‘©π’Š Γ— 𝟏. πŸ˜πŸ’

By comparing the applied electric power 𝑹𝒇 with the imposed thermal power π‘Ήπ’–π’Š, the energy losses is approximately 7 % and 10 % for uniform and non-uniform test section, respectively.

𝒓"π’Žπ’‘π’… = 𝒍 π’š βˆ†π‘Όπ‘¬

Non-Uniform test section Uniform test section 𝑼𝒙 = 𝑼𝑼𝑫 βˆ’ 𝒓"π’Žπ’‘π’…π’– 𝒍 βˆ†π‘Όπ‘¬

slide-10
SLIDE 10

18th IGORR Conference, Sydney, Australia, Dec 3rd -7th, 2017.

  • 4. Results and Discussion (1/7)

Transverse normalized heat flux distribution (at 3.9 kW)

(a) (b)

ONB incipience on the heated surface; (a) Non-uniformly heated surface, (b) Uniformly heated surface.

  • The local heat flux for the

uniform test section is similar at any location on the heated surface.

  • the local heat flux near the edges

is much higher than the middle part of the non-uniform heated section.

Γ—1.3 Γ—2.8

slide-11
SLIDE 11

18th IGORR Conference, Sydney, Australia, Dec 3rd -7th, 2017.

  • 4. Results and Discussion (2/7)

The ONB incipience for uniform and non-uniform heat flux, (Mass flow rate is 0.08 kg/s, Inlet temperature is 50 oC)

  • The ONB incipience is local

phenomenon that is highly depends on the local conditions such as the local heat flux rather than the total power deposited in the channel.

  • The ONB in the case of non-

uniformly heater occurs at power lower than the one for the case of uniformly heated surface due to high heat flux near the edges.

slide-12
SLIDE 12

18th IGORR Conference, Sydney, Australia, Dec 3rd -7th, 2017.

  • 4. Results and Discussion (3/7)

The ONB heat flux for the uniform and non-uniform test section (Mass flow rate is 0.08 kg/s, Inlet temperature is 50 oC).

  • The local heat flux at the ONB is

similar for the case of uniform and non-uniform heated test section, as well as the local wall temperature.

slide-13
SLIDE 13

18th IGORR Conference, Sydney, Australia, Dec 3rd -7th, 2017.

  • 4. Results and Discussion (4/7)

Thermal hydraulic parameters under constant power (Power 3.57 kW, inlet Temperature 50oC)

Non-Uniform test section Uniform test section

  • Pressure drop is different, the inlet pressure fluctuation conditions are same

OIPF OIPF

slide-14
SLIDE 14

18th IGORR Conference, Sydney, Australia, Dec 3rd -7th, 2017.

ሢ 𝑛 = 0.050 𝑙𝑕/𝑑 ሢ 𝑛 = 0.043 𝑙𝑕/𝑑 ሢ 𝑛 = 0.035 𝑙𝑕/𝑑 ሢ 𝑛 = 0.031 𝑙𝑕/𝑑 ሢ 𝑛 = 0.024 𝑙𝑕/𝑑 Center Edge t=0ms t=9 ms t=27ms t=18ms

  • In the case of non-uniform heating, the

pressure drop after OIPF is not increased due to low void fraction in the middle.

  • When the flow pattern changes to churn slug

flow near the edge, the pressure drop suddenly increases.

Comparison of void fraction under constant power (Power 3.00 kW, Inlet temperature 65 oC)

Before OIPF Constant dP Sudden increase of dP

Change of flow pattern under non-uniform heating

  • 4. Results and Discussion (5/7)
slide-15
SLIDE 15

18th IGORR Conference, Sydney, Australia, Dec 3rd -7th, 2017.

  • 4. Results and Discussion (6/7)

Thermal hydraulic parameters under constant mass flow rate (Mass flow rate 0.03 kg/s, inlet Temperature 50oC)

Non-Uniform test section Uniform test section

  • Pressure drop is different, the inlet pressure fluctuation conditions are same
slide-16
SLIDE 16

18th IGORR Conference, Sydney, Australia, Dec 3rd -7th, 2017.

Non-Uniform test section

  • 4. Results and Discussion (7/7)

Point (4) and (6) Point (5)

  • Different moments in bubbles

generation and condensation between edge and center lead to have different pressure fluctuation behavior

  • When the flow pattern changes

between edge and center, the pressure drop changes.

slide-17
SLIDE 17

18th IGORR Conference, Sydney, Australia, Dec 3rd -7th, 2017.

  • 5. Conclusion

(a) Effects of transversely heat flux distribution on the ONB and OFI are experimentally investigated through a narrow rectangular channel heated form one-side. (b) At the same total power, the local heat flux of the non-uniformly heated surface is much higher than the one in the uniform case. (c) ONB is local phenomena, it occurs at the same heat fluxes and wall temperature, even though the thermal power in the case of non-uniform heat flux is around 25 % less than the one in uniform case. (d) OFI is global phenomena. OFI occurs at similar thermal power and mass fluxes for the same operation conditions. (e) The differences in the heat flux distribution lead to different bubble behavior: the pressure drop behavior and void generation are different between uniform and non-uniform heat fluxes.

slide-18
SLIDE 18

18th IGORR Conference, Sydney, Australia, Dec 3rd -7th, 2017.

Thank you for your attention

slide-19
SLIDE 19

18th IGORR Conference, Sydney, Australia, Dec 3rd -7th, 2017.

References

1. Bergles, A.E., Rohsenow, W.M., β€œThe Determination of Forced Convection Surface Boiling Heat Transfer”, ASME J. Heat Transfer 86 (1964) 365–372. 2. Lee, J., Jo, D., Chae, H., Chang, S.H., Jeong, Y.H., Jeong, J.J., β€œThe Characteristics of Premature and Stable Critical Heat Flux for Downward Flow Boiling at Low Pressure in a Narrow Rectangular Channel”, Experimental Thermal and Fluid Science 69 (2015) 86-98. 3. Al-Yahia, O.S., Jo, D., β€œOnset of Nucleate Boiling for Subcooled Flow through a One-Side Heated Narrow Rectangular Channel”, Annals of Nuclear Energy 109 (2017) 30-40. 4. Al-Yahia, O.S., Jo, D., β€œONB, OSV, and OFI for Subcooled Flow Boiling Through a Narrow Rectangular Channel Heated on One-Side”, International Journal of Heat and Mass Transfer 116 (2018) 136-151. 5. Jo, D., Seo, C.G., β€œEffects of Transverse Power Distribution on Thermal Hydraulic Analysis”, Progress

  • f Nuclear Energy 81 (2015) 16-21.

6. Al-Yahia, O.S., Lee, Y.J., Jo, D., β€œEffect of Transverse Power Distribution on the ONB Location in the Subcooled Boiling Flow”, Annals of Nuclear Energy 100 (2017) 98-106. 7. Al-Yahia, O.S., Kim, T., Jo, D.,s, β€œExperimental Study Of Uniform And Non-Uniform Transverse Heat Flux Distribution Effect on The Onset of Nucleate Boiling”, Proceedings of the 25th International Conference on Nuclear Engineering ICONE25 May 14-18, 2017, Shanghai, China.