# Moisture Transfer Model and Simulation for Dehumidification of HTGR Core

Authors: Jun Li; Zhenzhong Zhang; Yingchao Meng; Huaqiang Yin; Shengchao Ma; Xuedong He; Xingtuan Yang; Shengyao Jiang
Journal: Science and Technology of Nuclear Installations
 (2018)

---

## Abstract

A large number of carbon materials are adopted in high-temperature gas-cooled reactor (HTGR). These carbon materials mainly include graphite IG-110 and boron-containing carbon material (BC), both of which are typical porous materials and normally absorb moisture. In order to inhibit the chemical corrosion reaction between core internals materials and moisture, the core needs to be strictly dehumidified before the reactor is put into operation. This paper mainly analyzed the moisture transfer mechanism in these carbon materials. Moisture transfer models were developed, and the dehumidification process of HTR-PM core was simulated. In addition, the influence of working temperature and system pressure on dehumidification was studied as well.


---

## Body

## 1. Introduction

Due to the excellent comprehensive performance of carbon materials, they have been extensively adopted in high-temperature gas-cooled reactor (HTGR), serving as structural material, fuel matrix material, etc. [1–3]. These carbon materials used in HTGR typically include isostatic pressure graphite IG-110 and boron-containing carbon material (BC) [4]. Graphite IG-110 mainly serves as fuel matrix material, neutron reflecting layer, and other reactor core internals, while for BC, since it contains boron, which has a big neutron absorption cross section, it is made into the shielding layer [5]. It has been reported that about 60 tons and 1000 tons of graphite have been applied in 10 MW high-temperature gas-cooled reactor (HTR-10) [6] and high-temperature gas-cooled reactor-pebblebed modules (HTR-PM) [7], respectively.As porous material, carbon materials have numerous pore structures and large specific surface area. The pore structures and molecule groups attached on the surface of pores result in strong moisture absorption capacity [8, 9]. Also, other impurities, such as H2O, O2, and CO2, can be absorbed by carbon materials as well. These impurities will react with the carbon materials at high temperature, thus, weakening the material strength and decreasing the service life of HTGR [10–13]. In addition, under some extreme nuclear accidents, e.g., air or water leakage accident, the corrosion can be even more serious [14, 15]. In order to attenuate the corrosion, the limited value of moisture content in the helium coolant of the primary circuit is set as 2 cm3/m3 [16]. Therefore, before power operation and test, the reactor internals and fuel balls need to be dehumidified to control the moisture content to a certain level.The dehumidification system [17] of HTR-PM is shown in Figure 1. Firstly, the primary circuit is evacuated and then filled with coolant helium at pressure of 7 MPa. The circulated helium is driven by helium circulator, whose mechanical energy is transformed into heat of helium. Thus, the temperature of helium can be heated to 250°C. Due to the high temperature, moisture absorbed in the carbon internals of reactor core will be released and then carried away by flowing helium. The helium flow containing relatively high moisture content will flow out of internals and then be cooled to about 20°C in two water-helium coolers, so the water vapor will condense to liquid water which can be separated in gas-water separator. Finally, dry helium will return to the reactor core. This circulation of helium flow described above will gradually decrease the moisture content in the carbon internals. After being dehumidified for several dozens of days, the moisture content of reactor core will meet the design requirement.Figure 1
Dehumidification system of HTR-PM.Dehumidification process mainly includes moisture transfer in carbon materials, moisture transfer in helium flow, and the condensation and separation of moisture. During dehumidification and purification process of reactor core, mass transfer in carbon material is considered as the dominating factor. Thus, mass transfer behavior in these carbon materials is worth investigating.In this paper, Section2 will present the analysis on moisture transfer mechanism in carbon materials, and then mathematical models will be developed in Section 3. In Section 4, the simulation on dehumidification of HTR-PM will be introduced, and simulation result will be analyzed.


## 2. Moisture Transfer Mechanism

### 2.1. Mass and Heat Transfer

Dehumidification process of porous material involves both moisture mass transfer and heat transfer, the relationship of which can be described by Lewis NumberLe:(1)Le=αDwhere α (m2/s) andD (m2/s) are heat diffusion coefficient and mass diffusion coefficient respectively.WhileLe> 60, it is thought that the heat transfer is faster much more than mass transfer, so that it can be neglected in the analysis [18]. For graphite, the heat diffusion coefficient is about 1.07 × 10−4 m2/s, and mass diffusion coefficient ranges 10−7 ~ 10−6 m2/s. Thus,Le of graphite ranges 102 ~ 103, which means that the heat transfer can be neglected, so it is thought that the mass transfer takes place at constant temperature in analysis.Considering the moisture in material is in gaseous state, the mass transfer mechanism of gas in porous material involves molecular diffusion and seepage flow.


### 2.2. Molecular Diffusion

Molecular diffusion is the result of molecular thermal motion and collision, the mass flux densityJ (kg·m−2·s−1) caused by which can be described by Fick’s law:(2)J=-De∇ρwhereρ (kg/m3) is moisture concentration andDe (m2/s) is effective diffusion coefficient.De is influenced by diffusion type, gas component, porosity of material and tortuosity of pore, etc.De can be expressed as(3)De=D0×εeτwhere D0 (m2/s) is the gas diffusion coefficient. εe is effective porosity of material, which is equal to open porosity in general.τ is tortuosity of pore, which is defined as the ratio of the actual length to the apparent length of pore. In addition, the tortuosity is in reciprocal relationship to the effective porosity; thus, relation τ=εe-1 has been adopted to the analysis of graphite [19, 20]. Table 1 summarizes the expressions of D0 according to diffusion type and gas component.Table 1
Expressions ofD0 according to diffusion type and gas component.
K n Diffusion type Gas component D 0 K n < 0.1 general diffusion [21] unitary gas D F = 1 3 λ 8 R T π M binary gas D A B = 4.36 × T 1.5 ( 1 / M A + 1 / M B ) 0.5 1 0 4 × P ( ν A 1 / 3 + ν B 1 / 3 ) 2  0.1 < K n < 10 Knudsen diffusion - D K = 1 3 d 8 R T π M  K n > 10 transition diffusion - D T = 1 1 / D F + 1 / D K or 11/DAB+1/DKIn Table1,P (Pa) is the gas pressure,T (K) is the temperature, andR is gas constant which is 8.314 J/(mol·K) in general. MA (g/mol) and MB (g/mol) are the molar mass of gas component A and B, respectively. νA (m3/kmol) and νB (m3/kmol) are the molar volume of component A and B respectively, which are in liquid state and at normal boiling point.Diffusion type is influenced by pore diameterd and molecular average free pathλ, which can be determined by Knudsen number Kn. Kn can be expressed as(4)Kn=λdwhered is also the most probable pore diameter of material. The value ofKn and corresponding diffusion type are presented in Table 1.λ can be expressed as [22](5)λ=3.2μPRT2πMwhere μ (Pa·s or kg/(m·s)) is the kinetic viscosity of gas,P (Pa) is the pressure, andM (kg/mol) is the molar mass. Relevant parameters of material IG-110 and BC, which are acquired by mercury intrusion porosimetry (MIP), are shown in Table 2.Table 2
Pore parameters of carbon materials.
Parameter IG-110 BC Volume density/ (g/cm3) 1.7435 1.6960 Total porosityεt/% 22.85 25.61 Open porosityεo/% 19.68 14.37 Closed porosityεc/% 3.17 11.24 Specific surface area/ (m2/g) 15.93 7.34 Most probable pore diameterd/nm 2512 30262


### 2.3. Seepage

Besides molecular diffusion, another possible transfer mechanism is convection, which can be defined as seepage in porous media. Seepage flow is a kind of viscous flow driven by pressure gradient. According to Darcy’s law, the seepage velocityu   (m/s) is(6)u=-Kμ∇PwhereK (m2) is the permeability of porous material, which is influenced by porosity, pore diameter, tortuosity, etc. According to Timur formula [23], K∝εe4.4. In addition, research has shown that, for graphite, whose open porosity is about 10%,K is around 10−16 m2 [24].


### 2.4. Comparison of Transfer Mechanisms at Different Pressure

Different drying methods involve different pressure conditions. For the dehumidification of HTR-PM, the pressure of its working atmospherePo (Pa) is 7 MPa, while for vacuum drying method and hot air drying method which are extensively adopted in experiment and industry,Po is vacuum and 0.1 MPa, respectively. The working pressure has a great influence on the moisture transfer mechanism in porous materials. A comparison is given in Table 3. In addition, the comparison is based on the following conditions: temperature is 250°C, the working atmosphere is helium, and the initial moisture content Cr,0 is 0.1%. Moisture content Cr is defined as(7)Cr=mH2OmC×100%where mC (kg) is the dry material mass and mH2O (kg) is the moisture mass. In addition, both of the density of IG-100 and BC are set to 1700 kg/m3 in analysis.Table 3
Comparison of transfer mechanisms at different pressure.
P o/MPa < 0.1 (vacuum) 0.1 (normal pressure) 7 (high pressure) Moisture transfer mechanism diffusion +seepage diffusion + seepage diffusion  Gas component H2O H2O + He H2O + He  Diffusion type general diffusion + Knudsen diffusion general diffusion general diffusion  Gas diffusion coefficientD0 variable variable constantThe comparison is discussed as follows. For hot air drying method, wherePo = 0.1 MPa, the gas mixture in pores of material is binary, which consists of H2O and He. The released moisture can generate a high pressurePi in pores, which can be 2.86 MPa in early stage as if all the absorbed moisture transforms to vapor. Then the pressure gradient betweenPi andPo drives the seepage flow in pores. Meanwhile, molecular diffusion also exists. According to (5), molecular free pathλ of H2O is about 114 nm (250°C, 0.1 MPa), while the most probable pore diametersd of IG-110 and BC are 2.5 μm and 30.3 μm, respectively, soKn < 0.1, which means that the general diffusion dominates in pores. As for vacuum drying method, wherePo < 0.1 MPa, only water vapor exists in pores and Knudsen diffusion also exists according toKn. As for the dehumidification of HTR-PM, wherePo = 7 MPa, seepage flow driven by the pressure gradient can be neglected andPi can be considered constant; thus the gas diffusion coefficientD0 influenced byPi is constant as well.


## 2.1. Mass and Heat Transfer

Dehumidification process of porous material involves both moisture mass transfer and heat transfer, the relationship of which can be described by Lewis NumberLe:(1)Le=αDwhere α (m2/s) andD (m2/s) are heat diffusion coefficient and mass diffusion coefficient respectively.WhileLe> 60, it is thought that the heat transfer is faster much more than mass transfer, so that it can be neglected in the analysis [18]. For graphite, the heat diffusion coefficient is about 1.07 × 10−4 m2/s, and mass diffusion coefficient ranges 10−7 ~ 10−6 m2/s. Thus,Le of graphite ranges 102 ~ 103, which means that the heat transfer can be neglected, so it is thought that the mass transfer takes place at constant temperature in analysis.Considering the moisture in material is in gaseous state, the mass transfer mechanism of gas in porous material involves molecular diffusion and seepage flow.


## 2.2. Molecular Diffusion

Molecular diffusion is the result of molecular thermal motion and collision, the mass flux densityJ (kg·m−2·s−1) caused by which can be described by Fick’s law:(2)J=-De∇ρwhereρ (kg/m3) is moisture concentration andDe (m2/s) is effective diffusion coefficient.De is influenced by diffusion type, gas component, porosity of material and tortuosity of pore, etc.De can be expressed as(3)De=D0×εeτwhere D0 (m2/s) is the gas diffusion coefficient. εe is effective porosity of material, which is equal to open porosity in general.τ is tortuosity of pore, which is defined as the ratio of the actual length to the apparent length of pore. In addition, the tortuosity is in reciprocal relationship to the effective porosity; thus, relation τ=εe-1 has been adopted to the analysis of graphite [19, 20]. Table 1 summarizes the expressions of D0 according to diffusion type and gas component.Table 1
Expressions ofD0 according to diffusion type and gas component.
K n Diffusion type Gas component D 0 K n < 0.1 general diffusion [21] unitary gas D F = 1 3 λ 8 R T π M binary gas D A B = 4.36 × T 1.5 ( 1 / M A + 1 / M B ) 0.5 1 0 4 × P ( ν A 1 / 3 + ν B 1 / 3 ) 2  0.1 < K n < 10 Knudsen diffusion - D K = 1 3 d 8 R T π M  K n > 10 transition diffusion - D T = 1 1 / D F + 1 / D K or 11/DAB+1/DKIn Table1,P (Pa) is the gas pressure,T (K) is the temperature, andR is gas constant which is 8.314 J/(mol·K) in general. MA (g/mol) and MB (g/mol) are the molar mass of gas component A and B, respectively. νA (m3/kmol) and νB (m3/kmol) are the molar volume of component A and B respectively, which are in liquid state and at normal boiling point.Diffusion type is influenced by pore diameterd and molecular average free pathλ, which can be determined by Knudsen number Kn. Kn can be expressed as(4)Kn=λdwhered is also the most probable pore diameter of material. The value ofKn and corresponding diffusion type are presented in Table 1.λ can be expressed as [22](5)λ=3.2μPRT2πMwhere μ (Pa·s or kg/(m·s)) is the kinetic viscosity of gas,P (Pa) is the pressure, andM (kg/mol) is the molar mass. Relevant parameters of material IG-110 and BC, which are acquired by mercury intrusion porosimetry (MIP), are shown in Table 2.Table 2
Pore parameters of carbon materials.
Parameter IG-110 BC Volume density/ (g/cm3) 1.7435 1.6960 Total porosityεt/% 22.85 25.61 Open porosityεo/% 19.68 14.37 Closed porosityεc/% 3.17 11.24 Specific surface area/ (m2/g) 15.93 7.34 Most probable pore diameterd/nm 2512 30262


## 2.3. Seepage

Besides molecular diffusion, another possible transfer mechanism is convection, which can be defined as seepage in porous media. Seepage flow is a kind of viscous flow driven by pressure gradient. According to Darcy’s law, the seepage velocityu   (m/s) is(6)u=-Kμ∇PwhereK (m2) is the permeability of porous material, which is influenced by porosity, pore diameter, tortuosity, etc. According to Timur formula [23], K∝εe4.4. In addition, research has shown that, for graphite, whose open porosity is about 10%,K is around 10−16 m2 [24].


## 2.4. Comparison of Transfer Mechanisms at Different Pressure

Different drying methods involve different pressure conditions. For the dehumidification of HTR-PM, the pressure of its working atmospherePo (Pa) is 7 MPa, while for vacuum drying method and hot air drying method which are extensively adopted in experiment and industry,Po is vacuum and 0.1 MPa, respectively. The working pressure has a great influence on the moisture transfer mechanism in porous materials. A comparison is given in Table 3. In addition, the comparison is based on the following conditions: temperature is 250°C, the working atmosphere is helium, and the initial moisture content Cr,0 is 0.1%. Moisture content Cr is defined as(7)Cr=mH2OmC×100%where mC (kg) is the dry material mass and mH2O (kg) is the moisture mass. In addition, both of the density of IG-100 and BC are set to 1700 kg/m3 in analysis.Table 3
Comparison of transfer mechanisms at different pressure.
P o/MPa < 0.1 (vacuum) 0.1 (normal pressure) 7 (high pressure) Moisture transfer mechanism diffusion +seepage diffusion + seepage diffusion  Gas component H2O H2O + He H2O + He  Diffusion type general diffusion + Knudsen diffusion general diffusion general diffusion  Gas diffusion coefficientD0 variable variable constantThe comparison is discussed as follows. For hot air drying method, wherePo = 0.1 MPa, the gas mixture in pores of material is binary, which consists of H2O and He. The released moisture can generate a high pressurePi in pores, which can be 2.86 MPa in early stage as if all the absorbed moisture transforms to vapor. Then the pressure gradient betweenPi andPo drives the seepage flow in pores. Meanwhile, molecular diffusion also exists. According to (5), molecular free pathλ of H2O is about 114 nm (250°C, 0.1 MPa), while the most probable pore diametersd of IG-110 and BC are 2.5 μm and 30.3 μm, respectively, soKn < 0.1, which means that the general diffusion dominates in pores. As for vacuum drying method, wherePo < 0.1 MPa, only water vapor exists in pores and Knudsen diffusion also exists according toKn. As for the dehumidification of HTR-PM, wherePo = 7 MPa, seepage flow driven by the pressure gradient can be neglected andPi can be considered constant; thus the gas diffusion coefficientD0 influenced byPi is constant as well.


## 3. Mathematical Model

Two mathematical models for moisture transfer, diffusion model and diffusion-seepage model, are developed and presented here. In addition, the gas components in the equations below include H2O and He.


### 3.1. Diffusion Model

Diffusion model is based on molecular diffusion and Fick’s law. The governing equation for moisture mass in porous material is(8)∂ρH2O∂τ=-∇J=De∇2ρH2Owhere ρH2O (kg/m3) is moisture concentration.


### 3.2. Diffusion-Seepage Model

Diffusion-seepage model is based on diffusion and seepage, and the corresponding governing equation of moisture is(9)∂ρH2O∂τ=∇KRT·ρH2OμHe+H2O·MHe+H2O·∇ρHe+H2O+De∇ρH2Owhere MHe+H2O (kg/mol), μHe+H2O (Pa·s or kg/(m·s)), and ρHe+H2O (kg/m3) are the molar mass, kinetic viscosity, and density of binary gas consists of H2O, and He, respectively.The kinetic viscosity of binary gas is [25](10)uHe+H2O=aH2OMH2O0.5uH2O+aHeMHe0.5uHeaH2OMH2O0.5+aHeMHe0.5where aH2O and aHe are the volume percentages of H2O and He, respectively.For the binary gas at 0.1 MPa, due to the relatively small proportion of He in pore, it has little effect on moisture transfer. Therefore, to simplify the governing Equation (9), ρHe can be considered as a constant and then the simplified equation is(11)∂ρH2O∂τ=∇KRT·ρH2OμHe+H2O·MHe+H2O+De∇ρH2O


## 3.1. Diffusion Model

Diffusion model is based on molecular diffusion and Fick’s law. The governing equation for moisture mass in porous material is(8)∂ρH2O∂τ=-∇J=De∇2ρH2Owhere ρH2O (kg/m3) is moisture concentration.


## 3.2. Diffusion-Seepage Model

Diffusion-seepage model is based on diffusion and seepage, and the corresponding governing equation of moisture is(9)∂ρH2O∂τ=∇KRT·ρH2OμHe+H2O·MHe+H2O·∇ρHe+H2O+De∇ρH2Owhere MHe+H2O (kg/mol), μHe+H2O (Pa·s or kg/(m·s)), and ρHe+H2O (kg/m3) are the molar mass, kinetic viscosity, and density of binary gas consists of H2O, and He, respectively.The kinetic viscosity of binary gas is [25](10)uHe+H2O=aH2OMH2O0.5uH2O+aHeMHe0.5uHeaH2OMH2O0.5+aHeMHe0.5where aH2O and aHe are the volume percentages of H2O and He, respectively.For the binary gas at 0.1 MPa, due to the relatively small proportion of He in pore, it has little effect on moisture transfer. Therefore, to simplify the governing Equation (9), ρHe can be considered as a constant and then the simplified equation is(11)∂ρH2O∂τ=∇KRT·ρH2OμHe+H2O·MHe+H2O+De∇ρH2O


## 4. Simulation and Result Analysis

### 4.1. Initial Parameters and Boundary Conditions

The dehumidification process of HTR-PM core has been simulated by using COMSOL Multiphysics software. Table3 and (8) demonstrate the corresponding physical hypothesis and mathematical model, respectively.The geometric model of HTR-PM core is shown in Figure2. The internals in HTR-PM core mainly include the reflecting layer and insulation layer. Reflecting layer, mainly made of graphite bricks (IG-110), can be approximately considered as a cylinder, the dimension of which is 15600 mm in height and 4500 mm in outside diameter. In addition, the reflecting layer creates an interior space to contain fuel balls, the dimension of which is 11800 mm in equivalent height and 3000 mm in inside diameter. Insulation layer, mainly made of carbon bricks (BC), is also an approximate cylinder located outside the reflecting layer, the dimension of which is 4500 mm in inside diameter and 5000 mm in outside diameter. In addition, the thickness of top and bottom insulation layers is 400 mm and 800 mm, respectively. In fact, the horizontal cross section is not a perfect round but a regular polygon with 30 edges. Moreover, diameter of fuel ball is 60 mm. During the dehumidification process, the helium flow (250°C, 7 MPa) takes the moisture away from internals. Helium flows up from core bottom to cold gas chamber through cold gas channel and then flows down to hot gas chamber through core chamber and leaves the core at last.Figure 2
Internals in HTR-PM core.
(a)
Part of internals (b)
A module of internalsSome simplifications of geometric model are adopted here. The neutron source channels, absorption ball channels, hot gas duct, dowels and keys, gap between graphite and carbon bricks, and so on are neglected. Thus, the internals model can be divided into 30 same trigonal prism modules, one of which is shown in Figure2(b). Each module is mirror symmetrical to the adjacent one. Therefore, we can adopt only one module to represent the moisture transfer process in the whole core internals. In addition, pore structure of materials features isotropy, and the material parameters of fuel ball are considered the same as graphite IG-110.The mass transfer process is considered as internal diffusion controlled, so the first kind of boundary condition is adopted to the model, which means moisture contentCr,b on the boundary, surface of which is exposed to helium atmosphere, is constant. Then Cr,b=0 is adopted in this simulation. In addition, the kinetic viscosity u of gas is considered constant due to the very small effect of pressure variation.Some key parameters adopted in simulation are presented in Table4.Table 4
Some parameters of model (250°C, 7 MPa).
Parameter Value Unit Density of heliumρHe 6.524 kg/m3 Molar mass of He 0.004 kg/mol Molar mass of H2O 0.018 kg/mol Kinetic viscosity of heliumuHe 2.925×10−5 kg/(m·s) Kinetic viscosity of water vaporuH2O 1.822×10−5 kg/(m·s) Most probable aperture of BC 30.3 μm Most probable aperture of IG-110 2.5 μm Open porosity of BCεo 14.37 % Open porosity of IG-110εo 19.68 % Gas diffusion coefficientDAB of binary gas (H2O + He) 1.226×10−6 m2/s Effective diffusion coefficientDe D e =D AB × ε o 2 m2/s PermeabilityK of BC 0.49×10−15 m2 PermeabilityK of IG-110 1.97×10−15 m2


### 4.2. Result Analysis of Simulation

According to the numerical simulation, the distribution of moisture contentCr (%) in HTR-PM core at timet = 50 d is shown in Figure 3, where d represents units day. It can be found from Figure 3(a) that moisture mainly concentrates at the top, bottom, and side of core. Figure 3(b) presents the moisture distribution at the axial cross section. It can be found that the maximum Cr locates at the top of core. In addition, Figure 3(c) and Figure 3(d) further demonstrate the moisture distribution at the top and bottom surface, respectively. It can be seen that Cr near the cold gas and control rod channels, fuel feed tube, and fuel discharge tube is relatively low, while it is much higher elsewhere.Figure 3
Distribution ofCr (%) in HTR-PM core (t= 50 d).
(a)
Horizontal cross sections (b)
Axial cross section (c)
Top surface (d)
Bottom surfaceThe curves of moisture content versus time are shown in Figure4. It can be found that the average moisture content Cr,a of internals and fuel ball decreases exponentially, in which it can be explained that moisture content decreases much faster in the early stage because of the high moisture concentration gradient and decreases slower in the later stage as a result of the decrease of concentration gradient. Cr,a of fuel ball decreases far faster than that of core internals. Thus, the internals, not fuel balls, are the major limiting factor for the dehumidification. It can also be found that the required time to reach Cr,a = 0.02%, 0.01%, and 0.005% of whole internals is 43 d, 72 d, and 108 d, respectively. In addition, Cr,a of reflecting layer (IG-110) decreases much faster than that of insulation layer (BC), which can be attributed to the larger open porosity εo of IG-110. Thus effective diffusion coefficient De of IG-110 is much larger, which promotes the moisture transfer.Figure 4
Curves of moisture versus time.


### 4.3. Sensitivity Analysis

In order to accelerate the moisture transfer to improve the dehumidification efficiency, it is reasonable to promote the molecular diffusion and seepage flow. According to the relations,De∝T1.5 and De∝P-1, molecular diffusion can be promoted by increasing the temperatureT and decreasing the system pressurePo, while seepage can be promoted by increasing the pressure gradient, which can be achieved by decreasing the system pressure Po.The influence of temperatureT on dehumidification of HTR-PM is shown in Figure 5(a), whereCr,a on Y-axis denotes the average moisture content of whole internals obtained by simulation. It can be found that the dehumidification efficiency increases as temperature increases. Figure 5(b) presents the influence of system pressure Po on dehumidification, and it shows that the dehumidification efficiency of 0.1 MPa is certainly much higher than that of 7 MPa. In addition, in the simulation where Po=0.1 MPa, the diffusion-seepage model is adopted and De is considered constant to simplify the computation.Figure 5
Curves of moisture versus time at different temperature and pressure.
(a)
Temperature (b)
PressureAs for the actual dehumidification for HTR-PM,T = 250°C andPo = 7 MPa. Since the temperature is already pretty high, the promotion by increasing temperature is limited, while as for the system pressure, if it is reduced greatly from 7 MPa to 0.1 MPa, this change will promote both the diffusion and seepage significantly. But on the other hand, the decrease of pressure is not advantageous for the core to be heated. It is because the heat of core is transformed from mechanical energy of helium circulator; thus the heating power is relatively low at low pressure.


### 4.4. Analytical Solution

In order to obtain an analytical solution of the mathematical model, the irregular geometric model of core internals can be represented by an infinite plate, whose characteristic lengthLc (m) is the same as that of irregular geometry. In addition, Lc=V/Ae, whereV (m3) and Ae (m2) are the volume and effective surface area of geometry, respectively.The infinite plate with2Lc thickness is shown in Figure 6. Diffusion model (8) is adopted. The moisture content is Cr(x,t), and effective diffusion coefficient De is constant. Thus, the governing equation can be rewritten as(12)∂Cr∂t=De∂2Cr∂x2The boundary condition is Cr(x=Lc,t)=Cr,b. The initial condition is Cr(x,t=0)=Cr,0. Thus, the analytical solution, which is an infinite series, can be obtained as(13)Cr=Cr,b+∑n=1∞an·e-bnt·coscnxwhere cn=(n-0.5)π/Lc (n≥1, and n is an integer), bn=De·cn2, and an=Cr,0-Cr,b.Figure 6
Infinite plate.In order to simplify the solution, only the first polynomial of (13), which means n=1, is adopted. Thus, the average moisture content Cr,a can be calculated, which is(14)Cr,a=∫0LcCrdxLc=Cr,b+2Cr,0-Cr,bπ·e-Deπ2/4Lc2tAs for the geometric model of internals of HTR-PM, according to the statistic obtained in Solidworks software, its Lc is 0.3066 m. In addition, considering Cr,b=0 and Cr,0=0.1%, the analytical solution for dehumidification of HTR-PM can be obtained.Figure7 shows the comparison between analytical and numerical simulation. It is found that curvesCr,a-t of analytical and numerical solution of infinite plate are in good agreement. While comparing curves of analytical and numerical solution of HTR-PM, it is found that the agreement is not that good. It is suggested that the simplified analytical solution is pretty accurate for an ideal regular geometry model, while it is difficult to obtain the characteristic lengthLc to accurately represent an irregular geometry. The uneven distribution of channels, chambers, and other irregular structures in reactor core results in the deviation of solution. It is also suggested that, in order to derive a more representativeLc, some optimization methods can be adopted to the calculation ofLc, e.g., introducing empirical correction coefficient, dividing the whole geometry into several sections, etc.Figure 7
Comparison of analytical and numerical solutions.


## 4.1. Initial Parameters and Boundary Conditions

The dehumidification process of HTR-PM core has been simulated by using COMSOL Multiphysics software. Table3 and (8) demonstrate the corresponding physical hypothesis and mathematical model, respectively.The geometric model of HTR-PM core is shown in Figure2. The internals in HTR-PM core mainly include the reflecting layer and insulation layer. Reflecting layer, mainly made of graphite bricks (IG-110), can be approximately considered as a cylinder, the dimension of which is 15600 mm in height and 4500 mm in outside diameter. In addition, the reflecting layer creates an interior space to contain fuel balls, the dimension of which is 11800 mm in equivalent height and 3000 mm in inside diameter. Insulation layer, mainly made of carbon bricks (BC), is also an approximate cylinder located outside the reflecting layer, the dimension of which is 4500 mm in inside diameter and 5000 mm in outside diameter. In addition, the thickness of top and bottom insulation layers is 400 mm and 800 mm, respectively. In fact, the horizontal cross section is not a perfect round but a regular polygon with 30 edges. Moreover, diameter of fuel ball is 60 mm. During the dehumidification process, the helium flow (250°C, 7 MPa) takes the moisture away from internals. Helium flows up from core bottom to cold gas chamber through cold gas channel and then flows down to hot gas chamber through core chamber and leaves the core at last.Figure 2
Internals in HTR-PM core.
(a)
Part of internals (b)
A module of internalsSome simplifications of geometric model are adopted here. The neutron source channels, absorption ball channels, hot gas duct, dowels and keys, gap between graphite and carbon bricks, and so on are neglected. Thus, the internals model can be divided into 30 same trigonal prism modules, one of which is shown in Figure2(b). Each module is mirror symmetrical to the adjacent one. Therefore, we can adopt only one module to represent the moisture transfer process in the whole core internals. In addition, pore structure of materials features isotropy, and the material parameters of fuel ball are considered the same as graphite IG-110.The mass transfer process is considered as internal diffusion controlled, so the first kind of boundary condition is adopted to the model, which means moisture contentCr,b on the boundary, surface of which is exposed to helium atmosphere, is constant. Then Cr,b=0 is adopted in this simulation. In addition, the kinetic viscosity u of gas is considered constant due to the very small effect of pressure variation.Some key parameters adopted in simulation are presented in Table4.Table 4
Some parameters of model (250°C, 7 MPa).
Parameter Value Unit Density of heliumρHe 6.524 kg/m3 Molar mass of He 0.004 kg/mol Molar mass of H2O 0.018 kg/mol Kinetic viscosity of heliumuHe 2.925×10−5 kg/(m·s) Kinetic viscosity of water vaporuH2O 1.822×10−5 kg/(m·s) Most probable aperture of BC 30.3 μm Most probable aperture of IG-110 2.5 μm Open porosity of BCεo 14.37 % Open porosity of IG-110εo 19.68 % Gas diffusion coefficientDAB of binary gas (H2O + He) 1.226×10−6 m2/s Effective diffusion coefficientDe D e =D AB × ε o 2 m2/s PermeabilityK of BC 0.49×10−15 m2 PermeabilityK of IG-110 1.97×10−15 m2


## 4.2. Result Analysis of Simulation

According to the numerical simulation, the distribution of moisture contentCr (%) in HTR-PM core at timet = 50 d is shown in Figure 3, where d represents units day. It can be found from Figure 3(a) that moisture mainly concentrates at the top, bottom, and side of core. Figure 3(b) presents the moisture distribution at the axial cross section. It can be found that the maximum Cr locates at the top of core. In addition, Figure 3(c) and Figure 3(d) further demonstrate the moisture distribution at the top and bottom surface, respectively. It can be seen that Cr near the cold gas and control rod channels, fuel feed tube, and fuel discharge tube is relatively low, while it is much higher elsewhere.Figure 3
Distribution ofCr (%) in HTR-PM core (t= 50 d).
(a)
Horizontal cross sections (b)
Axial cross section (c)
Top surface (d)
Bottom surfaceThe curves of moisture content versus time are shown in Figure4. It can be found that the average moisture content Cr,a of internals and fuel ball decreases exponentially, in which it can be explained that moisture content decreases much faster in the early stage because of the high moisture concentration gradient and decreases slower in the later stage as a result of the decrease of concentration gradient. Cr,a of fuel ball decreases far faster than that of core internals. Thus, the internals, not fuel balls, are the major limiting factor for the dehumidification. It can also be found that the required time to reach Cr,a = 0.02%, 0.01%, and 0.005% of whole internals is 43 d, 72 d, and 108 d, respectively. In addition, Cr,a of reflecting layer (IG-110) decreases much faster than that of insulation layer (BC), which can be attributed to the larger open porosity εo of IG-110. Thus effective diffusion coefficient De of IG-110 is much larger, which promotes the moisture transfer.Figure 4
Curves of moisture versus time.


## 4.3. Sensitivity Analysis

In order to accelerate the moisture transfer to improve the dehumidification efficiency, it is reasonable to promote the molecular diffusion and seepage flow. According to the relations,De∝T1.5 and De∝P-1, molecular diffusion can be promoted by increasing the temperatureT and decreasing the system pressurePo, while seepage can be promoted by increasing the pressure gradient, which can be achieved by decreasing the system pressure Po.The influence of temperatureT on dehumidification of HTR-PM is shown in Figure 5(a), whereCr,a on Y-axis denotes the average moisture content of whole internals obtained by simulation. It can be found that the dehumidification efficiency increases as temperature increases. Figure 5(b) presents the influence of system pressure Po on dehumidification, and it shows that the dehumidification efficiency of 0.1 MPa is certainly much higher than that of 7 MPa. In addition, in the simulation where Po=0.1 MPa, the diffusion-seepage model is adopted and De is considered constant to simplify the computation.Figure 5
Curves of moisture versus time at different temperature and pressure.
(a)
Temperature (b)
PressureAs for the actual dehumidification for HTR-PM,T = 250°C andPo = 7 MPa. Since the temperature is already pretty high, the promotion by increasing temperature is limited, while as for the system pressure, if it is reduced greatly from 7 MPa to 0.1 MPa, this change will promote both the diffusion and seepage significantly. But on the other hand, the decrease of pressure is not advantageous for the core to be heated. It is because the heat of core is transformed from mechanical energy of helium circulator; thus the heating power is relatively low at low pressure.


## 4.4. Analytical Solution

In order to obtain an analytical solution of the mathematical model, the irregular geometric model of core internals can be represented by an infinite plate, whose characteristic lengthLc (m) is the same as that of irregular geometry. In addition, Lc=V/Ae, whereV (m3) and Ae (m2) are the volume and effective surface area of geometry, respectively.The infinite plate with2Lc thickness is shown in Figure 6. Diffusion model (8) is adopted. The moisture content is Cr(x,t), and effective diffusion coefficient De is constant. Thus, the governing equation can be rewritten as(12)∂Cr∂t=De∂2Cr∂x2The boundary condition is Cr(x=Lc,t)=Cr,b. The initial condition is Cr(x,t=0)=Cr,0. Thus, the analytical solution, which is an infinite series, can be obtained as(13)Cr=Cr,b+∑n=1∞an·e-bnt·coscnxwhere cn=(n-0.5)π/Lc (n≥1, and n is an integer), bn=De·cn2, and an=Cr,0-Cr,b.Figure 6
Infinite plate.In order to simplify the solution, only the first polynomial of (13), which means n=1, is adopted. Thus, the average moisture content Cr,a can be calculated, which is(14)Cr,a=∫0LcCrdxLc=Cr,b+2Cr,0-Cr,bπ·e-Deπ2/4Lc2tAs for the geometric model of internals of HTR-PM, according to the statistic obtained in Solidworks software, its Lc is 0.3066 m. In addition, considering Cr,b=0 and Cr,0=0.1%, the analytical solution for dehumidification of HTR-PM can be obtained.Figure7 shows the comparison between analytical and numerical simulation. It is found that curvesCr,a-t of analytical and numerical solution of infinite plate are in good agreement. While comparing curves of analytical and numerical solution of HTR-PM, it is found that the agreement is not that good. It is suggested that the simplified analytical solution is pretty accurate for an ideal regular geometry model, while it is difficult to obtain the characteristic lengthLc to accurately represent an irregular geometry. The uneven distribution of channels, chambers, and other irregular structures in reactor core results in the deviation of solution. It is also suggested that, in order to derive a more representativeLc, some optimization methods can be adopted to the calculation ofLc, e.g., introducing empirical correction coefficient, dividing the whole geometry into several sections, etc.Figure 7
Comparison of analytical and numerical solutions.


## 5. Conclusion

In the dehumidification process of HTGR, the major moisture transfer mechanism in carbon materials includes molecular diffusion and seepage. Thus, diffusion model and diffusion-seepage model are developed to describe the dehumidification process of internals in reactor core. The simulation result of dehumidification of HTR-PM suggests that the average moisture content of internals decreases exponentially with time and the moisture mainly concentrates in the top and bottom of core. According to the sensitivity analysis, the dehumidification efficiency can be promoted by the increase of temperature and decrease of system pressure, especially by reducing the system pressure, and the dehumidification efficiency can be greatly improved. The analysis on analytical solution of dehumidification problem shows that the simplified analytical solution is pretty accurate for an ideal regular geometry model, while the calculation of characteristic length of irregular geometry needs to be improved.



---

*Source: 4163914-2018-08-05.xml*