_{1}

^{*}

In this technical paper, the oxidation mechanism and kinetics of aluminum powders are discussed in great details. The potential applications of spherical aluminum powders after oxidation to be part of the surging arresting materials are discussed. Theoretical calculations of oxidation of spherical aluminum powders in a typical gas fluidization bed are demonstrated. Computer software written by the author is used to carry out the basic calculations of important parameters of a gas fluidization bed at different temperatures. A mathematical model of the dynamic system in a gas fluidization bed is developed and the analytical solution is obtained. The mathematical model can be used to estimate aluminum oxide thickness at a defined temperature. The mathematical model created in this study is evaluated and confirmed consistently with the experimental results on a gas fluidization bed. Detail technical discussion of the oxidation mechanism of aluminum is carried out. The mathematical deviations of the mathematical modeling have demonstrated in great details. This mathematical model developed in this study and validated with experimental results can bring a great value for the quantitative analysis of a gas fluidization bed in general from a theoretical point of view. It can be applied for the oxidation not only for aluminum spherical powders, but also for other spherical metal powders. The mathematical model developed can further enhance the applications of gas fluidization technology. In addition to the development of mathematical modeling of a gas fluidization bed reactor, the formation of oxide film through diffusion on both planar and spherical aluminum surfaces is analyzed through a thorough mathematical deviation using diffusion theory and Laplace transformation. The dominant defects and their impact to oxidation of aluminum are also discussed in detail. The well-controlled oxidation film on spherical metal powders such as aluminum and other metal spherical powders can potentially become an important part of switch devices of surge arresting materials, in general.

Oxidation of aluminum to form aluminum oxide film has been playing significant roles in modern industry and receiving wide applications. Aluminum hard anodization has been widely using in semiconductor industry as chamber surface coating. The anodized aluminum through different oxidizing processes at low temperature in acidic solutions has demonstrated excellent performance in corrosion resistance for semiconductor plasma etching equipment from 200 mm silicon wafer fabrication to current 300 mm silicon wafer fabrication worldwide [

As an n-type of semiconductor, aluminum oxide (Al_{2}O_{3}), is a metal excess oxide. The oxidation mechanism and the formation of metal oxide on metal surface at high temperature are summarized by Jones in his book [_{ox}/V_{metal}) [_{2}O_{3} is less than aluminum metal, the oxide layer can cover the entire surface of aluminum, and the aluminum oxide is very protective. Typical PRB and the properties of metal oxides are summarized in

Ellingham diagram has been used to show the relationship of Gibbs free energy versus temperature for oxidation of metals [

In order to make the oxide film on metals as switch material, one has to control the growth rate and the uniformity of the oxide film on metals very well. It means that the thickness of oxide on metals and the uniformity of thickness of oxide layer have to be controlled very well during oxidation. The major properties are the thickness of oxide layer and the uniformity of the oxide film. Taking consideration of aluminum powder oxidation, the selection of aluminum powders and the methodology to form the aluminum oxide film are very important.

Metal | Oxide | PR Ratio | Protectiveness | Type of Oxide |
---|---|---|---|---|

Aluminum | Al_{2}O_{3} | 1.28 | P | n |

Calcium | CaO | 0.64 | NP | n |

Cadmium | CdO | 1.42 | NP | n |

Cobalt | Co_{2}O_{3} | 2.40 | P | p |

Copper | Cu_{2}O | 1.67 | P | p |

Chromium | Cr_{2}O_{3} | 2.02 | P | p |

Iron | FeO | 1.78 | P | p |

Magnesium | MgO | 0.81 | P | n/p |

Manganese | MnO_{2} | 2.37 | P? | n |

Molybdenum | MoO_{3} | 3.27 | NP | n |

Nickel | NiO | 1.70 | P | p |

Lead | PbO | 1.28 | NP | p |

Silicon | SiO_{2} | 2.15 | P | n |

Tantalum | Ta_{2}O_{5} | 2.47 | NP | n |

Titanium | Ti_{2}O_{3} | 1.76 | NP | n/p |

Uranium | UO_{2} | 1.97 | NP | p |

Tungsten | WO_{3} | 1.87? | NP | n |

Zinc | ZnO | 1.58 | NP | n |

Zirconium | ZrO_{2} | 1.57 | P | n |

P: Protective; NP: Not Protective; n: n-type semiconductor, metal excess oxide; p: p-type semiconductor, metal deficit oxide from Jones [

The spherical shape of aluminum powders in a gas phase fluidized bed can help to produce a uniform aluminum oxide film. A fluidized bed has been used for gas and solid particles interaction for many years. The term FLUIDIZATION is used to designate the gas-solid contacting process in which a bed of finely divided solid particles is lifted and agitated by a rising stream of process gas. At the lower end of the velocity range, the amount of lifting is slight, the bed behaving likes a boiling liquid. At the other extreme, the particles are fully suspended in the gas stream and are carried along with it. The terms of suspension, suspensoid, and entrainment contact have all been used to designate this action.

The fluidized technique was born from the pioneering work of Standard Oil Development Company and was introduced commercially in 1937. The development of the technique was summarized by Zenz [

1) Excellent temperature control to reach isothermal conditions.

2) Continuity of operation.

3) Heat transfer.

4) Catalytic gas reactions.

5) Aeratable powders have the most desirable fluidization properties.

6) Interparticle forces in the particles are significantly smaller than hydrodynamic forces and gas bubbles are limited in size.

The heat transfer, minimum fluidization velocity, fluidization technology, and packing properties have been discussed by various authors [

Although gas fluidization bed has been using in different applications for many years, there is no theoretical modeling available. It means that the oxide formation relies on experimental results. The author needs to develop the theoretical model based on the dynamic system with specific hydrodynamic flow conditions. The current work is to create and to build up a theoretical mathematic modeling based on the real dynamic system and system parameters to estimate the thickness of oxide layer in a gas fluidization bed. Therefore, one can carry on the theoretical calculation and compare the results with experiments. In fact, the mathematical model developed in this work can help to carry out quantitative calculation of oxide growth rate at different conditions and compare the calculation with experimental results. Taking aluminum ball powder as an example, the author successfully created and developed the mathematical modeling and applied the model to compare with experimental results on a gas fluidization bed [

1) Basic parameters used in the kinetic calculation of a gas fluidization bed reactor:

π = 3.1416.

S_{sur} = πD^{2} = surface area of aluminum ball with diameter, D.

D_{in} = Inner diameter of bed = 70 mm = 7.0 cm = 0.070 m.

D_{out} = Outer diameter of bed = 75 mm = 7.5 cm = 0.075 m.

Bed wall thickness = 2.5 mm = 0.25 cm = 0.0025 m.

Area of cross section of bed = π D i n 2 / 4 = 38.48 cm^{2}.

Height of bed = 9.0 inch = 229.4 mm = 22.94 cm = 0.229 m.

Total volume of bed = 883 cm^{3}.

Material of fluidization bed is quartz.

Thickness of the gas distributor = 5 mm = 0.5 cm = 0.005 m.

Hole size of the gas distributor = 10 μm - 15 μm.

Number of holes on the gas distributor ≈ 1000.

2) Parameters and bed operation conditions:

Volume flow rate of dry air = 8.333 × 10^{−}^{5} m^{3}/sec = 83.33 cm^{3}/sec.

Weight of aluminum powders = 1.5 pounds.

Average diameter of aluminum powders = 20.0 μm.

The velocity of dry air = 2.16 cm/sec.

The height of air flowing = 130.3 cm/min.

Particle size distribution of aluminum powder is between 10 μm and 100 μm.

Temperature is between 25˚C and 400˚C.

Density of aluminum particles = 2700 kg/m^{3}.

Bed bulk density = 1.414 kg/m^{3}.

Bed void before fluidization = 0.4764.

Thermal conductivity of air at 100˚C = 0.0286 W・m^{−1}・K^{−1}.

Thermal conductivity of air at 200˚C = 0.0370 W・m^{−1}・K^{−1}.

Specific heat of quartz = 0.80 kJ・kg^{−1}・K^{−1}.

Specific heat of air = 1.03 kJ・kg^{−1}・K^{−1}.

Specific heat of aluminum powders = 0.982 kJ・kg^{−1}・K^{−1}.

Air flow rate in kg/sec. = 5.80 × 10^{−5} kg/sec.

Density of air at 200˚C = 0.696 kg/m^{3}.

Density of air at 100˚C = 0.927 kg/m^{3}.

Height of gently settled aluminum powders = 12.5 cm = 0.125 m.

3) Basic calculation of a gas fluidized bed reactor

In order to carry out the calculation of the fluidization bed and the temperature-dependent relations of parameters in fluidization bed reactor (FBR), Shih wrote a general computer program 1 called “FLOW.EXE”. This computer program was used to carry on the complete calculation of all important parameters in a FBR. It generates a very large database which can be used as a handbook of a fixed FBR. It basically carries on nine calculations. The particle size is between 10 μm and 100 μm, and the FBR temperatures are between 25˚C and 480˚C in the calculation. In order to reduce the generated database, a software program 2 called “FLOW1.EXE” was written by Shih [

Details of calculations of “FLOW.EXE” and FLOW1.EXE” are listed below:

a) Average particle size vs. number of particles.

b) Average particle size vs. effective surface area.

c) Average viscosity of gases vs. temperature.

d) Effective gas density and specific heat of air vs. temperature.

e) The height of the gentle settled bed.

f) Relationship between diameter of aluminum powder and H_{max}.

g) Relations of particle size, Archimedes number, and minimum fluidization velocity.

h) Thermal Conductivity of Gases at Different Temperatures.

1) Diameter, particle number, and surface area of spherical aluminum powders

The relations of average particle diameter, particle number, and effective surface area of aluminum powders are calculated and shown in ^{3}. The average spherical aluminum particle size is between 10 μm and 100 μm. The effective surface area is in m^{2}. With the increase of average particle diameter (d_{v}), both particle number (N_{p}) and effective surface area (S_{p}) decrease (see the table of nomenclature for the definitions of parameters).

Diameter in mm (d_{v}) | Particle Number (N_{p})* | Surface Area in m^{2} (S_{p})^{#} |
---|---|---|

10 | 1.000E+15 | 314,160 |

20 | 1.250E+14 | 157,080 |

30 | 3.330E+13 | 94,154 |

40 | 1.560E+13 | 78,414 |

50 | 8.000E+12 | 62,832 |

60 | 4.700E+12 | 53,156 |

70 | 2.900E+12 | 44,642 |

80 | 1.900E+12 | 38,202 |

90 | 1.400E+12 | 33,081 |

100 | 1.000E+12 | 31,416 |

^{+}: total volume of aluminum powder = 1.0 m^{3}; *: N_{p} = 1.0 × 10^{18}/(d_{v})^{3}; ^{#}: S_{p} = πN_{p} × (d_{v})^{2} × 10^{−12}.

2) The density and specific heat of air at various temperatures

The density and specific heat of air at different temperature is calculated and summarized by Shih [

3) Viscosity of air, argon, and oxygen at various temperatures

The viscosity of air, argon, and oxygen at different temperatures is calculated and summarized in

4) Calculation of heat transfer

One of the principal attractions of fluidized bed is the high rate of heat removal and addition. There are three modes of heat transfer which are of interest in gas fluidized bed.

・ Gas to particle.

・ Particle to particle.

・ Particle to heat transfer surface.

For particles less than 1mm diameter, the overall heat transfer rate between a fluidized bed and gas entering through the distributor plater is so large that the gas attains the bed temperature within a few centimeters. For particles less than 100 μm diameter, the overall heat transfer rate between a fluidized bed and gas entering through the distributor plate is even larger that the gas attains the bed temperature within 1 centimeter. Since the high thermal conductivity of aluminum powders, the heat transfer among aluminum particles can be negligible. In fact, aluminum powders are pre-heated in argon environment and reach the temperature of the fluidization bed. Therefore, the surface-to-bed heat transfer coefficient can reach a maximum value α_{max} and it can maintain over a wide range of velocity above the U_{mf}. By considering the masked surface by bubbles, the Zabrodsky correction for α_{max} is shown below:

α max = 35.8 K g 0.6 ρ p 0.2 / d p 0.36 (1)

For air at 200˚C

T˚K | Air Density (kg/m^{3}) | Specific Heat (kJ/kg・K) |
---|---|---|

100 | 3.5560 | 1.0320 |

200 | 1.7460 | 1.0070 |

300 | 1.1610 | 1.0070 |

373 | 0.9270 | 1.0100 |

500 | 0.6960 | 1.0300 |

1000 | 0.3480 | 1.1410 |

Temperature (˚C) | μ_{air}* | μ_{argon} | μ_{oxygen} |
---|---|---|---|

25 | 18.0 | 22.00 | 20.00 |

50 | 19.5 | 23.50 | 21.50 |

100 | 21.8 | 26.50 | 24.00 |

150 | 24.0 | 30.00 | 26.50 |

100 | 26.0 | 33.00 | 29.00 |

250 | 28.0 | 35.50 | 31.50 |

300 | 30.0 | 38.00 | 33.50 |

350 | 32.0 | 40.50 | 35.00 |

400 | 34.0 | 44.00 | 36.50 |

*: The unit of μ is (kg/m.s.) ×10^{6} and μ_{air} = {17.4 + 0.0420 × T(˚C)}.

α max = 1182 W ⋅ m − 2 ⋅ K − 1 (2)

The overall heat transfer coefficient, α_{o}, is shown as

1 / α o = 1 / α max + x / k q z + 1 / α f (3)

and

α o = 145 . 41 W ⋅ m − 2 ⋅ K − 1 (4)

Assuming the heat height of aluminum powder is about 22.0 centimeter, the total heat through the wall to the fluid of the bed can be calculated as

Q = α o A m ( T b − T f ) (5)

For a given fluidized bed from 25˚C to 200˚C as an example, the total heat is

Q b = 1231 W (6)

Assuming that 1.5 pounds of aluminum powders have been pre-heated in argon environment to 200˚C, the total heat to raise air from 25˚C to 200˚C as an example at a constant flow rate is shown below:

Q a i r = C a i r ( T b − T f ) V a i r (7)

and

Q a i r = 10.45 W (8)

The specific heat of aluminum powders at 200˚C is calculated as

C A l = 20.7 + 12.310 T (9)

and

C A l = 0.982 kJ ⋅ kg − 1 ⋅ K (10)

If aluminum powders are injected to the bed in a feed rate f_{Al}, the heat needed for aluminum powers can be calculated as

Q A l = C A l ( T b − T f ) f A l (11)

Since the volume of flow rate of air is fixed as

J a i r = 8.33 × 10 − 5 m 3 / sec (12)

Therefore, one can calculate the velocity of air

( 1 / 4 ) ( π D i n 2 ) V = 8.33 × 10 − 5 m 3 / sec (13)

and

V = 2.16 cm / sec (14)

From Equation (14), the calculated air velocity equals to the actual air velocity.

Now we should consider the required surface area of the wall to raise air from 25˚C to 200˚C. Here we have

Q a i r = 10.45 W = α max A a i r Δ T = 1185 × 175 × A a i r (15)

and

A a i r = 5.052 × 10 − 5 m 2 (16)

Since the circumference of the reactor is

C b e d = 22.94 cm = 0.229 m (17)

The heat length of air equals

H a i r = A a i r / C b e d = 22.0 × 10 − 4 m = 0.022 cm (18)

The conclusion is that for particles with diameter less than 100 μm, air can attain the bed temperature as soon as it enters the bed. By increasing the flow rate by 100 times, A_{air} increases 100 times and the heat length equals 2.2 centimeters. Hair is function of D_{in}, T, d_{p}, V, and J_{air}.

5) Maximum bed height

In order to keep the fluidized bed in bubble condition without slugging, the maximum bed height (H_{max}) and the inner diameter of the bed should satisfy the following relation:

( H max / D i n ) < 1.9 / { ρ p d p } 0.3 (19)

For a given fluidized bed with 0.070 m inner diameter, the maximum allowable depth for 100 μm to 10 μm diameter aluminum powders, the maximum bed depth should be less than 19.7 cm and 39.3 cm, respectively (see _{max} is a function of

Diameter of aluminum powder in mm | H_{max} in meter* |
---|---|

10 | <0.39 |

20 | <0.32 |

30 | <0.28 |

40 | <0.26 |

50 | <0.24 |

60 | <0.23 |

70 | <0.22 |

80 | <0.21 |

90 | <0.20 |

100 | <0.20 |

*: Equation (19) is used to estimate H_{max}.

D_{in}, ρ_{p}, and d_{p}. The complete values of H_{max} are shown below:

The relation will be used for the comparison.

π D i n L = 0.0484 m 2 (20)

and

L = 22.0 cm (21)

The minimum heated length from Equation (21) is 22.0 cm and the maximum bed depth should be less or close to L to avoid slugging.

6) The Maximum bubble diameter

The maximum bubble diameter at the surface of a bed can be calculated as

D b d = ( 2 / g ) ( V f f ) 2 (22)

Assuming the free fall velocity of the particles equals 0.1 m/sec., the maximum bubble diameter is shown below:

D b d = 0.2 cm ≪ 3.5 cm (23)

7) The Average particle size

The average particle size can be calculated by Equation (24)

d p = 1 / Σ ( x i / d i ) (24)

d_{p} is taken as 20 μm in the calculation.

8) The Height of the gently settled bed

The gently settled bed is calculated as follows:

1.5 × 454 / 2.7 = 252.22 cm 3 (25)

π D i n 2 / 4 = 38.48 cm 2 (26)

H s e t = ( 252.22 / 38.48 ) × 1 / ( 1.0 − ε v o i d ) (27)

and

H s e t = 12.52 cm (28)

The pressure of the bed is

Δ P = 1.414 × 12.52 × 0.0981 = 1.736.7 N ⋅ m − 2 = 1.737 × 10 − 3 bar (29)

The mass of the bed is

M b = ρ b s H s e t π D i n 2 / 4 = 0.681 kg (30)

9) Archimedes number and minimum fluidization velocity

The basic equation to calculate Archimedes number, A_{r}, is as below:

A r = ρ a i r d v 3 ( ρ p − ρ f ) g / μ 2 (31)

and the minimum fluidization velocity is

U m f = [ μ / ( ρ a i r d v ) ] ⋅ { [ 1135.7 + 0.0408 A r ] 1 / 2 − 33.7 } (32)

For particles with 20 μm diameter and at 200˚C, A_{r} is given as

A r = 0.2464 (33)

and

U m f = 0.10 cm / sec (34)

There are many factors which effect Archimedes number and minimum fluidization velocity. _{v} increasing, both A_{r} and U_{mf} increase. For the same particle size with temperature increasing, both d_{v} and U_{mf} decrease. _{v}, A_{r}, and U_{mf} at 100˚C and 200˚C, respectively [_{r} and U_{mf} values for different d_{v}, respectively [

10) Thermal conductivity of air, argon, and oxygen

Particle Size (d_{v}) in mm | Archimedes Number (A_{r})^{+} | Minimum Velocity (U_{mf}*) |
---|---|---|

10 | 0.0516 | 0.05 |

20 | 0.4130 | 0.10 |

30 | 1.3937 | 0.23 |

40 | 3.3037 | 0.51 |

50 | 6.4500 | 0.95 |

60 | 11.1456 | 1.62 |

70 | 17.6988 | 2.56 |

80 | 26.4192 | 3.80 |

90 | 37.6164 | 5.39 |

100 | 51.6000 | 7.38 |

^{+}: A r = ρ g d v 3 ( ρ p − ρ f ) g / μ 2 . *: U m f = μ / ( ρ a i r d v ) ⋅ { [ 1135.7 + 0.0408 A r ] 1 / 2 − 33.7 } .

Particle Size (d_{v}) in mm | Archimedes Number (A_{r})^{+} | Minimum Velocity (U_{mf}*) |
---|---|---|

10 | 0.0308 | 0.05 |

20 | 0.2464 | 0.10 |

30 | 0.8314 | 0.21 |

40 | 1.9712 | 0.44 |

50 | 3.8500 | 0.81 |

60 | 6.6528 | 1.37 |

70 | 10.5644 | 2.14 |

80 | 15.7696 | 3.17 |

90 | 22.4532 | 4.49 |

100 | 30.8000 | 6.14 |

^{+}: A r = ρ g d v 3 ( ρ p − ρ f ) g / μ 2 . *: U m f = μ / ( ρ a i r d v ) ⋅ { [ 1135.7 + 0.0408 A r ] 1 / 2 − 33.7 } .

The thermal conductivity of air, argon, and oxygen is calculated in

After the basic calculation, we can start to develop the dynamic model of velocity of pressure in the defined gas fluidization bed. The geometry of the fluidized bed is shown in ^{3}/s. There are n holes on the gas distributor, the radius of the hole is R. Assuming the particles are gently settled and when the bed is fluidized, the gas flow from the distributor forms n flow pipes through the bed with radius of the pipe equals to R. We also assume that the flow is evenly distributed among the n pipes. Therefore, the flow rate in each pipe is

q = Q / n (35)

From Equation (12), one has

Q = 8.33 × 10 − 5 m 3 / s (36)

Gas at | 100^{˚}K | 200^{˚}K | 300^{˚}K | 400^{˚}K | 500^{˚}K | 600^{˚}K |
---|---|---|---|---|---|---|

Air | 0.0094 | 0.0184 | 0.0262 | 0.0333 | 0.0397 | 0.0457 |

Argon | 0.0062 | 0.0124 | 0.0179 | 0.0226 | 0.0268 | 0.0306 |

Oxygen | 0.0093 | 0.0184 | 0.0263 | 0.0337 | 0.0410 | 0.0481 |

Nitrogen | 0.0098 | 0.0187 | 0.0260 | 0.0323 | 0.0383 | 0.0440 |

^{+}: Thermal conductivity in unit W・m^{−1}・K^{−1}.

If n = 1000, then we have

q = 8.33 × 10 − 8 m 3 / s (37)

1) Reynolds Number

Since R ≈ 10^{−5} m, velocity of air = 0.022 m/s and assume that the flow has a similar viscosity as air, at 200˚C, μ_{air} = 26.0 × 10^{−6} kg/ms, the Reynolds number is

R e = ( v R ) / μ = ( 0.022 × 10 − 5 ) / ( 26.0 × 10 − 6 ) ≈ 0.01 (38)

Since R_{e} is very small, the flow is laminar flow.

2) Governing Equations and Boundary Conditions

We chose cylindrical and Spherical coordinate system as illustrated in _{1}). Assume the flow is steady and incompressible, i.e., change of gas density is negligible, the basic equations for gas flow in the cylindrical and spherical coordinates [

Continuity:

( ∂ V r / ∂ r ) + ( 1 / r ) ( ∂ V θ / ∂ θ ) + ( ∂ V z / ∂ z ) + ( V r / r ) = 0 (39)

R-direction momentum equation:

V r ( ∂ V r / ∂ r ) + ( V θ / r ) ( ∂ V r / ∂ θ ) + V z ( ∂ V z / ∂ z ) − V θ 2 / r = g r − ( 1 / ρ ) ( ∂ P / ∂ r ) + ( μ / ρ ) ∇ 2 V r − ( V r / r 2 ) − ( 2 / r 2 ) ( ∂ V θ / ∂ θ ) (40)

where

∇ 2 V r = ( 1 / r ) { ( ∂ / ∂ r ( r ( ∂ V r / ∂ r ) ) ) } + ( 1 / r 2 ) ( ∂ 2 V r / ∂ θ 2 ) + ( ∂ 2 V r / ∂ z 2 ) (41)

Equation (41) is the Laplace operator of V_{r}.

θ-direction momentum equation:

V r ( ∂ V θ / ∂ r ) + ( V θ / r ) ( ∂ V θ / ∂ θ ) + V z ( ∂ V z / ∂ z ) + ( V r V θ / r ) = g θ − { 1 / ( ρ r ) } ( ∂ P / ∂ θ ) + ( μ / ρ ) ∇ 2 V θ − ( V θ / r 2 ) − ( 2 / r 2 ) ( ∂ V r / ∂ θ ) (42)

z-direction momentum equation:

V z ( ∂ V r / ∂ z ) + ( V θ / r ) ( ∂ V θ / ∂ r ) + V z ( ∂ V r / ∂ z ) = g z + ( μ / ρ ) ∇ 2 V z − ( 1 / ρ ) ( ∂ P / ∂ z ) (43)

Ñ^{2}V_{θ} and Ñ^{2}V_{z} are the Laplace operators of V_{θ} and V_{z}, respectively.

Since the radius of the pipe is very small and the bed height is much larger than the pipe’s radius, it is reasonable to say that the dominant flow velocity component is along Z direction. The rest of two components is negligible. The flow in the pipe is also axisymmetric. There two conditions can be written as

V θ = V r = 0 (44)

( ∂ / ∂ θ ) = 0 (45)

Substitute Equations (44) and (45) into continuity Equation (39) and momentum Equations (40)-(43), we have the following relations obtained

∂ V r / ∂ θ = ∂ V r / ∂ z = 0 (46)

g r − ( 1 / ρ ) ( ∂ P / ∂ r ) = 0 (47)

g θ − { 1 / ( ρ r ) } ( ∂ P / ∂ θ ) = 0 (48)

g z − ( 1 / ρ ) ( ∂ P / ∂ z ) + ( μ / ρ ) ( 1 / r ) ( ∂ ( r ( ∂ V r / ∂ r ) ) / ∂ r ) = 0 (49)

From Equation (46)

V z = V z ( r ) (50)

g r = 0 (51)

g θ = 0 (52)

g z = 0 (53)

where g is the gravity.

Substitute Equations (51)-(53) into Equations (47)-(49), one obtains

∂ P / ∂ r = 0 (54)

∂ P / ∂ θ = 0 (55)

g z − ( 1 / ρ ) ( ∂ P / ∂ z ) + ( μ / ρ ) ( 1 / r ) ( ∂ ( r ( ∂ V r / ∂ r ) ) / ∂ r ) = 0 (56)

From Equations (54) and (55), we know that P is only a function of Z, therefore, Equation (56) can be written as

( 1 / r ) ( d ( r ( d V r / d r ) ) / d r ) = ( 1 / μ ) ( d P / d z ) (57)

The boundary condition at the pipe wall is

V z | r = R = 0 (58)

3) Solution of the Equations

Integration of Equation (57), on obtains

r ( d V z / d r ) = [ ( 1 / μ ) ( d P / d z ) − ρ g / μ ] ( r 2 / 2 ) + C 1 (59)

where C_{1} is an integration constant and the density ρ is assumed to be constant.

Integrate Equation (59) again, we have

V z = [ ( 1 / μ ) ( d P / d z ) − ρ g / μ ] ( r 2 / 4 ) + C 1 ln r + C 2 (60)

where C_{2} is another integration constant.

At r = 0 (pipe axis), V_{z} is finite, therefore, using Equation (59)

C 1 = 0 (61A)

0 = [ ( 1 / μ ) ( d P / d z ) − ρ g / μ ] ( R 2 / 4 ) + C 2 (61B)

Therefore, one has

C 2 = − ( R 2 / 4 ) [ ( 1 / m ) ( ∂ P / ∂ z ) − ρ g / μ ] (61C)

Substitute Equation (61C) into Equation (60), one has

V z = [ ( r 2 − R 2 ) / 4 μ ] ( ( d P / d z ) − ρ g ) (62)

In the pipe, maximum velocity occurs at r = 0

V z , maximum = − ( R 2 / 4 μ ) ( ( d P / d z ) − ρ g ) (63)

To calculate the mean velocity, we need to calculate flow rate q from the V_{z} distribution.

q = ∫ 0 R V z 2 π r d r = ∫ 0 R [ ( r 2 − R 2 ) / 4 μ ] ( ( d P / d z ) − ρ g ) 2 π r d r = − ( π R 4 / 8 μ ) [ ( d P / d z ) − ρ g ] (64)

Mean velocity is shown below:

V z , mean = q / π R 2 = − ( R 2 / 8 μ ) [ ( d P / d z ) − ρ g ] (65)

It is interesting to note that

V z , mean = 1 2 V z , maximum (66)

Equation (65) can be rewritten as

d P / d z = − ( 8 μ / R 2 ) V z , mean + ρ g (67)

Integrate Equation (67), we then have

P − P 1 = − ( 8 μ / R 2 ) V z , mean ( Z − Z 1 ) + g ∫ z 1 z ρ d z (68)

where P_{1} is the pressure at the inlet of the pipe Z_{1}.

We assume the variation of gas density inside the bed is negligible as discussed before, i.e.,

ρ = constant (69)

Thus, Equation (68) becomes

P 1 − P = ( 8 μ / R 2 ) V z , mean ( Z − Z 1 ) + g ρ ( Z − Z 1 ) (70)

From Equation (70), one can see that the pressure of air depends on many factors. The factors are T, μ, ρ, Z, q, R, and V_{z}_{,mean}. At a given temperature and for a fixed bed, T, μ, ρ, Z, q, and R are constants, the pressure difference is only a function of Z. The concentration of oxygen is proportional to the pressure at a Z position. Therefore, the bulk concentration of oxygen at different Z positions can be estimated using Equation (70).

1) Calculation of air pressure in FBR.

We have the following relation as shown in Equation (70).

P 1 − P = ( 8 μ / R 2 ) V z , mean ( Z − Z 1 ) + g ρ ( Z − Z 1 ) (70)

Assuming P 1 = 1 atm = 1.01325 × 10 5 N / m 2 = 1.01325 × 10 5 kg / m ⋅ s 2 ;

Z 1 = 0.0 m (inlet of the FBR).

where N = m・kg/s^{2} and

Equation (70) becomes

P 1 − P = ( 8 μ / R 2 ) V z , mean Z + g ρ Z (71)

Since ρ is the fluid density which equals to the density of gas phase, V_{z}_{,mean} is determined by the overall flow rate in FBR and the inner diameter of the bed, D_{in}, Z is limited by the maximum bed depth (H_{max}), dynamic viscosity μ and fluid density ρ can be calculated and obtained from Program “FLOW.EXE” and FLOW1.EXE, so as to V_{z}_{,mean} (V_{ave}) and H_{max}, P can be calculated using Equation (71).

Since ρ and μ are function of temperature, therefore, P is function of μ, ρ, V_{z}_{,mean}, Z, T, and R and R is the radius of flow pipe. Z is also a function of particle size. We can express these relations as follows (see Equation (72)):

P = P 1 − ( 8 μ ( T ) / R 2 ) × Z × V z , mean ( Q , D i n ) × Z ( R ) − g ρ ( T , C o m p ) × Z ( R ) (72)

where Q is the overall flow rate entering the FBR, T is the temperature, and Comp is the composition of gas phase in FBR.

From Equation (72), one can see that P decreases linearly with Z increasing and the typical calculation is shown below when the following parameters are selected in the calculation:

T = 200 ˚ C ; Q = 5.80 × 10 − 5 kg / m 3 ; D i n = 0.070 m ; r o 1 = 10 μm = 1.0 × 10 − 5 m ; r o 2 = 60 μm = 6.0 × 10 − 5 m ; r o 3 = 100 μm = 1.0 × 10 − 4 m ; Comp = air ( fixed ) ; D g a s = D a i r = 0.696 kg / m 3 ; V z , mean = 4 Q / ( π D i n 2 ) = 0.0216 m ; μ a i r = 26.0 × 10 − 6 kg / m ⋅ s ; Z 1 , max = H max = 0.39 m ; Z 2 , max = H max = 0.23 m ; Z 3 , max = H max = 0.20 m ; P 1 = 1.0 atm = 101325 N / m 2 = 101325 kg / m ⋅ s 2 .

Using Equation (72), one can obtain P at different Z positions.

Case 1: When Z = Z 1 , max = 0.39 m and r o = r o 1 = 10 μm , one obtains

P min = ( 101325 − 17521.92 − 2.66 ) kg / m ⋅ s 2 = 83800.4 kg / m ⋅ s 2 = 0.827 atm (73)

Δ P = 0.173 atm

The average pressure of air in the FBR bed is 0.9135 atm.

In this case, the concentration of oxygen decreases greatly in the FBR.

Case 2: When Z = Z 2 , max = 0.23 m and r o = r o 2 = 60 μm , one obtains

P min = ( 101325 − 287.04 − 1.57 ) kg / m ⋅ s 2 = 101036.4 kg / m ⋅ s 2 = 0.997 atm (74)

Δ P = 0.003 atm

The average pressure of air in the FBR bed is 0.9985 atm.

In this case, the concentration of oxygen is uniformly distributed in the FBR.

Case 3: When Z = Z 3 , max = 0.20 m and r o = r o 3 = 100 μm , one obtains

P min = ( 101325 − 89.86 − 1.37 ) kg / m ⋅ s 2 = 101233.8 kg / m ⋅ s 2 = 0.999 atm (75)

Δ P = 0.001 atm

The average pressure of air in the FBR bed is 0.9995 atm.

In this case, the concentration of oxygen is uniformly distributed in the FBR.

From the above calculation, one can get the following conclusions:

- Oxygen is uniformly distributed in a FBR with large size of particles.

- Oxygen is not uniformly distributed in a FBR with small size of particles.

- The effect of gravity to the pressure P is negligible.

2) Estimation of Diffusion Coefficient in FBR

From the theory of random walk analysis with a vacancy mechanism, one has the following relation:

D o = ( 1 / 6 ) n ω α 2 (76)

where D_{o} is specific diffusion coefficient in cm^{2}/s; n is the nearest position of jumping; ω is the frequency to jump from a position to a specific nearest neighbor and is about 10^{12} to 10^{13} s^{−1}; α is the jumping distance of lattice sites.

For B.C.C., α = ( 3 / 2 ) a o and n = 8.

For F.C.C., α = ( 2 / 2 ) a o and n = 4.

For interstitial problem, α = ( 1 / 2 ) a o and n = 4.

Therefore, one has

D o , b c c = ω a o 2 (77)

D o , f c c = ( 1 / 3 ) ω a o 2 (78)

D o , i n t = ( 1 / 6 ) ω a o 2 (79)

Considering the correction factor f_{cf}, Equations (77)-(79) can be written as follows:

D o , b c c = 0.72149 ω a o 2 (77F)

D o , f c c = 0.78145 × ( 1 / 3 ) ω a o 2 (78F)

D o , i n t = ( 1 / 6 ) ω a o 2 (79F)

Since aluminum has a f.c.c. crystal structure and a_{o} equals 2.86 Å, Equations (78F) holds and D o , f c c . A l equals

D o , f c c . A l = 0.78145 × ( 1 / 3 ) × ( 2.86 ) 2 × 10 − 16 × 10 12 = 2.13 × 10 − 4 cm 2 / s (80)

Considering the anion vacancy mechanism, the diffusion coefficient is shown as Equation (81).

D v , a n i o n = D o , f c c . A l exp ( − Δ E t o t / R T ) (81)

where D_{v,anion} is the volume diffusion coefficient of anion vacancies; ΔE_{tot} is the total activation energy including motion and vacancies in kcal/mol; R is gas constant in 1.9872 cal/K・mol; T is the temperature in K. D_{o,fcc.Al} is specific diffusion coefficient of f.c.c. aluminum

Generally, the activation energy is between 20 kcal/mol and 40 kcal/mol. By taking the average value of 30 kcal/mol, one obtains

D v , a n i o n ( 100 ˚ C ) = 5.636 × 10 − 22 cm 2 / s (82)

D v , a n i o n ( 200 ˚ C ) = 2.931 × 10 − 18 cm 2 / s (83)

D v , a n i o n ( 300 ˚ C ) = 7.965 × 10 − 16 cm 2 / s (84)

D v , a n i o n ( 400 ˚ C ) = 3.858 × 10 − 14 cm 2 / s (85)

Since at low temperature oxidation of aluminum powders, diffusion may also be affected by grain boundaries, the actual diffusion coefficient may be less than the values calculated in Equations (82)-(85).

3) The Effect of Oxygen Pressure to Dominant Oxide Defect

As D. D. Macdonald proposed a “dry” oxidation of aluminum [

The overall reaction of the oxidation is shown as below:

Al + ( 3 / 2 ) V o 2 + + ( 3 / 4 ) O 2 + 6 e − = AlO 3 / 2 + V Al 3 − (86)

The principal crystallographic defects are

1) Vacanies: V o 2 + , V Al 3 − for AlO_{3/2}.

2) Interstitials: O i 2 − , and Al i 3+ .

The growth of the oxide occurs via simultaneous vacancy and electron transport. Either of them can be rate-determining step.

The reaction constant of Equation (86), K AlO 3 / 2 , equals

K AlO 3 / 2 = Po 2 − 3 / 4 [ e − ] − 6 [ V o 2+ ] − 3 / 2 (87)

Because stoichiometrically, 6 [ e − ] = ( 3 / 2 ) [ V o 2 + ] , Equation (87) becomes

K AlO 3 / 2 = Po 2 − 3 / 4 [ V o 2+ ] − 3 (88)

and the oxygen anion vacancy equals

[ V o 2+ ] = K ′ AlO 3 / 2 Po 2 − 1 / 4 (89)

Considering the cation vacancy V Al 3 − ,

K AlO 3 / 2 = [ V Al 3 − ] / Po 2 − 3 / 4 (90)

and

[ V Al 3 − ] = K AlO 3 / 2 Po 2 3 / 4 (91)

The effect of oxygen pressure to the dominant oxide defect can be studied by using Equations (89) and (91).

Let’s first consider the effect of oxygen pressure to oxygen anion vacancy at the worst case when aluminum powders have an average diameter as of 10 μm (see Equation (73)).

[ V o 2 + ] ∝ Po 2 − 1 / 4 (92)

P air,max = 1 atm (93)

P air,min = 0.827 atm (94)

P air,max − 1 / 4 = 1 atm (95)

P air,min − 1 / 4 = 1.049 atm (96)

and

P oxygeninair,max − 1 / 4 = 0 .21 atm (97)

P oxygeninair,min − 1 / 4 = 0 .22 atm (98)

From the above calculation, we can conclude that oxygen pressure in FBR at current working condition has no effect to the concentration of oxygen anion vacancies. The maximum variation of oxygen anion vacancies is below 5% for aluminum spherical powders with at the 10 μm diameter as the worst case 1. The vacancy concentration increases about 5.0% at Z_{max} (H_{max}) position (0.39 m).

For aluminum powder with 50 μm diameter or larger, oxygen pressure has no effect to both oxygen anion vacancies and cation vacancies. The maximum variation on defect concentration of 50 μm diameter aluminum powder, for example at 200˚C for both cases, is less than 0.4% (see Case 2 and Case 3 for Equations (74) and (75)).

4) Oxidation on A Planar Aluminum Surface

From the diffusion theory on a planar surface, the dependence of rate of oxidation and the thickness of the oxide film formed through diffusion can be expressed as shown in Equations (99A) and (99B).

K = 1.0 / ( a A l ∗ P 02 3 / 2 ) (99A)

and

ξ = 2 r ( D o t ) 1 / 2 (99B)

where ξ is the thickness of the oxide, D_{o} is the diffusion coefficient of oxygen, and r is a constant.

From Fick’s 2^{nd} law, the diffusion of oxygen and metal is written as

( ∂ N o / ∂ t ) = − D o ( ∂ 2 N o / ∂ x 2 ) (100)

( ∂ N A l / ∂ t ) = D A l ( ∂ 2 N A l / ∂ x 2 ) (101)

where D_{Al} is the diffusion coefficient of aluminum species, N_{o} and N_{Al} are atomic concentrations of oxygen and aluminum species, respectively.

In order to solve the differential equations, the following boundary conditions are assumed:

N o = N o s for x = 0 at t ≥ 0 (102)

N o = 0 for x ≥ ξ at t > 0 (103)

N A l = 0 for x ≤ ξ at t > 0 (104)

N A l = N A l b for x ≥ 0 at t = 0 (105)

where N o s is the concentration of oxygen species at the surface and N A l b is the concentration of bulk aluminum.

By combining Equations (100)-(105), the concentrations of N_{o} and N_{Al} are given

N o = N o s { 1 − erf [ x / 2 ( D o t ) 1 / 2 ] / erf ( r ) } (106)

N A l = N A l b { 1 − erf [ x / 2 ( D A l t ) 1 / 2 ] / erf ( r ) } (107)

Let

λ = D o / D A l (108)

v = N o / N A l (109)

J o = v ∗ J A l (110)

− D o ( ∂ N o / ∂ x ) = v D A l ( ∂ N A l / ∂ x ) (111)

where J_{o} and J_{Al} are the flux of oxygen and aluminum species, respectively.

By combining Equations (106), (107), and (111), one obtains

N o s / ( v N A l b ) = { exp ( r 2 ) erf ( r ) } / { λ 2 exp ( r 2 λ ) erfc ( r λ 1 / 2 ) } (112)

From Equations (112), r may be evaluated and ξ may be calculated or compared with experimental results. For internal oxidation, the following relations may have

1 2 O 2 = O (113)

kP O 2 = [ O ] (114)

When r ≪ 1 or r λ 1 / 2 ≫ 1 , the error and exponential functions can be simplified as

erf ( r ) ≈ ( 2 / π 1 / 2 ) r (115)

exp ( r 2 ) ≈ 1 (116)

erfc ( r λ 1 / 2 ) ≈ { exp ( − r 2 v ) } / π 1 / 2 r λ 1 / 2 (117)

Introduction Equations (115), (116), and (117) to (112), one has

r ≈ { N o s / ( 2 v N A l b ) } 1 / 2 (118)

and

ξ = { 2 N o s D o t / ( v N A l b ) } 1 / 2 (119)

If { 2 N o s D o / ( v N A l b ) } is a constant and N o s is a function of position Z, one obtains

ξ 2 = f ( z , T , t ) (120)

at a fixed temperature,

ξ 2 = K c ( z ) t (121)

where K c = 2 N o s D o / ( v N A l b )

A parabolic rate law is applied. It is associated with Wagner’s mechanism [

Assuming

D o ≈ 10 − 8 cm 2 / s (122)

N o s / N A l b ≤ 10 − 7 (123)

υ = 3 / 2 (124)

The thickness of the oxide film after 60 minutes oxidation from Equation (121) is shown in

N o s / N A l b | ξ in Å |
---|---|

1 × 10^{−7} | 219 |

5 × 10^{−8} | 155 |

3 × 10^{−8} | 120 |

2 × 10^{−8} | 98 |

1 × 10^{−8} | 69 |

8 × 10^{−8} | 62 |

6 × 10^{−8} | 54 |

4 × 10^{−8} | 44 |

2 × 10^{−8} | 31 |

1 × 10^{−9} | 22 |

5) Oxidation on Spherical Aluminum Powders

General diffusion in sphere has been studied by Crank [

The growth rate of oxide layer can be estimated from the total amount of diffusing substance entering (oxygen lattice anions) or leaving (aluminum ions) the sphere in a gas fluidization bed considering the initial oxide thickness at time zero.

M t / ( M ∞ − M i n i t ) = 1 − ( 6 / π 2 ) f ( r o ) ∑ n = 1 ∞ [ ( 1 / n 2 ) exp ( − V t ) ] (125)

where M t is the amount of substance entering or leaving the sphere at time t; M ∞ is the amount of substance entering or leaving the sphere at time infinite and M ∞ = 1.0 ; M i n i t is the initial amount of substance entering or leaving the sphere at time t = 0 ; f ( r o ) is a correction function due to effect of oxygen pressure to the dominant oxide defect and equals to 1.0 for spherical particle with 50 μm diameter or larger; V = D v n 2 π 2 / r o 2 ; r o is radius of aluminum powder.

f ( r o ) is a correction function of oxygen pressure at a position Z in a gas fluidization bed and the estimation of f ( r o ) to the dominant oxide defect can be described in two cases. Laplace transforms is used to describe the calculation of growth of oxide film and Laplace Transforms and Laplace Transforms is shown in the standard tables [

Case 1: Oxygen anion vacancy is the dominant oxide defect.

The correction factor considering the oxygen vacancy concentration with oxygen pressure is shown as below:

f ( r o ) o x y = { P air,min − 1 / 4 + P air,max − 1 / 4 } / 2 (126)

Equation (126) is based on the linear increase of P^{-1/4} with Z increasing in the FBR.

For aluminum spherical powders with diameters between 10 μm and 100 μm, f ( r o ) o x y _{ }is very close to 1.0 and the average maximum variation of oxygen anion vacancy is +2.5%.

Case 2: Aluminum cation vacancy is the dominant oxide defect.

The correction factor considering the oxygen vacancy concentration with oxygen pressure is shown as below:

f ( r o ) c a t = { P air,min 3 / 4 + P air,max 3 / 4 } / 2 (127)

Equation (127) is based on the linear decrease of P^{3/4} with Z increasing in the FBR.

For aluminum spherical powders with diameters between 10 μm and 100 μm, f ( r o ) c a t is also very close to 1.0 and the average maximum variation of aluminum cation vacancy is −6.0%.

The typical values of f ( r o ) o x y and f ( r o ) c a t are 1.025 and 0.994 for 10 μm diamter spherical aluminum powders, respectively. Detailed calculation were generated through a computer program called “APRIL.EXE” written by Shih [_{init}, M_{t.cal} and D_{v}, respectively. Equations (70) and (125) are the base of the computer software for the calculation.

For spherical aluminum powders with diameter at 50 μm or larger, both f ( r o ) o x y and equal 1.0. It means that the oxygen pressure is so uniform in the FBR that oxygen pressure has no effect to the concentration of the dominant oxide defects at different Z positions.

Typical calculation of oxide thickness using software “APRIL.EXE” [

In _{t,exp} and M_{t,cal} for Oxygen Anion Vacancy Mechanism is calculated using computer software “APRIL.EXE” at various temperatures when r_{o} = 10 μm, f(r_{o})_{oxy} = 1.025, and M_{init} = 0.00075 (25 Å) are considered.

M_{∞} | M_{t}_{,exp}* | M_{t,cal} |
---|---|---|

1.0 | 0.0030* | 0.000865 (373^{ }K) |

0.001113 (473^{ }K) | ||

0.006475 (573^{ }K) | ||

0.040230 (673^{ }K) |

*: M_{t,exp} is the experimental thickness in a FBR obtained at the optimized FBR temperature based on the theoretical calculation between 473 K and 573 K.

From

In _{t,exp} and M_{t,cal} for Oxygen Anion Vacancy Mechanism is calculated using computer software “APRIL.EXE” at various temperatures when r_{o} = 60 μm, f(r_{o})_{oxy} = 1.0, and M_{init} = 0.000125 (25 Å) are considered.

From

In _{t,exp} and M_{t,cal} for Oxygen Anion Vacancy Mechanism is calculated using computer software “APRIL.EXE” at 373˚ K when r_{o} = 10 μm, f(r_{o})_{oxy} = 1.025, and M_{init} = 0.00075 (25 Å) are considered.

From

In _{t,exp} and M_{t,cal} for Oxygen Anion Vacancy

M_{∞} | M_{t}_{,exp}* | M_{t,cal} |
---|---|---|

1.0 | 0.0005 | 0.000240 (373 K) |

0.000253 (473 K) | ||

0.00107 (573 K) | ||

0.00676 (673 K) |

*: M_{t,exp} is the experimental thickness in a FBR obtained at the optimized FBR temperature based on the theoretical calculation between 473 K and 573 K.

Time in seconds | M_{t}_{,exp}* | M_{t,cal} |
---|---|---|

1000 | 0.0030 | 0.0008642 |

2000 | 0.0008644 | |

3600 | 0.0008648 | |

7200 | 0.0008650 | |

14400 | 0.0008652 |

*: M_{t,exp} is the experimental thickness in a FBR obtained at the optimized FBR temperature based on the theoretical calculation between 473 K and 573 K. But the target oxide layer thickness cannot be achieved in the theoretical calculation when FBR sets at 373^{ }K.

M_{t}_{,exp}* | M_{t,cal} | D_{v} in cm^{2}/sec | T K |
---|---|---|---|

0.00050 | 0.000240 | 5.636 × 10^{−22} | 373 |

0.000265 | 2.931 × 10^{−18} | 473 | |

0.000348 | 1.643 × 10^{−17} | 500 | |

0.000407 | 2.970 × 10^{−17} | 510 | |

0.000443 | 3.959 × 10^{−17} | 515 | |

0.000487 | 5.248 × 10^{−17} | 520 | |

0.000537 | 6.919 × 10^{−17} | 525 |

*: M_{t,exp} is the experimental thickness in a FBR obtained at the optimized FBR temperature based on the theoretical calculation between 473 K and 573 K. M_{t,exp} can be achieved between 520 K and 525 K.

Mechanism is calculated using computer software “APRIL.EXE” at various temperatures when r_{o} = 60 μm, f(r_{o})_{oxy} = 1.0, and M_{init} = 0.000125 (25 Å) are considered. The oxidizing time in the gas fluidization bed reactor is fixed as of 7200 seconds.

From

1) The mathematical modeling considering mass, momentum, energy conservation, effects of powder size, temperature, gas velocity, and other related factors has been developed for aluminum powder oxidation in a gas fluidized bed. This is the 1^{st} modeling for a quantitative calculation of thickness of aluminum oxide in a gas fluidization bed.

- At cross section of the pipe, the distribution of V_{z} is parabolic and V_{z} does not change along Z direction.

- Pressure is linearly decreasing along Z direction and it is proportional to the pressure of air.

2) In the mathematical model, Equations (65) and (70) as a solution can be used in the calculation of the growth rate of aluminum oxide quantitatively. V_{z}_{,mean} is the mean velocity of air. The average oxidation rate of aluminum powders can be estimated from Equations (65), (70), and (125) for spherical aluminum powders and Equations (92) and (93) for planar aluminum surface.

3) The oxidation kinetics of planar and spherical aluminum powders has been discussed in great details. A uniform temperature field and a nice distribution of gas phase provide the uniform oxide growth rate of aluminum powders in the gas fluidization bed.

4) The basic theories of oxidation and growth rates have been summarized. The most possible growth rate of aluminum powders at low temperature initially is linear due to the surface boundary process or reaction is rate-determining. After initial oxidation the growth rate will follow parabolic law due to the diffusion processes in the growth of aluminum oxide. In this case, Wagner’s mechanism holds.

5) The oxidation mechanisms of aluminum powders can be determined by interstitial cations diffusion process or anion vacancies diffusion process. As an n-type semiconductor, the rate-determining process of aluminum oxide is diffusion.

6) If the oxidation is controlled by a diffusion process, the growth rate of the oxide layer can be calculated by solving partial differential equations. Mathematical models have been developed and presented in the article. The modeling of internal oxidation process has also been discussed.

7) The growth rate from the experimental results is compared with the estimation of the models developed in this article. Very good correlation between theoretical calculation and experimental results are observed.

8) The optimization of gas fluidization of the oxide of aluminum powders involve many factors which have to be considered in the theoretical analysis of mathematical modeling. The detail information can be obtained through theoretical calculation and its comparison with experimental results. A computer program written by the author has been used in the calculation. An example calculation has been shown in the technical paper at a fixed temperature, particle size, gas velocity, bulk concentration of oxygen and aluminum, Archimedes number, minimum fluidization velocity, bed depth, and effective surface area.

9) The oxidation at high temperature generates stress in the outer oxide films on aluminum particles and leads to the formation of cracks. With the increase in the intensity of crack formation, the oxidation rate increases markedly, particularly for the γ-α-Al_{2}O_{3} transition. The oxidation film is not protective at high temperature oxidation. By controlling the oxidation at low temperature with the variations of oxygen concentration and relative humidity, the oxide layer is more protective and uniform.

The author would like to express thanks to Leviton Manufacturing Company for sponsoring the project through Henry Martínez and Steven Chen. The author also would like to express great thanks to Stanford Research Institute for conducting the gas fluidization bed tests to confirm the theoretical calculation and experimental results of aluminum spherical powder oxidation through Dr. Kai-Hung Lau and A. Sanjurjo. Helpful discussion with Dr. D. B. Macdonald, Dr. Richard Gottscho, and Dr. John Daugherty are greatly appreciated.

The author declares no conflicts of interest regarding the publication of this paper.

Shih, H. (2019) Oxidation Kinetics of Aluminum Powders in a Gas Fluidized Bed Reactor in the Potential Application of Surge Arresting Materials. Materials Sciences and Applications, 10, 253-292. https://doi.org/10.4236/msa.2019.103021

Diffusion in sphere was studied by Crank and others. The partial diffusion equation has a standard form, but the solution depends on the initial and boundary conditions for each special case. For aluminum spherical powders in FBR, Shih provided the solution through Laplace Transformation as shown below [

Since the formation of aluminum oxide involves diffusion of metal interstitial cations or anion vacancies at different ranges of temperature. The differential equations can be considered in the following ways.

∂ N o / ∂ t = D o { ∂ 2 N o / ∂ r 2 + ( 2 / r ) ( ∂ N o / ∂ r ) } (A1)

∂ N A l / ∂ t = D A l { ∂ 2 N A l / ∂ r 2 + ( 2 / r ) ( ∂ N A l / ∂ r ) } (A2)

Equations (A1) and (A2) are the diffusion equations for anion vacancies and metal interstitial cations, respectively.

Considering only oxygen anion vacancies diffusion process as rate-determining step and assuming that the concentration of N A l b at metal/oxide interface keeps constant, then at steady-state which indicates

∂ N o / ∂ t = D o { ∂ 2 N o / ∂ r 2 + ( 2 / r ) ( ∂ N o / ∂ r ) } = 0 (A3)

of which the general solution is

N o = B + A / r (A4)

where A and B are constants to be determined from the boundary conditions.

If the thickness of oxide is ξ and the radius of the sphere is r_{o}, and

( r o − ξ ) ≤ r ≤ r o (A5)

N o ( ∞ , t ) = 0 (A6)

N o ( r , 0 ) = 0 (A7)

N o ( r o , t ) = C s (A8)

N o ( r − ξ , t ) = C o (A9)

Using Laplace Transform, one has

N ˜ o ( r , p ) = ∫ o ∞ exp ( − p t ) N o ( r , t ) d t (A10)

∫ o ∞ ( ∂ N o / ∂ r ) exp ( − p t ) d t = N o exp ( − p t ) | o ∞ + p ∫ o ∞ N o exp ( − p t ) d t = − N o ( r , 0 ) + p N ˜ o ( r , p ) = p N ˜ o ( r , p ) (A11)

where N o ( r , 0 ) = 0

∫ o ∞ D o ( ∂ 2 N o / r 2 ) exp ( − p t ) d t = D o ( ∂ 2 / ∂ r 2 ) ∫ o ∞ N o exp ( − p t ) d t = D o ( ∂ 2 N ˜ o ( r , p ) / ∂ r 2 ) (A12)

∫ o ∞ ( 2 / r ) ( ∂ N o / ∂ r ) exp ( − p t ) d t = ( 2 ∂ / r ∂ r ) ∫ o ∞ N o exp ( − p t ) d t = 2 ∂ N ˜ o ( r , p ) / r ∂ r (A13)

Equation (A1) becomes

D o ( d 2 N ˜ o ( r , p ) / d r 2 ) + ( 2 / r ) d N ˜ o ( r , p ) / d r = p N ˜ o ( r , p ) (A14)

or

D o ( d 2 N ˜ o ( r , p ) / d r 2 ) + ( 2 / r ) d N ˜ o ( r , p ) / d r − p N ˜ o ( r , p ) = 0 (A15)

Let N ˜ o = x / r

d N ˜ o / d r = ( d x / d r ) ( 1 / r ) − x / r 2 (A16)

d 2 N ˜ o / d r 2 = ( 1 / r ) ( d 2 x / d r 2 ) − ( 2 / r 2 ) ( d x / d r ) + 2 x / r 3 (A17)

and

d 2 x / d r 2 − ( p / D o ) x = 0 (A18)

Therefore

x = k 1 exp ( − p / D o ) r + k 2 exp ( p / D o ) r (A19)

and

N ˜ o = ( k 1 / r ) exp ( − p / D o ) r + ( k 2 / r ) exp ( p / D o ) r (A20)

at r = r o

N o = C s and N ˜ o = C s / p (A21)

at

r = r o − ξ

N o = C o and N ˜ o = C o / p (A22)

We have

N ˜ o = C s / p = ( k 1 / r o ) exp ( − p / D o ) r o + ( k 2 / r o ) exp ( p / D o ) r o (A23)

and

N ˜ o = C o / p = ( k 1 / ( r o − ξ ) ) exp ( − p / D o ) ( r o − ξ ) + ( k 2 / ( r o − x ) ) exp ( p / D o ) ( r o − ξ ) (A24)

at r = ∞ , N o = 0 , we have

k 2 = 0 (A25)

N ˜ o = C s / p = ( k 1 / r o ) exp ( − p / D o ) r o (A26)

N ˜ o = C o / p = ( k 1 / ( r o − ξ ) ) exp ( − p / D o ) ( r o − ξ ) (A27)

with

k 1 = exp ( p / D o ) r o { A / 2 p + [ ( B − E ) / 2 p ] exp ( − p / D o ) ξ } (A28)

where

A = C s r o (A29)

B = C o r o (A30)

E = C o ξ (A31)

Using reverse Laplace Transformation, we have

N o ( r , t ) = [ a / ( 2 ( π t 3 ) 1 / 2 ) ] exp ( − a 2 / 4 t ) { [ A / 2 r + ( B − E ) / 2 r ] erfc ( − ξ / ( 2 D o t ) 1 / 2 ) } (A32)

where

a = ( r o − r ) / D o 1 / 2 (A33)

Equation (A32) can be used to estimate N_{o}(r, t) at a fixed time when ξ and r_{o} are selected.

Symbol Explanation Unit

A_{air} area needed for heating air m^{2}

A_{m} surface area of the wall in contact with fluidized bed m^{2}

A_{r} Archimedes number -

C_{air} specific heat of air kJ/kgK

C_{Al} specific heat of aluminum power kJ/kgK

C_{bed} circumference of the reactor m

D_{in} inner diameter of bed m

D_{out} outer diameter of bed m

d_{bd} diameter of sphere having same volume as a bubble m

d_{i} arithmetic mean of adjacent sieve apertures m

d_{p} mean sieve size of Al powder m

d_{v} diameter of sphere particles m

f_{Al} feed rate of aluminum powder kg/s

Q total heat rate W

Q_{air} total heat to heat air W

Q_{Al} total heat rate to heat powders W

Q_{b} total heat rate at fixed temperature W

J_{air} volume flow rate of air m^{3}/s_{ }

H_{air}_{ }heat length of air m

H_{max}_{ } maximum bed height m

H_{set} height of gently settled bed m

M_{b}_{ }mass of powder in bed kg

N_{p} total particle number of powder -

S_{p} total surface area of powder m^{2}

K_{qz}_{ }thermal conductivity of quartz W/mK

V_{ff}_{ }free fall velocity of Al powder m/s

T_{b}_{ }bed temperature K

T_{f}_{ }fluid temperature K

α_{max}_{ } maximum heat wall-to-bed transfer coefficient W/m^{2}K

α_{o}_{ }overall heat transfer coefficient W/m^{2}K

α_{f}_{ }heat transfer coefficient between tube and fluid W/m^{2}K

ρ_{p} Al particle density kg/m^{3 }

ρ_{f} density of fluid kg/m^{3}

ρ_{air}_{ } air density kg/m^{3 }

V total superficial gas velocity m/s

V_{air}_{ }air flow rate kg/s

U_{mf} minimum fluidization velocity m/s

x thickness of reactor wall (quartz) m

x_{i} mass fraction of particles of size d_{i}_{ } -

μ viscosity kg/ms

ε_{void} bed voidage -

ΔP pressure drop N/m^{2}

Symbol Explanation Unit

r, θ, z components of cylindrical and spherical coordinates -

V_{r}, V_{θ}, V_{z} velocity components corresponding to r, θ, and z -

ρ gas density kg/m^{3}

μ gas viscosity kg/ms

ν kinematic viscosity of gas, ν = μ/ρ m^{2}/s

P pressure N/m^{2 }

g acceleration due to gravity 9.81 m/s^{2}

Q total flow rate entering the bed m^{3}/s

q flow rate of the small flow pipes and and q = Q/n m^{3}/s

R radius of flow pipe m

R_{e} Reynolds number -

n Number of holes of gas distributors -

M_{t} the diffusion substance at time t kg

M_{¥} the diffusion substance kg

M_{init} the diffusion substance at time zero kg

M_{t,exp}_{ }the diffusion substance in fluidized bed kg

M_{t,cal} the diffusion substance at time t calculated kg

r_{o} radium of aluminum spherical powders μm

f(r_{o}) the correction function of oxygen pressure -

f(r_{o})_{oxy}_{ }correction factor considering the oxygen vacancy concentration -

f(r_{o})_{cat}_{ }correction factor considering aluminum cation vacancy concentration -

D_{ν} specific volume diffusion coefficient m^{2}/s

T time in second s

C_{1} integration constant -

C_{2} integration constant -

ξ the thickness of oxide layer m

ν ratio of oxygen atom and aluminum atom -

V o 2 + oxygen anion vacancy -

V Al 3 − ^{ }aluminum cation vacancy -

O i 2 − oxygen interstitial -

Al i 3 + aluminum interstitial -

N o s atomic concentration of oxygen species atoms/m^{3}

N A l b atomic concentration of bulk aluminum atoms/m^{3}

K reaction rate of oxide growth on planar Al -

K_{c} reaction constant at a fixed z position -

D_{o} specific diffusion coefficient cm^{2}/s

n the nearest position of jumping in random walk -

ω the frequency to jump from a position a specific nearest neighbor. 10^{12} to 10^{13} s^{−1}

α the jumping distance of lattice sites cm

B.C.C. Body-centered cubic structure -

F.C.C. Face-centered cubic structure -

D_{o,bcc} specific diffusion coefficient of B.C.C. cm^{2}/s

D_{o,fcc} specific diffusion coefficient of F.C.C. cm^{2}/s

D_{o,int} specific diffusion coefficient of interstitial cm^{2}/s

D_{o,fcc.Al}_{ }_{ }specific diffusion coefficient of aluminum of anion vacancy mechanism cm^{2}/s

D_{v,anion} volume diffusion coefficient of aluminum of anion vacancy mechanism cm^{2}/s

DE_{tot} the total activation energy including motion and vacancies kcal/mol

R gas constant 1.9872 cal/K・mol

erf(x) error function erf ( x ) = 2 π ∫ 0 x e − t 2 d t

erfc(x) the complementary error function erfc ( x ) = 1 − erf (x)

erfc ( x ) = 1 π ∫ x ∞ e − t 2 d t