3 research outputs found

    Analysis of Conductance Probes for Two-Phase Flow and Holdup Applications

    Full text link
    [EN] In this paper we perform an analysis of the conductance probes used in two-phase flow applications especially for two-phase flow tomography of annular flow, to measure the waves produced in the interface with different boundary conditions without perturbing the flow, and in addition we examine the holdup applications as measuring the average void fraction in a given region. The method used to obtain the detector conductance between the electrodes is to solve analytically the generalized Laplace equation in 3D with the boundary conditions of the problem, and then to obtain the average potential difference between the detector electrodes. Then, dividing the current intensity circulating between the emitter and the receiver electrodes by the average potential difference yields the probe conductance, which depends on the geometric and physical characteristics of the measured system and the probe. This conductance is then non-dimensionalized by dividing by the conductance of the pipe full of water. In this way a set of analytical expression have been obtained for the conductance of two-plate sensors with different geometries and locations. We have performed an exhaustive comparison of the results obtained using the equations deduced in this paper with the experimental data from several authors in different cases with very good agreement. In some cases when the distribution of bubbles is not homogeneous, we have explored the different alternatives of the effective medium theory (EMT) in terms of the self-consistent EMT and the non-consistent EMT.This research was funded by Spain Ministry of Science and Technology previously dependent on MINECO (Ministry of Economy and Competitiveness), grant number ENE2016-79489-C2-1-P under Plan Nacional de I+D.Muรฑoz-Cobo Gonzรกlez, JL.; Rivera-Durรกn, Y.; Berna Escriche, C.; Escrivรก Castells, FA. (2020). Analysis of Conductance Probes for Two-Phase Flow and Holdup Applications. Sensors. 20(24):1-29. https://doi.org/10.3390/s20247042S129202

    Feasibility Test of a Liquid Film Thickness Sensor on a Flexible Printed Circuit Board Using a Three-Electrode Conductance Method

    No full text
    Liquid film thickness measurements under temperature-varying conditions in a two-phase flow are of great importance to refining our understanding of two-phase flows. In order to overcome the limitations of the conventional electrical means of measuring the thickness of a liquid film, this study proposes a three-electrode conductance method, with the device fabricated on a flexible printed circuit board (FPCB). The three-electrode conductance method offers the advantage of applicability under conditions with varying temperatures in principle, while the FPCB has the advantage of usability on curved surfaces and in relatively high-temperature conditions in comparison with sensors based on a printed circuit board (PCB). Two types of prototype sensors were fabricated on an FPCB and the feasibility of both was confirmed in a calibration test conducted at different temperatures. With the calibrated sensor, liquid film thickness measurements were conducted via a falling liquid film flow experiment, and the working performance was tested

    ๊ณ ์ •๋ฐ€ ์‹คํ—˜ ๋ฐ ํ•ด์„์— ๊ธฐ๋ฐ˜ํ•œ ์›์ž๋กœ ๊ฐ•์ˆ˜๋ถ€ ์ƒ๋ถ€์—์„œ์˜ ์•ก๋ง‰ ๊ฒฌ์ธ ๋ชจ๋ธ๋ง

    Get PDF
    ํ•™์œ„๋…ผ๋ฌธ(๋ฐ•์‚ฌ) -- ์„œ์šธ๋Œ€ํ•™๊ต๋Œ€ํ•™์› : ๊ณต๊ณผ๋Œ€ํ•™ ์›์žํ•ต๊ณตํ•™๊ณผ, 2021.8. ์กฐํ˜•๊ทœ.๋ณธ ์—ฐ๊ตฌ์˜ ๋ชฉ์ ์€ ๊ณ ์ •๋ฐ€ ์‹คํ—˜๊ณผ CFD ํ•ด์„์œผ๋กœ๋ถ€ํ„ฐ ์–ป์€ ๊ตญ์†Œ ์œ ๋™ ๋ณ€์ˆ˜๋“ค์„ ๊ธฐ๋ฐ˜์œผ๋กœ ์›์ž๋กœ ๊ฐ•์ˆ˜๋ถ€ ์ƒ๋ถ€์—์„œ ๋ฐœ์ƒํ•  ์ˆ˜ ์žˆ๋Š” ์•ก๋ง‰ ๊ฒฌ์ธ ํ˜„์ƒ์„ ๋ชจ๋ธ๋งํ•˜๋Š” ๊ฒƒ์ด๋‹ค. ๋จผ์ € ์•ˆ์ „์ฃผ์ž…์ˆ˜ ์šฐํšŒ ํ˜„์ƒ์—์„œ์˜ ์•ก๋ง‰ ๊ฑฐ๋™์— ๋Œ€ํ•œ ์ดํ•ด๋ฅผ ์ œ๊ณ ํ•˜๊ณ , CFD ํ•ด์„ ๊ฒ€์ฆ์— ์ ํ•ฉํ•œ ๊ณ ์ •๋ฐ€ ์‹คํ—˜ ๋ฐ์ดํ„ฐ๋ฅผ ํ™•๋ณดํ•˜๊ธฐ ์œ„ํ•ด 1/10 ์Šค์ผ€์ผ๋กœ ์ถ•์†Œ๋œ ๊ฐ•์ˆ˜๋ถ€์—์„œ ์‹คํ—˜ ์—ฐ๊ตฌ๊ฐ€ ์ˆ˜ํ–‰๋˜์—ˆ๋‹ค. ์‹คํ—˜์—์„œ์˜ ์ฃผ์š” ์ธก์ • ๋ณ€์ˆ˜๋Š” ์•ก๋ง‰ ๋‘๊ป˜์˜€์œผ๋ฉฐ, ์•ก๋ง‰ ๋‘๊ป˜ ๋ถ„ํฌ๋ฅผ ํšจ์œจ์ ์œผ๋กœ ์ธก์ •ํ•  ์ˆ˜ ์žˆ๋Š” ์ „๊ธฐ์  ๊ธฐ๋ฒ•์˜ ์•ก๋ง‰ ์„ผ์„œ๊ฐ€ ๋ณธ ์—ฐ๊ตฌ์—์„œ ๊ฐœ๋ฐœ๋˜์—ˆ๋‹ค. ๊ฐœ๋ฐœ ์„ผ์„œ๋Š” ํด๋ฆฌ์ด๋ฏธ๋“œ(polyimide) ํ•„๋ฆ„ ๊ธฐ๋ฐ˜์˜ ์—ฐ์„ฑํšŒ๋กœ๊ธฐํŒ์„ ์ด์šฉํ•˜์—ฌ ์ œ์ž‘๋˜์—ˆ๊ธฐ ๋•Œ๋ฌธ์— ๊ฐ•์ˆ˜๋ถ€์™€ ๊ฐ™์€ ๊ณก๋ฉด๋ถ€์—๋„ ์‰ฝ๊ฒŒ ๋ถ€์ฐฉ๋  ์ˆ˜ ์žˆ์—ˆ๋‹ค. ๊ณต๊ธฐ ์œ ๋™์— ๋”ฐ๋ผ์„œ ๋ณ€ํ•˜๋Š” ์•ก๋ง‰ ๊ฑฐ๋™์— ๋Œ€ํ•ด์„œ ์‹œ๊ฐ„ ํ‰๊ท  ์•ก๋ง‰ ๋‘๊ป˜๋ฅผ ์„ผ์„œ๋กœ ์ธก์ •ํ•˜๊ณ  ์ด๋ฅผ ๊ด€์ธก ๊ฒฐ๊ณผ์™€ ๋น„๊ตํ•˜์˜€๋‹ค. ๊ณต๊ธฐ ์œ ์†์ด ์ฆ๊ฐ€ํ• ์ˆ˜๋ก ์•ก๋ง‰์ด ํผ์ง€๋Š” ํญ์€ ํŒŒ๋‹จ๋ถ€ ์ค‘์‹ฌ์œผ๋กœ ์ข์•„์กŒ์œผ๋ฉฐ, ๋‘๊บผ์šด ์•ก๋ง‰ ๊ฒฝ๊ณ„์—์„œ ์•ก์  ์ดํƒˆ์ด ๋ฐœ์ƒํ•˜๋Š” ๊ฒƒ์ด ๊ด€์ธก๋˜์—ˆ๋‹ค. ๋˜ํ•œ ๊ณต๊ธฐ ์œ ์†์ด ํฐ ๊ฒฝ์šฐ์— ํŒŒ๋‹จ๋ถ€ ์ฃผ์œ„์—์„œ์˜ ๋‘๊บผ์šด ์•ก๋ง‰์œผ๋กœ๋ถ€ํ„ฐ wisp๊ฐ€ ์ƒ์„ฑ๋˜์–ด ํŒŒ๋‹จ๋ถ€๋กœ ๊ฒฌ์ธ๋˜๋Š” ๊ฒƒ์ด ํ™•์ธ๋˜์—ˆ๋‹ค. ์ด์™€ ๊ฐ™์€ ๊ด€์ธก์„ ํ†ตํ•ด ์•ก๋ง‰ ์„ผ์„œ ์ธก์ •์— ๋Œ€ํ•œ ์–‘์ , ์งˆ์  ํ‰๊ฐ€๋ฅผ ์ˆ˜ํ–‰ํ•  ์ˆ˜ ์žˆ์—ˆ๋‹ค. CFD ํ•ด์„์€ 3์ฐจ์› ํ˜•์ƒ ํšจ๊ณผ๊ฐ€ ์ˆ˜๋ฐ˜๋˜๋Š” ์ด์ƒ์œ ๋™์˜ ๊ฑฐ๋™์„ ์ •ํ™•ํžˆ ์˜ˆ์ธกํ•˜๊ณ ์ž ํ•  ๋•Œ ์œ ์šฉํ•˜๊ฒŒ ํ™œ์šฉ๋  ์ˆ˜ ์žˆ๋‹ค. ๊ทธ๋Ÿฌ๋‚˜ ์ „ํ†ต์ ์ธ ์ด์ƒ์œ ๋™ ๋ชจ๋ธ์€ ๋ถ„์‚ฐ์ƒ๊ณผ ์—ฐ์†์ƒ์ด ํ˜ผ์žฌํ•ด ์žˆ๋Š” ์ด์ƒ์œ ๋™์˜ ๊ฑฐ๋™์„ ์ ์ ˆํ•˜๊ฒŒ ๋‹ค๋ฃฐ ์ˆ˜ ์—†๊ธฐ ๋•Œ๋ฌธ์— CFD ํ•ด์„์˜ ํ™œ์šฉ์— ์ œํ•œ์ด ์žˆ์—ˆ๋‹ค. ๋ณธ ์—ฐ๊ตฌ์—์„œ๋Š” VOF ๋ชจ๋ธ๊ณผ Mixture ๋ชจ๋ธ์ด ๊ฒฐํ•ฉ๋œ ํ˜•ํƒœ์˜ VOF-slip ๋ชจ๋ธ์„ ํ™œ์šฉํ•˜์—ฌ ๊ฐ•์ˆ˜๋ถ€์—์„œ์˜ ์•ก๋ง‰๊ฒฌ์ธ ํ˜„์ƒ์„ ์ƒ์šฉ CFD ์ฝ”๋“œ์ธ STAR-CCM+๋กœ ํ•ด์„ํ•˜์˜€๋‹ค. VOF-slip ๊ธฐ๋ฒ•์—์„œ๋Š” ์•ก๋ง‰๊ณผ wisp์™€ ๊ฐ™์€ large-scale interface๋Š” VOF๋กœ ํ•ด์„ํ•˜๊ณ , ์•ก์ ๊ณผ ๊ฐ™์€ subgrid-scale interface๋Š” Mixture ๋ชจ๋ธ๋กœ ํ•ด์„ํ•˜๋Š” ๊ฒƒ์ด ๊ฐ€๋Šฅํ•˜๋‹ค. ๋ณธ ์—ฐ๊ตฌ์—์„œ๋Š” VOF-slip ๋ชจ๋ธ์„ ํ™œ์šฉํ•  ๋•Œ, ์•ก์  ์ง๊ฒฝ๊ฐ’๊ณผ ๊ณ„๋ฉด๋‚œ๋ฅ˜๋Œํ•‘ ๊ณ„์ˆ˜๊ฐ€ ํ•ด์„ ๊ฒฐ๊ณผ์— ์ฃผ์š”ํ•œ ์˜ํ–ฅ์„ ์ฃผ๋Š” ํŒŒ๋ผ๋ฏธํ„ฐ์ž„์„ ํ™•์ธํ•  ์ˆ˜ ์žˆ์—ˆ๋‹ค. ๋ฏผ๊ฐ๋„ ๋ถ„์„์„ ํ†ตํ•ด ํ•ฉ๋ฆฌ์ ์ธ ์•ก๋ง‰ ํญ ๋ฐ ์šฐํšŒ์œจ์„ ๋„์ถœํ•  ์ˆ˜ ์žˆ๋Š” ์•ก์  ์ง๊ฒฝ๊ณผ ๊ณ„๋ฉด๋‚œ๋ฅ˜๋Œํ•‘ ๊ณ„์ˆ˜๋ฅผ ๊ฒฐ์ •ํ•˜์˜€๋‹ค. ๊ณ ์ •๋ฐ€ ์‹คํ—˜ ๋ฐ ํ•ด์„์—ฐ๊ตฌ๋กœ๋ถ€ํ„ฐ ์•ก๋ง‰ ๊ฒฌ์ธ ํ˜„์ƒ์€ ๊ธฐ์ƒ ๋ฐ ์•ก์ƒ์˜ ์œ ๋Ÿ‰๊ณผ ์•ก๋ง‰ ๊ฒฝ๊ณ„์˜ ์œ„์น˜์— ์˜ํ•ด ๊ฒฐ์ •๋œ๋‹ค๋Š” ๊ฒƒ์„ ํ™•์ธํ•  ์ˆ˜ ์žˆ์—ˆ๋‹ค. ๋”ฐ๋ผ์„œ ์ด๋“ค์˜ ํŠน์„ฑ์„ ๋Œ€๋ณ€ํ•  ์ˆ˜ ์žˆ๋Š” ํŒŒ๋ผ๋ฏธํ„ฐ๋“ค์„ Re_la, We, ๊ทธ๋ฆฌ๊ณ  R_fb์™€ ๊ฐ™์ด ์ •์˜ํ•˜๊ณ  ์ •๊ทœํ™”๋œ ์•ก๋ง‰๊ฒฌ์ธ๋ฅ ์„ ์˜ˆ์ธกํ•  ์ˆ˜ ์žˆ๋Š” ์ƒ๊ด€์‹์„ ๊ฐœ๋ฐœํ–ˆ๋‹ค. ๊ฐœ๋ฐœ ๋ชจ๋ธ์—์„œ๋Š” ํŒŒ๋‹จ ์ €์˜จ๊ด€์„ ์ค‘์‹ฌ์œผ๋กœ off-take ๋ณผ๋ฅจ์ด ๋‘ ์„น์…˜(LEFT, RIGHT)์œผ๋กœ ๊ตฌ๋ถ„๋œ๋‹ค. ๊ทธ๋ฆฌ๊ณ  ๊ฐ ์„น์…˜์—์„œ ์˜ˆ์ธก๋œ ์šฐํšŒ์œจ์„ ํ•ฉ์‚ฐํ•˜์—ฌ ์ „์ฒด ์šฐํšŒ์œจ์ด ๋„์ถœ๋œ๋‹ค. ๊ฐœ๋ฐœ๋œ ์•ก๋ง‰๊ฒฌ์ธ๋ฅ  ๋ชจ๋ธ์€ MARS-multiD ์ฝ”๋“œ์— ์‚ฝ์ž…๋˜์—ˆ์œผ๋ฉฐ, SNU ์‹คํ—˜๊ณผ KAERI์—์„œ ์ˆ˜ํ–‰๋œ DIVA ์‹คํ—˜ (1/5 ์Šค์ผ€์ผ) ํ•ด์„์„ ํ†ตํ•ด ๊ฒ€์ฆ๋˜์—ˆ๋‹ค. ๊ฒ€์ฆ ๊ฒฐ๊ณผ, ๊ฐœ๋ฐœ ๋ชจ๋ธ์ด ๊ธฐ์กด MARS ํ•ด์„๋ณด๋‹ค ์ „์ฒด์ ์œผ๋กœ ๊ฐœ์„ ๋œ ์šฐํšŒ์œจ ์˜ˆ์ธก ์„ฑ๋Šฅ์„ ๋ณด์˜€๋‹ค. ์ถ”๊ฐ€์ ์œผ๋กœ ๊ฐœ๋ฐœ ๋ชจ๋ธ์˜ ํ™œ์šฉ์ด ์ฆ๊ธฐ-๋ฌผ ์กฐ๊ฑด์œผ๋กœ ํ™•์žฅ๋  ์ˆ˜ ์žˆ๋Š”์ง€ ํ™•์ธํ•ด๋ณด๊ธฐ ์œ„ํ•ด MIDAS ์‹คํ—˜์— ๋Œ€ํ•œ ํ•ด์„์„ ์ˆ˜ํ–‰ํ•˜์˜€๋‹ค. ์ฆ๊ธฐ-๋ฌผ ์กฐ๊ฑด์—์„œ๋Š” ๋ชจ๋ธ์ด ๊ฐœ์„ ๋  ์—ฌ์ง€๊ฐ€ ์žˆ์Œ์ด ํ™•์ธ๋˜์—ˆ์œผ๋ฉฐ, ๋ชจ๋ธ ๊ฐœ์„ ์„ ์œ„ํ•ด์„œ๋Š” ์‘์ถ• ํ˜„์ƒ์ด ์ˆ˜๋ฐ˜๋œ ์‹คํ—˜ ์—ฐ๊ตฌ๊ฐ€ ์ถ”ํ›„์— ์ˆ˜ํ–‰๋˜์–ด์•ผ ํ•  ๊ฒƒ์ด๋‹ค. ๋ณธ ์—ฐ๊ตฌ์—์„œ ์ˆ˜ํ–‰๋œ ์›์ž๋กœ ๊ฐ•์ˆ˜๋ถ€์—์„œ์˜ ์•ก๋ง‰ ๊ฒฌ์ธ ํ˜„์ƒ์— ๋Œ€ํ•œ ๊ณ ์ •๋ฐ€ ์‹คํ—˜ ๋ฐ ํ•ด์„, ๊ทธ๋ฆฌ๊ณ  ์ƒ๊ด€์‹ ๊ฐœ๋ฐœ ๋ฐ MARS-multiD ์ฝ”๋“œ ์ ์šฉ์€ ํ–ฅํ›„ ์›์ž๋กœ ๋ƒ‰๊ฐ์žฌ์ƒ์‹ค์‚ฌ๊ณ  ์กฐ๊ฑด ์‹œ, ์•ˆ์ „์ฃผ์ž…์ˆ˜ ์šฐํšŒ๋Ÿ‰ ์˜ˆ์ธก์— ์œ ์šฉํ•˜๊ฒŒ ์‚ฌ์šฉ๋  ์ˆ˜ ์žˆ๋‹ค. ํŠนํžˆ, ๊ฐœ๋ฐœ ๋ชจ๋ธ์€ ์žฌ๊ด€์ˆ˜ ๊ธฐ๊ฐ„ ๋™์•ˆ ์šฐํšŒ์œจ์— ๋”ฐ๋ฅธ ๋…ธ์‹ฌ์˜ ์žฌ๊ฐ€์—ด์„ ์˜ˆ์ธกํ•˜๋Š”๋ฐ ์ค‘์š”ํ•œ ์—ญํ• ์„ ํ•  ๊ฒƒ์œผ๋กœ ๊ธฐ๋Œ€๋œ๋‹ค.This study focuses on the modelling of the film off-take phenomenon in a reactor vessel downcomer based on local flow parameters obtained from experiment and computational fluid dynamics (CFD) analysis. Experiments are conducted in the reduced-scale downcomer annulus of a nuclear reactor pressure vessel to investigate the liquid film behaviors under emergency core coolant (ECC) bypass conditions and to obtain high-fidelity data for the validation of two-phase flow CFD codes. The main instrumentation is an electrical conductance sensor for measuring the local liquid film thickness, which is developed in this study. The fabrication of the electrodes on a flexible printed circuit board enabled the installation of the sensor on the curved surface. The developed sensor is used to measure the time-averaged liquid film thickness, which shows the influence of the lateral air flow on the liquid film flow, and the results are compared with visual observations. As the air velocity increased, a droplet that was created in the thick part of the liquid film appeared, and the wisps generated near the broken cold leg could be observed. In the experiment, qualitative and quantitative analyses of the measurement results showed the reliability of the developed sensor, and helped to understand the liquid film behavior in the ECC bypass phenomenon. Furthermore, the measured film thickness could contribute to film off-take modelling and to validating the CFD codes, which have not been validated sufficiently because of the absence of local measurement data. Recent advances in computational power have resulted in the application of CFD to nuclear reactor safety analyses, which require accurate predictability for two-phase flow with three-dimensional (3D) geometrical effects. Even though the different flow regimes can exist simultaneously in the real flow, the traditional two-phase CFD models have a disadvantage with respect to regime dependency. Therefore, the CFD study used VOF-slip, which is a hybrid model combining volume of fluid (VOF) and mixture model offered by STAR-CCM+ 15.04 was used. This approach enables the large-scale interface to be treated using the VOF method and the subgrid-scale interface to be treated with a mixture model that accounts for a phase slip via the drag law. The key parameters of the VOF-slip model for the film off-take phenomenon were the droplet diameter and the interface turbulence damping coefficient. Therefore, the sensitivity analyses are conducted by varying droplet diameter and damping coefficient and a suitable value was determined based on the film spreading width and ECC bypass fraction. The droplet diameter was determined to be 150 ฮผm for all simulation cases. The interface turbulence damping coefficients ranged from 0 to 30 and mesh-independent damping term ranged from 2.7ร—10-5m to 5.7ร—10-5m. From experiment and CFD analysis studies, it was confirmed that the liquid film off-take phenomenon is governed by the air flow rate, water flow rate, and the film boundary position. Considering these three parameters, the normalized film off take rate was correlated with Rela, We, Rfb, and Bo. The concept of the model was to divide the off-take volume into two sections (LEFT and RIGHT) through a virtual boundary so that the model could evaluate the film off-take rate in each section differently. The developed film off-take model was implemented in MARS-multiD, and it was validated with the SNU experiment (1/10 scale) and DIVA test (1/5 scale). The validation results showed that the newly developed film off-take model could improve the predictability of the bypass fraction. In addition, the MIDAS test with steam-water flow was simulated using the developed model, and the results was confirmed that the phenomenon accompanied by condensation should be experimentally investigated in future study to accurately predict the film off-take in the steam-water condition.Chapter 1 Introduction 1 1.1 Background and Motivation 1 1.1.1 Liquid film off-take in reactor vessel downcomer 1 1.1.2 Challenges with multi-dimensional system code simulation 3 1.2 Literature Review 5 1.2.1 ECC bypass experiment 5 1.2.2 Liquid film thickness sensor 7 1.2.3 CFD analysis 9 1.2.4 Modelling 11 1.3 Objectives and Scopes 12 Chapter 2 Liquid Film Thickness Sensor 18 2.1 Sensor Design 18 2.1.1 Features with flush-mounted electrode 18 2.1.2 Electrode design 19 2.1.3 Circuitry design 22 2.2 Sensor Calibration 23 2.2.1 Calibration method 23 2.2.2 Calibration result 24 Chapter 3 Experiment for Two-phase Film Flow 40 3.1 Scaling for ECC bypass phenomenon 40 3.2 Experimental Setup and Conditions 42 3.2.1 Experiment facility 42 3.2.2 Test matrix 43 3.3 Experimental Results 45 3.3.1 Time-averaged film thickness 45 3.3.2 Fluctuation of film thickness 51 3.3.3 ECC bypass fraction 53 Chapter 4 CFD Analysis 71 4.1 Two-phase CFD Models 71 4.1.1 VOF model 71 4.1.2 Mixture model 73 4.1.3 Two-fluid model 74 4.2 CFD Modelling 75 4.2.1 VOF-slip model 75 4.2.2 Interface turbulence damping 77 4.2.3 Computational domain and simulation cases 79 4.3 Simulation Results 79 4.3.1 No air flow conditions 79 4.3.2 Determination of droplet diameter 80 4.3.3 Effect of interface turbulence damping 83 Chapter 5 Modelling of Film Off-take 104 5.1 Difficulties Associated with Simulating Film Off-take Phenomenon 104 5.2 Development of Film Off-take Model 106 5.2.1 Strategy for model development 106 5.2.2 Definition of modelling parameters 108 5.2.3 Development of film off-take model 114 5.3 Validation of Developed Film Off-take Model 115 5.3.1 SNU experiment 115 5.3.2 DIVA experiment 116 5.3.3 MIDAS experiment 119 5.4 Applicability of Developed Film Off-take Model 121 Chapter 6 Conclusions 147 6.1 Summary 147 6.2 Recommendations 149 Nomenclature 150 References 154 Appendix A Uncertainty Analysis 162 Appendix B Implementation of Model in MARS 166 ๊ตญ๋ฌธ ์ดˆ๋ก 169๋ฐ•
    corecore