Site icon Turkish Forum

Refractory Failure Elimination in Fludised Bed Reactors

İşletmelerin 2020 Yılında Faydalanabileceği Bazı Kilit Teknolojiler
Spread the love

Refractory Failure Elimination in Fludised Bed Reactors

Technical paper IRE – Refractory Failure Elimination (Download .doc)

Dr Cengiz Erkan

Principle Engineer BHB Billiton

PhD.MSc.BEng.FIEAustCPEng.MIREng.DBus

Dr Tom Honeyands

BHP Billiton – Newcastle Technology Centre

 

Dr Andrew Shook

BHP Billiton – Newcastle Technology Centre

 

Dr Philip Clausen

School of Engineering, University of Newcastle

  1. Executive Summary

 

Spalling of fluidised bed reactor refractories and subsequent feed leg blockages have been a significant cause of lost production in Iron, alumina and nickel processing plants. This work describes the Strategic Refractory Management Plan developed to address high refractory failure rates during the first years of operation of the Port Hedland HBI plant.

 

An example is then given of a mathematical modelling study of the heat transfer and stresses in the reactor vessel to see whether any conclusions can be drawn about the maximum allowable heating or cooling rates.

 

The reactor hot face refractory and insulation layers, restrained by a stainless steel shell, were modelled using “strand7” finite element software. Here only the dome of a simplified reactor, not including expansion joints, anchors, gas inlets and outlet, was modelled to determine the stress profile through the wall when subject to thermal loading.

 

Several alternative refractory configurations were considered for this application. The compressive thermal stresses predicted on heat-up were found to always be lower than the crushing strength of the materials in this investigation, and no conclusions can be made regarding maximum heating rate.

 

When the internal surface of the reactor is cooled by cold purge gasses, tensile stresses can potentially be generated in the hot face refractory, along with shear stresses between the hot face and the insulation. Tensile stresses were only predicted in one case – 150mm hot face / 150mm insulation – during rapid cooling. This configuration was therefore rejected.

 

Shear stresses in the other refractory configurations were found to be relatively low (<1.1MPa). No indications of the reasons for previous refractory failures were found using this simple model. It is concluded that differences in the dry-out and curing of the 60% Alumina and 80 % Alumina based castables hot face refractories is responsible for their differing performance.

 

  1. Introduction

 

In this paper we consider the example of BHP Billiton’s hot Briquetted iron (HBI), processing plant in Port Hedland due to problematic processing characters. In this plant the high alumina refractories must contain abrasive solids at elevated temperatures and pressures under reducing gas conditions. A combination of reasons accelerated the deterioration and shortened the life expectancy of the refractories during this period. The following examples are typical of the problems experienced:

 

 

The objectives of this paper are two fold:

  1. To outline the overall refractory management plan that was put in place and to demonstrate the economic benefits

  2. To describe a modelling study used to assist in troubleshooting the refractory failures and in assessing alternative refractory thickness options for one of the vessels including:

 

 

  1. The Port Hedland Example

 

Cracking, spalling and or even complete failure of reactor refractories causes significant lost of production, increased maintenance and increased operational costs. These refractory failures can be caused by:

 

 

To address this problem, a technical team supervised by Genghis Erkan developed a long term “Strategic Refractory Management Plan in BHP Billiton”, to systematically solve refractory and related problems. This plan provided systematic material selection and evaluation process and enabled us to manage and monitor such things as uniformity of heat transfer during start-up, operation and shut down.

 

This resulted in longer campaign periods, less equipment wear, shorter shut-downs and decreased overall operating costs. The resulting benefits can be seen in the later years of operation shown in Table 1. The corresponding increase in refractory life is shown in Table 2.

 

Table 2: Refractory Operational Life (Briquetting Feed Bin)

 

Year

Days

1999

70

2000

75

2001

60

2002

140

2003

220

 

To achieve these improvements, the following actions were taken as part of the Strategic Refractory Management Plan:

 

 

Figure 1 show refractory failures, which have decreased significantly since the full implementation of this strategy.

Figure 1.  Refractory failures from 1999 to 2003

 

  1. Mathematical Modelling

 

To assist with the modification to the refractory design, and to guide changes to start up and operating strategy, a number of mathematical models were used.  In conjunction with this, the failures of the refractory lining were extensively studied during post mortem investigations, Historical plant data were analysed to determine refractory performance relative to operating conditions, and in service tests were carried out.

 

4.1 Heat Transfer Analysis

The first step was to carry out an analysis of the heat transfer environment. A set of data were examined to give a good range of heat up and cool down rates actually experienced in the plant.

 

The temperature profiles through the refractory were calculated using a commercial heat transfer package, Heating, using the thermal properties of the refractories supplied by the plant technology department and assuming that the thermocouple present in the gas space of the reactor vessel was a good approximation to the hot face temperature, and applying external convective heat transfer. An example of the temperatures calculated through the refractory thickness is given in Figure 2.

Figure 2. Temperatures calculated for the 100mm hot face / 200mm insulation reactor configuration

 

When the reactor is cooled due to process upsets, the internal surface of the hot face can be cooled below the temperature at the hot face / insulation boundary (see Figure 3). This sub-cooling has the potential to generate tensile stresses in the hot face refractory and shear between the hot face and insulation. These stresses are calculated in the next section.

 

4.2 Stress Analysis using Finite Element Modelling

 

A finite element model of the dome of the reactor was created to determine the stress profile through the reactor wall when subject to the thermal loading calculated in Section 4.1. The reactor consists of a 20mm thick, 3740mm internal diameter stainless steel cylinder with a hemispherical dome cap. The internal surface of the shell is coated with insulation refractory with this coated with hot face refractory.  The “normal” configuration has 200mm of insulation and 100mm of hot face.

 

To analyse the thermally induced stress in the upper wall section of the reactor requires a full 3D model to be built and solved.  Ignoring the two nozzles in the reactor’s dome section, however, simplifies this structure to a solid of revolution.  Here a planar section through the wall can be modelled to represent the complete stress field in the structure as long as this section can be extruded 360 around an axis to create the structure.

 

Modelling the dome as a solid of revolutions also requires the material properties, loading and constraints to be axisymmetric in nature, that is, they are circumferentially independent.  Here we assumed all material properties to be isotropic with the lower edge of this section of the bin restrained from moving in the axial direction only.

 

Figure 4 shows the finite element model of the simplified dome of the reactor vessel. The model consists of just over 5000 low-order plate elements with the red, blue and green coloured elements modelling the hot-face refractory, insulation refractory and stainless steel respectively.  Nodes on the lower edge of the cylindrical section were constrained from moving in the axial (y) direction. The material properties used for the stainless steel, insulation and hot-face are given in Table 3.

Early modelling work showed that the internal surface of the hot face was always in compression for the entire temperature cycle. Here it was assumed that the zero stress state occurred at 40C, the ambient air temperature. This assumption, however, may be in error as the hot face refractory was cured at a temperature significantly above ambient.  As such, the temperature profile through the wall during the curing process is more likely to coincide with the zero stress state.  At this zero stress state, the modelling has indicated that tensile stress will occur on the internal surface of the hot face refractory.

 

The model was also solved where the zero stress states were set at both 550°C and 700°C to see if the stress profile through the wall was sensitive to its initial temperature state.

 

Table 3.  Material properties used in finite element model

 

Material

Young’s Modulus of Elasticity (GPa)

Poisson’s ratio

Density (kg/m3)

Thermal expansion coefficient mm/mm°C

Stainless steel

199.95

0.32

7, 806

1.72×10-5

Insulation

4.865

0.3

800

8.5×10-6

Hot Face

4.865

0.3

2100

8.5×10-6

  1. Finite element Analysis Results

 

  1. Compressive and Shear Stress During the Early Part of the Operating Cycle

 

Table 4 below shows the maximum compressive stress and the shear stress at the interface of the insulation refractory and hot face refractory during the heat-up period and into the early cool down part of the bin’s operating cycle.  As expected, the maximum stress for all times considered occurred on the plate element adjacent to the inside surface of the hot face refractory.  The shear stress rz, shown in Table 4 was the value at the centroid of the hot face element adjacent to the boundary.  As can seen, the value of the shear stress is small in comparison to the normal stresses.

 

Table 4.  Maximum compressive stress and shear stress at the hot face – insulation junction during the early part of the operating cycle.

 

Time (seconds)

r  MPa (max)

z MPa (max)

  MPa (max)

rz, Mpa

465

-26.1

3.9

-26.1

-0.003

10000

-34.5

0.49

-34.5

0.15

20000

-34.9

0.12

-34.9

0.18

80000

-31.7

-0.08

-31.7

0.17

166800

-30.0

-0.19

-30.0

0.16

200000

-20.28

-0.23

-20.28

0.12

300000

-8.58

-0.02

-8.58

0.05

 

  1. Thermal Stresses when the Zero Stress State is above Ambient

 

The following two zero stress state temperatures were considered: a wall temperature of 550C and a wall temperature of 700C.

Zero Stress State at a Wall temperature of 550C

 

Figure 3 and 4 show contours of circumferential stress for the 200,000s and 350,000s respectively with the wall temperature at 550C. As can be seen, figure 5 shows compressive stresses in the element adjacent to the inside surface of the hot-face with the other two recording tensile stress. The centroid stress is shown in Table 4.

Figure 3. Contours of stress at 200,000s with an internal temperature of 550C.

Figure 4. Contours of stress at 300,000s with an internal temperature of 550C.

 

Table 5.  Stress at the centroid of the plate indicated in figure 3 for the 550C zero stress state.

 

Time (sec)

r  MPa

z Mpa

  MPa

rz, MPa

167000

-6.08

-0.0023

-6.08

0.035

200000

3.56

-0.017

3.56

-0.042

300000

15.8

-0.073

15.8

-0.107

Zero Stress State at a Wall temperature of 700C

 

Table 6 documents the stress for the case 167,000s and 200,000s. Results for the 300,000s case was not considered, as these stresses will simply be higher than for the case at 200,000s.

 

Table 6.  Stress at the centroid of the plate indicated in figure 6 for the 700C zero stress state.

 

Time (sec)

r  MPa

z MPa

  MPa

rz, MPa

167000

1.41

-0.077

1.41

-0.018

200000

10.96

-0.137

10.97

-0.086

 

  1. Stresses in Different Refractory Thickness

 

As previously stated, the normal wall had 200mm of insulation refractory and 100mm of hot face.  Here two other configurations are considered: 150 mm of insulation and 150mm of hot face and 50mm insulation and 150mm of hot face. The zero stress state was assumed to be at 40C.  Temperature profiles through the wall were supplied.

 

The results show for the 150mm hot face/ 50mm insulation case the refractory is in compression at both 200,000s and 300, 000s as indicated in Table 7. The shear stress at the hot face/insulation junction is shown in the last column in the table.

 

For the 150mm hot face /150mm insulation case, however, the plate element adjacent to the inside surface is in tension as shown in Table 8. Here the shear stress at the hot face /insulation junction is of similar magnitude to that in the 150mm hot face/50mm insulation case. The tensile stress at the surface of the hot face suggests that this insulation/hot face configuration is undesirable.

 

Table 7.  Normal stresses at the centroid of the plate adjacent to the inside surface of the hot face and the shear stress at the hot face insulation boundary for 150mm hot face /50mm insulation case.

 

Time (sec)

r  MPa

z MPa

  MPa

rz, MPa (interface)

200000

-16.0

-0.095

-16.0

0.081

300000

-7.55

-0.045

-7.55

0.038

 

Table 8.  Normal stresses at the centroid of the plate adjacent to the inside surface of the hot face and the shear stress at the hot face insulation boundary for 150mm hot face /150mm insulation case.

 

Time (sec)

r  MPa

z MPa

  MPa

rz, MPa (interface)

200000

13.54

-0.017

13.55

-0.052

300000

7.4

-0.018

7.5

-0.028

 

  1. Discussion

 

The heat-up rate used in Figure 2 is extremely high, as the gas thermocouple has essentially been set instantaneously to 700C. The state of stress in the bin is generally compression in the refractory with the stainless steel shell in tension. Even at this extreme heating rate, the compressive stresses in the refractory do not exceed the crushing strength of either 60% Al2O3 or 80% Al2O3 refractories where the crushing strength is between 68 and 72MPa.

 

These stresses were calculated for the three refractory configurations:

 

 

During the thermal analysis a simple comparison was made between 60% Al2O3 and 80% Al2O3. Analysis of historical data by plant suggests that the performance of 60% Al2O3 has been superior to that of 80% Al2O3. A simple measure of thermal shock resistance, based on avoiding the initaition of fracture, can be calculated using the following equation.

The results of the comparison between 60% Al2O3 and 80% Al2O3 using Equation (1) are given in Table 9.

The termal shock resistance, R, is given by:

R’ = k..(1-) / E.a (1)

Where;

R is Thermal shock resistance

k is the thermal conductivity,

 is the failure stress,

E is Young’s modulus (value of  4.8 Gpa used),

a is the coefficient of thermal expansion,

and  is Poisson’s ratio (value of 0.3 used).

 

Table 9. Simple comparison of thermal shock resistance of 60% Al2O3 and 80% Al2O3.

 

60% Al2O3

80% Al2O3

K (W/m.K)

0.92

1.56

 (MPa)

6

16

a (mm/mmC)

8.5 x 10-6

7.0 x 10-6

R’

95

520

 

This comparison suggests that 80% Al2O3 should have better thermal shock resistance than 60% Al2O3 refractories, which is presumably why it was used. There is clearly some other reason for the differences in performance. One possibility is that differences may exist in the behaviour of the two refractories during dry out and curing. 80% Al2O3 is a wet gunning mix, whereas 60% Al2O3 is a dry gunning mix. Both are reportedly dried out using the same temperature ramps and holds (up to 550C).

 

This was explored using special stress model, by changing the temperature at which the refractory was at a “zero stress” state from room temperature to 550C. This assumes that the castable refractories form their bonds at the curing temperature of 550C. Heating above this temperature puts the refractory into compression as before, however, cooling below this temperature puts the internal surface of the refractory into tension and could lead to failure.

 

This hypothesis is likely to be generally correct, in that large thermal cycles are bad for the refractories and that the reactor should be kept at it’s operating temperature whenever possible. Other more quantitative predictions cannot be made with our current level of knowledge.

 

  1. Conclusions

 

The establishment of a Strategic Refractory Management Plan resulted in a significant increase in refractory life at the Port Hedland HBI plant, and a corresponding decrease in refractory management costs.

 

A mathematical modelling study was carried out as part of this plan. The finite element modelling indicated that if the zero stress state was assumed to be at an ambient temperature of 40C, then the stress in the refractory wall was compressive at all times during a typical operating cycle.  If, however, the zero stress state was set at the refractory curing temperature, then the stress within the refractory wall maybe tensile.

 

Furthermore, raising the refractory curing temperature can lead to tensile stress in a greater percentage of the operating cycle.  This led the refractory engineers to closely examine the drying and curing of refractories, and led to improved practices. The refractory wall with a 150 mm hot face/150 mm insulation configuration produces tensile stress on the inside surface of the hot face irrespective of the temperature of the zero stress state.  These results suggest that this thickness combination is inappropriate for a refractory vessel.

 

To achieve the best results for the plant a proper evaluation is best accomplished by complete cooperation and communication between refractory manufacturer, installer, and internal stakeholders. To continue to operate the plant without a “Strategic Refractory Management Plan” will prove to be very costly and unsustainable.

 

References

 

Genghis Erkan, BChEng, MscEnv, MIEAust, CPEng, MIREng “Fludised Bed Reactor Refractory Implementation Strategy, 2003”.

 

Dr P.D Clausen, BE, PhD, “A Finite Element Analysis of the Refractory Wall, 2001”.

 

Dr Andrew Shook, BHP Billiton Technology Internal report, “Thermal Stress Analysis of BFB Refractories, 2001”.

Dr Tom Honeyands, BHP Billiton Technology Internal Report, “Thermal Stress Analysis of BFB Refractories, 2001.


Spread the love
Exit mobile version