Numerical study of the dynamic characteristics of a single-layer graphite core in a thorium molten salt reactor

Numerical study of the dynamic characteristics of a single-layer graphite core in a thorium molten salt reactor 1 2 first-author 1 2 1 2 1 1 corresp Corresponding author. E-mail address: derektsang@sinap.ac.cn Corresponding author. E-mail address: derektsang@sinap.ac.cn derektsang@sinap.ac.cn Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Shanghai 201800, China University of Chinese Academy of Sciences, Beijing 100049, China University of Chinese Academy of Sciences, Beijing 100049, China A reactor core in a thorium molten salt reactor uses graphite as a moderator and reflector. The graphite core is a multi-layered arrangement of graphite bricks that are loosely connected to each other using a system of keys and dowels. Consequently, the graphite core is a type of discrete stack structure with highly nonlinear dynamic behavior. Hence, it is important to investigate the dynamic characteristics of the graphite core. In this study, a three-dimensional single-layer graphite core model, which is a part of the thorium molten salt reactor side reflector structure, was analyzed using the explicit method in ABAQUS 2016 to study the core dynamic behavior when subjected to different excitations. The design parameters, such as the diameter of the dowel, the gap between key and keyway and the bypass flow gap between two adjacent bricks, were also considered in this model. To reduce excessive demands on available computational resources considering the effect of molten salt, the spring-dashpot model was applied to model the interaction forces between the molten salt and graphite bricks. Numerical simulation results show that the effect of molten salt is a reduction in the peak maximal principal stress, and a larger gap between two bricks is beneficial to maintain the integrity of the graphite core under earthquake loading. The results obtained by the simulation can be used as a reference for future designs of a molten salt graphite core. Graphite core Dynamic behavior ABAQUS 1. Introduction 1 The Chinese Academy of Sciences (CAS) has launched a strategic priority research project to utilize thorium energy via the development of a thorium molten salt reactor (TMSR) to secure a long-term nuclear fuel supply by diversifying the sources of nuclear fuel [ 1 1 , 2 2 ]. The reactor core in a TMSR contains graphite components used as moderator and reflector. The graphite core is a multi-layered arrangement of graphite bricks that are loosely connected to each other using a system of keys and dowels. The graphite core is a type of discrete stack structure with highly nonlinear dynamic behavior and its behave completely differently in terms of dynamic response compared to a welded structure or bolted structure when subjected to seismic loading[ 3 3 , 4 4 ]. Hence, it is important to investigate the dynamic characteristic of the TMSR graphite core assembly. According to Institute of Electrical and Electronics Engineers Standard 344[ 5 5 ], during an earthquake the reactor core requires structural integrity to be maintained and safe shutdown of the reactor. The core is required to withstand two types of earthquakes, namely an operating basis earthquake (OBE) and a safe shutdown earthquake (SSE). An OBE is an earthquake that could reasonably be expected to occur at the site during the operating life of the reactor and the reactor is designed to remain functional. An SSE is an earthquake that is based upon an evaluation of the maximal earthquake potential considering the regional and local geology and seismology and the specific characteristics of the local subsurface material. An SSE produces a maximal vibratory ground motion for which the capability to shut down the reactor and maintain it in a safe shutdown condition is required. The fundamental safety requirements of the graphite core are to allow free movement of control rods and to direct the flow of coolant to ensure adequate cooling of the fuel and the core structure, both in normal and fault/hazard conditions. Consequently, a safety analysis demonstrating the tolerability of the core safety functions to geometrical changes and earthquakes is required. An understanding of the dynamic response of the graphite core and the internal stress of individual graphite components under normal operation and transient conditions is essential in assessing core geometry; structure integrity; inspection strategy; and, consequently, core safety. A numerical method plays a principal role in the dynamic analysis of the graphite core. One major advantage of this method is that it can provide detailed displacement and stress distribution information for individual components. It also can allow for more complex analysis and inaccessible experimental conditions. Reactors including a gas-cooled reactor, advanced gas-cooled reactor (AGR), high-temperature gas-cooled reactor (HTGR) and molten salt reactor have a graphite core. The dynamic analysis of a graphite core using a numerical method has occurred quite often during the past decades. For an HTGR, Merson and Bennett developed a small number of degrees-of-freedom computer program to study the seismic response of an HTGR core[ 6 6 ]. Ikushima and Honma contrived a simplified model for HTGR core seismic analysis to reduce the considerable computational cost. In this code, one degree-of-freedom mass with nonlinear column characteristics has been used instead of a stacked column[ 7 7 ]. Considering the pebble-bed HTGR has various graphite shapes, a generalized three-dimensional (3D) model with the associated computer code was developed to treat these interconnected graphite bricks with arbitrary shapes by Zeng and Wang [ 8 8 ]. In this model, each brick is treated as a rigid body with six degrees-of-freedom: three translational displacements and three rotations around the brick center of gravity. A nonlinear spring dashpot model has been applied to model the collision between adjacent bricks and the interaction forces through the key systems. Tsuji et al. developed a vertical two-dimensional analytical model using multi-purpose commercial finite element method (FEM) code, which considers the multiple impacts and friction between block interfaces and rocking motion in contact with dowel pins of the HTGR core using contact elements[ 9 9 ]. For an accurate prediction response of an AGR core to an earthquake, Ahmed et al. [ 10 10 ] proposed an analytical method with several hundred thousand degrees of freedom for an AGR reactor. The graphite bricks which comprise the core were idealized as rigid masses, whilst the contact spring elements were used to represent the load transmissions or impacts that could occur between the bricks. The graphite core of the AGR core has been seismically qualified using the non-linear computer program AGRCOR[ 11 11 , 12 12 ]. Graphite bricks have been treated as rigid bodies and each brick has three degrees of freedom: two translational displacements and one rotation. However, few works have addressed the dynamic behavior of the graphite core in a TMSR. The difference between a TMSR and other types of graphite-moderated reactors is the existence of molten salt, which is a vital factor influencing the dynamic response of the graphite core, because the mass density of molten salt is much higher than that of a gas coolant. In addition, the aforementioned literature used a rigid body and limited dimensions in numerical analysis. For a better understanding of the dynamic response of graphite cores and the internal stress of individual core components in a TMSR, graphite bricks cannot be treated as a rigid body and the model should be three-dimensional with full stress analysis. In this study, analysis was conducted using the FEM program ABAQUS 2016. A single layer of graphite core was built. Each brick was treated as a deformable body. Because the design of the TMSR graphite core brick is different from that of the other graphite-moderated reactors, several design parameters such as the diameter of the dowel, the gap between key and keyway and bypass flow gaps were also considered in this model to investigate how these parameters influence the dynamic responses and stresses within the graphite bricks. 2. Numerical study 1 2.1 Explicit dynamic analysis 2 The equation of equilibrium governing the dynamic response of a system in the FEM program at time <math id="M1"><mrow><msub><mi>t</mi><mi>n</mi></msub></mrow></math> is generally expressed as follows: <math display="block" id="M2"><mrow><mi>M</mi><msub><mover accent="true"><mi>u</mi><mo>¨</mo></mover><mi>n</mi></msub><mo>+</mo><mi>C</mi><msub><mover accent="true"><mi>u</mi><mo>¨</mo></mover><mi>n</mi></msub><mo>+</mo><msubsup><mi>F</mi><mi>n</mi><mrow><mtext>int</mtext></mrow></msubsup><mo stretchy="false">(</mo><msub><mi>u</mi><mi>n</mi></msub><mo stretchy="false">)</mo><mo>=</mo><msub><mi>R</mi><mi>n</mi></msub><mo>,</mo></mrow></math> where, <math id="M3"><mstyle mathvariant="bold" mathsize="normal"><mi>M</mi></mstyle></math> is the mass matrix and <math id="M4"><mrow><mstyle mathvariant="bold" mathsize="normal"><mi>C</mi></mstyle><mtext> </mtext></mrow></math> is the damping matrix; <math id="M5"><mrow><mtext> </mtext><mstyle mathvariant="bold" mathsize="normal"><mi>R</mi></mstyle></mrow></math> n is the vector of the externally applied loads; <math id="M6"><mrow><msubsup><mstyle mathvariant="bold" mathsize="normal"><mi>F</mi></mstyle><mi>n</mi><mrow><mtext>int</mtext></mrow></msubsup></mrow></math> is the internal force vector; and ü n , <math id="M7"><mover accent="true"><mstyle mathvariant="bold" mathsize="normal"><mtext>u</mtext></mstyle><mo>˙</mo></mover></math> n , u n are the acceleration, velocity and displacement vectors, respectively. Eq. (1) is integrated by an explicit central difference integration rule and can be written as follows: <math display="block" id="M8"><mrow><msub><mover accent="true"><mi>u</mi><mo>¨</mo></mover><mi>n</mi></msub><mo>=</mo><msup><mstyle mathvariant="bold" mathsize="normal"><mi>M</mi></mstyle><mrow><mo>−</mo><mn>1</mn></mrow></msup><mrow><mo>[</mo><mrow><msub><mstyle mathvariant="bold" mathsize="normal"><mi>R</mi></mstyle><mi>n</mi></msub><mo>−</mo><mi>C</mi><msub><mstyle mathvariant="bold" mathsize="normal"><mover accent="true"><mi>u</mi><mo>˙</mo></mover></mstyle><mi>n</mi></msub><mo>−</mo><msubsup><mstyle mathvariant="bold" mathsize="normal"><mi>F</mi></mstyle><mi>n</mi><mrow><mtext>int</mtext></mrow></msubsup><mrow><mo>(</mo><mrow><msub><mstyle mathvariant="bold" mathsize="normal"><mi>u</mi></mstyle><mi>n</mi></msub></mrow><mo>)</mo></mrow></mrow><mo>]</mo></mrow><mo>.</mo></mrow></math> Velocities and displacements are updated in each time step as follows: <math display="block" id="M9"><mrow><msub><mstyle mathvariant="bold" mathsize="normal"><mover accent="true"><mi>u</mi><mo>˙</mo></mover></mstyle><mrow><mi>n</mi><mo>+</mo><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msub><mo>=</mo><msub><mstyle mathvariant="bold" mathsize="normal"><mover accent="true"><mi>u</mi><mo>˙</mo></mover></mstyle><mrow><mi>n</mi><mo>−</mo><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msub><mo>+</mo><mi>Δ</mi><msub><mi>t</mi><mi>n</mi></msub><msub><mover accent="true"><mi>u</mi><mo>¨</mo></mover><mi>n</mi></msub><mo>,</mo></mrow></math> <math display="block" id="M10"><mrow><msub><mstyle mathvariant="bold" mathsize="normal"><mi>u</mi></mstyle><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><msub><mstyle mathvariant="bold" mathsize="normal"><mi>u</mi></mstyle><mi>n</mi></msub><mo>+</mo><mi>Δ</mi><msub><mi>t</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msub><msub><mstyle mathvariant="bold" mathsize="normal"><mover accent="true"><mi>u</mi><mo>˙</mo></mover></mstyle><mrow><mi>n</mi><mo>+</mo><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msub><mo>,</mo></mrow></math> where, <math id="M11"><mrow><mi>Δ</mi><msub><mi>t</mi><mi>n</mi></msub><mo>=</mo><mo stretchy="false">(</mo><mi>Δ</mi><msub><mi>t</mi><mi>n</mi></msub><mo>+</mo><mi>Δ</mi><msub><mi>t</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub><mo stretchy="false">)</mo><mo>/</mo><mn>2</mn></mrow></math> . The explicit integration scheme has a very important advantage in that it improves the computational efficiency by using a diagonal mass matrix because the inversion of the mass matrix used in Eq. (2) is trivial. Among the disadvantages of the explicit integration procedure is the stability of the numerical integration. To achieve stability, the calculation time step size must be smaller than the critical time step ∆ t cr , which is determined by the character length of the elements and their material properties. For a finite element mesh of constant strain and rate independent materials, the critical time step can be calculated as follows: <math display="block" id="M12"><mrow><mo> </mo><mi>Δ</mi><msub><mi>t</mi><mrow><mtext>cr</mtext></mrow></msub><mo>≤</mo><msub><mrow><mi>min</mi></mrow><mtext>e</mtext></msub><mo stretchy="false">(</mo><msub><mi>l</mi><mtext>e</mtext></msub><mo>/</mo><msub><mi>c</mi><mtext>e</mtext></msub><mo stretchy="false">)</mo><mo>,</mo></mrow></math> where, l e is the character length of the element and <math id="M13"><mrow><msub><mi>c</mi><mtext>e</mtext></msub></mrow></math> is the wave speed in the element[ 13 13 ]. 2.2 Fluid-structure coupling method 2 In TMSR, all the graphite components with narrow gaps between them are immersed in molten salt. It is difficult to construct the finite element model of the totally assembled reactor internals with molten salt using the coupled Eulerian-Lagrangian method because of the complex multibody interaction, boundary conditions and limitation on computational time and physical memory of the computer. The effects of fluid on the brick can be regarded as the stiffness <math id="M14"><mrow><mtext> k</mtext></mrow></math> and damping <math id="M15"><mrow><mtext>c </mtext></mrow></math> of the spring element[ 14 14 ]. Spring-like elastic behavior and dashpot-like (damping) behavior are used in the model to consider fluid–structure interaction (FSI). The graphite brick is subjected to sinusoidal harmonic <math id="M16"><mrow><mtext>P</mtext><mrow><mo>(</mo><mtext>t</mtext><mo>)</mo></mrow></mrow></math> . The effect of the fluid–structure interaction between the two bricks can be simplified as an added mass, added damping and added stiffness coefficient as shown in Fig. 1 Fig. 1 . In this model, the graphite brick is supported by six parallel springs. If the effect of the molten salt is negligible, the free vibration equation of the brick is as follows[ 15 15 , 16 16 ]: <math display="block" id="M17"><mrow><mi>m</mi><mover accent="true"><mi>y</mi><mo>¨</mo></mover><mo>+</mo><mi>K</mi><mi>y</mi><mo>=</mo><mn>0</mn></mrow></math> where, m and K are the mass of the brick and the stiffness of the spring, respectively; <math id="M18"><mi>y</mi></math> is the displacement of the brick; and <math id="M19"><mover accent="true"><mi>y</mi><mo>¨</mo></mover></math> is the acceleration of the brick. When the brick vibrates in the fluid, the added damping and resistance force are produced. The vibration equation of the brick in the fluid is written as follows: <math display="block" id="M20"><mrow><mo> </mo><mrow><mo>(</mo><mrow><mi>m</mi><mo>+</mo><msub><mi>m</mi><mrow><mtext>ad</mtext></mrow></msub></mrow><mo>)</mo></mrow><mover accent="true"><mi>y</mi><mo>¨</mo></mover><mo>+</mo><msub><mi>C</mi><mrow><mtext>ad</mtext></mrow></msub><mover accent="true"><mi>y</mi><mo>˙</mo></mover><mo>+</mo><mrow><mo>(</mo><mrow><mi>K</mi><mo>+</mo><msub><mi>k</mi><mrow><mtext>ad</mtext></mrow></msub></mrow><mo>)</mo></mrow><mi>y</mi><mo>=</mo><mn>0</mn><mo>,</mo></mrow></math> where, C ad the added damping coefficient caused by the fluid; m ad is the added mass caused by the fluid; and k ad is the added stiffness coefficient caused by the fluid. Thus, the spring-dashpot is applied to present the interaction forces between the molten salt and graphite to reduce excessive demands on available computational resources. 3. Structural details and finite element model 1 3.1 TMSR graphite core 2 The design of the TMSR graphite core is shown in Fig. 2 Fig. 2 . Graphite serves as an internal structural material as well as a moderator in the TMSR. A large number of graphite bricks are involved in the TMSR core, forming reflectors. The reflectors, which consist of a top, side, and bottom reflector, together with the core support structure, form an inner space for fuel elements and at the same time serve as a flow guidance structure for the molten salt. 3.2 3D finite element model 2 A single-layer side reflector model of the TMSR core is shown in Fig. 3 Fig. 3 . The single-layer model is a typical unit of the side reflector core structure of the TMSR reactor, in which the dynamic response of an individual component and the interaction between bricks can be investigated. The model is composed of 16 graphite bricks, 16 keys and a plate. The bricks are all horizontally constrained by dowels on the plate. Each brick is interconnected with a key to restrain the circumferential relative motion. The whole model is mounted on a plate that generates the dynamic motion. As a simplification, a 3D finite element model without side metallic component support, which is part of the TMSR side reflector structure, was developed using the ABAQUS 2016 software according to the parameters of the TMSR. The main parameters of the graphite material are listed in Table 1 Table 1 . To improve the accuracy of the model, the mesh on the corner and channel designs was refined as shown in Fig. 3 Fig. 3 . Some descriptions regarding this model are shown in Table 2 Table 2 . 1 1 Parameters 1 1 Values 1 1 Density( g / cm 3 ) 1 1 1.85 1 1 Young’s modulus (GPa) 1 1 10 1 1 Poisson’s ratio 1 1 0.21 1 1 Compressive strength (MPa) 1 1 85.3 1 1 Tensile strength (MPa) 1 1 41.5 a The C3D8R element is 3D 8-node element, with reduced integration (one integration point). 1 1   2 1 Description 1 1 Components 2 1 16 bricks (173312 elements) 1 1   1 1   1 1 16 keys (16560 elements) 1 1   1 1   1 1 1 plate (50624 elements) 1 1 Type of element 2 1 C3D8R a 1 1 Contact 2 1 All with self 1 1 Procedure 2 1 Dynamic, explicit 1 1 Output 2 1 Whole model 1 1 Time 2 1 20 periods of every job 3.3 Boundary conditions 2 In this study, the aim was to investigate the dynamic behavior and response characteristics of the integrated graphite components. The seismic excitation applied on this model included the following: (case a) the excitation amplitude was from 3 to 20 mm in the X direction; (case b) the frequency ranges were from 2 to 10 Hz. During TMSR operation, the dimension of the graphite bricks will change under high temperature and fast neutron irradiation[ 17 17 , 18 18 ]. When the bricks contract, the gaps between the graphite bricks will increase, which will affect the dynamic response of the whole core structure. Consequently, the following cases were considered in the analyses: (case c) the gap between the dowel and brick, (case d) the gap between the key and keyway and (case e) the bypass flow gap between the two bricks. Another crucial factor that also needed to be considered was the effect of the molten salt in the TMSR (case f). All load and condition details are listed in Table 3 Table 3 . 1 1 Cases 1 1 Values 1 1 (a) Different input amplitude (mm) 1 1 3,5,10,15,20 1 1 (b) Different input frequency (Hz) 1 1 2,4,6,8,10 1 1 (c) Gap between dowel and brick (mm) 1 1 0.5,1,1.5,2 1 1 (d) Gap between key and keyway (mm) 1 1 0.5,1,1.5,3 1 1 (e) Gap between two bricks (mm) 1 1 0.5,1,3,5 1 1 (f) Effect of molten salt 1 1 - 3.4 Mesh sensitivity study 2 It is well known that the accuracy of numerical results is strongly dependent on the mesh density used. Better results would have been obtained by increasing the total number of elements, but the computational time and required computer memory would increase substantially as well. The 3D model with full degrees-of-freedom requires larger computational power and memory. Because of the limitations of the computer and software, it was impossible to use the largest mesh density in the analyses. Thus, an acceptable mesh density needed to be identified. Based on the numerical results, the peak maximal principal stress error of the total element numbers 140K, 240K and 340K (as shown in Fig. 4 Fig. 4 ) are 28.84 <math id="M21"><mtext>%</mtext></math> and 108.48 <math id="M22"><mtext>%</mtext></math> at 0.5s, 8.42 <math id="M23"><mtext>%</mtext></math> and 128.14 <math id="M24"><mtext>%</mtext></math> at 1.0s, and 7.05 <math id="M25"><mtext>%</mtext></math> and 149.17 <math id="M26"><mtext>%</mtext></math> at 1.5s. The errors in the peak maximal principal stress of different mesh numbers are shown in Table 4 Table 4 . Therefore, the total number of elements of 240K was used for all the analyses in this study. 1 1   1 1 Total element number 3 1 Peak maximal principal stress of model 1 1   1 1   1 1 T1 1 1 T2 1 1 T3 1 1 Error (%) 1 1 140 K 1 1 108.48 (%) 1 1 128.14(%) 1 1 149.17(%) 1 1   1 1 240 K 1 1 28.84 (%) 1 1 8.42(%) 1 1 7.05(%) 1 1   1 1 340 K 1 1 - 1 1 - 1 1 - 4. Results and discussion 1 4.1 Brick maximal principal stress location 2 The graphite bricks were numbered from 01 to 16 in an anticlockwise direction with the 01 brick at 0 degrees (+ Y -axis) as shown in Fig. 5a Fig. 5a . The motion of the plate was assumed to be a sine wave with a forcing frequency 10Hz and an amplitude of excitation of 5mm. The maximal principal stress of the brick and its location were recorded as shown in Fig. 5b Fig. 5b . It can be seen from the results that the location of the maximal principal stress always appears on 01, 08, 09 and 16. Compared to all the maximal principal stresses during the analysis process, the peak maximal principal stress (10.78 MPa) appears on number 08. Thus, further research objectives should be concentrated on the Number 08 brick. 4.2 Effect of excitation 2 4.2.1 Different amplitude 3 The displacements of the graphite bricks may affect the insertion of the control rods and shutdown material and consequently the reactor’s shut down. Therefore, it is of vital importance, from the viewpoint of seismic design, to investigate the effect of different amplitudes of excitation. The vibration was excited by a sinusoidal wave whose frequency was a fixed value in the X direction. The value of the amplitude was a variable parameter. It was set to be 3 mm, 5 mm, 10 mm, 15 mm and 20 mm respectively. Different amplitude time histories of brick: 08-point P in the Z direction were recorded as shown in Fig. 6 Fig. 6 . The obtained results indicated that the displacement of the point in the Z direction during a certain frequency was sensitive to the oscillation amplitude, and collapse of the layer occurred when the excitation amplitude was greater than 10 mm. As shown in Fig. 6 Fig. 6 , the largest amplitude of the total structure of the graphite core is less than 10 mm, which has the smallest peak displacement response as compared to the other amplitude situations. These load and boundary conditions are extremely conservative compared to the actual seismic motion. For the graphite core, side metallic component support was ignored, which reduced the total dynamic response of the model. Therefore, such a collapse is not expected in the following seismic experiments and in the actual reactor core. According to the figure, the whole single-layer model maintained its basic integrity and no graphite bricks were found to lose their initial position when the excitations were 3 mm and 5 mm. The No.08 brick time histories of a maximal principal stress of 3 mm and 5 mm were recorded as shown in Fig. 7 Fig. 7 . It is clearly observed that the greater the excitation amplitude, the greater the maximal principal stress in brick No.08. This result means the maximal principal stress also has a strong relationship with excitation amplitude. 4.2.2 Different frequencies 3 The excitation frequency also has a significant influence on the dynamic analysis results. Thus, the effect of different values of frequency was investigated as described in this section. The frequency of an earthquake is usually less than 10 Hz; thus, five different excitation frequencies (2 Hz, 4 Hz, 6 Hz, 8 Hz and 10 Hz) on the dynamic responses of the model are represented in Fig. 8 Fig. 8 . It is observed from the figure that the greater the excitation frequency, the greater the value of the maximal principal stress in brick No.08 and the gap between No.08 and the whole model bricks increases as the excitation frequency increases. In summary, the dynamic response of the model has a significant relationship with the excitation frequency. 4.3 Effect of design parameters 2 During TMSR operation, the dimensions of the graphite bricks will change under fast neutron irradiation. The graphite material will contract then turn around and return to the zero point, which is considered as the life end of the components. When the bricks contract, the gap between brick and dowel, the gap between the key and keyway and the bypass flow gap between two bricks will increase, which will affect the dynamic response of the whole core structure. 4.3.1 Gap between key and brick 3 The graphite brick nonlinear dynamic response is mainly caused by the interaction between the bricks and dowel. Different values of dowel diameter, namely the gap between the dowel and brick, are a vital factor. It was set to be 0.5 mm, 1 mm, 1.5 mm and 2 mm respectively in this section. The motion of the plate was assumed to be a sine wave with a forcing frequency of 10 Hz and an amplitude of excitation of 5 mm. The maximal principal stress of No.08 and the maximal displacement of point P in the X direction were recorded as shown in Fig. 9 Fig. 9 . As is shown, the larger the gap between dowel and brick, the greater the maximal principal stress and absolute displacement of P. The reason for this phenomenon is because the gap affects the displacement of the brick, which makes the impact forces between the graphite bricks increase. Consequently, the gap should be as small as possible to reduce the collision between two bricks to ensure the structural integrity of the core. 4.3.2 Gap between the key and keyway 3 The whole graphite core consists of a large number of graphite components interconnected by keys. Keys are employed to restrain the circumferential relative motion to maintain the integrity of the graphite core. Thus, the gap between the key and keyway is also a vital factor in dynamic analysis. The plate was excited by a sinusoidal wave at a frequency of 10 Hz and an amplitude of 5 mm. The gap between the key and keyway was a study parameter. It was set at 0.5 mm, 1 mm, 1.5 mm and 3mm. The key under consideration was that between brick:08 and brick:09. The maximal principal stress and absolute displacement in the X , Y and Z directions were recorded as shown in Fig. 10 Fig. 10 . It can be seen from Fig. 10 Fig. 10 that the peak maximal displacement in the X , Y and Z directions at P increases along with the gap between the key and keyway and the maximal principal stress of key also showed this trend. The reason for these phenomena may be because the larger gap provides more room for the displacement of the key, which means there is less restraint on the key and the impact forces between the graphite brick and key increase. We also found that the displacement in the Z direction increased rapidly compared to that of the X and Y directions. This was because the movement of the brick in the X and Y directions was not only confined by the gap between the key and keyway, but also the gap between the two bricks. When the gap between the key and keyway was 3 mm, the 1-mm gap between the two bricks had less of an effect on movement of the key. In the Z direction, the larger gap confined movement of the key less. Thus, the displacement in the Z direction increased rapidly. 4.3.3 Bypass flow gap 3 In this section, four different gaps were considered and compared to each other to obtain the dynamic response considering the gap effect between the two bricks. The plate was excited by a sinusoidal wave at a frequency of 10 Hz and amplitude of 5 mm. The gap between the two bricks was a study parameter. It was set to be 0.5 mm, 1 mm, 3 mm and 5 mm, respectively. The maximal principal stress and displacement in the X , Y and Z directions of P were recorded as shown as Fig. 11 Fig. 11 . The results showed that the gap between the two bricks had little influence on the displacement of point P in the three directions. The responsible factor for the slight variation in the displacement response at different gaps may be that the movement of the brick was also restricted by the gaps between the key and keyway. As previously discussed, the movement of the brick was confined not only by the bypass flow gap but also the gap between the key and keyway. The gap between the key and keyway was 0.5mm in this section. When the bypass flow gap was larger than 0.5 mm, the movement of the brick was restrained by the gap between the key and keyway. As for the peak maximal principal stress of brick: 08, the movement of the brick had strong a relationship with the stress. The bypass flow gap had little influence on the movement of the brick. Thus, the peak maximal principal stress also changed slightly. 4.4 Effect of molten salt 2 When a structural component is subjected to dynamic motion, the viscous fluid produces reaction forces that are generally referred to as added mass and added damping forces[ 19 19 - 21 21 19-21 ]. The added mass is associated with a mass of fluid that is accelerated by the structure and added damping is related to the shape of the structure and the viscosity of the fluid [ 22 22 ]. All the graphite components with a narrow gap between them are immersed in a TMSR and its mass density is much greater than that of gas coolant as in gas-cooled reactors. Therefore, the FSI between the molten salt and graphite components cannot be ignored when subjected to an excitation load. To reduce excessive demands on computational resources, such as computational time and computer physical memory, the spring-dashpot method was applied to represent the interaction forces between the molten salt and graphite bricks [ 8 8 ]. The dynamic characteristics of added mass, added damping and stiffness can be obtained by fitting the time history of vibration experiments with the analytical solution. Thus, the proper stiffness of the spring and the damping coefficient in the equivalent model can be obtained as shown in Table 5 Table 5 [ 23 23 ]. Then, the equivalent model can be used in the seismic analysis of the TMSR graphite core as shown as Fig. 12 Fig. 12 . The vibration was excited by a sinusoidal wave at a frequency of 10Hz and an amplitude of 5mm in this section. The results of the displacement of point P in the three directions ( X , Y and Z ) and the peak maximal principal stress were recorded around brick:08 as shown in Table 6 Table 6 . As indicated, the existence of molten salt can reduce the peak maximum principal stress of brick: 08. In addition, the displacement of X and Y at point P is different from that of Z . The peak maximal displacement of P in the X and Y directions increased when the effect of molten salt was added, but the peak maximal displacement in the Z direction decreased considering the effect of the liquid. The phenomenon of the response of the graphite considering the FSI effects as concluded using the spring-dashpot model is that it only acts in the X and Y directions when the structure is subjected to excitation. In the horizontal Y direction, there is no restriction and the there is a stronger effect in the X and Y directions making movement of the brick more possible. Thus, the larger bypass flow gap is the best choice for a design reducing the dynamic behavior under earthquake loading. This is because there is more molten salt in this gap, thus there is a stronger damping effect for a dynamic response of graphite bricks and maintain the integrity of the graphite core. 1 1 Group 1 1 Added mass (kg) 1 1 Added damping 1 1 Stiffness  ( N/m )  1 1 3 mm 1 1 9.76 1 1 833.40 1 1 39492 1 1 5 mm 1 1 8.41 1 1 755.59 1 1 37776 1 1 10 mm 1 1 6.65 1 1 540.25 1 1 32140 1 1 15 mm 1 1 4.87 1 1 455.38 1 1 29839 1 1 20 mm 1 1 4.51 1 1 432.80 1 1 29147 1 1   1 1 Air 1 1 Liquid 1 1 X (mm) 1 1 7.12 1 1 6.80 1 1 Y (mm) 1 1 1.25 1 1 1.05 1 1 Z (mm) 1 1 3.80 1 1 4.00 1 1 Max peak stress (MPa) 1 1 10.14 1 1 9.607 5. Conclusion 1 A single-layer model was built based on a TMSR graphite core structure side reflector using commercial software ABAQUS 2016. The overall objectives of the present study were to numerically investigate the dynamic response and determine the seismic response of the graphite core. The dynamic responses under different excitations were obtained. Some conclusions can be drawn based on the numerical results of the maximal principal stress and displacement: (1) The location of the maximal principal stress always appears on bricks 01, 08, 09 and 16 when subjected to excitation in the X direction. (2) The whole single-layer model cannot maintain its basic structural integrity when the excitation amplitude is greater than 10mm without side metallic component support and the maximal principal stress on brick No.08 increases with an increase in amplitude and excitation frequency. (3) The movement of a brick is restrained by two factors namely the bypass flow gap and the gap between the key and keyway. For a certain gap of key and keyway, the bypass flow gap has a slight influence on displacement in the three directions and the maximal principal stress. (4) The effect of molten salt can reduce the peak maximal principal stress, thus a larger bypass flow gap is beneficial to maintain the integrity of the graphite core under earthquake loading. References J. Serp , M. Allibert , O. Beneš et al ., The molten salt reactor (MSR) in generation IV: overview and perspectives . Pro. Nucl. Energ . 77 , 308 - 19 ( 2014 ). doi: 10.1016/j.pnucene.2014.02.014 http://doi.org/10.1016/j.pnucene.2014.02.014 The molten salt reactor (MSR) in generation IV: overview and perspectives Z. Dai , Thorium molten salt reactor nuclear energy system (TMSR) . Molten Salt Reactors and Thorium Energy . p, 531 - 40 ( 2017 ). doi: 10.1016/b978-0-08-101126-3.00017-8 http://doi.org/10.1016/b978-0-08-101126-3.00017-8 Thorium molten salt reactor nuclear energy system (TMSR) L. Sun , L. Shi , H. G. Wan , et al ., Seismic study on double-layer model of HTR graphite core structure . At. Energy. Sci. Technol . 46 , 839 - 44 ( 2012 ). doi: 10.1115/icone20-power2012-54961 http://doi.org/10.1115/icone20-power2012-54961 Seismic study on double-layer model of HTR graphite core structure W. Laframboise , T. Desmond , Nonlinear seismic analysis of a graphite reactor core . Nucl. Eng. Des . 110 , 131 - 46 ( 1988 ). doi: 10.1016/0029-5493(88)90013-1 http://doi.org/10.1016/0029-5493(88)90013-1 Nonlinear seismic analysis of a graphite reactor core Stations NPG . IEEE Recommended Practice for Seismic Qualification of Class 1E Equipment for Nuclear Power Generating Stations . doi: 10.1109/ieeestd.2005.96099 http://doi.org/10.1109/ieeestd.2005.96099 IEEE Recommended Practice for Seismic Qualification of Class 1E Equipment for Nuclear Power Generating Stations J. Merson , J. Bennett , Computer method for analyzing HTGR core block response to seismic excitation . Los Alamos Scientific Lab ., NM (USA) , ( 1976 ). doi: 10.2172/7334839 http://doi.org/10.2172/7334839 Computer method for analyzing HTGR core block response to seismic excitation T. Ikushima , T. Honma , Simplified Analytical Model for HTGR Core Seismic Response . J. Nucl. Sci. Technol . 16 , 605 - 12 ( 1979 ). doi: 10.1080/18811248.1979.9730952 http://doi.org/10.1080/18811248.1979.9730952 Simplified Analytical Model for HTGR Core Seismic Response C. Zeng , H. Wang , A Three-Dimensional Rigid-Body Model for Seismic Analysis of the Pebble-Bed HTR Graphite Core Structure . Publisher ( 2013 ). doi: 10.1115/icone21-15896 http://doi.org/10.1115/icone21-15896 A Three-Dimensional Rigid-Body Model for Seismic Analysis of the Pebble-Bed HTR Graphite Core Structure N. Tsuji , K. Ohashi , Development of seismic analysis model for HTGR core on commercial FEM code . ICONE . ( 2015 ). doi: 10.1299/jsmeicone.2015.23._icone23-1_105 http://doi.org/10.1299/jsmeicone.2015.23._icone23-1_105 Development of seismic analysis model for HTGR core on commercial FEM code K. Ahmed , J. Parker , D. Proffitt , Seismic response of the advanced gas cooled reactor core . Nucl. Eng. Des . 94 , 67 - 92 ( 1986 ). doi: 10.1016/0029-5493(86)90154-8 http://doi.org/10.1016/0029-5493(86)90154-8 Seismic response of the advanced gas cooled reactor core K. Ahmed , The dynamic response of multi-layers AGR core brick arrays . Nucl. Eng. Des . 104 , 1 - 66 ( 1987 ). doi: 10.1016/0029-5493(87)90303-7 http://doi.org/10.1016/0029-5493(87)90303-7 The dynamic response of multi-layers AGR core brick arrays K. Ahmed , S. Stojko , The non-linear seismic response of AGR core graphite brick slices correlation of experimental and analytical results . Earthq. Eng. Struc. D . 15 , 159 - 88 ( 1987 ). doi: 10.1002/eqe.4290150203 http://doi.org/10.1002/eqe.4290150203 The non-linear seismic response of AGR core graphite brick slices correlation of experimental and analytical results W. Rust , Non-linear finite element analysis in structural mechanics : Springer ; 2016 . doi: 10.1007/978-3-319-13380-5 http://doi.org/10.1007/978-3-319-13380-5 Non-linear finite element analysis in structural mechanics Abaqus V. 2016 Documentation. Dassault Systemes Simulia Corporation . 2016 . http://abaqus.software.polimi.it/v2016/books/hhp/default.htm?startat=pt02ch02s02.html http://abaqus.software.polimi.it/v2016/books/hhp/default.htm?startat=pt02ch02s02.html 2016 Documentation. Dassault Systemes Simulia Corporation L. Lu , Y. R. Yang , P. Li , et al ., Added mass, added stiffness and added damping coefficients for a parallel plate-type structure . Appl. Mech. Mater . 66 - 68 ( 2011 ). doi: 10.4028/www.scientific.net/AMM.66-68.1738 http://doi.org/10.4028/www.scientific.net/AMM.66-68.1738 Added mass, added stiffness and added damping coefficients for a parallel plate-type structure L. Sargentini , B. Cariteau , M. Angelucci , Experimental and Numerical Analysis for Fluid-Structure Interaction for an Enclosed Hexagonal Assembly . AMSE ( 2014 ). doi: 10.1115/pvp2014-28053 http://doi.org/10.1115/pvp2014-28053 Experimental and Numerical Analysis for Fluid-Structure Interaction for an Enclosed Hexagonal Assembly R. Nightingale , Radiation Effects in Graphite . Nucl. Sci. Eng . 26 , 286 - 7 ( 1966 ). doi: 10.2172/4785225 http://doi.org/10.2172/4785225 Radiation Effects in Graphite D. Tsang , B. Marsden , J. Vreeling et al ., Analyses of a restrained growth graphite irradiation creep experiment . Nucl. Eng. Des . 238 , 3026 - 30 ( 2008 ). doi: 10.1016/j.nucengdes.2007.12.017 http://doi.org/10.1016/j.nucengdes.2007.12.017 Analyses of a restrained growth graphite irradiation creep experiment T. Belytschko , Fluid-structure interaction . Comput. Struct . 12 , 459 - 69 ( 1980 ). doi: 10.1016/0045-7949(80)90121-2 http://doi.org/10.1016/0045-7949(80)90121-2 Fluid-structure interaction J. A. Jendrzejczyk , Added mass and damping of a vibrating rod in confined viscous fluids . J. Appl. Mech . 325 ( 1976 ). doi: 10.1115/1.3423833 http://doi.org/10.1115/1.3423833 Added mass and damping of a vibrating rod in confined viscous fluids W. D. Wilson , Added mass and damping coefficients fora hexagonal cylinder . J. Fluid. Struc . 5 , 503 - 19 ( 1991 ). doi: 10.1016/s0889-9746(05)80003-9 http://doi.org/10.1016/s0889-9746(05)80003-9 Added mass and damping coefficients fora hexagonal cylinder R. Fritz , The effect of liquids on the dynamic motions of immersed solids . J. Eng. Ind . 94 , 167 - 73 ( 1972 ). doi: 10.1115/1.3428107 http://doi.org/10.1115/1.3428107 The effect of liquids on the dynamic motions of immersed solids Y. Zhong , X. Yang , D. Ding , et al . Dynamic characteristics identification of two graphite bricks in molten salt reactor considering fluid-structure interaction . Nucl. Eng. Des . 335 , 409 - 416 ( 2018 ). doi: 10.1016/j.nucengdes.2018.06.015 http://doi.org/10.1016/j.nucengdes.2018.06.015 Dynamic characteristics identification of two graphite bricks in molten salt reactor considering fluid-structure interaction 10.1007/s41365-018-0488-8 This work was supported by the "Hundred Talent Program" of the Chinese Academy of Sciences, the Ministry of Human Resources and Social Security (No. Y419016031), and the Strategic Priority Research Program of Chinese Academy of Science (No. XDA02040100). 2018-01-04 2018-04-10 2018-04-11 2018-08-29 2018-10-01 2021-04-11 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2018 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-018-0488-8 http://dx.doi.org/10.1007/s41365-018-0488-8 . 2018 Nuclear Science and Techniques 10 29