Abstract

Increased interest in compact heat exchangers (CHXs) to serve as intermediate heat exchangers of very high temperature reactors resulted in significant research and development on their design, analysis, and construction. Printed circuit heat exchangers are a type of CHXs with high thermal efficiency and compactness achieved through diffusion bonding a stack of etched plates with millimeter scaled channels. The diffusion bonding process changes the microstructural and mechanical properties of the wrought metal plates. The current nonnuclear design code ASME section VIII, division 1 captures the material property change through a “joint efficiency factor.” However, the current nuclear design code ASME section III, division 5 does not address or support the diffusion bonded material properties. Hence, there is a need to develop allowable stresses, isochronous curves, and fatigue life curves for various diffusion bonded alloys. In this study, Alloy 800H material was selected to establish the diffusion bonded material properties under tension, creep, fatigue, and creep-fatigue loads at elevated temperatures in the range 550–760 °C. A set of tests on diffusion bonded Alloy 800H (DB 800H) were performed and the acquired data are used in developing allowable stresses Sy, Su, Sr, Sm, St, Smt, So, isochronous curves and fatigue life curves according to the ASME section III, division 5 requirements. This paper also presents detailed procedures used in developing the ASME code section III division 5 design provisions for diffusion bonded Alloy 800H.

1 Introduction

The operation temperature of generation IV (gen. IV) very high temperature reactors (VHTRs) ranges from 550 °C–950 °C [13]. Such elevated operation temperatures impose challenges in designing intermediate heat exchangers, which transfers heat from primary loop to secondary loop [1]. Traditional heat exchangers, such as shell and tube heat exchangers, are inefficient and bulky and hence are uneconomical. Printed circuit heat exchangers (PCHEs) are compact heat exchangers (CHXs) with a very high surface area for heat transfer. PCHEs are highly efficient and compact, making them an ideal option for Gen. IV power plants [4,5]. Typical PCHE consists of a channeled core attached to a set of headers for hot and cold fluid inlets and outlets. The steps to construct a PCHE are discussed below.

The first step is the photochemical etching of the thin plates. The photochemical etching process results in two sets of channeled plates with hot and cold fluid flow channel patterns. The typical channel dimensions range from 0.5 mm to 2 mm [4,6,7]. The next fabrication step stacks alternate hot and cold fluid plates inside a vacuum furnace and subjected to specified compressive stress at a specified temperature for a specified period. The compressive stress at elevated temperature promotes the atomic diffusion and grain growth across plate interface, consequently bonds all plates in a single monolithic PCHE core with thousands of millimeter scaled channels. Usually, multiple PCHE channeled cores are welded together to form the core of a PCHE unit. In the final fabrication step, headers are welded to the PCHE core for the inlet and outlet of hot and cold fluids.

The diffusion bonding process used in fabricating PCHE cores is an innovative technique that changes the plate alloy's microstructure, resulting in changes in its physical and mechanical properties. The microstructural and property changes of diffusion bonded Alloy 617, 316H, and Haynes 230 have been demonstrated [6,810]. The authors have determined the mechanical properties of diffusion bonded 800H (DB 800H) through a series of tension tests at elevated temperatures [11]. A summary of the tension test results are shown in Table 1. The reduction in strengths of the DB 800H can be observed by comparing the strengths of the wrought base metal Alloy 800H (BM 800H) as shown in Fig. 1. The current ASME section VIII nonnuclear design code captures such material property change by implementing the “joint efficiency factor” of 0.7 on all allowable stresses. The ASME section III design code does not specify any guidelines for incorporating material property change due to diffusion bonding. Also, in the ASME code there are no provisions to capture the diffusion bonded material performance at elevated temperatures, especially in the creep regime. Hence, developing allowable stresses, isochronous curves, and fatigue life curves of various diffusion bonded alloys is essential. Including such diffusion bonded material properties in the ASME code section III division 5 will allow analysis and design of PCHEs for nuclear service.

Fig. 1
Comparison of the yield and ultimate tensile stresses of the diffusion bonded 800H from tension tests at different temperatures [11] and the wrought metal given in the ASME section III, division 5 design code [12]
Fig. 1
Comparison of the yield and ultimate tensile stresses of the diffusion bonded 800H from tension tests at different temperatures [11] and the wrought metal given in the ASME section III, division 5 design code [12]
Close modal
Table 1

Elastic modulus, yield, and ultimate tensile strength, and % elongations of diffusion bonded 800H from tests performed at various temperatures under axial load perpendicular to the diffusion bonding planes

T (°C)E (GPa)σy (MPa)σu (MPa)% elongation
2019718752851.2
2019316050755.6
2019517852751.5
2019918553155.1
45016812345357.9
45017912545155.1
50017312541137.5
55016811632635.8
55015411633035.9
60015812131722.7
65015911029520.4
65015911329020.6
70014811122713.5
7601449514516.3
76013412017015.0
76013410316315.7
7601329715314.3
T (°C)E (GPa)σy (MPa)σu (MPa)% elongation
2019718752851.2
2019316050755.6
2019517852751.5
2019918553155.1
45016812345357.9
45017912545155.1
50017312541137.5
55016811632635.8
55015411633035.9
60015812131722.7
65015911029520.4
65015911329020.6
70014811122713.5
7601449514516.3
76013412017015.0
76013410316315.7
7601329715314.3

Keating et al. [13] investigated the current state of ASME code section III for PCHE design and analysis and provided a list of ASME code gaps to be addressed. The gaps relevant to this study, listed below, indicate the information needed to be developed for diffusion bonded alloys:

  1. Allowable stress and material property data of diffusion bonded materials.

  2. Creep–fatigue curves of diffusion bonded materials.

  3. Isochronous curves of diffusion bonded materials.

This study addresses these gaps for diffusion bonded Alloy 800H (DB 800H). A set of tensile, fatigue, creep, and creep-fatigue tests were performed on diffusion bonded 800H at elevated temperatures. The data developed are presented and discussed in Ref. [11]. This dataset is used to develop the allowable stresses, isochronous curves, fatigue life curves, and creep-fatigue damage envelop. The methodology used and material properties of DB 800H developed is discussed below.

2 Experimental Data Developed

A diffusion bonded solid block (200 × 200 × 200 mm dimensions) of Alloy 800H fabricated by Vacuum Process Engineering Inc. (VPE, Sacramento, CA) was acquired. Plates of 1.6 mm thickness were diffusion bonded together following ASME section VIII, division 1, code case 2621-1. Diffusion bonding methodologies and parameters are the proprietary information of the company and hence cannot be presented. Various ASTM standard specimens were fabricated from the diffusion bonded block using wire electrical discharge machining. A set of tension, fatigue, creep, and creep-fatigue tests were performed. The readers are referred to Ref. [11] for the mechanical tests and scanning electron microscope examinations performed and corresponding results with critical discussions on the performance of diffusion bonded Alloy 800H (DB 800H). A summary of the tensile test performed and results are shown in Table 1. Creep test results under prescribed stresses at different temperatures are presented in Table 2. Readers are referred to Ref. [11] for DB800H creep tests, which showed primary, secondary, and small tertiary creep regime. Few tests did not show significant tertiary creep regime. All tests failed through bond delamination mechanism with small reduction in area. Strain controlled fatigue and creep–fatigue experimental results are presented in Tables 3 and 4, respectively. Table 4 shows the positive strain peak dwell period, tdwell, prescribed in each test to determine the influence of creep–fatigue interaction on rupture life.

Table 2

Time to reach 1% strain, onset of tertiary creep, and rupture of DB 800H specimens

T (°C)σt1%ttertrup
5501802311773
5501978.9887.85
550197231.92441a
6001380.11306a
6501001522a
650114816.31306a
65011639301
6501602057b
700681261a
76040494.55931002.7
76040217017502624
76046190.9275706.2
76046212.7278960.8
7606221.2579129
7606835.593180
760687.66981.0
760689.84680.9
T (°C)σt1%ttertrup
5501802311773
5501978.9887.85
550197231.92441a
6001380.11306a
6501001522a
650114816.31306a
65011639301
6501602057b
700681261a
76040494.55931002.7
76040217017502624
76046190.9275706.2
76046212.7278960.8
7606221.2579129
7606835.593180
760687.66981.0
760689.84680.9
a

Creep tests discontinued before rupture at this time.

b

Specimen failed in primary creep regime.

Table 3

Fatigue life cycles for failure Nf of three DB 800H specimens tested at each temperatures, where failure is defined as 20% peak force drop

T (°C)ε˙ (/s)εrNf1Nf2Nf3
5500.0010.006473262543413
6500.0010.006221318932147
7600.0010.006476452323
T (°C)ε˙ (/s)εrNf1Nf2Nf3
5500.0010.006473262543413
6500.0010.006221318932147
7600.0010.006476452323
Table 4

Creep-fatigue life cycles for failure Nf of DB 800H specimens, where failure is defined as 20% peak force drop

T (°C)ε˙rεrtdwellNf1Nf2Nf3
5500.0010.006600 s282297258
6500.0010.006300 s87127133
6500.0010.006600 s113
7600.0010.006120 s7711176
7600.0010.006600 s60
T (°C)ε˙rεrtdwellNf1Nf2Nf3
5500.0010.006600 s282297258
6500.0010.006300 s87127133
6500.0010.006600 s113
7600.0010.006120 s7711176
7600.0010.006600 s60

3 Allowable Stresses of DB 800H

The procedure of ASME section III HBB-3221(b) is used for determining allowable stresses of diffusion bonded Alloy 800H (DB 800H). The tension tests performed by the authors [11] are used to determine the allowable stresses of DB 800H. In diffusion bonding, no filler material is used like in conventional welding; hence, the allowable stress determination procedure used for solid wrought metal is followed. Welding of components, such as header, and pipes, to DB 800H core is considered as weldments [13]. As mentioned before, the current nonnuclear design code captures the mechanical property change due to welding through a join efficiency factor of 0.7 as per the section VIII, division 1, code case 2621-1 [14]. This reduction factor is used as the upper bound of the allowable stress. Hence, the general format of allowable stress calculation as per Eq. (1) is adopted
(1)

[13]

where Si,DB is the allowable stress of DB 800H, Si is the calculated lower bound allowable stress from DB 800H tests (see below Secs. 3.1 and 3.2, Si,BM is allowable stress of wrought base metal Alloy 800H (BM 800H), and i is the selected allowable stress type. The detailed procedure in the allowable stress calculations and related discussions are presented below.

3.1 Yield Stress, Sy.

The yield stress is determined using a stress–strain curve with a 0.2% strain offset. A line parallel to the initial linear region of the stress–strain curve, used for calculating elastic modulus, is drawn from the total strain of 0.2%. The yield stress is determined by the intersection point of the offset line and the stress–strain curve. These calculated yield stress values are rounded off to the nearest MPa per section II, part D, nonmandatory Appendix D-130 [15]. The determined yield stresses from different tests are presented in Table 1. A lower bound of the calculated yield stresses is determined based on Eq. (2) [16]
(2)

where SDB and Ci are model parameters determined by fitting experimental yield stresses through a gradient descent method, and optimum parameters are listed in Table 5. The allowable yield stresses are then determined based on Eq. (1). The allowable yield stress values of DB 800H determined are listed in Table 6. A comparison of the test yield stresses to the allowable yield stresses determined is shown in Fig. 2.

Fig. 2
Comparison of yield and ultimate stresses test data against the allowable yield and ultimate stresses determined for DB 800H
Fig. 2
Comparison of yield and ultimate stresses test data against the allowable yield and ultimate stresses determined for DB 800H
Close modal
Table 5

Parameters of Eqs. (2) and (3) used for determining allowable yield and ultimate tensile stresses of diffusion bonded 800H

ParametersSySu
SDB155375
C00.0060270.002156
C1−0.00022−1.4 × 10−5
C2−2.2 × 10−6−2.6 × 10−6
C38.22 × 10−91.42 × 10−8
C4−1 × 10−11−2.5 × 10−11
C54.39 × 10−151.24 × 10−14
ParametersSySu
SDB155375
C00.0060270.002156
C1−0.00022−1.4 × 10−5
C2−2.2 × 10−6−2.6 × 10−6
C38.22 × 10−91.42 × 10−8
C4−1 × 10−11−2.5 × 10−11
C54.39 × 10−151.24 × 10−14
Table 6

Allowable stresses Sy, Su, Sm, and So determined for DB 800H

T (°C)SySuSmSo
20126.0315.0
45086.0280.469.766.7
47585.6280.069.366.1
50085.2278.069.066.4
52584.7274.168.959.3
55084.0268.168.052.9
57583.0259.867.247.2
60081.8249.166.342.2
62580.2236.165.037.7
65078.1221.263.332.2
67575.7204.761.326.8
70072.8187.056.122.3
72569.4165.049.518.6
75065.6144.343.315.5
T (°C)SySuSmSo
20126.0315.0
45086.0280.469.766.7
47585.6280.069.366.1
50085.2278.069.066.4
52584.7274.168.959.3
55084.0268.168.052.9
57583.0259.867.247.2
60081.8249.166.342.2
62580.2236.165.037.7
65078.1221.263.332.2
67575.7204.761.326.8
70072.8187.056.122.3
72569.4165.049.518.6
75065.6144.343.315.5

3.2 Ultimate Tensile Stress, Su.

The maximum stress values from the recorded tensile stress–strain curves are collected as the ultimate tensile strengths of DB 800H at different temperatures. A lower bound of the tensile strength values are determined based on Eq. (3) [16]
(3)

where SDB and Ci are model parameters determined by fitting the experimental ultimate tensile stresses, and model parameters are listed in Table 5. The allowable ultimate stresses are then determined based on Eq. (1). The allowable ultimate stress values of DB 800H determined are listed in Table 6. A comparison of ultimate stresses from the tests to the allowable ultimate stresses determined is shown in Fig. 2.

3.3 Time-Independent Allowable Stress, Sm.

Sm is the time-independent allowable stress determined per HBB-2160(d)(3), and is the minimum of

  1. The product of one-third of the minimum tensile strength at room temperature and the tensile strength reduction factor [Table HBB-3225-2, ASME BPVC.III.5-2017];

  2. The product of one-third of the minimum tensile strength at specified temperature and the tensile strength reduction factor [Table HBB-3225-2, ASME BPVC.III.5-2017];

  3. The product of two-thirds of the minimum yield strength at room temperature and the yield strength reduction factor [Table HBB-3225-2, ASME III.5-2017];

  4. The product of 90% of yield strength at specified temperature and the yield strength reduction factor [Table HBB-3225-2, ASME III.5-2017]

The yield strength and ultimate strength reduction factor for long time elevated temperature service for DB 800H are not available; hence, wrought base metal Alloy 800H reduction factors from Table HBB-3225-2 are used for DB 800H.

3.4 Primary Membrane Allowable Stress So.

So is the maximum allowable limit of the primary membrane stress to be used as a reference for stress calculations under design loading [17]. The provisions per HBB-3221(b)(1) These stresses are determined as follows:

  1. (1)

    For temperatures below the range where creep and stress rupture strength governs, the So is:

    1. the specified minimum tensile strength at room temperature divided by 3.5

    2. the tensile strength at specified temperature divided by 3.5

    3. two-thirds of the specified minimum yield strength at room temperature

    4. two-thirds of the yield strength at specified temperature

      As per ASME code section III, creep effects in Alloy 800H are neglected below 425 °C (800 °F); hence, below 425 °C, tension test data are used to calculated So as follows:

  1. (2)

    For temperatures in the range where creep and stress rupture strength governs, So is calculated as follows:

    1. 100% of the average stress to produce a creep rate of 0.01%/1000 h

    2. 100Favg% of the average stress to cause rupture at the end of 100,000 h, where for temperatures below 815 °C (1500 °F) Favg is 0.67.

    3. 80% of the minimum stress to cause rupture at the end of 100,000 h

The power law function is fitted to the test data as shown in Fig. 3. The power law function was fitted by using the wrought metal data. By keeping the stress exponent same, linear fitting parameter was determined. Ideally, a large creep dataset would be needed to determine the power term, which would capture the dislocation mechanisms at lower stresses accurately. Due to limited creep tests on DB 800H, the power exponent of wrought metal was used. By extrapolation, the stress to produce creep rate of 0.01%/1000 h is calculated. The final values of S0 after accounting the maximum bound of 0.7 reduction factor are presented in Table 6.

Fig. 3
Minimum creep rate against stress of wrought and diffusion bonded 800H at 760 °C
Fig. 3
Minimum creep rate against stress of wrought and diffusion bonded 800H at 760 °C
Close modal

3.5 Time and Temperature-Dependent Allowable Stress Sr.

The creep rupture data of DB 800H is used to calculate the Larson Miller Parameter using Eq. 4. In this equation, T is absolute temperature in Kelvin(K), t is rupture time in h, and C is a constant set to 15.21, which is determined for BM 800H by Swindeman [18]. The linear regression equation (Eq. (5)) is used to fit the calculated Larson Miller parameter (LMP) data from creep rupture tests shown with cross sign “X” in Fig. 4. LMP calculated based on data shown in cross sign “X” are used to estimate the rupture life of few tests. These tests were discontinued once test time exceeded the expected rupture life. These discontinued tests DB 800H (D) are shown with plus sign “+” in Fig. 4. Note that DB 800H (D) was not used to calculate the fitting line shown in Fig. 4. The X-ray scans of specimens showed microvoid distributed at diffusion bonds, which may change creep performance from specimen to specimen [11]. Hence, these tests served as validation for interpolated creep life and statistical study to ensure the LMP fit rupture predictions are reasonable. More creep tests will change LMP parameters, and improve rupture life predictions and allowable stress reliability.

Fig. 4
Larson Miller Parameter of diffusion bonded 800H determined based on creep rupture data shown in cross sign “X.” This LMP was used to estimate the creep rupture life for few tests. These tests are shown in plus sign “+” DB 800H (D), and were discontinued when test life exceeded the estimated rupture life. These tests validate the interpolation of rupture life to intermediate temperatures.
Fig. 4
Larson Miller Parameter of diffusion bonded 800H determined based on creep rupture data shown in cross sign “X.” This LMP was used to estimate the creep rupture life for few tests. These tests are shown in plus sign “+” DB 800H (D), and were discontinued when test life exceeded the estimated rupture life. These tests validate the interpolation of rupture life to intermediate temperatures.
Close modal
The linear regression equation (Eq. (5)) is used to fit the calculated LMP data with respective creep stress values of DB 800H as shown in Fig. 4 
(4)
(5)
Regression coefficients a0 and a1 are calculated for the LMP values for the time of 1% strain, time of onset of tertiary creep, and time of rupture, and presented in Table 7. The standard error of estimate is calculated per Eq. (6) [17]
(6)
Table 7

Constant C of LMP (Eq. (4)), coefficients of linear regression fit for calculated LMP of DB 800H (Eq. (5)) and SEE calculated from Eq. (6)

Data typeCa0a1SEE
1% strain 18.12 43325 −12962 0.543357 
Tertiary creep 14.35 27539 −6001.7 0.235697 
Creep rupture 15.21 29070 −6210.4 0.201313 
Data typeCa0a1SEE
1% strain 18.12 43325 −12962 0.543357 
Tertiary creep 14.35 27539 −6001.7 0.235697 
Creep rupture 15.21 29070 −6210.4 0.201313 

where Nd is the number of data points and Df is the degree-of-freedom calculated as a total constraint number in Larson Miller parameter and degree of the regression equation. The LMP function has one constant C, and first-order regression equation is used (Eq. (5)); hence, the degree-of-freedom in above equation is two [19]. An adjustment of 1.65 multiple of standard error of estimate (SEE) is applied to logarithmic time to represent the 95% lower bound of stress data. This procedure results in the calculation of the minimum stress to cause 1% strain for a given time and temperature. The same procedure is used with creep data of tertiary creep and for creep rupture life by using Eqs. (4) and (5) to calculated minimum stress for tertiary creep and creep rupture. The calculated regression coefficients of DB800H are presented in Table 7.

The maximum value of minimum stress to rupture is bounded by the upper stress limit of Su/1.1. The proposed minimum stress to rupture values of diffusion bonded 800H material is presented in Table 8.

Table 8

Time and temperature-dependent allowable stresses, Sr for diffusion bonded Alloy 800H

Sr
Time
T (°C)1 h3 h10 h30 h100 h300 h1000 h3000 h10,000 h
500252.7252.7252.7252.7252.7252.7252.7224.7192.6
525249.2249.2249.2249.2249.2249.2216.3187.4159.3
550243.7243.7243.7243.7243.7212.1180.6155.9131.8
575236.2236.2236.2236.2207.9178.5151.2129.5109.1
600226.4226.4226.4207.9175.0149.8126.0107.890.3
625214.7214.7207.2176.4147.7126.0105.088.974.2
650201.1201.1176.4149.8124.6105.087.573.560.9
675186.1181.3150.5126.7105.088.272.860.950.4
700170.0154.7128.1107.188.273.560.250.441.3
725150.0132.3108.590.374.261.649.741.333.6
750131.2113.492.476.362.351.141.334.327.3
Sr
Time
T (°C)1 h3 h10 h30 h100 h300 h1000 h3000 h10,000 h
500252.7252.7252.7252.7252.7252.7252.7224.7192.6
525249.2249.2249.2249.2249.2249.2216.3187.4159.3
550243.7243.7243.7243.7243.7212.1180.6155.9131.8
575236.2236.2236.2236.2207.9178.5151.2129.5109.1
600226.4226.4226.4207.9175.0149.8126.0107.890.3
625214.7214.7207.2176.4147.7126.0105.088.974.2
650201.1201.1176.4149.8124.6105.087.573.560.9
675186.1181.3150.5126.7105.088.272.860.950.4
700170.0154.7128.1107.188.273.560.250.441.3
725150.0132.3108.590.374.261.649.741.333.6
750131.2113.492.476.362.351.141.334.327.3

3.6 Temperature and Time-Dependent Allowable Stress St.

St are the temperature and time-dependent allowable stresses, which are determined using data from long-term, constant load, uniaxial tests. Per HBB-3221(b)(1), for each specific time t, St is the smallest of the following values:

  1. 100% of the average stress required to obtain a total (elastic, plastic, primary, and secondary creep) strain of 1%;

  2. 80% of the minimum stress to cause initiation of tertiary creep;

  3. 67% of the minimum stress to cause rupture

Similar to the Sr determination in Sec. 3.5, the Larson Miller parameter (Eq. (4)) and liner fitting function (Eq. (5)) with respective SEE values are used to calculate the minimum stresses of tertiary creep and rupture. For calculating average stress for 1% strain, no SEE values are used, but the rest of the calculation procedure remains the same as for Sr. The final calculated St values are tabulated in Table 9.

Table 9

Time and temperature-dependent allowable stresses, St for diffusion bonded Alloy 800H

St
Time
T (°C)1 h10 h30 h100 h300 h1000 h3000 h10,000 h
50089.689.689.687.381.876.171.366.4
52588.288.284.778.773.568.363.859.3
55086.882.076.570.966.161.257.152.9
57586.174.269.163.959.454.951.147.2
60078.567.262.457.553.449.345.842.2
62571.360.856.351.848.044.241.037.7
65064.855.050.946.743.239.636.733.6
67558.949.845.942.138.835.532.830.0
70053.645.141.537.934.931.929.426.8
72548.740.837.534.131.428.626.323.4
75044.236.933.830.828.225.623.519.4
St
Time
T (°C)1 h10 h30 h100 h300 h1000 h3000 h10,000 h
50089.689.689.687.381.876.171.366.4
52588.288.284.778.773.568.363.859.3
55086.882.076.570.966.161.257.152.9
57586.174.269.163.959.454.951.147.2
60078.567.262.457.553.449.345.842.2
62571.360.856.351.848.044.241.037.7
65064.855.050.946.743.239.636.733.6
67558.949.845.942.138.835.532.830.0
70053.645.141.537.934.931.929.426.8
72548.740.837.534.131.428.626.323.4
75044.236.933.830.828.225.623.519.4

3.7 Primary Membrane Allowable Stress Smt.

Smt is the allowable limit of general primary membrane stress, to be used as a reference for stress calculations for the actual service life under the level A and B service loading. Smt is determined as the lower value of Sm (time-independent) and St (time-dependent), and is documented in Table 10.

Table 10

Time and temperature-dependent allowable stresses, Smt for diffusion bonded Alloy 800H

Smt
Time
T (°C)1 h3 h10 h30 h100 h300 h1000 h3000 h10,000 h
50069.069.069.069.069.069.069.069.066.4
52568.668.668.668.668.668.668.363.859.3
55068.068.068.068.068.066.161.257.152.9
57567.267.267.267.263.959.454.951.147.2
60066.266.266.262.457.553.449.345.842.2
62564.964.960.856.351.848.044.241.037.7
65063.359.955.050.946.743.239.636.733.6
67558.954.449.845.942.138.835.532.830.0
70053.649.345.141.537.934.931.929.426.8
72548.744.740.837.534.131.428.626.323.4
75043.340.636.933.830.828.225.623.519.4
Smt
Time
T (°C)1 h3 h10 h30 h100 h300 h1000 h3000 h10,000 h
50069.069.069.069.069.069.069.069.066.4
52568.668.668.668.668.668.668.363.859.3
55068.068.068.068.068.066.161.257.152.9
57567.267.267.267.263.959.454.951.147.2
60066.266.266.262.457.553.449.345.842.2
62564.964.960.856.351.848.044.241.037.7
65063.359.955.050.946.743.239.636.733.6
67558.954.449.845.942.138.835.532.830.0
70053.649.345.141.537.934.931.929.426.8
72548.744.740.837.534.131.428.626.323.4
75043.340.636.933.830.828.225.623.519.4

4 Isochronous Curves

The isochronous curves are a series of curves that present the different creep data points at different times for different stresses at a given temperature. The construction of isochronous curves uses the following steps: the total strain is decomposed into elastic, plastic, and creep strain as shown in Eq. (7). Elastic strain is calculated using the Hook's law, Eq. (8), where E is elastic modulus at a given temperature
(7)
(8)
Creep strain is calculated using the Messner et al. [20] equation, Eq. (9), where C1, C2, and C3 are calculated using Eq. (10). These parameters are calculated using the strain versus time history from creep tests. A gradient optimization algorithm available in matlab is implemented to determine an optimum parameter set. Equation (9) can simulate the creep strains for a selected time and stress value
(9)
(10)
Plastic strains are calculated with an empirical equation, Eq. (11), where B̂i are temperature-dependent parameters at a specified T in °C calculated using Eq. (12)
(11)
(12)

While calculating the isochronous curves, a stress value is selected, and creep strains are calculated at different times. Eq. (11) is a stress function in terms of plastic strain; hence, this equation is solved numerically to determine the plastic strain values corresponding to a stress. By fitting the experimental responses from tension tests, parameters of Eq. (11) are evaluated for a specified temperature as shown in Table 11.

Table 11

Parameters of plastic stress–strain simulation for hot tensile curves in the isochronous stress–strain curves for all temperatures

i1234
Bi(1)3.1280.172290.0375670.004503
Bi(2)−0.3547−0.0006540.0132890.00498
Bi(3)0.55510.026628−0.0296220.048257
Bi(4)−0.4807−0.0936650.006146−0.101054
Bi(5)00.058690.0051200.054705
i1234
Bi(1)3.1280.172290.0375670.004503
Bi(2)−0.3547−0.0006540.0132890.00498
Bi(3)0.55510.026628−0.0296220.048257
Bi(4)−0.4807−0.0936650.006146−0.101054
Bi(5)00.058690.0051200.054705

The creep test provides a time history data, from which strain value at desired time is extracted from experiments and used to fit the isochronous curves simulation. The proposed isochronous curves are compared against the test data from tension tests (hot tension curve) and from the creep tests at different times in Fig. 5. In this figure, a good comparison between simulated isochronous curves and test data at 550 °C, 650 °C, and 760 °C is observed. Hence, the equation set presented above with the optimized parameter set is used to simulate the isochronous curves at intermediate temperatures for diffusion bonded Alloy 800H as shown in Figs. 6(a)6(c).

Fig. 5
Comparison of the isochronous curve simulation of the DB800H against experiments from the tension test and creep data
Fig. 5
Comparison of the isochronous curve simulation of the DB800H against experiments from the tension test and creep data
Close modal
Fig. 6
Isochronous curves of diffusion bonded alloy 800H at (a) 538 °C, (b) 649 °C, and (c) 760 °C
Fig. 6
Isochronous curves of diffusion bonded alloy 800H at (a) 538 °C, (b) 649 °C, and (c) 760 °C
Close modal

5 Fatigue Curves

The design fatigue curves in section III, division 5, figure HBB-T-1420-1C correlate the strain range versus the fatigue life, Nf. These curves are developed from hundreds of fatigue tests. A reduction factor of “2” on stress in the high-cycle fatigue regime and “20” on fatigue life in the low-cycle fatigue regime are usually used to develop the fatigue curves of metals.

Only a few fatigue tests have been conducted on DB 800H at three temperatures. Ideally, these few tests are not enough to generate a new set of design fatigue curves. Hence, a reduction factor of 10 is used on the life of wrought base metal 800H fatigue curves and compared to the DB 800H fatigue test data in Fig. 7. In this figure, a reduction factor of 10 on wrought base metal 800H (BM 800H) life bounds the lower values of the DB 800H test life. Hence, for the S–N curves of DB 800H, it is proposed to use the figure HBB-T-1420-1C and a reduction factor, Rf = 10, to develop S-N curves of DB 800H. Table 12 shows the S–N curve values of DB 800H. These proposed fatigue curves of DB 800H can be revised as a large set of test data becomes available.

Fig. 7
DB 800H fatigue test data and proposed DB 800H low cycle fatigue curves developed based on the BM 800H curves in the ASME code
Fig. 7
DB 800H fatigue test data and proposed DB 800H low cycle fatigue curves developed based on the BM 800H curves in the ASME code
Close modal
Table 12

Proposed design fatigue cycles for diffusion-bonded Alloy 800H

Number of cycle538 °C649 °C760 °C
100.011600.009320.00774
200.008490.006780.00562
400.006600.005330.00469
1000.005150.004170.00388
2000.004540.003660.00349
3000.004330.003470.00309
4000.004090.003270.00270
10000.002930.002340.00212
20000.002430.001970.00183
40000.002120.001750.00164
10,0000.001940.001550.00149
Number of cycle538 °C649 °C760 °C
100.011600.009320.00774
200.008490.006780.00562
400.006600.005330.00469
1000.005150.004170.00388
2000.004540.003660.00349
3000.004330.003470.00309
4000.004090.003270.00270
10000.002930.002340.00212
20000.002430.001970.00183
40000.002120.001750.00164
10,0000.001940.001550.00149

6 Creep–Fatigue Damage

Creep–fatigue damage in test specimens is calculated as per ASME section III, division 5. The first term in the damage equation (Eq. (13)) is the fatigue damage, a ratio of experienced fatigue cycle number (n), and estimated fatigue life (Nr) summation for multiple strain ranges. Estimated fatigue life for DB 800H is obtained from Fig. 7. The second term calculates the creep damage, a ratio of Δt, time increment at a stress σ and tr, estimated rupture life at the stress from Table 8 
(13)
The rupture life for a given stress (tr|σ) is calculated by equating LMP in Eqs. (4) and (5), resulting in Eq. (14)
(14)
This rupture equation replaces the denominator in the creep damage equation for a given stress. The total creep damage during strain dwell, when stress relaxation occurs, is calculated by integrating the incremental time damage over the strain dwell period, giving Eq. (15) to calculate the creep damage during strain dwell
(15)

The above equation's stress values can either be imported from creep–fatigue tests or a material model that can simulate the stress relaxation history for the dwell period. The following discussion presents the material model and associated calculations.

During the dwell period, the total strain rate is zero; however, the inelastic and elastic strain rates are not zero as per Eq. (16). The inelastic strain rate during dwell period is governed by the creep deformation, and can be calculated using the Power law, Eq. (17). Hook's law calculates the elastic strain rate. Thus, Eq. (18) can be obtained from Eq. (16)
(16)
(17)
(18)
Upon rearranging and integrating Eq. (18) from σmax at zero dwell time to σ at time tdwell, stress at the end of strain dwell is obtained as
(19)
The creep damage shown in Eq. (15) can be rewritten in terms of stress as
(20)
The stress term in Eq. (19) is plugged in Eq. (20), and integrated to get the final creep damage during strain dwell of a cycle as follows:
(21)

where M=Na1T+1 and σend is stress at the end of strain dwell per Eq. (19). Note that the C*=C+1.65×SEE for creep rupture per Table 7 to calculated the 95% lower bound of stress rupture. This equation calculates the creep damage during stress relaxation of strain dwell in a cycle. Multiplying this equation with the cycle number in the creep–fatigue test gives the total creep damage accumulated in the test specimen till rupture.

The model parameters N and K in Eq. (17) are determined for three temperatures (Table 13) based on relaxation histories obtained from experiments. As the N and K parameters are determined, the creep damage accumulation in creep–fatigue test specimens can be calculated using Eq. (21). The fatigue damage in creep–fatigue test specimens can be calculated by using the first term in Eq. (13) as stated earlier. As sufficient test data are not available to develop a creep–fatigue envelope for DB 800H, the calculated creep–fatigue damages of DB 800H specimens at three temperatures are compared against the creep–fatigue envelope of BM 800H from ASME section III, division 5, HBB-T-1420-2 in Fig. 8. This figure shows that the DB 800H test specimens creep–fatigue damage accumulations up to rupture fall outside the BM 800H ASME code creep–fatigue damage envelop. This indicates that the BM 800H envelope can be used for design and analysis of PCHE constructed of diffusion bonded 800H.

Fig. 8
Creep-fatigue damage calculation of the diffusion bonded 800H specimens at 550 °C, 650 °C, and 760 °C compared against creep-fatigue damage envelop of wrought metal 800H per HBB-T-1420-2
Fig. 8
Creep-fatigue damage calculation of the diffusion bonded 800H specimens at 550 °C, 650 °C, and 760 °C compared against creep-fatigue damage envelop of wrought metal 800H per HBB-T-1420-2
Close modal
Table 13

Material model parameter for creep damage calculations

Parameters550 °C650 °C760 °C
K1.00×10292.34×10243.981×1021
N9.28.17.81
Parameters550 °C650 °C760 °C
K1.00×10292.34×10243.981×1021
N9.28.17.81

7 Discussion, Conclusions, and Future Work

This study developed material properties and design curves of diffusion bonded Alloy 800H (DB 800H) to design and analyze PCHE following the ASME code section III division 5 provision. To facilitate this development, a set of tension, creep, fatigue, and fatigue–creep tests on DB 800H, and few fatigue and creep–fatigue tests on wrought 800H ASTM standard specimens were conducted by the authors [11]. The material properties and design curves development were motivated by the ASME code gap analysis by Keating et al. [13]. Hence, this study addressed the material property-related gaps toward developing ASEM code section III division 5 design and analysis methodology for PCHE.

The allowable stresses of DB800H are developed following the procedure of ASME section III HBB-3221(b). First, the allowable stresses Sy, Su, Sm, and So are developed based on the tension test data. Following the allowable stresses, Sr, St, and Smt are developed based on the creep test data. The LMPs are determined and fitted with a linear function for these three time and temperature-dependent allowable stresses. A SEE calculates the lower bound of rupture stresses from the LMP plot. A joint efficiency factor of 0.7 bounds the calculated allowable stresses to maintain consistency with the existing joint efficiency factor for nonnuclear components provided in section VIII, division 1. The calculated allowable stresses are tabulated similarly as in ASME code.

The tension tests at various temperatures and creep strain time histories are then used to determine the isochronous curves for DB 800H for the temperature range 510 °C–760 °C up to 10,000 h. With the small number of fatigue and creep–fatigue test data of DB 800H developed, it was not possible to develop fatigue curves and creep–fatigue envelopes for DB 800H. An extensive set of data are needed for the development of these curves and envelope. Hence, fatigue design (S–N) curves of DB 800H are proposed based on the BM 800H ASME curves in section III, division 5, figure HBB-T-1420-1C with a reduction factor of 10. These S–N curves are observed to bind the lower values of the DB 800H fatigue test data. When the creep–fatigue damage accumulations in the DB 800H test specimens are calculated and compared to the ASME BM 800H creep–fatigue envelope (Fig. HBB-T-1420-2), it is observed that the BM 800H envelope can be used PCHE design unlit sufficient data are available to develop a DB 800H envelop.

It is also noted here that the material properties and design curves proposed in this study for DB 800H are based on a limited set of tests. As more test data become available, the proposed allowable stresses and design curves can be reevaluated and revised as needed. The creep test data for developing the time and temperature-dependent allowable stresses and isochronous curves are limited to 1000 h of rupture life. The LMPs are determined and fitted with a linear equation to interpolate for intermediate temperatures and extrapolate allowable stresses up to 10,000 h. The equations developed should not be used to extrapolate the creep life of DB800H beyond 10,000 h. Creep tests with a target rupture life of 10,000–30,000 h will be needed to extend the allowable stresses to longer creep life.

It is demonstrated that a reduction factor of 10 on BM 800H fatigue curves bounds the DB 800H fatigue life data. Only a few fatigue tests, all tested with a strain range of 0.6%, were performed. Hence, additional fatigue tests on DB 800H with smaller and larger strain ranges would be needed to evaluate the proposed fatigue curves. Similarly, additional creep–fatigue tests can facilitate the development and evaluation of the creep-fatigue envelope of DB 800H. Note that all experiments conducted in this study were focused on solid ASTM specimens. This study does not focus on the influence of channels. The presence of channels may change the local stress flow during the diffusion bonding process, which may result in different bond properties than the solid diffusion bonded regions. A series of tests at elevated temperatures are needed on diffusion bonded specimens with channels to address the effect of stress concentrations near channel corners on the proposed allowable stresses and design life presented in this study.

Acknowledgment

This research is being performed using funding received from the DOE Office of Nuclear Energy's Nuclear Energy University Program under an Integrated Research Program entitled “ASME Code Application of the Compact Heat Exchanger for High Temperature Nuclear Service” (NEUP-16-10714) with the award number DE-NE0008576, and “Advancements toward ASME Nuclear code case for compact heat exchangers” (IRP-17-14227) with the award number DE-NE0008714. Authors would like to acknowledge Mr. Aaron Wildberger from Vacuum Process Engineering for his help with the diffusion bonded Alloy 800H block fabrication.

Funding Data

  • Office of Nuclear Energy (Grant Nos. DE-NE0008576 and DE-NE0008714; Funder ID: 10.13039/100006147).

References

1.
U.S. DOE OSTI,
2002
, “
A Technology Roadmap for Generation IV Nuclear Energy Systems
,” U.S. Department of Energy Office of Scientific and Technical Information, Oak Ridge, TN, accessed Mar. 21, 2022,
https://www.osti.gov/servlets/purl/859029
2.
Collins
,
J.
,
2009
, “
Next Generation Nuclear Plant Project Technology Development Roadmaps: The Technical Path Forward
,” Idaho National Laboratory, Idaho Falls, ID, Report No. INL/EXT-08-15148.
3.
Sabharwall
,
P.
,
Kim
,
E. S.
,
Siahpush
,
A.
, and
Patterson
,
M.
,
2012
, “
Preliminary Design for Conventional and Compact Secondary Heat Exchanger in a Molten Salt Reactor
,”
ASME
Paper No. HT2012-58007.10.1115/HT2012-58007
4.
Ravindran
,
P.
,
Sabharwall
,
P.
, and
Anderson
,
N. A.
,
2010
, “
Modeling a Printed Circuit Heat Exchanger With REPAL5-3d for the Next Generation Nuclear Plant
,” Idaho National Laboratory, Idaho Falls, ID, Report No.
INL/EXT-10-20619
.https://inldigitallibrary.inl.gov/sites/sti/sti/4731815.pdf
5.
Oh
,
C. H.
, and
Kim
,
E. S.
,
2008
, “
Design Option of Heat Exchanger for the Next Generation Nuclear Plant
,”
ASME
Paper No. HTR2008-58175.10.1115/HTR2008-58175
6.
Mylavarapu
,
S. K.
,
Sun
,
X.
,
Christensen
,
R. N.
,
Unocic
,
R. R.
,
Glosup
,
R. E.
, and
Patterson
,
M. W.
,
2012
,
Fabrication and Design Aspects of High-Temperature Compact Diffusion Bonded Heat Exchangers
, Vol.
249
,
Elsevier
, Amsterdam, The Netherlands, pp.
49
56
.
7.
Sabharwall
,
P.
,
Clark
,
D. E.
,
Mizia
,
R. E.
,
Glazoff
,
M. V.
, and
McKellar
,
M. G.
,
2013
, “
Diffusion-Welded Microchannel Heat Exchanger for Industrial Processes
,”
ASME J. Therm. Sci. Eng. Appl.
,
5
(
1
), p.
011009
.10.1115/1.4007578
8.
Li
,
X.
,
Kininmont
,
D.
,
Le Pierres
,
R.
, and
Dewson
,
S. J.
,
2008
, “
Alloy 617 for the High Temperature Diffusion-Bonded Compact Heat Exchangers
,”
Proceedings of ICAPP
, Vol.
8
, Anaheim, CA, June 8–12, pp.
282
288
.https://www.heatric.com/app/uploads/2018/04/Alloy-617-for-the-high-temperature-diffusionbonded-compact-heat-exchangers.pdf
9.
Sah
,
I.
,
Kim
,
D.
,
Lee
,
H. J.
, and
Jang
,
C.
,
2013
,
The Recovery of Tensile Ductility in Diffusion-Bonded ni-Base Alloys by Post-Bond Heat Treatments
, Vol.
47
,
Elsevier
, Amsterdam, The Netherlands, pp.
581
589
.
10.
Rozman
,
K.
,
Carl
,
M.
,
Kapoor
,
M.
,
Doğan
,
Ö.
, and
Hawk
,
J.
,
2019
,
Creep Performance of Transient Liquid Phase Bonded Haynes 230 Alloy
, Vol.
768
,
Elsevier
, Amsterdam, The Netherlands, p.
138477
.
11.
Mahajan
,
H. P.
,
Lima
,
L. M. A.
, and
Hassan
,
T.
,
2022
, “
Mechanical and Microstructural Performance Evaluation of Diffusion Bonded Alloy 800H for Very High Temperature Nuclear Service
,”
J. Eng. Mater. Technol.
,
144
(
2
), p.
11
.10.1115/1.4052825
12.
ASME
,
2017
, “
ASME Boiler & Pressure Vessel Code Section III Division 5
,”
ASME
, New York, Standard.
13.
Keating
,
R.
,
Nestell
,
J.
,
McKillop
,
S.
,
Allen
,
T.
, and
Anderson
,
M.
,
2019
, “
Potential ASME Code Case for Construction of Compact Heat Exchangers in High Temperature Reactors
,” Proceedings of the International Conference on Nuclear Engineering (
ICONE
), The Japan Society of Mechanical Engineers, Tsukuba, Japan, May 19–24, p.
1138
.10.1299/jsmeicone.2019.27.1138
14.
ASME
,
2017
, “
ASME Boiler & Pressure Vessel Code, Section VIII Division 1
,”
ASME
, New York, Standard.
15.
ASME
,
2017
, “
ASME Boiler & Pressure Vessel Code, Section II Part D
,”
ASME
, New York, Standard.
16.
ASME
,
2016
,
Fitness-For-Service
,
ASME
, New York, Standard No. API 579-1/ASME FFS-1.
17.
Swindeman
,
R. W.
,
Swindeman
,
M. J.
,
Roberts
,
B. W.
,
Thurgood
,
B. E.
, and
Marriott
,
D. L
,
2007
, “
Verification of Allowable Stresses in ASME Section III Subsection NH For Grade 91 Steel & Alloy 800H
,” ASME, New York.10.2172/974278
18.
Swindeman
,
R. W.
,
Zamrik
,
S. Y.
, and
Maziasz
,
P. J.
,
2007
, “
Effects of Long-Term Service on the Microstructure and Tensile Properties of Alloy 800H
,”
ASME
Paper No. PVP2007-26167.https://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.1079.7874&rep=rep1&type=pdf
19.
Abe
,
F.
,
Tabuchi
,
M.
, and
Hayakawa
,
M.
,
2017
, “
Influence of Data Scattering on Estimation of 100,000 Hrs Creep Rupture Strength of Alloy 617 at 700 c by Larson–Miller Method
,”
ASME J. Pressure Vessel Technol.
, 139(1), p.
011403
.10.1115/1.4033290
20.
Messner
,
M.
,
Jetter
,
R.
, and
Sham
,
T.-L.
,
2020
, “
A High Temperature Primary Load Design Method Based on Elastic Perfectly-Plastic and Simplified Inelastic Analysis
,”
ASME
Paper No. PVP2020-21470.10.1115/PVP2020-21470