THOR(Tsing Hua Open-pool Reactor)為一裝載小型爐心的大型水池式反應器，由於其內部仍含有燃料，故以大量的液態水覆蓋，帶走其所釋放之大量熱能。由於台灣位處環太平洋地震帶，時常發生地震等天災事故，恐易導致池中池水潑灑出，使池水靜止水位高度下降，造成燃料裸露情況。因此，本研究決定採用計算流體力學(CFD)技術模擬水池在地震力下的池水晃蕩情況，並以數值方法中的流體體積法Volume of fluid (VOF)追蹤波動的水面和晃蕩所引起的碎波，藉以檢測出燃料池在未加有額外檔板和附加檔板後的條件下所損失的水量、地震過後的水位高度、地震中水位的波高、壁面壓力變化和每一秒的最大壓力。
由於計算時間過長，本研究先採用結構較為簡單的沸水式用過燃料池(SFP)作為初步模擬的物理模型，藉以檢測出此方法的可行性，待地震模式之方法論成立後便進階下一步，比較未附加檔板的THOR原模型及另兩項修改方案的晃蕩後結果，其中，此兩項方案分別為加有垂直檔板的模型A及加有水平檔版的模型B。
模擬結果顯示壁面壓力與地震加速度呈正相關性，並非受到自由液面垂直高度的影響。另外，垂直檔板的增加能減少水量的損失，但也會增加水波的最大高度；水平檔板的增加則能有效阻擋水量流失且降低平均波高，但於特殊情況時亦有增加最大波高的風險。

THOR(Tsing Hua Open-pool Reactor) is a large Open-pool Reactor with a small core inside. Since there are still fuels inside THOR, it needs a large amount of water to cover the fuels so as to bring out a great deal of heat. Taiwan situates in Circum-Pacific Seismic Belt, so earthquakes happen in Taiwan very often. The accidents may cause water spill from the pool and hence the exposure of fuels, which results in the damage of nuclear safety. Hence, this study provides a CFD methodology to simulate the liquid sloshing in the pool, and uses Volume of fluid (VOF) to track the free surface of the liquid and small water waves. The methodology can receive the results of the amount of water loss, water level after earthquakes, wave height, and maximum wall pressure per second.
This study took the spent fuel pool model to be the first simulation in order to test the feasibility of the methodology because the structure of SFP is simpler than THOR. After the methodology was verified, this study use it to simulate the liquid sloshing of THOR and other two revised models, Model A and Model B. Model A is the THOR model with vertical baffles and Model B is the one with horizontal baffles.
Finally, the results validate that the vertical baffles and horizontal baffles could reduce the water loss, although they may probably induce higher water wave. Other than that, the results also indicate that wall pressure can be strongly affected by earthquake acceleration.

摘要 I
Abstract II
致謝 III
表目錄 VII
圖目錄 VIII
符號對照表 XI
第一章 緒論 1
1.1 研究動機和目的 1
1.2 文獻回顧 3
第二章 理論基礎及計算模型 7
2.1 THOR的物理模型 7
2.1.1 未加檔板模型 7
2.1.2 增加二公尺垂直檔板之模型 10
2.1.3 增加水平檔板之模型 11
2.2 數學模型 13
2.3 LES模式 (Turbulence Model) 15
2.4 流體體積法(VOF) 16
2.5 地震力之加速度 17
2.6 案例設定與數值方法 18
2.6.1 案例設定 18
2.6.2 數值方法 19
2.7 研究流程 19
第三章 初步模擬 23
3.1 SFP的物理模型 23
3.2 SFP的模擬結果 24
3.2.1 晃蕩後的損失水量和靜止水位 24
3.2.2 用過燃料池的自由液面 25
3.2.3用過燃料池的最高水位、振幅和波高 30
3.2.4 用過燃料池的壁面壓力結果 33
第四章 THOR的結果與討論 37
4.1 晃蕩後的損失水量和靜止水位 37
4.2 THOR模型的自由液面 39
4.2.1 原模型的自由液面 39
4.2.2 模型A-10.53 m的自由液面 45
4.2.3 模型B-8.83m的自由液面 51
4.3 THOR模型的最高水位、振幅和波高 57
4.3.1 最高水位評比 57
4.3.2 振幅和波高 58
4.4 THOR模型的的壁面壓力結果 60
4.5 THOR模型的的速度結果 64
第五章 結論 65
參考文獻 67

國. 原科中心, 清華大學水池式反應器(THOR), 國立清華大學 原科中心, 2017，取自:http://thor.web.nthu.edu.tw/files/132-1028-19,r4-1.php?Lang=zh-tw.
國. 原科中心, 清華大學水池式反應器(THOR), 原科中心, 國立清華大學, 2017，取 自:http://thor.web.nthu.edu.tw/files/11-1028-1583.php?Lang=zh-tw.
D. M. Liu., “A numerical study of three-dimensional liquid sloshing in tanks. Journal of Computational Physics,” Journal of Computational Physics 227.8, pp. 3921-3939, 2008.
O. M. Faltinsen, “A numerical nonlinear method of sloshing in tanks with two-dimensional flow,” Journal of Ship Research 22.3, 1978.
O. M. e. a. Faltinsen, “Multidimensional modal analysis of nonlinear sloshing in a rectangular tank with finite water depth,” Journal of Fluid Mechanics 407 , pp. 201-234, 2000.
Faltinsen, Odd M., and Alexander N. Timokha, “An adaptive multimodal approach to nonlinear sloshing in a rectangular tank,” Journal of Fluid Mechanics 432, pp. 167-200, 2001.
V. e. a. Singal, “CFD analysis of a kerosene fuel tank to reduce liquid sloshing,” Procedia Engineering 69, pp. 1365-1371, 2014.
Kandasamy, T., S. Rakheja, and A. K. W. Ahmed, “An analysis of baffles designs for limiting fluid slosh in partly filled tank trucks,” Open Transportation Journal 4, pp. 23-32, 2010.
Hou, Ling, Fangcheng Li, and Chunliang Wu, “A numerical study of liquid sloshing in a two-dimensional tank under external excitations,” Journal of Marine Science and Application 11.3, pp. 305-310, 2012.
Singal, Vaibhav, et al., “CFD analysis of a kerosene fuel tank to reduce liquid sloshing.,” Procedia Engineering 69 , pp. 1365-1371, 2014.
Craig, K. J., and T. C. Kingsley., “Design optimization of containers for sloshing and impact.,” Structural and Multidisciplinary Optimization 33.1, pp. 71-87, 2007.
國立清華大學原子科學技術發展中心, “國立清華大學水池式反應器運轉執照更新安全分析報告,” 2012.
王福军, 计算流体动力学分析: CFD 软件原理与应用, 清华大学出版社有限公司, 2004.
J. Smagorinsky, “General Circulation Experiments with the Primitive Equations,” Monthly Weather Review. 91 (3), pp. 99-164, March 1963.
J. Deardorff, “A numerical study of three-dimensional turbulent channel flow at large Reynolds numbers,” Journal of Fluid Mechanics. 41 (2), pp. 453-480, 1970.
C. Y. e. a. Wu, “Investigating the advantages and disadvantages of realistic approach and porous approach for closely packed pebbles in CFD simulation.,” Nuclear Engineering and design 240.5, pp. 1151-1159, 2010.
H. a. C. N. S. Ayhan, “CFD modeling of thermal mixing in a T-junction geometry using LES model.,” Nuclear Engineering and Design 253, pp. 183-191, 2012.
D. e. a. Gueyffier, “Volume-of-fluid interface tracking with smoothed surface stress methods for three-dimensional flows.,” Journal of Computational physics 152.2, pp. 423-456, 1999.
行. 核能管制處, “「核能電廠耐震安全再評估精進作業」之地震危害度分析與設計地震檢討原能會安全評估報告,” 行政院原子能委員會 核能管制處 , 中華民國 103年8月.
S. V. a. D. B. S. Patankar, “A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows,” International journal of heat and mass transfer 15.10, pp. 1787-1806, 1972.