Design of Laminated Seal for Triple Offset Butterfly Valve (350 °C) Used in Combined Cycle Power Plants

In combined cycle power plants (CCPPs), the bypass butterfly valve is a key component to facilitate regulation of exhaust gas energy available at the turbine and to not produce too much boost pressure. The conventional damper valve causes leakage, back flow into the turbine, and damage of the blade, and the existing dual-layered seal with polytetrafluoroethylene (PTFE) and metal should be frequently replaced owing to its low durability and deterioration of mechanical properties under a high temperature. This study devised a triple offset butterfly valve with a new type of seal by alternatively laminating stainless steel and graphite to improve valve performance at the high temperature (350 °C). The slope angles of the seal contact surface to prevent friction were calculated using the mathematical models of the triple offset. Thermal-structure coupled analyses by varying the number of graphite and thickness were conducted, and the seven-layer model with the graphite thickness of 0.8 mm, which shows airtightness and smooth operation, was chosen. The contact stresses behaviors of the graphite at 350 °C and at −196 °C were investigated, and it was found that the graphite is in charge of improving driving performance of the disc at the high temperature and sealing performance at the cryogenic temperature. The performance tests and the field tests of the suggested model verified its performance at the working temperature.


Introduction
As global warming and the depletion of energy resources have intensified, the growth in worldwide energy demand currently faces the difficulty of installing new power generation facilities as a result of limited funding and the strengthening of environmental regulations. Thus, combined cycle power plants (CCPPs), which consist of gas turbines and steam turbines, are becoming more important as the global demand for electrical power increases [1,2]. CCPPs manipulate discarded heat energy and store in hot gas into rotational energy, and the gas turbine compresses air and mixes it with fuel that is heated to a high temperature. The hot air-fuel mixture moves through the gas turbine blades to make them spin, and then the fast-spinning turbine drives the generator that converts a portion of the spinning energy into electricity [3][4][5]. With the stoichiometric burn, the power output of turbine drops sharply for flow rates above a critical value as a result of aerodynamic losses produced by the rotor blade stalling, and the turbine performs poorly in excessive load and exhaust gas pressure-states, which brings about a decrease in energy production.
In order to prevent the above problems, the bypass system facilitates the regulation of the exhaust gas energy available at the turbine wheel, and thus allows the turbocharger to be controlled and to not produce too much boost pressure. In this system, the bypass valve, which is mounted between the chamber and the atmosphere in parallel with the turbine, limits the maximum pressure in the chamber to prevent the instantaneous air flow rate through the turbine from exceeding the values above which aerodynamic stalling at the rotor blades would produce a severe fall in power output [6][7][8]. As the modern CCPPs require an increase of exhaust pressure in order to fulfill the restrictions of pollutant emission and to improve energy productivity, the following complementary functions have to be taken into account while designing the bypass valve: resistance to the harsh engine boundary conditions (temperature of 350 • C and working pressure of 1 MPa); smooth operation of the disc during opening and closing; and maintenance of sealing in closed valve operation to maximize gas energy directed to the turbine [9,10].
In the damper valve, which had been used in the existing bypass system, thermal deformations and stress concentrations occurred at a high temperature because there is no offset of the disc rotary axis from the seal. Also, a gap between the body and the disc due to absence of seal caused leakage, back flow into the turbine, and damage of the blade, as illustrated in Figure 1.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 2 of 19 in the chamber to prevent the instantaneous air flow rate through the turbine from exceeding the values above which aerodynamic stalling at the rotor blades would produce a severe fall in power output [6][7][8]. As the modern CCPPs require an increase of exhaust pressure in order to fulfill the restrictions of pollutant emission and to improve energy productivity, the following complementary functions have to be taken into account while designing the bypass valve: resistance to the harsh engine boundary conditions (temperature of 350 °C and working pressure of 1 MPa); smooth operation of the disc during opening and closing; and maintenance of sealing in closed valve operation to maximize gas energy directed to the turbine [9,10].
In the damper valve, which had been used in the existing bypass system, thermal deformations and stress concentrations occurred at a high temperature because there is no offset of the disc rotary axis from the seal. Also, a gap between the body and the disc due to absence of seal caused leakage, back flow into the turbine, and damage of the blade, as illustrated in Figure 1. The dual-layered seal with polytetrafluoroethylene (PTFE) and metal had been adopted to solve the above-mentioned problems, but PTFE should be frequently replaced owing to its low durability. Metal has a worse sealing performance than PTFE [11,12], and its mechanical properties are deteriorated at the high temperature [13], so new sealing material, particularly composite material, is needed to improve the sealing ability and durability of other original seals.
The existing single and double offset butterfly valves (SOBV and DOVB) have simple geometric structures to be easily manufactured, but rubbing and friction between the seal and the disc cause life shortening and wear of the seal during opening and closing, as shown in Figure 2. Therefore, this study applies the triple offset (radial eccentricity, axis eccentricity, and cone axis angle eccentricity) butterfly valve (TOBV), whose surface of the seal is in full contact when the valve is fully closed. This eliminates galling and minimizing seal wear and rubbing for long life and tight shut-off. In the previous studies related to the TOBV, Chen determined suitable triple eccentricities to improve the seal performance, and checked the static and dynamic interference through motion simulation [14]. Liang developed a computer programs for calculate the feasible the three eccentricities and the value of the cone's half angle [15]. Kan provided a criterion for non-interference for the metal seal pair considering the radial offset, the diameter of valve disc, the sealing thickness, and the cone angle of the sealing surface [16]. Bodhayan validated the PTFE seal design for gas lift valves through checking leakage amount and performed finite element analysis (FEA) on the seal design under a high pressure to predict the sealing contact pressures [17]. Ping implemented FEA with different sealing ring surface using the numerical method of symmetrical penalty function. The sealing ring surface and the sealing ratio pressure were suggested to enhance the sealing performance The dual-layered seal with polytetrafluoroethylene (PTFE) and metal had been adopted to solve the above-mentioned problems, but PTFE should be frequently replaced owing to its low durability. Metal has a worse sealing performance than PTFE [11,12], and its mechanical properties are deteriorated at the high temperature [13], so new sealing material, particularly composite material, is needed to improve the sealing ability and durability of other original seals.
The existing single and double offset butterfly valves (SOBV and DOVB) have simple geometric structures to be easily manufactured, but rubbing and friction between the seal and the disc cause life shortening and wear of the seal during opening and closing, as shown in Figure 2. Therefore, this study applies the triple offset (radial eccentricity, axis eccentricity, and cone axis angle eccentricity) butterfly valve (TOBV), whose surface of the seal is in full contact when the valve is fully closed. This eliminates galling and minimizing seal wear and rubbing for long life and tight shut-off.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 2 of 19 in the chamber to prevent the instantaneous air flow rate through the turbine from exceeding the values above which aerodynamic stalling at the rotor blades would produce a severe fall in power output [6][7][8]. As the modern CCPPs require an increase of exhaust pressure in order to fulfill the restrictions of pollutant emission and to improve energy productivity, the following complementary functions have to be taken into account while designing the bypass valve: resistance to the harsh engine boundary conditions (temperature of 350 °C and working pressure of 1 MPa); smooth operation of the disc during opening and closing; and maintenance of sealing in closed valve operation to maximize gas energy directed to the turbine [9,10].
In the damper valve, which had been used in the existing bypass system, thermal deformations and stress concentrations occurred at a high temperature because there is no offset of the disc rotary axis from the seal. Also, a gap between the body and the disc due to absence of seal caused leakage, back flow into the turbine, and damage of the blade, as illustrated in Figure 1. The dual-layered seal with polytetrafluoroethylene (PTFE) and metal had been adopted to solve the above-mentioned problems, but PTFE should be frequently replaced owing to its low durability. Metal has a worse sealing performance than PTFE [11,12], and its mechanical properties are deteriorated at the high temperature [13], so new sealing material, particularly composite material, is needed to improve the sealing ability and durability of other original seals.
The existing single and double offset butterfly valves (SOBV and DOVB) have simple geometric structures to be easily manufactured, but rubbing and friction between the seal and the disc cause life shortening and wear of the seal during opening and closing, as shown in Figure 2. Therefore, this study applies the triple offset (radial eccentricity, axis eccentricity, and cone axis angle eccentricity) butterfly valve (TOBV), whose surface of the seal is in full contact when the valve is fully closed. This eliminates galling and minimizing seal wear and rubbing for long life and tight shut-off. In the previous studies related to the TOBV, Chen determined suitable triple eccentricities to improve the seal performance, and checked the static and dynamic interference through motion simulation [14]. Liang developed a computer programs for calculate the feasible the three eccentricities and the value of the cone's half angle [15]. Kan provided a criterion for non-interference for the metal seal pair considering the radial offset, the diameter of valve disc, the sealing thickness, and the cone angle of the sealing surface [16]. Bodhayan validated the PTFE seal design for gas lift valves through checking leakage amount and performed finite element analysis (FEA) on the seal design under a high pressure to predict the sealing contact pressures [17]. Ping implemented FEA with different sealing ring surface using the numerical method of symmetrical penalty function. The sealing ring surface and the sealing ratio pressure were suggested to enhance the sealing performance In the previous studies related to the TOBV, Chen determined suitable triple eccentricities to improve the seal performance, and checked the static and dynamic interference through motion simulation [14]. Liang developed a computer programs for calculate the feasible the three eccentricities and the value of the cone's half angle [15]. Kan provided a criterion for non-interference for the metal seal pair considering the radial offset, the diameter of valve disc, the sealing thickness, and the cone angle of the sealing surface [16]. Bodhayan validated the PTFE seal design for gas lift valves through checking leakage amount and performed finite element analysis (FEA) on the seal design under a high pressure to predict the sealing contact pressures [17]. Ping implemented FEA with different sealing ring surface using the numerical method of symmetrical penalty function. The sealing ring surface and the sealing ratio pressure were suggested to enhance the sealing performance of the rotating ball valve with double direction metal sealing [18]. Ahn proposed a design criterion to ensure the seat tightness in which the contact pressure between the metal seal and the valve disc would be compared with the fluid pressure, and mechanical behavior of a flexible solid metal seal for a cryogenic butterfly valve was investigated [19]. However, there were deficiencies of previous works in relation with design of laminated seal in TOBV to possess good sealing property, operability of the disc, and capability of withstanding high temperatures.
In this paper, a study on performance enhancement of the triple-offset butterfly valve used in the bypass system of a CCPP's gas turbine was conducted. After calculating the slope angles of the contact surface to eliminate rubbing on the seal contact surface, with them, the mathematic model of triple-offset was generated. In order to solve the problems of the existing valve, a laminated seal was devised by substituting PTFE for graphite, which is known to excel in extreme conditions, withstanding heat and pressure. The technique for the thermal-structural coupled analysis was established to consider thermal deformation because of exposure to the high temperature (350 • C), and the maximum contact stresses at the contact surface of the seal were derived by varying the number of graphite layers and thicknesses, slope angle of the contact surface (66.1 • and 86.1 • ) calculated and temperature (22 • C and 350 • C). On the basis of the analysis results, the seven-layer model with the graphite thickness of 0.8 mm was suggested to guarantee not only airtightness, but also smooth operation of the disc, and the structural reliability was evaluated. When the contact stresses' behavior of the graphite at 350 • C was compared with that at −196 • C, it was found that the graphite is in charge of improving driving performance of the disc at high temperature and sealing performance at cryogenic temperature. The internal pressure, leakage tests were implemented based on the international regulation, and the actual verification experiment for three months proved the durability and suitability for the actual working environment.

Structure of Triple Offset Butterfly Valve
The geometry of the triple offset butterfly valve and its component is shown in Figure 3. The disc rotates by the actuator to control amount of the exhaust gas, and the shaft connects the valve body to the actuator. The sealing part consists of the laminated seal and the retainer. of the rotating ball valve with double direction metal sealing [18]. Ahn proposed a design criterion to ensure the seat tightness in which the contact pressure between the metal seal and the valve disc would be compared with the fluid pressure, and mechanical behavior of a flexible solid metal seal for a cryogenic butterfly valve was investigated [19]. However, there were deficiencies of previous works in relation with design of laminated seal in TOBV to possess good sealing property, operability of the disc, and capability of withstanding high temperatures.
In this paper, a study on performance enhancement of the triple-offset butterfly valve used in the bypass system of a CCPP's gas turbine was conducted. After calculating the slope angles of the contact surface to eliminate rubbing on the seal contact surface, with them, the mathematic model of triple-offset was generated. In order to solve the problems of the existing valve, a laminated seal was devised by substituting PTFE for graphite, which is known to excel in extreme conditions, withstanding heat and pressure. The technique for the thermal-structural coupled analysis was established to consider thermal deformation because of exposure to the high temperature (350 °C), and the maximum contact stresses at the contact surface of the seal were derived by varying the number of graphite layers and thicknesses, slope angle of the contact surface (66.1° and 86.1°) calculated and temperature (22 °C and 350 °C). On the basis of the analysis results, the seven-layer model with the graphite thickness of 0.8 mm was suggested to guarantee not only airtightness, but also smooth operation of the disc, and the structural reliability was evaluated. When the contact stresses' behavior of the graphite at 350 °C was compared with that at −196 °C, it was found that the graphite is in charge of improving driving performance of the disc at high temperature and sealing performance at cryogenic temperature. The internal pressure, leakage tests were implemented based on the international regulation, and the actual verification experiment for three months proved the durability and suitability for the actual working environment.

Structure of Triple Offset Butterfly Valve
The geometry of the triple offset butterfly valve and its component is shown in Figure 3. The disc rotates by the actuator to control amount of the exhaust gas, and the shaft connects the valve body to the actuator. The sealing part consists of the laminated seal and the retainer.

Sealing Mechanism of Laminated Seal
This study devised a laminated seal by insertion of graphite with metal (A240-316) to aide in strength and stability for extremely high temperature applications, as shown in Figure 4. To achieve tight shut-off, the disc is driven into the valve body using larger and more expensive actuator than that used with other valve designs, thereby it is necessary to design a shape of laminated seal to achieve the high sealing property and smooth driving with low torque [20].

Sealing Mechanism of Laminated Seal
This study devised a laminated seal by insertion of graphite with metal (A240-316) to aide in strength and stability for extremely high temperature applications, as shown in Figure 4. To achieve tight shut-off, the disc is driven into the valve body using larger and more expensive actuator than that used with other valve designs, thereby it is necessary to design a shape of laminated seal to achieve the high sealing property and smooth driving with low torque [20].  A working pressure of 1 MPa is imposed inside the valve and then produces contact stress on the contact surface between the seal and the body at the final shut-off position. When the working pressure (Pworking) is higher than the contact stress (Pcontact), the fluid is squeezed into the sealing clearance; otherwise, sealing performance is guaranteed as expressed.

Mathematical Modeling of TOBV Seal Pair
Conventional single offset valve with center shaft, which penetrates a seal, is suitable for low temperature and low pressure services only, and double offset valve with eccentric shaft results in an uninterrupted seal that can be used at higher pressures and temperatures, but the friction and seal wear deteriorate sealing performance.
In the triple offset valve, both the seal and the disc are surfaces of a cone that is sectioned at an angle. The valve shaft is located slightly to one side of the seal center and above the plane of the seat, and its rotation center is also offset from the axis of the imaginary cone that extends from the surface of the seat. When the valve is closed, the surfaces of the seal and the disc are in full contact at all points, and opening the valve results in the disc moving away from the seal at all points, eliminating galling and minimizing seal wear due to non-rubbing sealing surfaces for long life and tight shut-off. The three separate offsets are designed in the butterfly valve, as shown in Figure 5.
Offset 1 (H): The shaft is offset from the seal plane providing an uninterrupted sealing surface. Offset 2 (E): The centerline of disc is offset from the centerline of the shaft allowing the disc to freely lift off and away from the seal during opening.
Offset 3 (γ): The centerline of seal cone with conical angle of 2α is offset by γ° from the centerline of the body to eliminate rubbing on the contact surface during opening and closing, thereby rotating the disc without interference [21]. A coordinate system is established by taking the apex of the cone, and the geometry of the triple offset seal pair is shown in Figure 6. A working pressure of 1 MPa is imposed inside the valve and then produces contact stress on the contact surface between the seal and the body at the final shut-off position. When the working pressure (P working ) is higher than the contact stress (P contact ), the fluid is squeezed into the sealing clearance; otherwise, sealing performance is guaranteed as expressed.

Mathematical Modeling of TOBV Seal Pair
Conventional single offset valve with center shaft, which penetrates a seal, is suitable for low temperature and low pressure services only, and double offset valve with eccentric shaft results in an uninterrupted seal that can be used at higher pressures and temperatures, but the friction and seal wear deteriorate sealing performance.
In the triple offset valve, both the seal and the disc are surfaces of a cone that is sectioned at an angle. The valve shaft is located slightly to one side of the seal center and above the plane of the seat, and its rotation center is also offset from the axis of the imaginary cone that extends from the surface of the seat. When the valve is closed, the surfaces of the seal and the disc are in full contact at all points, and opening the valve results in the disc moving away from the seal at all points, eliminating galling and minimizing seal wear due to non-rubbing sealing surfaces for long life and tight shut-off. The three separate offsets are designed in the butterfly valve, as shown in Figure 5.
Offset 1 (H): The shaft is offset from the seal plane providing an uninterrupted sealing surface. Offset 2 (E): The centerline of disc is offset from the centerline of the shaft allowing the disc to freely lift off and away from the seal during opening.
Offset 3 (γ): The centerline of seal cone with conical angle of 2α is offset by γ • from the centerline of the body to eliminate rubbing on the contact surface during opening and closing, thereby rotating the disc without interference [21].  A working pressure of 1 MPa is imposed inside the valve and then produces contact stress on the contact surface between the seal and the body at the final shut-off position. When the working pressure (Pworking) is higher than the contact stress (Pcontact), the fluid is squeezed into the sealing clearance; otherwise, sealing performance is guaranteed as expressed.

Mathematical Modeling of TOBV Seal Pair
Conventional single offset valve with center shaft, which penetrates a seal, is suitable for low temperature and low pressure services only, and double offset valve with eccentric shaft results in an uninterrupted seal that can be used at higher pressures and temperatures, but the friction and seal wear deteriorate sealing performance.
In the triple offset valve, both the seal and the disc are surfaces of a cone that is sectioned at an angle. The valve shaft is located slightly to one side of the seal center and above the plane of the seat, and its rotation center is also offset from the axis of the imaginary cone that extends from the surface of the seat. When the valve is closed, the surfaces of the seal and the disc are in full contact at all points, and opening the valve results in the disc moving away from the seal at all points, eliminating galling and minimizing seal wear due to non-rubbing sealing surfaces for long life and tight shut-off. The three separate offsets are designed in the butterfly valve, as shown in Figure 5.
Offset 1 (H): The shaft is offset from the seal plane providing an uninterrupted sealing surface. Offset 2 (E): The centerline of disc is offset from the centerline of the shaft allowing the disc to freely lift off and away from the seal during opening.
Offset 3 (γ): The centerline of seal cone with conical angle of 2α is offset by γ° from the centerline of the body to eliminate rubbing on the contact surface during opening and closing, thereby rotating the disc without interference [21]. A coordinate system is established by taking the apex of the cone, and the geometry of the triple offset seal pair is shown in Figure 6. A coordinate system is established by taking the apex of the cone, and the geometry of the triple offset seal pair is shown in Figure 6. The functions of cone where the disc is situated about y-axis are expressed in Equation (1), and the sections of the disc (Si) in Equation (2) depend on Hi, which are the distances between the coordinate origin (O) and Ni, given by Equation (3) (i = 1, 2, and 3).
where k = tan ( ) and k0 = tan ( ). When the sealing surface is truncated by an arbitrary plane that is in parallel with yOz plane with its x-coordinate as x0, the equation of the hyperbola (ajbckd) shown in Figure 7 is obtained by Equation (4), and the coordinate of two intersection points (x, yi, zi) between hyperbola and disc surface (Si) is obtained by Equation (5).
: Intermediate diameter of disc H0, H1, and H2 : Distances from the coordinate origin (O) to N0, N1 and N2 N0, N1 and N2 : Intersections between cone center and S0, S1 and S2 : Centerline of valve body O0 : Rotation center of the disc S1 , S2 : Upper and lower sections of disc S0 : Section of disc center T : Thickness of disc α : Half conical angle The functions of cone where the disc is situated about y-axis are expressed in Equation (1), and the sections of the disc (S i ) in Equation (2) depend on H i , which are the distances between the coordinate origin (O) and N i , given by Equation (3) (i = 1, 2, and 3).
where k = tan (α) and k 0 = tan (γ). When the sealing surface is truncated by an arbitrary plane that is in parallel with yOz plane with its x-coordinate as x 0 , the equation of the hyperbola (ajbckd) shown in Figure 7 is obtained by Equation (4), and the coordinate of two intersection points (x, y i , z i ) between hyperbola and disc surface (S i ) is obtained by Equation (5). The functions of cone where the disc is situated about y-axis are expressed in Equation (1), and the sections of the disc (Si) in Equation (2) depend on Hi, which are the distances between the coordinate origin (O) and Ni, given by Equation (3) (i = 1, 2, and 3).
where k = tan ( ) and k0 = tan ( ). When the sealing surface is truncated by an arbitrary plane that is in parallel with yOz plane with its x-coordinate as x0, the equation of the hyperbola (ajbckd) shown in Figure 7 is obtained by Equation (4), and the coordinate of two intersection points (x, yi, zi) between hyperbola and disc surface (Si) is obtained by Equation (5).
: Intermediate diameter of disc H0, H1, and H2 : Distances from the coordinate origin (O) to N0, N1 and N2 N0, N1 and N2 : Intersections between cone center and S0, S1 and S2 : Centerline of valve body O0 : Rotation center of the disc S1 , S2 : Upper and lower sections of disc S0 : Section of disc center T : Thickness of disc α : Half conical angle z 2 The elliptical equations for the intersections between the disc cutting planes and the curved surface of cone are expressed by Equation (6).
The intersections between the vertical lines of N 0 , N 1 , and N 2 with the plane x = x 0 are the center point of jk, bc, and ad, and apexes N i of minor axis of the intersection constitutes the line N i O. Friction behavior of TOBV seal pair depends on the following five parameters: offset 1 (E), offset 2 (H), offset 3 (γ), half conical angle (α), and disc thickness (T). During the valve opening and closing, point P(y, z), which is on the hyperbola ab, rotates about O 0 with speed of → PN perpendicular to O 0 P, as shown in Figure 8. The angle of θ in Equation (7) is the friction angle between the seal pair (PM and PN), and the angle of β in Equation (8) is between PM and the axial line crossing the seal point and parallel with y-axis. The angle of θ 1 in Equation (9) is between O 0 P and the vertical line crossing the rotary center (O 0 ). The higher tical line crossing the rotary center (allel with fiv, which means that the disc can rapidly be disengaged from the valve seal, and the friction moment decreases [15].
Appl. Sci. 2019, 9, x FOR PEER REVIEW 6 of 19 The elliptical equations for the intersections between the disc cutting planes and the curved surface of cone are expressed by Equation (6).
The intersections between the vertical lines of N0, N1, and N2 with the plane x = x0 are the center point of jk, bc, and ad, and apexes Ni of minor axis of the intersection constitutes the line . Friction behavior of TOBV seal pair depends on the following five parameters: offset 1 (E), offset 2 (H), offset 3 (γ), half conical angle (α), and disc thickness (T). During the valve opening and closing, point P(y, z), which is on the hyperbola ab, rotates about 0 O with speed of ⃗ perpendicular to , as shown in Figure 8. The angle of θ in Equation (7) is the friction angle between the seal pair ( and ), and the angle of β in Equation (8) is between and the axial line crossing the seal point and parallel with y-axis. The angle of θ1 in Equation (9) is between and the vertical line crossing the rotary center (O0). The higher tical line crossing the rotary center (allel with fiv, which means that the disc can rapidly be disengaged from the valve seal, and the friction moment decreases [15].

Calculation of Slope Angle at Contact Surface
When the TOBV with the inner diameter of 200 mm is subjected to the temperature of 350 °C and the internal pressure of 1 MPa, the diameter and the thickness of the disc, the outer diameter and the thickness of the seal, the three offsets (E, H, and γ), and the half conical angle (α) are given from the international regulations (API609 and ASME B16.10). On the basis of the dimensions shown in Table 1, the geometry of the disc surface was created using commercial software, GEOGEBRA, as follows. (Step 1) Cone surface is generated by inputting α (13.9°) using Equation (3), as shown in Figure 9a. (Step 2) When γ (10°), disc thickness (15 mm) is given, H0 = 369.9 mm, H1 = 377.5 mm, H2 = 362.3 mm and rotary center of disc is obtained from Equation (5), and then S0, S1, and S2 are generated using Equation (4), as shown in Figure 9b.

Calculation of Slope Angle at Contact Surface
When the TOBV with the inner diameter of 200 mm is subjected to the temperature of 350 • C and the internal pressure of 1 MPa, the diameter and the thickness of the disc, the outer diameter and the thickness of the seal, the three offsets (E, H, and γ), and the half conical angle (α) are given from the international regulations (API609 and ASME B16.10). On the basis of the dimensions shown in Table 1, the geometry of the disc surface was created using commercial software, GEOGEBRA, as follows.

(Step 3)
The cone surface is cut by x = x 0 (0.05, green plane) to obtain hyperbola in Equation (6), which enables the outer diameter of disc to be 190.1 mm, and then the disc surface is obtained as shown in Figure 9c. (Step 4) Rotary center of the disc moves by the eccentricity H (23 mm) and E (2.5 mm) to avoid interference between the seal and the disc during opening and closing, as shown in Figure 9d. Step 3) The cone surface is cut by x = x0 (0.05, green plane) to obtain hyperbola in Equation (6), which enables the outer diameter of disc to be 190.1 mm, and then the disc surface is obtained as shown in Figure 9c. (Step 4) Rotary center of the disc moves by the eccentricity H (23 mm) and E (2.5 mm) to avoid interference between the seal and the disc during opening and closing, as shown in Figure  9 (d).

Geometry and Analysis Conditions
On the basis of the above generated geometry, 2D modeling of the sealing part at final shut-off position was conducted using a commercial software, Inventor 2019, as shown in Figure 10. The contact stress rises with the increase of the slope angle based on the Hertzan theory shown in Equation (10) and Figure 11, so simulations were carried out only in the case of minimum (66.1°) and maximum angle (86.1°).

Geometry and Analysis Conditions
On the basis of the above generated geometry, 2D modeling of the sealing part at final shut-off position was conducted using a commercial software, Inventor 2019, as shown in Figure 10. The contact stress rises with the increase of the slope angle based on the Hertzan theory shown in Equation (10) and Figure 11, so simulations were carried out only in the case of minimum (66.1 • ) and maximum angle (86.1 • ).   Figure 12. Full-thickness (7 mm) of the seal is fixed, so thickness of 1-layer stainless steel depending on that of the graphite was calculated according to the number of laminated layers, as listed in Table 2. The analysis models were simulated according to the temperatures (the room temperature of 22 °C and the high temperature of 350 °C) using commercial finite element software, ANSYS workbench ver. 18.   In order to investigate contact stress behavior according to the seal design, a series of finite element (FE) simulations were performed with the different total number of seal layers by laminating alternately the American Society for Testing and Materials (ASTM) A240 Grade 316 stainless steel (SS) and the graphite (G) (1-layer (SS), 3-layer (SS/G/SS), 5-layer (SS/G/SS/G/SS), 7-layer (SS/G/SS/G/SS/G/SS)) and thickness of graphite (0.6 mm~1.2 mm at the interval of 0.2 mm), as shown in Figure 12. Full-thickness (7 mm) of the seal is fixed, so thickness of 1-layer stainless steel depending on that of the graphite was calculated according to the number of laminated layers, as listed in Table 2. The analysis models were simulated according to the temperatures (the room temperature of 22 °C and the high temperature of 350 °C) using commercial finite element software, ANSYS workbench ver. 18.  In order to investigate contact stress behavior according to the seal design, a series of finite element (FE) simulations were performed with the different total number of seal layers by laminating alternately the American Society for Testing and Materials (ASTM) A240 Grade 316 stainless steel (SS) and the graphite (G) (1-layer (SS), 3-layer (SS/G/SS), 5-layer (SS/G/SS/G/SS), 7-layer (SS/G/SS/G/SS/G/SS)) and thickness of graphite (0.6 mm~1.2 mm at the interval of 0.2 mm), as shown in Figure 12. Full-thickness (7 mm) of the seal is fixed, so thickness of 1-layer stainless steel depending on that of the graphite was calculated according to the number of laminated layers, as listed in Table 2. The analysis models were simulated according to the temperatures (the room temperature of 22 • C and the high temperature of 350 • C) using commercial finite element software, ANSYS workbench ver. 18.  In order to investigate contact stress behavior according to the seal design, a series of finite element (FE) simulations were performed with the different total number of seal layers by laminating alternately the American Society for Testing and Materials (ASTM) A240 Grade 316 stainless steel (SS) and the graphite (G) (1-layer (SS), 3-layer (SS/G/SS), 5-layer (SS/G/SS/G/SS), 7-layer (SS/G/SS/G/SS/G/SS)) and thickness of graphite (0.6 mm~1.2 mm at the interval of 0.2 mm), as shown in Figure 12. Full-thickness (7 mm) of the seal is fixed, so thickness of 1-layer stainless steel depending on that of the graphite was calculated according to the number of laminated layers, as listed in Table 2. The analysis models were simulated according to the temperatures (the room temperature of 22 °C and the high temperature of 350 °C) using commercial finite element software, ANSYS workbench ver. 18.   The chemical composition and mechanical properties of the ASTM A240 Grade 316 stainless steel are listed in Table 3. Thermal and mechanical properties of A240-316 (retainer and steel layer), A216 Wrought Carbon with Grade B (WCB) (body), and graphite according to the temperatures (25 • C~350 • C) are listed in Table 4 [16,22,23]. Thermal-structural coupled analysis was conducted to calculate contact stress considering deformation caused by the thermal loads and the working pressure (1 MPa). In the thermal analysis, the convective heat transfer coefficient, 3.675 × 10 −3 W/mm 2 · • C, was employed to the edges of the body, the retainer and the seal, which was exposed to the room temperature, and the thermal load of 350 • C was imposed on the edges, where contact with the working fluid occurred, as shown in Figure 13a. Friction coefficients between the disc and the seal were 0.3 for stainless steel and 0.21 for graphite, respectively, as shown in Figure 13b [19].
Appl. Sci. 2019, 9, x FOR PEER REVIEW 11 of 19 Thermal-structural coupled analysis was conducted to calculate contact stress considering deformation caused by the thermal loads and the working pressure (1 MPa). In the thermal analysis, the convective heat transfer coefficient, 3.675 × 10 −3 W/mm 2 ·°C, was employed to the edges of the body, the retainer and the seal, which was exposed to the room temperature, and the thermal load of 350 °C was imposed on the edges, where contact with the working fluid occurred, as shown in Figure  13a. Friction coefficients between the disc and the seal were 0.3 for stainless steel and 0.21 for graphite, respectively, as shown in Figure 13b [19]. The global hexagonal grid size of the disc, the body, and the retainer was 1 mm, and that of the seal was 0.5 mm. The denser sizing of 0.1 mm was set up to the edges of the seal and the body at which contact stress occurs as shown in Figure 14a, and the number of nodes was 555,953 and that of elements was 182,554. The orthogonal quality of 98.9% was close to 1, and skewness of 92.8% was close to 0 as shown in Figure 14b, which verifies the FE model.   The global hexagonal grid size of the disc, the body, and the retainer was 1 mm, and that of the seal was 0.5 mm. The denser sizing of 0.1 mm was set up to the edges of the seal and the body at which contact stress occurs as shown in Figure 14a, and the number of nodes was 555,953 and that of elements was 182,554. The orthogonal quality of 98.9% was close to 1, and skewness of 92.8% was close to 0 as shown in Figure 14b, which verifies the FE model. Thermal-structural coupled analysis was conducted to calculate contact stress considering deformation caused by the thermal loads and the working pressure (1 MPa). In the thermal analysis, the convective heat transfer coefficient, 3.675 × 10 −3 W/mm 2 ·°C, was employed to the edges of the body, the retainer and the seal, which was exposed to the room temperature, and the thermal load of 350 °C was imposed on the edges, where contact with the working fluid occurred, as shown in Figure  13a. Friction coefficients between the disc and the seal were 0.3 for stainless steel and 0.21 for graphite, respectively, as shown in Figure 13b [19]. The global hexagonal grid size of the disc, the body, and the retainer was 1 mm, and that of the seal was 0.5 mm. The denser sizing of 0.1 mm was set up to the edges of the seal and the body at which contact stress occurs as shown in Figure 14a, and the number of nodes was 555,953 and that of elements was 182,554. The orthogonal quality of 98.9% was close to 1, and skewness of 92.8% was close to 0 as shown in Figure 14b, which verifies the FE model.   The modeling and outputs of the thermal analysis were coupled to the static structural analysis, as shown in Figure 15a. In this step, deformations at the disc, which were caused by torque applied on the shaft to rotate the disc and to maintain it at the desired position, should be considered [24]. Moment load condition is usable for only axisymmetric geometry in 2D analysis, but is inapplicable because the 2D model is asymmetric because of its eccentricities, thus the maximum deformation in the y-axis (−5.5 × 10 −3 mm) at the disc, which was obtained from the 3D simulation, was considered in the 2D simulation, as shown in Figure 15b. Also, the temperature distribution with the maximum value of 350 • C and the minimum of 36.5 • C was imported to the boundary condition, and the working pressure (1 MPa) was applied as shown in Figure 15c. The body and the retainer were fixed, and the x-displacement, which was a relatively very small value (1.7 × 10 −4 mm), was not considered. The modeling and outputs of the thermal analysis were coupled to the static structural analysis, as shown in Figure 15a. In this step, deformations at the disc, which were caused by torque applied on the shaft to rotate the disc and to maintain it at the desired position, should be considered [24]. Moment load condition is usable for only axisymmetric geometry in 2D analysis, but is inapplicable because the 2D model is asymmetric because of its eccentricities, thus the maximum deformation in the y-axis (−5.5 × 10 −3 mm) at the disc, which was obtained from the 3D simulation, was considered in the 2D simulation, as shown in Figure 15b. Also, the temperature distribution with the maximum value of 350 °C and the minimum of 36.5 °C was imported to the boundary condition, and the working pressure (1 MPa) was applied as shown in Figure 15c. The body and the retainer were fixed, and the x-displacement, which was a relatively very small value (1.7 × 10 −4 mm), was not considered.

Results of the Thermal-Structural Analysis
High contact stress improves sealing performance, while it does not allow the smooth driving of the disc. It is thus necessary to find an effective trade-off between needs for good operation and prevention of leakage. Contact stress, which is a little higher than the working pressure, could guarantee not only airtightness, but also smooth rotation of the disc. Tables 5 and 6 represent the maximum contact stresses for the different seal designs, temperatures (25 °C and 350 °C) and slope angles (66.1° and 86.1°). SS and G indicate stainless steel and graphite, and min and max mean the minimum (66.1°) and maximum (86.1°) slope angles.

Results of the Thermal-Structural Analysis
High contact stress improves sealing performance, while it does not allow the smooth driving of the disc. It is thus necessary to find an effective trade-off between needs for good operation and prevention of leakage. Contact stress, which is a little higher than the working pressure, could guarantee not only airtightness, but also smooth rotation of the disc. Tables 5 and 6 represent the maximum contact stresses for the different seal designs, temperatures (25 • C and 350 • C) and slope angles (66.1 • and 86.1 • ). SS and G indicate stainless steel and graphite, and min and max mean the minimum (66.1 • ) and maximum (86.1 • ) slope angles. For both stainless steel and graphite, the models with the slope angle of 66.1 • showed lower contact stresses than those of 86.1 • , which proves that the increases of the contact surface led to a reduction of the peak contact stress, and sealing performance was improved with increasing temperature as a result of thermal expansion.
In respect to the graphite layers, leakages were expected for the all models at 25 • C (G min : 0.19~0.28 MPa, G max : 0.14~0.18 MPa), and there were small changes in contact stress by varying the amount of the graphite, whereas at 350 • C, the larger contact stresses (G min : 12.33~4.73 MPa, G max : 5.03~0.79 MPa) beyond the working pressure were observed. This means that the graphite displayed a sealing effect at high temperature rather than at room temperature, and the contact stress decreases with increment of the number of laminated layer. In respect to the stainless steel layers, the graphite did not withstand the contact force at 25 • C, so the contact stress at the stainless steel increased with more amount of the graphite (SS min : 2.0 MPa → 6.3 MPa, SS max : 1.2 MPa → 4.9 MPa) because the contact force was intensively applied on the stainless steel. On the other hand, at 350 • C, the contact force was distributed to the graphite, so more graphite leads to a reduction of the peak contact stress at the stainless steel (SS min : 169.  17 MPa,SSmax: 199.71 MPa) at 350 °C, but they were still too large. The 7-layer model with the graphite thickness of 1.2 mm showed the lowest contact stresses of the stainless steel (SSmin: 74.13 MPa, SSmax: 63.02 MPa) at 350 °C, but the maximum contact stresses of the graphite were below the working pressure (Gmax: 0.79 MPa). Therefore, the 7-layer model with the graphite thickness of 0.8 mm was suggested to guarantee not only airtightness, but also smooth rotation of the disc, and its contact stress distributions are shown in Figures 16 and 17.
Analyses results of the new laminated seal at 350 °C were compared to those at the 1-layer, 3layer, and 5-layer models in Table 7. The maximum contact stresses in case of the slope angle of 66.1° were reduced by 65.5%, 42.6%~15.8%, and 18.3%~7.3%, and those in case of the slope angle of 86.1° were reduced by 66.5%, 48.8%~29.9%, and 35.8%~12.7%, respectively, which means that the shape design of the laminated seal allows to mitigate the excessive contact stress.

Discussions-Effect of Graphite on Sealing Ability
On the basis of the previous study in relation with the cryogenic triple-offset butterfly valve with the laminated seal, which is installed in liquefied natural gas (LNG) marine engine, to control the flow of liquid nitrogen (−196 °C) to liquefy natural gas [25], the contact stress behaviors under the three different temperatures (−196 °C, 25 °C, and 350 °C) were analyzed to discover the graphite effect 282.17 MPa, SSmax: 199.71 MPa) at 350 °C, but they were still too large. The 7-layer model with the graphite thickness of 1.2 mm showed the lowest contact stresses of the stainless steel (SSmin: 74.13 MPa, SSmax: 63.02 MPa) at 350 °C, but the maximum contact stresses of the graphite were below the working pressure (Gmax: 0.79 MPa). Therefore, the 7-layer model with the graphite thickness of 0.8 mm was suggested to guarantee not only airtightness, but also smooth rotation of the disc, and its contact stress distributions are shown in Figures 16 and 17.
Analyses results of the new laminated seal at 350 °C were compared to those at the 1-layer, 3layer, and 5-layer models in Table 7. The maximum contact stresses in case of the slope angle of 66.1° were reduced by 65.5%, 42.6%~15.8%, and 18.3%~7.3%, and those in case of the slope angle of 86.1° were reduced by 66.5%, 48.8%~29.9%, and 35.8%~12.7%, respectively, which means that the shape design of the laminated seal allows to mitigate the excessive contact stress.  On the basis of the previous study in relation with the cryogenic triple-offset butterfly valve with the laminated seal, which is installed in liquefied natural gas (LNG) marine engine, to control the flow of liquid nitrogen (−196 °C) to liquefy natural gas [25], the contact stress behaviors under the three different temperatures (−196 °C, 25 °C, and 350 °C) were analyzed to discover the graphite effect Analyses results of the new laminated seal at 350 • C were compared to those at the 1-layer, 3-layer, and 5-layer models in Table 7. The maximum contact stresses in case of the slope angle of 66.1 • were reduced by 65.5%, 42.6%~15.8%, and 18.3%~7.3%, and those in case of the slope angle of 86.1 • were reduced by 66.5%, 48.8%~29.9%, and 35.8%~12.7%, respectively, which means that the shape design of the laminated seal allows to mitigate the excessive contact stress. Table 7. Comparison of peak contact stress at the new seal with those at the 1-layer, 3-layer, and 5-layer models (350 • C).

Discussions-Effect of Graphite on Sealing Ability
On the basis of the previous study in relation with the cryogenic triple-offset butterfly valve with the laminated seal, which is installed in liquefied natural gas (LNG) marine engine, to control the flow of liquid nitrogen (−196 • C) to liquefy natural gas [25], the contact stress behaviors under the three different temperatures (−196 • C, 25 • C, and 350 • C) were analyzed to discover the graphite effect at each temperature. In the stainless steel layer, the contact stress was increased at 350 • C as a result of undergoing thermal expansion, while in the graphite layer, the larger contact stresses were observed at 350 • C and −196 • C rather than 25 • C, which indicates that the graphite is a favorable for sealing property at severe temperature. As more laminated layers of the graphite, the contact stresses at both graphite and stainless steel decrease at 350 • C, but opposite behavior was observed at −196 • C. It is worth noticing that the graphite contributes improving driving performance of the disc at the high temperature and sealing ability at the cryogenic temperature, as shown in Table 8. In order to verify the reliability of the TOBV with the proposed laminated seal, the 3D thermal-structural coupled analysis was conducted as shown in Figure 18. The simulation conditions (convective heat transfer coefficient, thermal load, friction coefficients, and torque) were the same as in the 2D analysis. Table 9 compares the maximum equivalent stress at each part with the allowable strength [26], and it is expected that all components are not yielded at 350 • C. Also, the contact surfaces of the final model showed the sticking status as shown in Figure 19, so deboning of the seal layers due to the different expansion coefficients does not occur.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 15 of 19 at each temperature. In the stainless steel layer, the contact stress was increased at 350 °C as a result of undergoing thermal expansion, while in the graphite layer, the larger contact stresses were observed at 350 °C and −196 °C rather than 25 °C, which indicates that the graphite is a favorable for sealing property at severe temperature. As more laminated layers of the graphite, the contact stresses at both graphite and stainless steel decrease at 350 °C, but opposite behavior was observed at −196 °C. It is worth noticing that the graphite contributes improving driving performance of the disc at the high temperature and sealing ability at the cryogenic temperature, as shown in Table 8. In order to verify the reliability of the TOBV with the proposed laminated seal, the 3D thermalstructural coupled analysis was conducted as shown in Figure 18. The simulation conditions (convective heat transfer coefficient, thermal load, friction coefficients, and torque) were the same as in the 2D analysis. Table 9 compares the maximum equivalent stress at each part with the allowable strength [26], and it is expected that all components are not yielded at 350 °C. Also, the contact surfaces of the final model showed the sticking status as shown in Figure 19, so deboning of the seal layers due to the different expansion coefficients does not occur.     at each temperature. In the stainless steel layer, the contact stress was increased at 350 °C as a result of undergoing thermal expansion, while in the graphite layer, the larger contact stresses were observed at 350 °C and −196 °C rather than 25 °C, which indicates that the graphite is a favorable for sealing property at severe temperature. As more laminated layers of the graphite, the contact stresses at both graphite and stainless steel decrease at 350 °C, but opposite behavior was observed at −196 °C. It is worth noticing that the graphite contributes improving driving performance of the disc at the high temperature and sealing ability at the cryogenic temperature, as shown in Table 8. In order to verify the reliability of the TOBV with the proposed laminated seal, the 3D thermalstructural coupled analysis was conducted as shown in Figure 18. The simulation conditions (convective heat transfer coefficient, thermal load, friction coefficients, and torque) were the same as in the 2D analysis. Table 9 compares the maximum equivalent stress at each part with the allowable strength [26], and it is expected that all components are not yielded at 350 °C. Also, the contact surfaces of the final model showed the sticking status as shown in Figure 19, so deboning of the seal layers due to the different expansion coefficients does not occur.   . Contact status between stainless and graphite layers. Figure 19. Contact status between stainless and graphite layers.

Internal Pressure and Leakage Tests
The performance tests of the proposed TOBV were conducted using the experimental equipment, based on the regulation of KS B 2304:2001 (General rules for inspection of valves). The boiler and the temperature sensor controlled the temperature of 350 • C, and the pressure gauge checked the internal pressure. In the internal pressure test, the valve body was filled with the water at the half-shut position, and then the internal pressure of 1.7 MPa, which exceeded the working pressure (1 MPa) × 1.5, was maintained for 182 s. In the leakage test, the valve body was filled with the water at the fully shut-off position, and the internal pressure of 1.15 MPa, which exceeded the working pressure (1 MPa) × 1.1, was maintained for 125 s. The results of the flowrate measurement exhibit that there were no leakage and exudation in both tests, as shown in Table 10 and Figure 20.

Internal Pressure and Leakage Tests
The performance tests of the proposed TOBV were conducted using the experimental equipment, based on the regulation of KS B 2304:2001 (General rules for inspection of valves). The boiler and the temperature sensor controlled the temperature of 350 °C, and the pressure gauge checked the internal pressure. In the internal pressure test, the valve body was filled with the water at the half-shut position, and then the internal pressure of 1.7 MPa, which exceeded the working pressure (1 MPa) × 1.5, was maintained for 182 s. In the leakage test, the valve body was filled with the water at the fully shut-off position, and the internal pressure of 1.15 MPa, which exceeded the working pressure (1 MPa) × 1.1, was maintained for 125 s. The results of the flowrate measurement exhibit that there were no leakage and exudation in both tests, as shown in Table 10 and Figure 20.

Field Test of Proposed TOBV in CCPPs
The prototype of the proposed TOBV was installed at the gas turbine in the CCPPs and was operated for three months to evaluate durability and suitability for the actual working environment. Leakage could not be measured at the rear end of the turbine, so sealing performance of the new

Field Test of Proposed TOBV in CCPPs
The prototype of the proposed TOBV was installed at the gas turbine in the CCPPs and was operated for three months to evaluate durability and suitability for the actual working environment. Leakage could not be measured at the rear end of the turbine, so sealing performance of the new model was predicted by comparing its temperature with that of the conventional valve, as shown in Figure 21. Perfect sealing blocks the hot gas flow into the outlet pipe, which allows the temperature in the rear end to decrease. The measurement results show that the gas (350 • C), which entered into the inlet pipe during the disc opening, elevated sharply the temperature at the rear end from 27 • C to 100 • C. After then, the temperature decreased at fully shut-off position, and it kept in an almost steady state. The time to be required for opening and closing the disc was 4.25 s, and its cycle was one day. Temperature of the new model (32 • C) decreased by 20%, compared with that of the conventional one (40 • C), which demonstrates that the proposed laminated seal improved the sealing ability of the TOBV.
Appl. Sci. 2019, 9, x FOR PEER REVIEW 17 of 19 model was predicted by comparing its temperature with that of the conventional valve, as shown in Figure 21. Perfect sealing blocks the hot gas flow into the outlet pipe, which allows the temperature in the rear end to decrease. The measurement results show that the gas (350 °C), which entered into the inlet pipe during the disc opening, elevated sharply the temperature at the rear end from 27 °C to 100 °C. After then, the temperature decreased at fully shut-off position, and it kept in an almost steady state. The time to be required for opening and closing the disc was 4.25 s, and its cycle was one day. Temperature of the new model (32 °C) decreased by 20%, compared with that of the conventional one (40 °C), which demonstrates that the proposed laminated seal improved the sealing ability of the TOBV.

Conclusions
This study developed the laminated seal to enhance sealing performance and operability of the TOBV. Our conclusions are summarized as follows.
(1) When the TOBV has the inner diameter of the pipe: 200 mm, the diameter of disc: 190.1 mm, the outer diameter of seal: 218 mm, and the thickness of disc: 15 mm, the three eccentricities (offset 1: 2.5 mm, offset 2: 23 mm, offset 3: 10°) were specified based on the international regulation, and the slope angles of contact surface (66.1° (minimum angle)~86.1° (maximum angle)) for eliminating rubbing were calculated. (2) Thermal-structural coupled analyses were conducted for the different seal designs (the number of seal layers and graphite thicknesses), temperatures (25 °C and 350 °C), and slope angles (66.1° and 86.1°), and the maximum contact stresses' behaviors at the contact surfaces were investigated as follows: -It was observed that the increases of contact surface led to a reduction of the peak contact stress, so that sealing performance was improved with increasing temperature as a result of thermal expansion. -In respect to the graphite layers, leakages were expected for the all models at 25 °C (Gmin: 0.19~0.28 MPa, Gmax: 0.14~0.18 MPa), and there were small changes in contact stress by varying the amount of the graphite, whereas at 350 °C, the larger contact stresses (Gmin: 12.33~4.73 MPa, Gmax: 5.03~0.79 MPa) beyond the working pressure were observed, which means that the graphite showed a sealing ability at the high temperature rather than at the room temperature. -In respect to the stainless steel layers, the graphite did not withstand the contact force at 25 °C, so it was intensively applied on the stainless steel whose surface was decreased as the amount of the graphite was increased, and then the contact stresses at the stainless steel were increased (SSmin: 2.0 MPa → 6.3 MPa, SSmax: 1.2 MPa → 4.9 MPa). While at 350 °C, the contact force was distributed to the graphite, and the peak contact stress at the stainless steel was reduced (SSmin: 169.2 MPa → 74.13 MPa, SSmax: 129.9 MPa → 63.0 MPa) as the amount of graphite was increased.

Conclusions
This study developed the laminated seal to enhance sealing performance and operability of the TOBV. Our conclusions are summarized as follows.
(1) When the TOBV has the inner diameter of the pipe: 200 mm, the diameter of disc: 190.1 mm, the outer diameter of seal: 218 mm, and the thickness of disc: 15 mm, the three eccentricities (offset 1: 2.5 mm, offset 2: 23 mm, offset 3: 10 • ) were specified based on the international regulation, and the slope angles of contact surface (66.1 • (minimum angle)~86.1 • (maximum angle)) for eliminating rubbing were calculated. (2) Thermal-structural coupled analyses were conducted for the different seal designs (the number of seal layers and graphite thicknesses), temperatures (25 • C and 350 • C), and slope angles (66.1 • and 86.1 • ), and the maximum contact stresses' behaviors at the contact surfaces were investigated as follows: -It was observed that the increases of contact surface led to a reduction of the peak contact stress, so that sealing performance was improved with increasing temperature as a result of thermal expansion. -In respect to the graphite layers, leakages were expected for the all models at 25 • C (G min : 0.19~0.28 MPa, G max : 0.14~0.18 MPa), and there were small changes in contact stress by varying the amount of the graphite, whereas at 350 • C, the larger contact stresses (G min : 12.33~4.73 MPa, G max : 5.03~0.79 MPa) beyond the working pressure were observed, which means that the graphite showed a sealing ability at the high temperature rather than at the room temperature. -In respect to the stainless steel layers, the graphite did not withstand the contact force at 25 • C, so it was intensively applied on the stainless steel whose surface was decreased as the amount of the graphite was increased, and then the contact stresses at the stainless steel were increased (SSmin: 2.0 MPa → 6.3 MPa, SSmax: 1.2 MPa → 4.9 MPa). While at 350 • C, the contact force was distributed to the graphite, and the peak contact stress at