Mechanical Integrity Assessment of Two-Side Etched Type Printed Circuit Heat Exchanger with Additional Elliptical Channel

Printed circuit heat exchangers (PCHEs) are often subject to high pressure and temperature difference between the hot and cold channels which may cause a mechanical integrity problem. A conventional plate heat exchanger where the channel geometries are semi-circular and etched at one side of the stacked plate is a common design in the market. However, the sharp edge tip channel may cause high stress intensity. Double-faced type PCHE appears with the promising ability to reduce the stress intensity and stress concentration factor. Finite element analysis simulation has been conducted to observe the mechanical integrity of double-etched printed circuit heat exchanger design. The application of an additional ellipse upper channel helps the stress intensity decrease in the proposed PCHE channel. Five different cases were simulated in this study. The simulation shows that the stress intensity was reduced up to 24% with the increase in additional elliptical channel radius. Besides that, the horizontal offset channels configuration was also investigated in this study. Simulation results show that the maximum stress intensity of 2.5 mm offset configuration is 9% lower compared to the maximum stress intensity of 0 mm offset. This work proposed an additional elliptical upper channel with a 2.5 mm offset configuration as an optimum design.


Introduction
A heat exchanger is a device that can exchange heat between two fluids at different temperatures without mixing [1,2]. In industry, the heat exchanger plays a critical role especially in heat-based power such as nuclear power. Practically, shell-and-tube heat exchanger is the most popular in industry. However, the printed circuit heat exchanger (PCHE) is now widely developed due to several promising aspects, such as a large area density and good pressure and temperature capabilities [3][4][5]. The development of highly efficient PCHE is important, which is strongly related to economy and safety aspects. In many cases, PCHE may be working at a high temperature and pressure [6]. This means that safety is an important feature of PCHE because some critical issues may occur under high temperature and pressure conditions. For example, Rakesh and Anand found thermal stress to be more dominant than pressure cycle stress due to the higher temperature gradient [7].
Many studies have been conducted to investigate thermal-hydraulic performance of PCHE. Mylavarapu et al. studied the thermal-hydraulic performance of PCHE which contains helium as the fluid [8]. The thermal hydraulic performance of a PCHE which has a longitudinal corrugation flow channel was studied by Kim et al. [9]. Sung and Lee conducted a study about the tangled channel heat PCHE for thermoelectric power generation [10]. These studies show that the development of a cross-section shape may be able to enhance the thermal-hydraulic performance.
Besides the thermal-hydraulic performance study, mechanical integrity study is also important in PCHE development. Lee and Lee [11] conducted a study about the structural assessment of a sodium and SCO2 zigzag channel by calculating the temperature and stress value. They found that the maximum stresses occurred at the channel tips which are the sharp corner of the semi-circular flow cross-section shape. Song et al. also conducted a study about the structural integrity evaluation of a lab-scale intermediate PCHE of very-high temperature reactor (VHTR). Under the test condition, the maximum Tresca stress was far below the allowable stress limit [12]. Armanto and Lee conducted a study about the mechanical integrity analysis of a PCHE with channel misalignment. This study also shows high-stress intensity at the tip edge of PCHE channels [13]. The phenomenon where the stress is shown high at the tip has also been observed at the nanoscale. Winter et al. [14] conducted a study which shows that the different pore shape and pattern will cause stress accumulation inside the structure and may lead to failure. Vo et al. [15] investigated the impact of nanopores and porosity on the mechanical properties of amorphous silica (a-SiO 2 ). Simulation results shows that high stress concentration occurs at the top and bottom of the pore.
A double-faced heat exchanger comes with a promising ability. Unlike the conventional plate heat exchanger where the channel geometries are semi-circle and etched at one side of the stacked plate, this double-faced PCHE has two faces on one plate. During the fabrication process, the metal stacked layer will be etched on two faces and after that, they will be stacked together through the diffusion bonding process. The shape of two-etched sides can be different or like each other. Usually, in double-faced PCHE, one side is the main channel while another face can be an additional upper channel. One of the reasons for double-faced application is to enhance the heat transfer by enlarging the fluid flow cross-section area. By considering an additional upper channel by etching two sides of one stacking plate, the flow cross-section area can be larger. Another reason is that the double-faced PCHE can reduce the friction as mentioned by Lee et al. [16].
There are several studies conducted about double-etched circular channel PCHE. Ma et al. [17] considered a double-etched PCHE at the high temperature of 900 • C and a high pressure above 7 MPa. One side of the stacked plate is considered as a semi-circle flow channel while another side is as an elliptical additional upper channel. A thermo-hydraulic performance of the proposed PCHE channel design was studied numerically. Based on the study, it was found that the additional elliptical upper channel can improve the heat transfer of PCHE and increase the pressure drop simultaneously. It was found that the Nusselt number was increased by 7−17% with a 0.234 mm additional upper channel, while the fanning factor increased by 3−5%. Lee and Kim [18] investigated the performance of a zigzag printed circuit heat exchanger with various channel shapes. Four shapes of channel cross-section were investigated such as semi-circle, rectangle, trapezoid, and circle. The investigation aims to see the global Nusselt number, Colburn j-factor, effectiveness, and friction factor. This study showed that a rectangular channel provides the best thermal performance while the circular channel has the worst thermal performance.
Lee and Kim [19] conducted a study about the thermal performance of a double-faced PCHE with thin plates. This study considered 2 mm thin plates inserted between the hot and cold channel. A thin plate is used to increase the overlapped area between the cold and hot channels which can improve the thermal performance of PCHE. This proposed geometry was studied under Reynold's number in cold channels ranging from 67,000 to 280,000. The result was compared to a reference PCHE design channel which has only one etched side. It was claimed that the new PCHE design can maximize the heat transfer by minimizing the distance of heat conduction and maximizing the overlapping area between the hot and cold channels. Hence, the proposed design shows a remarkably higher thermal performance than the reference design.
The present study examined the capability of double-faced PCHE design from a mechanical integrity point of view since all double-faced PCHE previous studies claimed that the proposed design can enhance the thermal-hydraulic performance. However, considering PCHE for nuclear power plant applications should include a mechanical integrity assessment which is related to the safety concern.
Energies 2020, 13, 4711 3 of 17 This study examined the additional elliptical upper channel effect in PCHE design and compared it to the design specification by considering the simplified design by rule assessment. Since the experiment was almost impractical to be conducted due to the tiny scale of geometry, the finite element method study was chosen to be conducted. The design condition used for the assessment was the PHCE design condition of sodium fast-cooled reactor (SFR) which has a 20 MPa and 0.5 MPa pressure of the cold and hot channel, respectively as shown in Table 1 [20]. The system will work under 525 • C maximum working environment. Two-dimensional finite element method (FEM) simulation was conducted to observe the stress distribution, the weakest area based on deformation analysis and construction code compliance including primary stress, primary and secondary stress, and yield stress assessment. Even though this model of geometry is not usual and still needs many supporting technology developments, especially for the assembly and etching processes, this study is important and innovative to prepare the more advanced PCHE design and manufacturing process in the future. In this study, the American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code was used. The obtained stress distributions were compared with the design criteria as stated in the construction code. ASME boiler and pressure vessel code Section III is the rule for the construction of nuclear facility components [21]. It consists of Division 1, Division 2, Division 3, and Division 5. These four divisions are broken down into several subsections and several subsections will be broken down into several subparts. Since the design condition used in this study is the sodium-fast cooled reactor of which the maximum temperature reached 538 • C, ASME Section III Division 5 Subsection HB, Subpart B, was chosen as the criteria. ASME Section III Division 5 Subsection HB Subpart B contains rules for the construction of nuclear facility components of high temperature reactors specified for class A metallic pressure boundary components under elevated temperature service.

Design and Approach
PCHE design is mostly influenced by the common semi-circular channel cross-section. For instance, PCHEs which have been commercially manufactured by Heatric typically have a semi-circular flow cross-section. A new cross-section design based on combined a semi-circular and elliptical shape is proposed. Even though this model of geometry is not usual and still needs many supporting technology developments, especially for the assembly and etching processes, this study is important and innovative to prepare the more advanced PCHE design and manufacturing processes in the future. This study aims to observe the effect of an additional elliptical shape channel from a mechanical integrity point of view after the analysis of thermal-hydraulic performance has been conducted for many studies. The common flow cross-section of the PCHE is shown in Figure 1a.
As per the investigation result of Lee and Lee [11], the maximum stresses occurred at the channel tips. These channel tips are formed during the stacking process of the semi-circular side of the plate and the flat side of the plate. To reduce this sharp corner, an additional elliptical upper channel is proposed to reduce the stress intensity of this flow cross-section design. This elliptical is located above the semi-circular channel. This change in design causes one layer of the stacked plate to have two sides: one semi-circular side, and an additional elliptical side. The proposed flow channel design is shown in Figure 1b. In this design, there will be two etched sides for one layer. One side is a semi-circle side and the other side is a semi-ellipse. The geometric parameters used in this study are shown in Figure 2 and Table 2. As per the investigation result of Lee and Lee [11], the maximum stresses occurred at the channel tips. These channel tips are formed during the stacking process of the semi-circular side of the plate and the flat side of the plate. To reduce this sharp corner, an additional elliptical upper channel is proposed to reduce the stress intensity of this flow cross-section design. This elliptical is located above the semi-circular channel. This change in design causes one layer of the stacked plate to have two sides: one semi-circular side, and an additional elliptical side. The proposed flow channel design is shown in Figure 1b. In this design, there will be two etched sides for one layer. One side is a semicircle side and the other side is a semi-ellipse. The geometric parameters used in this study are shown in Figure 2 and Table 2.     As per the investigation result of Lee and Lee [11], the maximum stresses occurred at the channel tips. These channel tips are formed during the stacking process of the semi-circular side of the plate and the flat side of the plate. To reduce this sharp corner, an additional elliptical upper channel is proposed to reduce the stress intensity of this flow cross-section design. This elliptical is located above the semi-circular channel. This change in design causes one layer of the stacked plate to have two sides: one semi-circular side, and an additional elliptical side. The proposed flow channel design is shown in Figure 1b. In this design, there will be two etched sides for one layer. One side is a semicircle side and the other side is a semi-ellipse. The geometric parameters used in this study are shown in Figure 2 and Table 2.

Simplified Design Assessment by Rule Assessment
Simplified design assessment is an analytical approach that can be used in PCHE mechanical design. This approach is based on a conservative model where the channels are modeled as a stayed plate on a flat plate. This design method has been used by Heatric in mechanical design and applied by Nestell and Sham [22]. In this method, PCHE channel design is simplified by using a model proposed by ASME Section VIII [23] codes for non-circular design vessels. As shown in Figure 3, there are six variables considered in this simplified design. Stayed thickness (t 3 ) and thin wall thickness (t 2 ) are Energies 2020, 13, 4711 5 of 17 critical in the stress analysis of the rectangular pressure vessel. By considering those design parameters, the membrane stress at a specific area can be estimated. The membrane stress at the stayed plate t 3 can be found by using Equation (1): where P is the design pressure of PCHE, h is the channel width, and t 3 is the ligament between the channel. Membrane stress at the thin wall thickness (t 2 ) is calculated by using Equation (2) below: where H is the channel height and t 2 is the thin wall thickness. Design by rule also considers the bending stress at the thin wall thickness but stayed plate. Bending stress is not considered the stayed plate because all channels in the same plane are considered at the same pressure. The bending stress can be calculated by using the below Equation (3): where I is a moment of inertia and c is the distance from the neutral axis to the extreme fiber. Finally, the stress value can be used in the design criteria which is: a.
The membrane stress at the stayed plate shall not exceed the maximum allowable design stress: The membrane stress at the thin wall thickness shall not exceed the maximum allowable design stress: c.
The total membrane and bending stress shall not exceed 1.5 times the design stress, where F is a joint factor which is recommended to be 0.7 for diffusion bonding block [24]: Energies 2020, 13, x FOR PEER REVIEW 5 of 18

Simplified Design Assessment by Rule Assessment
Simplified design assessment is an analytical approach that can be used in PCHE mechanical design. This approach is based on a conservative model where the channels are modeled as a stayed plate on a flat plate. This design method has been used by Heatric in mechanical design and applied by Nestell and Sham [22]. In this method, PCHE channel design is simplified by using a model proposed by ASME Section VIII [23] codes for non-circular design vessels. As shown in Figure 3, there are six variables considered in this simplified design. Stayed thickness (t3) and thin wall thickness (t2) are critical in the stress analysis of the rectangular pressure vessel. By considering those design parameters, the membrane stress at a specific area can be estimated. The membrane stress at the stayed plate t3 can be found by using Equation (1): where P is the design pressure of PCHE, h is the channel width, and t3 is the ligament between the channel. Membrane stress at the thin wall thickness (t2) is calculated by using Equation (2) below: where H is the channel height and t2 is the thin wall thickness. Design by rule also considers the bending stress at the thin wall thickness but stayed plate. Bending stress is not considered the stayed Maximum allowable design stress is given by ASME Section III and is a function of the working temperature. Since PCHEs are constructed through the diffusion bonding process, the common understanding is that the bonded structure's strength is below the material's value even though the purpose of welding is to achieve the same structural strength as the material. In diffusion bonding technology, there are two materials commonly used for construction: stainless steel 314, 316/316 L and Energies 2020, 13, 4711 6 of 17 alloy 617. However, in ASME Section III Division 5 at mandatory appendix HBB I−14, the maximum allowable stress intensity for design loading calculation is provided. There are five types of material for elevated temperature nuclear component construction listed as shown in Table 3. For the application of a sodium fast-cooled reactor with the working temperature around 550 • C, five listed materials are usable. However, pressure working conditions should be considered further to decide whether this is acceptable or not. In this study, SS316 was used as a material. The maximum allowable stress intensity for the working environment no more than 550 • C is 88 MPa for SS316.  154  500  99  107  101  81  133  525  86  101  99  64  117  550  74  88  89  48  102  575  69  77  74  35  81  600  65  76  68  26  62  625  51  62  62  -46  650  42  51  51  -29  675  34  39  41  --700  27  30  34  --725  21  23  28  --750  17  18  23  --775  14  13  ---800  11  11 ---By considering the geometrical parameter in Table 2, a simplified calculation by considering Equations (1)-(5) was conducted. The calculation result and comparison to the maximum allowable stress intensity is shown in Figure 4. The calculation result shows that the geometrical design can satisfy the design limit of the diffusion bonded printed circuit heat exchanger since the result of both the primary membrane and the total combined primary membrane and bending stress is below the maximum allowable stress intensity for the design pressure up to 35 MPa. The proposed geometry is aimed to be operated under a working pressure and temperature of 16.7 MPa to 22 MPa and 240 • C to 550 • C, respectively.

Mechanical Integrity Assessment by Finite Element Method Simulation
A two-dimensional finite element method simulation by using commercial COMSOL Multiphysics ® multiphysics was conducted to see the mechanical integrity of the proposed PCHE

Mechanical Integrity Assessment by Finite Element Method Simulation
A two-dimensional finite element method simulation by using commercial COMSOL Multiphysics ® multiphysics was conducted to see the mechanical integrity of the proposed PCHE design. Comsol multiphysics ® is a simulation software for modeling designs, devices, and processes in all fields of engineering, manufacturing, and scientific research [25]. The software was run on ASUS Intel ® Xeon ® CPU 2.0 GHz (32 CPUs) 64 GB RAM. The pressure boundary condition was used in the simulation by considering the SFR steam generator's working pressure. The cold side will be subjected to 21 MPa pressure while the hot channel will be subjected to 0.5 MPa pressure. This high difference in pressure conditions can create problems in terms of the mechanical integrity of the system. Lee and Lee's study has a different boundary condition. Lee and Lee [11] considered unconstrained boundary conditions for the top, bottom, left, and sides of the simulated geometry. However, in this study, the top and the bottom sides were considered as periodic boundary conditions while the left and the right side were symmetric boundary conditions on SUS316 material. The boundary condition used in the simulation is shown in Figure 5.

Mechanical Integrity Assessment by Finite Element Method Simulation
A two-dimensional finite element method simulation by using commercial COMSOL Multiphysics ® multiphysics was conducted to see the mechanical integrity of the proposed PCHE design. Comsol multiphysics ® is a simulation software for modeling designs, devices, and processes in all fields of engineering, manufacturing, and scientific research [25]. The software was run on ASUS Intel ® Xeon ® CPU 2.0 GHz (32 CPUs) 64 GB RAM. The pressure boundary condition was used in the simulation by considering the SFR steam generator's working pressure. The cold side will be subjected to 21 MPa pressure while the hot channel will be subjected to 0.5 MPa pressure. This high difference in pressure conditions can create problems in terms of the mechanical integrity of the system. Lee and Lee's study has a different boundary condition. Lee and Lee [11] considered unconstrained boundary conditions for the top, bottom, left, and sides of the simulated geometry. However, in this study, the top and the bottom sides were considered as periodic boundary conditions while the left and the right side were symmetric boundary conditions on SUS316 material. The boundary condition used in the simulation is shown in Figure 5. A mesh independent test was conducted to validate the geometry and mesh construction. The number of elements was varied, and the maximum stress intensity was observed. The error was determined by subtracting the stress intensity from the highest stress intensity. The mesh sensitivity showed that the results converged for different numbers of elements. The triangular element type was used in this numerical simulation study. The mesh sensitivity result is shown in Table 4 for the different additional elliptical channel radius conditions and in Table 5 for the different channel offset arrangement conditions. Figure 6 shows the offset channel arrangement (δ). A mesh independent test was conducted to validate the geometry and mesh construction. The number of elements was varied, and the maximum stress intensity was observed. The error was determined by subtracting the stress intensity from the highest stress intensity. The mesh sensitivity showed that the results converged for different numbers of elements. The triangular element type was used in this numerical simulation study. The mesh sensitivity result is shown in Table 4 for the different additional elliptical channel radius conditions and in Table 5 for the different channel offset arrangement conditions. Figure 6 shows the offset channel arrangement (δ).    Stainless steel 316 is used as the material of PCHE. Material properties used in this study is shown in Figure 7. Two-dimensional and three-dimensional stress distribution under pressure load is shown in Figure 8. It is shown that the maximum stress intensity is 160 MPa and located at the tip channel of the 21 MPa channel. This is an obvious phenomenon since the upper channel has a higher pressure compared to the bottom one. Stress intensity along with line (a) and (b) mentioned in Figure  5 is shown in Figure 9. Line (a) represents the vertical ligament between the upper and lower channels while line (b) represents the horizontal ligament between channels. Stainless steel 316 is used as the material of PCHE. Material properties used in this study is shown in Figure 7. Two-dimensional and three-dimensional stress distribution under pressure load is shown in Figure 8. It is shown that the maximum stress intensity is 160 MPa and located at the tip channel of the 21 MPa channel. This is an obvious phenomenon since the upper channel has a higher pressure compared to the bottom one. Stress intensity along with line (a) and (b) mentioned in Figure 5 is shown in Figure 9. Line (a) represents the vertical ligament between the upper and lower channels while line (b) represents the horizontal ligament between channels.
The proposed system is modeled under sodium fast-cooled reactor heat exchanger pressure and temperature condition. Hence, the thermal stresses of the individual channel are also important to analyze. To accomplish this, the constant temperature boundary condition of 525 • C and 505 • C were applied for the hot and cold channel, respectively. These boundary conditions were considered as the highest temperature condition in the sodium fast-cooled reactor heat exchanger operation. The periodic boundary condition was also considered and applied at the top and bottom side of the modeled simulation. Symmetric boundary conditions were at the left side and right side. A multilinear isotropic hardening model was used to see the thermal stress of the proposed model. The average isochronous stress-strain curve of SS31600 at various temperatures was obtained from the ASME boiler and pressure vessel code (BPVC) Section III Figure NH-T-1800-B-10. Figure 10a shows the yield stress intensity of SS31600 as a function of temperature. Figure 10b shows the average isochronous stress-strain curve of SS31600 at 631 • C.   The proposed system is modeled under sodium fast-cooled reactor heat exchanger pressure and temperature condition. Hence, the thermal stresses of the individual channel are also important to analyze. To accomplish this, the constant temperature boundary condition of 525 °C and 505 °C were applied for the hot and cold channel, respectively. These boundary conditions were considered as the highest temperature condition in the sodium fast-cooled reactor heat exchanger operation. The periodic boundary condition was also considered and applied at the top and bottom side of the modeled simulation. Symmetric boundary conditions were at the left side and right side. A multilinear isotropic hardening model was used to see the thermal stress of the proposed model. The average isochronous stress-strain curve of SS31600 at various temperatures was obtained from the ASME boiler and pressure vessel code (BPVC) Section III Figure NH-T-1800-B-10. Figure 10a shows the yield stress intensity of SS31600 as a function of temperature. Figure 10b shows the average isochronous stress-strain curve of SS31600 at 631 °C. The simulation result under the high temperature and pressure which is shown in Figure 11 indicates that the highest stress intensity around 238.5 MPa occurs at the tip area of the hot channel due to the high temperature there. A comparison between the geometry without the additional elliptical upper channel and with a 0.25 mm elliptical upper channel shows that 0.25 mm elliptical upper channel stress intensity is 9% lower It again shows that the additional elliptical channel can lower the stress intensity. This is due to the reduced tip channel angle with the addition of a 0.25 mm elliptical upper channel.  The proposed system is modeled under sodium fast-cooled reactor heat exchanger pressure and temperature condition. Hence, the thermal stresses of the individual channel are also important to analyze. To accomplish this, the constant temperature boundary condition of 525 °C and 505 °C were applied for the hot and cold channel, respectively. These boundary conditions were considered as the highest temperature condition in the sodium fast-cooled reactor heat exchanger operation. The periodic boundary condition was also considered and applied at the top and bottom side of the modeled simulation. Symmetric boundary conditions were at the left side and right side. A multilinear isotropic hardening model was used to see the thermal stress of the proposed model. The average isochronous stress-strain curve of SS31600 at various temperatures was obtained from the ASME boiler and pressure vessel code (BPVC) Section III Figure NH-T-1800-B-10. Figure 10a shows the yield stress intensity of SS31600 as a function of temperature. Figure 10b shows the average isochronous stress-strain curve of SS31600 at 631 °C. The simulation result under the high temperature and pressure which is shown in Figure 11 indicates that the highest stress intensity around 238.5 MPa occurs at the tip area of the hot channel due to the high temperature there. A comparison between the geometry without the additional elliptical upper channel and with a 0.25 mm elliptical upper channel shows that 0.25 mm elliptical upper channel stress intensity is 9% lower It again shows that the additional elliptical channel can lower the stress intensity. This is due to the reduced tip channel angle with the addition of a 0.25 mm elliptical upper channel. The simulation result under the high temperature and pressure which is shown in Figure 11 indicates that the highest stress intensity around 238.5 MPa occurs at the tip area of the hot channel due to the high temperature there. A comparison between the geometry without the additional elliptical upper channel and with a 0.25 mm elliptical upper channel shows that 0.25 mm elliptical upper channel stress intensity is 9% lower It again shows that the additional elliptical channel can lower the stress intensity. This is due to the reduced tip channel angle with the addition of a 0.25 mm elliptical upper channel. ASME BPVC suggests the different classifications of stress for the judgment of compliance with the construction specification code. Pm and Pb are used based on Equations (6)-(9) as per ASME Section III. In this study, Pm and Pb are obtained by using the stress linearization function provided by COMSOL. Stress linearization method is a method to obtain the primary membrane and primary bending stress for the judgment purpose according to ASME Section VIII Division II. Pm and Pb were obtained by using the linear elastic model simulated under mechanical forces and internal pressure only. Five stress classification lines (SCLs) were defined on the surface to obtain the stress linearization result. These five SCLs represent the two most important observation areas, between two hot channels, between two cold channels, and between the hot and cold channel. Five SCLs were shown in Figure 12a. SCL 1 represents the area between the cold channels where the high pressure occurs. SCL 2, SCL 3, and SCL 4 represent the area between the cold and hot channel. This ligament is important to be observed due to the high-pressure difference between the hot and cold channel. SCL 5 represents the area between two hot channels. ASME BPVC suggests the different classifications of stress for the judgment of compliance with the construction specification code. Pm and Pb are used based on Equations (6)-(9) as per ASME Section III. In this study, Pm and Pb are obtained by using the stress linearization function provided by COMSOL. Stress linearization method is a method to obtain the primary membrane and primary bending stress for the judgment purpose according to ASME Section VIII Division II. Pm and Pb were obtained by using the linear elastic model simulated under mechanical forces and internal pressure only. Five stress classification lines (SCLs) were defined on the surface to obtain the stress linearization result. These five SCLs represent the two most important observation areas, between two hot channels, between two cold channels, and between the hot and cold channel. Five SCLs were shown in Figure 12a. SCL 1 represents the area between the cold channels where the high pressure occurs. SCL 2, SCL 3, and SCL 4 represent the area between the cold and hot channel. This ligament is important to be observed due to the high-pressure difference between the hot and cold channel. SCL 5 represents the area between two hot channels.
The stress linearization result is shown in Figure 12. It shows that the proposed geometry can satisfy the design loading limit evaluation of ASME Section III. The primary membrane stress levels of all SCLs are lower than the maximum allowable stress intensity (So and 1.5 So) at SCL 1 and 5, and the primary membrane stress and combined primary membrane stress are the same due to the identical pressure condition between the subjected channels. It causes no bending stress to occur in that area. If the linearization result and simplified calculation based on the ASME code result (in Figure 4) are compared, then an almost identical stress intensity can be seen between them. This shows that the simplified calculation and linearization method can be used in mechanical integrity calculation.

Effect of Additional Elliptical Channel Radius
An additional elliptical channel radius is an important geometry condition in double-faced PCHE design. This parameter determines the fluid path surface area. To observe the effect of the elliptical channel radius, five different cases were conducted for 0 < r1/r < 0.5 where r1 is the additional elliptical short radius and r is the semi-circular channel radius. One important parameter to be The stress linearization result is shown in Figure 12. It shows that the proposed geometry can satisfy the design loading limit evaluation of ASME Section III. The primary membrane stress levels of all SCLs are lower than the maximum allowable stress intensity (So and 1.5 So) at SCL 1 and 5, and the primary membrane stress and combined primary membrane stress are the same due to the identical pressure condition between the subjected channels. It causes no bending stress to occur in that area. If the linearization result and simplified calculation based on the ASME code result (in Figure 4) are compared, then an almost identical stress intensity can be seen between them. This shows that the simplified calculation and linearization method can be used in mechanical integrity calculation.

Effect of Additional Elliptical Channel Radius
An additional elliptical channel radius is an important geometry condition in double-faced PCHE design. This parameter determines the fluid path surface area. To observe the effect of the elliptical channel radius, five different cases were conducted for 0 < r 1 /r < 0.5 where r 1 is the additional elliptical short radius and r is the semi-circular channel radius. One important parameter to be observed is the stress concentration factor. Since the geometry is classified as mini dimension plates with a hole, an accurate calculation of stress distribution and stress concentration factor (SCF) is critical. Besides that, a deeper variety of geometric discontinuities or effects with different sizes and shapes is also important to get the best design based on the standard limit given. The stress concentration factor was defined as the ratio of maximum stress intensity to the value of nominal stress. Figure 13 shows the effect of an additional elliptical channel radius to the stress concentration factor of the system. The figure clearly shows that the addition of the elliptical channel radius can reduce the stress concentration factor by double-etched PCHE. Increasing an additional elliptical shape radius up to 0.65 mm can reduce the stress concentration factor by more than 100%. The condition near the edge of the channel has been changed due to the addition of the elliptical upper channel. Stress intensity also slightly changes with the increase in additional elliptical channel radius. Figure 14 shows the stress intensity between two cold channels (a) and between the hot and cold channels (b). It shows that stress intensity along the line between two cold channels slightly decreases with the increase in elliptical shape factor. Stress intensity also slightly changes with the increase in additional elliptical channel radius. Figure 14 shows the stress intensity between two cold channels (a) and between the hot and cold channels (b). It shows that stress intensity along the line between two cold channels slightly decreases with the increase in elliptical shape factor.
Stress linearization result of the membrane stress intensity as shown in Figure 15 shows that the increase in elliptical channel radius can increase the membrane stress at SCL 1 which is a ligament between two cold channels with a pressure of 21 MPa. However, this reduces the membrane stress at SCLs 3 and 4 which is a ligament between the hot and cold channels with a pressure difference of 20.5 MPa. Another SCL shows a slightly similar stress intensity. The combined primary + bending stress intensity of the linearization method shows that the stress intensity of SCL 1 increases by around 33% with the increase in the elliptical upper channel radius. However, it is interesting that there was a decrease in the stress intensity of SCL 2, SCL 3, and SCL 4 with the increase in elliptical upper channel Energies 2020, 13, 4711 13 of 17 radius around 5%. This means that the increase in elliptical channel radius can reduce the bending stress that occurs on the surface. Stress intensity also slightly changes with the increase in additional elliptical channel radius. Figure 14 shows the stress intensity between two cold channels (a) and between the hot and cold channels (b). It shows that stress intensity along the line between two cold channels slightly decreases with the increase in elliptical shape factor. Stress linearization result of the membrane stress intensity as shown in Figure 15 shows that the increase in elliptical channel radius can increase the membrane stress at SCL 1 which is a ligament between two cold channels with a pressure of 21 MPa. However, this reduces the membrane stress at SCLs 3 and 4 which is a ligament between the hot and cold channels with a pressure difference of 20.5 MPa. Another SCL shows a slightly similar stress intensity. The combined primary + bending stress intensity of the linearization method shows that the stress intensity of SCL 1 increases by around 33% with the increase in the elliptical upper channel radius. However, it is interesting that there was a decrease in the stress intensity of SCL 2, SCL 3, and SCL 4 with the increase in elliptical upper channel radius around 5%. This means that the increase in elliptical channel radius can reduce the bending stress that occurs on the surface.

Effect of Shifted Channel Arrangement
A PCHE is a device that consists of a series of plates with numbers of cut-out channels. They are stacked on top of each other to assemble one PCHE unit. Shifted channel arrangement, also known as an offset model, is a model where the plates are arranged in offset horizontally. In the PCHE construction side, this model can be obtained by shifting plates before the diffusion bonding process.
The simulation by using COMSOL Multiphysics was conducted for the stress analysis of this shifted channel configuration. Then, 21 MPa and 0.5 MPa were applied for the cold and hot channel boundary conditions, respectively. Figure 16 shows the von Mises stress distribution for 0.25 elliptical radius channels under non-offset configuration and offset configurations. The maximum stress still occurs on the corner of the channel because of the corner strain of pressure boundary conditions. The maximum von Mises stress is 152.7 MPa and 139.8 MPa of non-offset and 2.5 mm offset configuration, respectively. This means that the maximum stress intensity of 1.5 mm offset configuration is 2% lower compared to the non-offset configuration. A 2.5 mm offset configuration lowers 6% of THE maximum stress intensity while 3% of the maximum stress intensity is lowered under a 3.5 mm offset configuration. The simulation result also shows that the stress concentration factor (Kt) decreases with the offset condition. Figure 17 shows that the stress concentration factor of 2.5 mm offset configuration is the lowest compared to 0.1 and 3.5 mm offset configuration. Simulated cases show that 2.5 mm can be an optimum offset length for a double-etched printed circuit heat exchanger.

Effect of Shifted Channel Arrangement
A PCHE is a device that consists of a series of plates with numbers of cut-out channels. They are stacked on top of each other to assemble one PCHE unit. Shifted channel arrangement, also known as an offset model, is a model where the plates are arranged in offset horizontally. In the PCHE construction side, this model can be obtained by shifting plates before the diffusion bonding process.
The simulation by using COMSOL Multiphysics was conducted for the stress analysis of this shifted channel configuration. Then, 21 MPa and 0.5 MPa were applied for the cold and hot channel boundary conditions, respectively. Figure 16 shows the von Mises stress distribution for 0.25 elliptical radius channels under non-offset configuration and offset configurations. The maximum stress still occurs on the corner of the channel because of the corner strain of pressure boundary conditions. The maximum von Mises stress is 152.7 MPa and 139.8 MPa of non-offset and 2.5 mm offset configuration, respectively. This means that the maximum stress intensity of 1.5 mm offset configuration is 2% lower compared to the non-offset configuration. A 2.5 mm offset configuration lowers 6% of THE maximum stress intensity while 3% of the maximum stress intensity is lowered under a 3.5 mm offset configuration. The simulation result also shows that the stress concentration factor (K t ) decreases with the offset condition. Figure 17 shows that the stress concentration factor of 2.5 mm offset configuration is the lowest compared to 0.1 and 3.5 mm offset configuration. Simulated cases show that 2.5 mm can be an optimum offset length for a double-etched printed circuit heat exchanger.

Conclusions
Finite element analysis simulation has been conducted to observe the mechanical integrity of double-etched printed circuit heat exchanger design. Even this design still needs more efforts to be applied in the future, but the design of double-etched design shows a promising and safer design in terms of a mechanical integrity point of view after some previous study mentioned double-faced PCHE from a thermal-hydraulic point of view. In this case, a sodium fast-cooled reactor steam

Conclusions
Finite element analysis simulation has been conducted to observe the mechanical integrity of double-etched printed circuit heat exchanger design. Even this design still needs more efforts to be applied in the future, but the design of double-etched design shows a promising and safer design in terms of a mechanical integrity point of view after some previous study mentioned double-faced PCHE from a thermal-hydraulic point of view. In this case, a sodium fast-cooled reactor steam

Conclusions
Finite element analysis simulation has been conducted to observe the mechanical integrity of double-etched printed circuit heat exchanger design. Even this design still needs more efforts to be applied in the future, but the design of double-etched design shows a promising and safer design in terms of a mechanical integrity point of view after some previous study mentioned double-faced PCHE from a thermal-hydraulic point of view. In this case, a sodium fast-cooled reactor steam generator heat exchanger is taken as the model for the proposed PCHE design. By considering the ASME Section 3Section construction rule, the proposed design can comply with the design criteria where the primary membrane stress (Pm) and combined with primary bending stress (Pm + Pb) are still below maximum requirement criteria. The application of an additional ellipse upper channel helps the stress intensity decrease in the proposed PCHE channel. Five different cases were conducted for 0 < r 1 /r < 0.5 where r 1 is the additional elliptical short radius and r is the semi-circular channel radius. The simulation shows that the stress intensity was reduced by up to 24%. Apart from that, the stress concentration factor also reduces with the increase in additional elliptical channel radius.
Besides that, the offset channel arrangement is also one of the possible arrangements that can be used in this type of PCHE. This offset arrangement can be done by shifting one plate to another during the construction process. This offset arrangement comes with the promising stress intensity decrease. Simulation results show that the 2.5 mm offset configuration can decrease 9% of the maximum stress intensity compared to the non-offset configuration. This work proposed an additional elliptical upper channel with a 2.5 mm offset configuration is an optimum design that can be applied in PCHE design.

Conflicts of Interest:
The authors declare no conflict of interest.