WO2011097961A1

WO2011097961A1 - Method for production of ultrapure phosphorus by zone melting and apparatus thereof

Info

Publication number: WO2011097961A1
Application number: PCT/CN2011/000218
Authority: WO
Inventors: Yongli Mi; Xiaoping Wang; Zhen Zhang; Xiaojin He
Original assignee: Dongying Jieda Chemical Technology Co., Ltd.
Priority date: 2010-02-12
Filing date: 2011-02-12
Publication date: 2011-08-18
Also published as: US20110212405A1

Abstract

Disclosed are a method for the production of ultrapure phosphorus by zone melting in a non-flammable environment and an apparatus thereof. The method comprises the steps of: a. fixing a heater in a non-flammable environment, b. fixing a phosphorus bar through the heater, c. setting up a motor to control the heater moving along the phosphorus bar, d. heating a part of the bar to form a heat zone, e. moving the heater along the bar to have the heat zone moving with the heater, and f. repeating step e if necessary. The apparatus comprises a heater, a phosphorus bar fixed through the heater and between fixing units, and a motor set up to control the heater moving along the phosphorus. The method is clean, chemical free, fast and energy efficient.

Description

Method for the Production of Ultrapure Phosphorus b Zone Melting in a Non-Flammable Environment, and the apparatus Used in Such Method

Technical field

Generally, the present invention relates to the production of ultrapure phosphorus. In particular, the present invention relates to the method for the production of ultrapure phosphorus by zone melting in a Non-flammable environment and the apparatus used in such method. Background art

Phosphorus acid is commonly used as food additive. It is commonly added in almost all kinds of soft drinks, juice, wine etc. As an important element for life forms, the phosphorus is mainly supplied for human being and animal in the form of the food intake.

For food industries, industrial phosphorus must be purified to lower the impurities concentration, especially the arsenic concentration. The high arsenic concentration can increase the risk of several types of cancer, such as skin and lung cancer. Industrial grade phosphorus contains arsenic of the concentration range from 50 to 140 ppm (some to even 300 ppm). Phosphorus and arsenic co-exist because they are both group V elements and share similar chemical properties, so it is difficult to separate them by simple chemical method. The arsenic concentration for the food grade phosphorus should be below 10 ppm, which can produce 85% phosphorus acid with the arsenic concentration of 2.8 ppm. According to both the Food and Drug Administration standard (FDA)^la and the Food Chemicals Codex standard (FCC)^lb, the arsenic concentration in food grade phosphoric acid (85%) should be below 3 ppm.

Therefore, there is a continued need in the art to provide a process of phosphorus purification which is not energy costly and as clean as possible.

Contents of the invention:

In one aspect of the present invention, there is provided a method of production of ultrapure phosphorus by zone melting in a non-flammable environment, comprising the steps of:

a. Fixing a heater in a non-flammable environment;

b. Fixing a phosphorus bar through the heater;

c. Setting up a motor to control the heater moving along the phosphorus bar;

d. Heating a part of the said bar to become a heat zone;

e. Moving the heater along the said bar in order to have the heat zone moving with the heater;

f. Repeat step e if necessary.

In a preferred aspect of the present invention, the length of the heater D, the radius of the phosphorus bar R, the temperature of the heater Th, and the moving velocity of the heater v should control the solutions of the following governing equations:

In r<R and 0<z<L,

c, dt ₊ W¾ dt 2¾w dt = v^.[*_Wv_n+^ dt v

The boundary conditions are:

kW)? = h(T_h - T) , for v - t≤z≤v - t + D and r=R, dn

k( )— = h(T_Q— T) , for z < v - t , z > v - t + D and r=R, z=0, z=L,

dn

^ = 0 , for r=R, z=0, z=L.

dn

In these equations, φ(χ, ί) is an continuous function which labels the two phases at the position x and the time /, with φ = 0 in the solid phase, φ = 1 in the liquid phase and 0 < φ < 1 in the interfacial region. M is the mobility coefficient, K_E and K_s are interfacial energy coefficients, W_E and W_s are barrier heights for double well potential function g(tf>) = φ² (\ - φ)² . ρ(φ) = φ² (3 - 2φ) . L_y is the latent heat per unit volume and T_M is the melting point. T(x, t) is the temperature distribution at the position Λ: and the time t. c_p is the heat capacity per unit volume, r(<#) = ρ(φ) . k(<fi) permits different thermal conductivities k_L = k(0) and k_s = k(l) in the bulk liquid and solid phases, respectively. These governing equations are solved on a cylindrical domain x = (z, r) e [0, L] x [0, R] . L is the length of phosphorus cylinder and R is the radius. D is the length of the heater, z is the axial direction coordinate, r is the radial direction coordinate, n is the outer normal vector of the boundary, h is the heat transfer coefficient, T₀ is the room temperature and T_h is the temperature of the heater.

In another aspect of the present invention, there is provided an apparatus for the method of the production of ultrapure phosphorus as mentioned above, which comprises in a non-flammable environment:

a) . a heater;

b) . a phosphorus bar fixed through the heater; c). a motor set up to control the heater moving along the phosphorus bar.

Description of figures

Figure 1 : A schematic drawing illustrating the apparatus for the zone-melting purification of phosphorus according to one embodiment of the present invention.

Figure 2: A schematic drawing showing the detail of the heater used in the apparatus according to one embodiment of the present invention. Figure 3 : A schematic drawing of the experimental apparatus for sample analysis used in one embodiment of the present invention.

Figure 4: A figure showing the moving characteristic of the fusion region and the heater according to one embodiment of the present invention, wherein temperature of the heater is 50°C, and the velocity of the heater is set at 0. lmm/s.

Figure 5: A figure showing the temperature over the whole phosphorus cylinder at r=0, r=R/2 and r=R while the heater moves as described in figure 4. Figure 6: A figure showing the temperature over the whole phosphorus cylinder of figure 9 in 2D wherein t =400s

Figure 7: A figure showing the moving characteristic of the fusion region and the heater according to one embodiment of the present invention, wherein temperature of the heater is 50°C, and the velocity of the heater is set at 0. lmm/s.

Figure 8: A figure showing the moving characteristic of the fusion region and the heater according to one embodiment of the present invention, wherein temperature of the heater is 50°C, and the velocity of the heater is set at 0.2mm/s. Figure 8(a) shows the left interface of the fusion region lags behind the left boundary of the heater; the temperature at different radii at t=250s are shown in (b).

Figure 9: A figure showing the moving characteristic of the fusion region and the heater according to one embodiment of the present invention, wherein temperature of the heater is 50°C, and the velocity of the heater is set at 0.05mm/s. Figure 9(a) shows the fusion region stays in the middle of the heater; the temperature at different radii at t=600s are shown in (b).

Figure 10: A figure showing the moving characteristic of the fusion region and the heater according to one embodiment of the present invention, wherein temperature of the heater is 50°C, and the velocity of the heater is set at 0.02mm/s. Figure 10(a) shows the fusion region stays in the middle of the heater; the temperature at different radii at t=2000s are shown in (b). Figure 1 1 : A figure showing the moving characteristic of the fusion region and the heater according to one embodiment of the present invention, wherein temperature of the heater is 50°C, and the velocity of the heater is set at 0.5mm/s. Figure 1 1(a) shows the left interface of the fusion region lags behind the left boundary of the heater; the temperature at different radii at t=150s are shown in (b).

Figure 12: A figure showing the moving characteristic of the fusion region and the heater at different time, wherein the temperature of the heater is 45°C, and the velocity of the heater is set at 0.05mm/s. The width of the melting region is much narrower than the length of the heater.

Figure 13 : A figure showing the temperature over the whole phosphorus cylinder at r=0, r=R/2 and r=R while the heater moves as described in figure 12.

Figure 14: A figure showing the moving characteristic of the fusion region and the heater at different time, wherein the temperature of the heater is 45°C, and the velocity of the heater is set at 0.02mm/s.

Figure 15: A figure showing the moving characteristic of the fusion region and the heater according to one embodiment of the present invention, wherein the temperature of the heater is 55°C, and the velocity of the heater is set at 0.05mm/s. The temperature at different radii at t=400s are shown in (b).

Figure 16: A figure showing the moving characteristic of the fusion region and the heater according to one embodiment of the present invention, wherein the temperature of the heater is 60°C, and the velocity of the heater is set at 0.05mm/s. The temperature at different radii at t=600s are shown in (b)

Figure 17: A figure showing the moving characteristic of the fusion region and the heater according to one embodiment of the present invention, wherein the temperature of the heater is 70°C, and the velocity of the heater is set at 0.02mm/s. The temperature at different radii at t= 1600s are shown in (b).

Figure 18: A figure showing the moving characteristic of the fusion region and the heater according to one embodiment of the present invention, wherein the temperature of the heater is 80°C, and the velocity of the heater is set at 0.02mm/s. The temperature at different radii at t=1600s are shown in (b).

Figure 19: A figure showing the moving characteristic of the fusion region and the heater according to one embodiment of the present invention, wherein the temperature of the heater is 50°C, and the velocity of the heater is set at O. lmm/s, the length d of the heater is 1.5cm, and the radius of phosphorus cylinder is 5mm. The temperature at different radii at t=500s are shown in (b).

Figure 20: A figure showing the moving characteristic of the fusion region and the heater according to one embodiment of the present invention, wherein the temperature of the heater is 50°C, and the velocity of the heater is set at O. lmm/s, the length d of the heater is 1.0cm, and the radius of phosphorus cylinder is 5mm. The temperature at different radii at t=800s are shown in (b). The core of the cylinder cannot be melted. It will keep this shape in (a) after t=800s.

Figure 21 : A figure showing the moving characteristic of the fusion region and the heater according to one embodiment of the present invention, wherein the temperature of the heater is 50°C, and the velocity of the heater is set at 0.05mm/s, the length d of the heater is 3.0cm, and the radius of phosphorus cylinder is 1cm. The temperature at different radii at t=1300s are shown in (b).

Figure 22: A figure showing the moving characteristic of the fusion region and the heater according to one embodiment of the present invention, wherein the temperature of the heater is 50°C, and the velocity of the heater is set at 0.05mm/s, the length d of the heater is 2.0cm, and the radius of phosphorus cylinder is 1cm. The temperature at different radii at t=2000s are shown in (b). The core of the cylinder cannot be melted. It will keep this shape in (a) after t=2000s.

Figure 23 : A figure showing the moving characteristic of the fusion region and the heater according to one embodiment of the present invention, wherein the temperature of the heater is 50°C, and the velocity of the heater is set at 0.05mm/s, the length d of the heater is 4.0cm, and the radius of phosphorus cylinder is 1.5cm. The temperature at different radii at t= 1900s are shown in (b). Figure 24: A figure showing the moving characteristic of the fusion region and the heater according to one embodiment of the present invention, wherein the temperature of the heater is 50°C, and the velocity of the heater is set at 0.05mm/s, the length d of the heater is 3.0cm, and the radius of phosphorus cylinder is 1.5cm. The temperature at different radii at t=2000s are shown in (b). The core of the cylinder cannot be melted. It will keep this shape in (a) after t=2000s. Figure 25: A figure showing the moving characteristic of the fusion region and the heater according to one embodiment of the present invention, wherein the temperature of the heater is 60°C, and the velocity of the heater is set at 0.05mm/s, the length d of the heater is 10.0cm, and the radius of phosphorus cylinder is 5.0cm. The fusion region has been separated to two parts so that it cannot wipe over the whole cylinder.

Figure 26: A figure showing the moving characteristic of the fusion region and the heater according to one embodiment of the present invention, wherein the temperature of the heater is 60°C, and the velocity of the heater is set at 0.02mm/s, the length d of the heater is 10.0cm, and the radius of phosphorus cylinder is 5.0cm. The fusion region can wipe over the whole cylinder.

Figure 27: A figure showing the moving characteristic of the fusion region and the heater according to one embodiment of the present invention, wherein the temperature of the heater is 60°C, and the velocity of the heater is set at O. lmm/s, the length d of the heater is 15.0cm, and the radius of phosphorus cylinder is 5.0cm. The fusion region has been separated to two parts so that it cannot wipe over the whole cylinder. Figure 28: A figure showing the moving characteristic of the fusion region and the heater according to one embodiment of the present invention, wherein the temperature of the heater is 60°C, and the velocity of the heater is set at 0.05mm/s, the length d of the heater is 15.0cm, and the radius of phosphorus cylinder is 5.0cm. The fusion region can wipe over the whole cylinder.

Figure 29: A figure showing the moving characteristic of the fusion region and the heater according to one embodiment of the present invention, wherein the temperature of the heater is 60°C, and the velocity of the heater is set at O. lmm/s, the length d of the heater is 25.0cm, and the radius of phosphorus cylinder is 5.0cm. The fusion region can wipe over the whole cylinder.

Mode of carrying out the invention

Definitions

Throughout the present description, the term "non-flammable environment" refers to any environment which could prevent the flame of phosphors. Preferably the method of the present invention is carried out in water.

The purification methods of phosphorus in the prior art Phosphorus has a melting point of 44.2°C and is extremely flammable when exposed to air. Because it does not have conventional solvent, purification of phosphorus is very difficult. In summary, there are the following purification methods in the prior art.

(1) Vacuum Distillation⁷

Excessive iodine (around 10,000 times in excess) is added to the white phosphorus by the stoichiometric amount. Iodine would then react with arsenic to form arsenic iodide. Arsenic iodide has a much higher boiling point than phosphorus and, thus, arsenic can be separated from phosphorus by this vacuum distillation method. Most of the arsenic can be separated from phosphorus through the vacuum distillation by the difference in boiling point. The disadvantage is that large amounts of iodine, which is corrosive and poisonous, would be produced by this process. This method is, therefore, very dangerous to operate.

(2) Liquid-liquid Extraction⁸

Crude yellow phosphorus is first mixed with the oxidant, nitric acid. Most of the arsenic inside the phosphorus would react to form arsenic acid but only part of the phosphorus would react to form phosphoric acid. When the mixture is heated to melting point or more phosphorus is melted, liquid phosphorus and nitric acid are separated into two layers by the difference of specific gravity. Most of the arsenic acid is dissolved into nitric acid due to different selectivity and thus separation of the arsenic can be done by extracting the liquid phosphorus. This method is not widely used because it uses nitric acid as a solvent, which is a strong acid and highly corrosive and highly oxidant.

(3) Fractional Distillation⁹

Impure phosphorus is continually passed through a heat exchanger and heated to vaporize a portion of the phosphorus at a temperature not above 200°C. The materials with large boiling point difference to vaporize phosphorus will be separated in a disengaging zone, while arsenic will be separated in a fractionating zone due to the vapour-liquid equilibrium. Distillation is capital and labour intensive and also leaves a distillation residue product that is much higher in arsenic and other heavy metals due to the concentration caused by the distillation.

(4) Mixing and Extraction¹⁰

The purification of yellow phosphorus can be done in this process by mixing liquid yellow phosphorus with active carbon and then separating the purified phosphorus from the spent active carbon containing adhering phosphorus. The spent active carbon is mixed with sufficient liquid phosphorus to form slurry which is pumped to a disposal zone. The disposal of the spent active carbon can cause several problems. For example, the active carbon spent gives rise to a useless and insignificant consumption of phosphorus.

(5) Phase Separation¹¹

Impure white phosphorus mixed with an aqueous suspension of a purifying amount of active charcoal can be used to separate phosphorus from the separation funnel due to different phases. The mixture of aqueous impurities is first separated into a settled phase consisting essentially of a spent active charcoal and phosphorus-containing phase and an aqueous phase. Then the phosphorus-containing phase and thus, the purified liquid white phosphorus can be separated from another separation. However, the production of the active charcoal causes much pollution as the acid is involved and heavy metal is produced.

(6) Oxidation¹²

The purification of crude yellow phosphorus can be treated by oxidizing those with mixed acid of sulphuric acid and nitric acid. The crude yellow phosphorus is first pretreated by alkali leaching with ammonium carbonate. Then, the reactant is oxidized with mixed acid of sulphuric acid and nitric acid. The purified phosphorus can be found after washing with hot water. The strong acid, sulphuric acid and nitric acid, which is highly corrosive, and strongly oxidizing agent is used as the oxidant. This can result in water pollution when disposing them into the sea. The loss of phosphorus in reaction forming phosphorus acid would also occur.

It can be seen that the methods summarized above all have significant disadvantages such as high energy consumption with hazardous operation-conditions (toxic solvent and vapor) and low yield. Therefore, in one aspect of the present invention, the objects thereof are to find an improved process for the production of ultrapure phosphorus which is clean, chemical free, fast, and energy efficient. It is absolutely a zero pollution process¹. In particular, the present invention provides a method for the production of ultrapure phosphorus by zone melting in non-flammable environment.

The apparatus and methods used in the present invention

According to one embodiment of the present invention, the apparatus used for zone melting includes a row of ring-shaped heater that moves slowly along a tube containing the solid to be purified. The heater heats a narrow band of the material, above the melting point, form the heat zone. Then the heat zone moves along with the heater. After heater passes, the heat zone cools and freezes. The impurities tend to stay in the heat zone and are carried to one end of the tube.

The impurities which could be purified by the zone melting purification should satisfy two important facts, which are (1) the diffusion coefficients of the impurities in the liquid zone are several orders higher in magnitude than that in the solid region, and (2) the solubility of the impurities in the amorphous liquid is significantly higher than that in the crystal solid. Therefore, the term "impurities" used in the present disclosure refers to any materials satisfying the above two requirements, preferably Arsenic.

The apparatus for the zone-melting purification of phosphorus according to one embodiment of the present invention is illustrated in Figure 1

Because phosphorus is very flammable in the air, all the operation must be in non-flammable environment, preferably in water. The apparatus used in the present invention could be put in a stainless steel tank which is a water container with certain volume, for example 600 mm length, 420 mm width, and 300 mm height.

Figure 1 shows the apparatus for zone-melting purification that includes three fixed pulleys, a step motor, and three slide rails. There is a plate which could be plastic plate or iron plate, with one side connected with the rails and its opposite side connected with the heater by the junction, as well as the top side fixed on the string of the pulleys by a junction. When the motor starts to work, the string draws the plate to move along rails, and then the heater will move along the phosphorus bar at the same time. There are two plastic end blocks for fixing the phosphorus bar by two screws on a rail. The phosphorus bar should be fixed through the heater and between the two plastic end blocks. The tank should be filled with material which provides non-flammable environment. Preferably the tank is filled with water. The water level must be kept above the holder of the heater in order to make sure all the conduction can be done under the water safely.

Although the apparatus illustrated in figure 1 is assembled in horizontal direction, the persons skilled in the art might understand that they could also be in vertical direction.

Figure 2 shows the design of the heater. According to one embodiment, the heater is made by an Aluminum cylinder wrapped by the copper wire which is embedded in a plastic block. The length of the heater should be longer than the diameter of the phosphorus bar. According to one embodiment of the present invention, the length of the heater is at least 2 times, preferably at least 3 times, more preferably at least 4 times, longer than the diameter of the phosphorus bar. In the middle of the Aluminum cylinder, there is a thermal couple which can control heater's temperature. The thermal couple could be any suitable type, for example T type, and the sensitivity of the temperature can be set till 0.1 °C.

The phosphorus bar can be molded by a plastic mold. The plastic mold was made by high density polyethylene or Teflon tube with a round hole in the middle. The diameter of the hole is 10 mm.

In this invention, the phosphorus bar is purified as it was molded in a cylindrical bar. But, for the convenience of future production, the phosphorus can also be confined in a case which can be plastic (polyethylene or Teflon, or glass). The detailed information of all the chemicals and machines which were used in the present disclosure are listed below in Table 1 and Table 2:

Table 1 : Chemicals

Examples

Determine of the As concentration

The As concentration in the solid phosphorus were determined before and after the purification. Any test methods could be used in the present invention for the determine of As concentration. Preferably ICP-OES (Inductively Coupled Plasma Optical Emission Spectrometry) was used for analyzing the trace concentration of Arsenic.

Generally, the procedure of ICP-OES includes the following steps. 1. Sample preparation:

1) 15 ml of HN0₃ (70%) was added in 250 ml of flat bottom boiling flask and weighed Wi. A schematic drawing showing the apparatus could be seen in figure 3. 2) A piece of phosphorus (around Ig) was added into the flat bottom boiling flask afterwards and weighed W₂.

The weight difference between above procedures is the weight of phosphorus. Thus, the weight of phosphorus W_p=W₂-Wi. The weight of phosphorus should be around lg.

3) Then the fiat bottom boiling flask was connected with the Jacket condenser quickly. During the reaction, the phosphorus should be kept below the HN0₃ liquid layer in order to keep the phosphorus from reacting with the air. The cool water should be kept to flow in the jacketed space in the condenser during the reaction. Finally, the phosphorus will react with HNO3 (70%) completely.

4) After the reaction, the Jacketed Condenser was rinsed with distilled water and then removed. 5) The sample solution was transferred from the flat bottom boiling flask to 100 ml ("Svoi") volumetric flask. Then the sample solution was diluted by distilled water till the mark of the flask.

2. Making Standard Curve

The concentration of the standard arsenic solution is lOOOppm. 0.2ml and 0.5 ml of the standard arsenic solution were taken accordingly and diluted to 100ml separately. So the concentrations of the standard solution that determines the standard curve should be 2 and 5 ppm. 3. Measure

The concentration of As(CAs)in the diluted sample solution can be directly measured by ICP-OES with the standard curve.

4. Calculation:

The reading of ICP result is CA_S (mg/1)

The weight of phosphorus is W_p (g)

Sample volume is S_voi (1)

The concentration of As in P is A (ppm) (mg/kg)

Recovery test

In order to know if there are some problems with the method or machine, it is necessary to do the recovery test. The procedure is described below: Prepare the sample as discussed before, before HN0₃ (70%) was added in 250 ml of FB boiling flask, 0.1 ml of 506 ppm standard Arsenic solution was added in the FB boiling flask first.

As can be seen from Table 3-7, all the recovery values are above 80%, so the concentration of Arsenic tested by this method are reliable. It is also able to get the concentration of Arsenic in this industrial phosphorus. From above results, it is clear that the initial concentration of Arsenic in the industrial solid phosphorus is between 52-56 ppm.

Table 3 : Concentration of Arsenic in solid P

Table 4: Results of recovery experiments

Recovery 86.91%

%

Table 5 : Results of recovery experiments

As from solid in the solution 42.6 50.4 from spike 43.3

Recovery 85.6%

%

Table 7: Results of recovery experiments

Sample Sample + Spike

Sample weight 1.1175 0.9855

/g

Sample volume 100 100

/ml

As spiked 0 50.6

Solution cone. 0.588 0.964

/ mg/1

As in solution 58.8 96.4

As in solid 52.6

As from solid in the solution 58.8 51.8

As from spike 44.6

g

Recovery 88.1%

%

Purification

In the experiments of the present invention, the step motor moving speed was controlled from 0.02 mm/min to 0. Imm/min, and the heater temperature range was set from 45°C to 80°C.

Table 8: Concentration of Arsenic in solid P

0.7798 100 0.426 54.6

1.1175 100 0.588 52.6

A phosphorus bar (Bar No.1) was put inside the stainless steel tube. The temperature of the heater is set at 48°C), and speed was set at 0.375mm/minute. Run three times at this condition. Then the two samples from the middle of the bar were taken and test by ICR The results are shown in table 9.

Table 9: Concentration of Arsenic after Purification. (Bar No.1)

In order to improve the purification, more running times have been tried. A new phosphorus bar (Bar No.2) was put inside the stainless steel tube. The temperature of the heater is set at 46°C this time, and speed was set at 0.375mm/minute as last time. Run eight times at this condition. Then the two samples from the middle of the bar was taken and tested by ICP. The results are shown in table 10.

Table 10: Concentration of Arsenic after Purification. (Bar No.2)

A new phosphorus bar (Bar No.3) was put inside the stainless steel tube. The speed was set at 0.375mm/minute. Run three times at 46°C, two times at 48°C, one time at 47°C. Then, the Arsenic concentration was successfully reduced below lOppm. See Table 11.

Table 11 : Concentration of Arsenic after Purification. (Bar No.3)

Mathematical Model and Numerical Simulation:

The inventor of the present invention has developed a mathematical model to simulate the above zone melting purification process. The mathematical model is based on a phase field model for solid-liquid phase transition¹⁴' ¹⁵' ¹⁶ and is described in the following: (1)

c^ + Lyr \φ) ^+ -^l W_Eg \φ) = V^■ [*(0 V Γ] + ^ ^ W (2) ot dt 2 (¾ 3t

The first equation describes the evolution of the two phases: solid and liquid phase. φ(χ,ί) is an unknown continuous function which labels the two phases at the position x and the time /, with φ = 0 in the solid phase, φ = 1 in the liquid phase and

0 < φ < 1 in the interfacial region. The left side describes the time evolution of phase function φ(χ,ί) : M is the mobility coefficient. The force leading to the evolution is given on the right side. The first term describes the effect of the interfacial energy on the dynamics of φ(χ, ί) : K_E and K_s are interfacial energy coefficients; the Laplace ν²φ indicates the variation of φ(χ,ί) in the interface. The second term describes the effect of the potential energy: W_E and W_s are barrier heights for double well potential function g(φ) = ²(\ - φ)² . K_E , K_s , i^ and W_s are phenomenological parameters and can be determined later. The third term describes the force due to the latent heat released in the phase transition: ρ(φ) = φ² (3 - 2φ) is an interpolating function that records the proportion of the latent heat released in an intermediate state, L_y is the latent heat per unit volume and T_M is the melting point.

The second equation, which is just the generalization of a regular heat equation, describes the evolution of the temperature T(x, t) at the position ; and the time t. The first term on the left side describes the time evolution of the temperature: c_p is the heat capacity per unit volume. The second term describes the effect of phase transition on the dynamics of temperature: τ(φ) = ρ(φ) has the same meaning of ρ(φ) . The third term describes the effect of potential energy. The first term on the right side describes the heat conduction by Fourier's Law: permits different thermal conductivities k_L = k(0) and k_s = k(Y) in the bulk liquid and solid phases, respectively, and varies continuously in the interfacial region. The second term on the right side describes the effect of interfacial energy.

All the parameters interpreted here are fixed physical parameters of phosphorus except the unknown functions φ(χ,ί) and T(x, t) .

This system of equations can be derived from the conservation laws for energy and entropy and the fact that the entropy production is positive¹⁶' ¹⁷' ¹⁸.

Let J = t / R , t = t / (R² / tc) , u = (T - T_M) I (L_y I Cp ) , by rescaling the above equations, we obtain the dimensionless form:

where the macroscopic length scale is R , the thermal timescale is R²1 κ , the temperature scale is Ly I c_p . K

is the average of the thermal diffusivities K_L=k_Ll c_p and K_s=k_s I c_p in the liquid and solid, respectively.

In the above equations, the following new parameters were used:

K -K_E + T_MK_S , W = W_E + T_MW_S,

It is easy to find that Q(0) = Q_L =K_L I κ and Q(\) ~Q_S = K_S I κ satisfy Q_L+Q_S=2.

By the thin interface asymptotic analysis¹⁶, special forms of the above parameters could be chosen as follows.

where we require 0.0689 < K_S I K_L < 14.5, c, =19/5, c₂ =173/525.

/ is the interface thickness which can be prescribed by ourselves. Therefore the undetermined parameters K_E, K_s, W_E&nd W_s , can be related to /:

The boundary condition is given by

dT

k(<)— = h(T_h-T) for v-i≤z≤v-t + D and r=R, (5)

dn k(<t>)— = h(T₀ - T) for z < v - t , z > v - l + D and r=R, z=0, z=L, (6) dn

= 0 for z=0, z=L, and r=R. (7)

The first two boundary conditions indicate that the two different directions of heat transfer: In the region covered by the heater, i.e., v - i≤ z≤v - t + D _^ heat transfer occurs between the heater (with temperature T_h ) and phosphorus (with temperature T(z, R, t) at position (z, R) on the boundary); in the other boundary region (z, R) with z < v t or z > v- t + D and z = 0 or L, heat transfer occurs between the water (with temperature T₀ ) and phosphorus (with temperature T(z, R, t) at position (z, R) and T(z, r, t) at position z = 0 or I on the boundary). The last boundary condition — = 0 indicates mass conservation: n is the normal vector of the dn

boundary, h is the heat transfer coefficient, T₀ is the room temperature and T_h is the temperature of the heater.

The four parameters, the length of the heater D, the radius of the phosphorus bar R, the temperature of the heater is T_h , and the moving velocity of the heater v, determine the temperature distribution and the effect of phase separation and therefore can be adjusted to obtain good results of phase separation. By solving the governing equations with different sets of adjustable parameters numerically, both good and bad results can be observed.

In the dimensionless form, by noticing that x = t I R and u

boundary condition can be reduced to be

in the heating region,

hR(T₀ - uL_y I c_p T_M ) otherwise, (9) ^ = 0 (10) dn

By the radial symmetry, it is only needed to consider a cross section along a certain radius. The 3-d problem can be reduced to a 2-d problem with a rectangular computational domain (z, r) s [0, L] x [0, R] . At r = 0 , it is needed = 0 , a„d U= 0.

dr dr

In cylindrical coordinate system, the dimensionless equations are dr ' r δθ' dz dr' 'dz '

1 5 . δφ. 1 a ^■ + 5 1 d</> &Φ δ²Φ

r dr^{^} dr' r² dO² dz² r dr dr² dz²

The explicit finite difference scheme could be applied to solve the differential equations. Suppose 1+1 grid points and J+l grid points are distributed along r-direction and z-direction, respectively. The dimensionless mesh size is uniform, i.e., r=RI I IR = \I I and Δζ = L 111 R = L I (RI) , where L is the length of the phosphorus cylinder in the computational domain. The functions u and φ evaluated at the point (η, z .) = Ar Az) at time t_n = nAt are simply written in the form w V and φ"_ί .. For =/,..., 7-7, y , J- J, the desired function values in the inner points are updated as follows:

ε²τ

The boundary values are updated by introducing ghost points out of the computational domain but next to the boundary:

.

other fs

hR(To - ^H,jLr / c_p - T_M) duⁿ \ _ «' ij+l ^{U 1}, -1 for j=0, J,

Here /?os describes the left end position of the heater and D = d I Az is the number of grid points along z-direction inside the region under heating.

When the temperature at the boundary reaches the melting point, the phase transition from solid to liquid starts to happen.

Now only one parameter remains to be determined. It could be simply assumed that the interface thickness is /=0.1mm. In order to make the computation numerically stable, the time step At is chosen to be about 10^"4 in dimensionless case and 10^"3 second in the dimensional case.

All the parameters used in the above simulations are listed in Table 12

Table 12: List of parameters

Numerical Results:

Starting with room temperature T₀ = 25°C, numerically several issues related to the purification process were investigated and compared with the experimental results.

Quantities of concern numerically are: (1) fusion time: the time for the region under the heater to be fused,

(2) the time for the fusion region to reach stable state ("stable" means the fusion region keeps the same shape and a constant velocity to move),

(3) the velocity of both interface,

(4) the width of fusion region at different radii, and

(5) the lagging distance of both interfaces to both boundaries of the heater.

It is important to make clear how these quantities depend on temperature and moving velocity of the heater. In particular, the quantities that should be concern the most are the fusion time and lagging distance of left interface to left end of heater, since they have a great influence on how to choose the temperature of the heater and moving velocity of the heater.

In the following tests, the length of the heater D is set as 4cm, and the Radius of the phosphors bar is set as 5mm, unless specifically indicated.

When temperature of heater is T_h =50°C and velocity is v= 0. lmm/s, the fusion time is

77.8s. It can be seen that the melting region is entirely under the heater as the heater moves. The melting region almost keeps the same shape while moving (Figure 4). The temperature is shown in Figure 5. Figure 6 shows the temperature in 2D at time t=400s. Figure 7 shows the moving velocity of the fusion region is exactly equal to the velocity of the heater.

As shown in Figure 8, when v=0.2mm/s, the fusion region will lag behind the heater. The characteristics are shown in Table 13.

Table 13 : Characteristics when T_h = 50°C and v = 0.2mm/s

As shown in Figure 9, when v=0.05mm/s, the fusion region can stay in the middle of the heater. The characteristics are shown in Table 14.

Table 14: Characteristics when T_h = 50°C and v = 0.05mm/s Characteristics Quantity Value

Fusion time 77.8 s

Time for attaining stable state 146 s

Velocity of left interface 0.05 mm/s

Velocity of right interface 0.05 mm/s

Width of fusion region at r=R 3.32 cm

Width of fusion region at r=0 3.06 cm

Lagging distance of left interface -2.7 mm

Lagging distance of right interface 6.7 mm

As shown in Figure 10, when v=0.02mm/s, the situation is even better. The characteristics are shown in Table 15.

Table 15: Characteristics when T_h = 50°C and v = 0.02mm/s

Figure 1 1 shows that when v=0.5mm/s, the situation is not very good. Then, as to the case that T_h = 45°C and v=0.05mm/s, the fusion time is longer, i.e., 400.6s. This can be seen in Table 16.

Table 16: Characteristics when T_h = 45°C and v = 0.05mm/s

The width of the melting region is much narrower than the length of the heater. According to Figure 12, the result seems to be good. The temperature is shown is Figure 13. The result is much better if the velocity is slowed down to be 0.02mm/s (Figure 14).

The result is expected to be better if the temperature of the heater is rising. Since for each fixed temperature of the heater, the situation is of course better if the velocity is slowed down. In each case, the numerical result was given when the velocity is as large as possible in consideration of the energy cost:

T_h = 55°C and v = 0.05mm/s (Figure 15, Table 18);

T_h = 60°C and v = 0.05mm/s (Figure 16, Table 19);

T_h = 70°C and v = 0.02mm/s (Figure 17, Table 20);

T_h = 80°C and v = 0.02mm/s (Figure 18, Table 21).

Table 17: Characteristics when T_h = 55°C and v = 0.05mm/s

Table 18: Characteristics when T_h = 60°C and v = 0.05mm/s

Table 19: Characteristics when T_h = 70°C and v = 0.02mm/s Characteristics Quantity Value

Fusion time 25.9 s

Time for attaining stable state 80 s

Velocity of left interface 0.02 mm/s

Velocity of right interface 0.02 mm/s

Width of fusion region at r=R 3.86 cm

Width of fusion region at r=0 3.82 cm

Lagging distance of left interface -0.5 mm

Lagging distance of right interface 1.3 mm

Table 20: Characteristics when T_h = 80°C and v = 0.02mm/s

It can be seen from above that although the fusion time is reduced as the temperature of the heater rises, the velocity has to be slowed down. The reduction in the fusion time is very slight, but the slowdown in velocity is remarkable. This will lead to more energy cost. Therefore, it is suggested that the temperature of the heater is between 45°C and 50°C. This is consistent with the experiments. In Table 21 the numerical results were listed when the temperature of the heater is between 45°C and 50°C and the velocity is fixed at lmm/min.

Table 21 : Numerical results when T_h is between 45°C and 50°C

Of course this heating system can work if the length of heater is enlarged. But the purpose of the present invention is to save the energy. Therefore, it is worth to see whether it works if the length d decreases. The temperature of heater was kept to be 50 degree Celsius and the situation that R=5mm, d=1.5cm and v=0.1mm/s was tested. The result could be seen in Figure 19. The core of the cylinder cannot be melted when d is reduced to be 1cm (see Figure 20).

Figure 21 shows the case wherein R=lcm, d=3cm and v=0.05mm/s. The core of the cylinder cannot be melted when d is reduced to be 2cm (see Figure 22).

Figure 23 shows the case wherein R=1.5cm, d=4cm and v=0.05mm/s. The core of the cylinder cannot be melted when d is reduced to be 3 cm (see Figure 24). Referring to the above results and Table 22, it could be roughly spoken that the ratio of the width of heater to the radius should be greater than 2, preferably greater than 3.

Table 22: Fusion time for different widths of heater and different radii at 50°C

In the above experiments, the radius of cylinder was almost fixed to be 5mm. It is obvious that if the radius is enlarged, for instance, to be 5cm, the volume of phosphorus in every unit length is 100 times more and the efficiency of purification may be increased. But increasing in radius may lead to decreasing in moving velocity since the heat conduction in the cylinder will cost more time. Some numerical results are listed here to illustrate this phenomenon. First, the temperature of the heater is fixed to be 60 degree Celsius and the length of the heater is 10cm. Two different velocities were used to move the heater. The results could be seen in Figure 25 and Figure 26.

When increasing the length of the heater to be 15cm, the results could be seen in Figure 27, Figure 28, and the case wherein the length of the heater is 25cm could be seen in Figure 29. The fusion time and admitted velocities for each length of the heater are listed in Table 23.

Table 23 : Fusion time and velocities for different lengths of heater at a fixed radius

R=5cm and a fixed temperature T_h = 60°C

The effects for R=5mm and R=5cm in the above numerical experiments were mainly tested. The parameters in the listed tables can be used to realize the purification in numerical sense. For some specified cases, for example, the case that the radius is 5mm and the length of heater is 4cm, both numerical results and chemical experiments give the same claim that the moving velocity should be around 0. lmm/s.

In sum, if the radius of the phosphorus bar is enlarged, it should be applied a little higher temperature of the heater to be between 50°C and 60°C, and also slowed down the velocity to be smaller than 0. lmm/s, in order that the melting region can move together with the heater.

The experiments agree with the numerical simulations very well. As shown in Table 24, the inventors performed a series of experiments with different sets of parameters. "Pass" indicates that the melting zone formed. These parameters support the numerical results.

Table 24: Experimental results with different sets of parameters

In fact, the velocity can be raised up to 0. lmm/s (=6mm/min). The most optimized result is: (1) the range of speed of heater is 0.02 mm per second to 0.1 mm per second. (2) the diameter of the phosphorus is 10 mm, (3) the length of heater is 40 mm, (4) the temperature of heater is between 45°C - 60°C. The mathematical model and numerical simulations provide a guideline on how to choose the parameters so that the experiments can be carried out successfully. Although the purpose of the mathematical methods is to determine the parameter for successful heat conduction rather than to predict the concentration of As in phosphorus, it establishes a fundament for this invention. For obtaining a desired As concentration in phosphorus, the persons skilled in the art are able to repeat the running of the heater. All of these fall within the scope of the present invention.

Although the above description exemplifies few embodiments of the present invention, it should be understood that various omissions, substitutions, and changes in the form of the detail of the apparatus, system, and/or method as illustrated as well as the uses thereof, may be made by those skilled in the art, without departing from the spirit of the present invention.

Cited references:

[la] Code and Federal Regulations Title 21, U. S. Food and Drug Administration, Vol. 1, Sec. 73.85, 2009.

[lb] Food Chemicals Codex, (5th edition), The National Academies Press, Washington DC, 2004, pp. 331-332.

[I] Yongli Mi, and Ning Li, "Removal of arsenic from the Industrial phosphorus by Zone Melting Technique", Patent, CHN, ZL03107677.7, 2007

[2] Bruce, N., and Hunt, J. 1993. The Science Corner, University of Guelph, Ontario, Canada.

[3] The Phosphate File, Centre Europeen d'Etudes des Polyphosphates 1997. CEFIC publications.

[4] Phosphorus (properties, occurrence, and uses)— Encyclopedia Britannica online, available at http://search.eb.com/eb/article-5697

[5] Van Wazer, J.R. 1961. Phosphorus and its compounds: Volumes 1 & 2,

Interscience Publisher, Library of Congress No. 58-10100.

[6] Arsenic Poisoning — Encyclopedia Britannica online, available at http ://search. eb. com/eb/article-9009639

[7] Uchida Hiroto, Nishinaka Keiji, and Ogi Katsumi, "Production of high purity phosphorus", Patent, JP 05-043210, 1993

[8] Satou Rikito, and Masagaki Takeshi, "Production of high-purity phosphorus" , JP

05-116916, 1993

[9] H. Rosenhouse, and J. C. Dore, "Process for phosphorus purification" , Patent, US 4,483,746, 1984

[10] R. E. Hall, "Process for purifying yellow phosphorus, US 4,664,896, 1987

[I I] J. Legrand, P. Bourdauducq, and T. Bulinge, "Process for the purification of white phosphorus", Patent, US 5,283,042, 1994

[12] Z. He, and J. Hu, "Dearsenification method of yellow phosphorus", Patent, CN 1035810, 1992

[13] http://www.state.nj.us/dep/dsr/arsenic/home-guide.htm

[14] J.B. Collins, H. Levine, "Diffuse interface model of diffusion-limited crystal growth", Phys. Rev. B31 (1985) 6119.

[15] J.S. Langer, in: G. Grinstein, G. Mazenko (Eds.), "Directions in Condensed

Matter Physics", World Scientific, Singapore, 1986 p. 165.

[16] G.B. McFadden, A.A. Wheeler, D.M. Anderson, "Thin interface asymptotics for an energy/entropy approach to phase-field models with unequal conductivities",

Physica D 144(2000) 154-168.

[17] R.F. Sekerka, "Morphology: from sharp interface to phase field models", Journal of Crystal Growth 264 (2004) 530-540.

[18] S-L. Wang, RF. Sekerka, A.A. Wheeler, B.T. Murray, S.R. Coriell, RJ. Braun,

GB. McFadden, "Thermodynamically-consistent phase-field models for solidification", Physica D 69 (1993) 189.

[19] H. Kramer, K. Schmeiser, "Thermal conductivity of yellow phosphorus in the temperature range + 16 to +80 X" ^', Zeitschrift fur Physikalische Chemie Neue

Folge Bd.35 (1962) 1-9.

[20] C.Y. Ho, R.W. Powell, P.E. Liley, "Thermal conductivity of the elements", (1972)

1, 279-421.

[21] O. Penrose, P.C. Fife, Phsica D 43 (1990) 44.

[22] J.W. Cahn, J.E. Hilliard, J. Chem. Phys 28(1958) 258.

[23] G.W.C. Keye, T.H. Laby, "Tables of Physical and Chemical Constants", (1995) p.90.

[24] I.S. Grigoriev, E.Z. Meilikhov, "Handbook of Physical Quantities", (1997) p.419-448.

[25] http://en.wikipedia.org/wiki/Phosphorus

Claims

2011/097961 The specification PCT/CN2011/000218 The claims

1. A method of production of ultrapure phosphorus by zone melting in a

non-flammable environment, comprising the steps of:

a. Fixing a heater in a non-flammable environment;

b. Fixing a phosphorus bar through the heater;

c. Setting up a motor to control the heater moving along the phosphorus bar;

d. Heating a part of the said bar to become a heat zone

e. Moving the heater along the said bar in order to have the meat zone moving with the heater;

f. Repeat step e if necessary.

2. The method of production of ultrapure phosphorus as mentioned in claim 1,

wherein the length of the heater D, the radius of the phosphorus bar R, the temperature of the heater T_h, and the moving velocity of the heater v shall control the solutions of the following governing equations:

In r<R and 0<z<L,

-L^ = (K_{E +} TK_S ) v - (w_E + TW_S )_G χ ) + ^Lz _{p w}

M dt 2 T_M

_£ρ ^₊^ φ)^₊^χ_Φ)^ = ν .[ φ)ντ]₊κ_Ε ^φ

^P dt dt 2 ^B dt ^E dt

The boundary conditions are:

Ιι(φ)— = h(T_h - T) , for v - t≤z≤v - t + D and r=R, dn

k^)— = h(T₀ - T) , for z < v - t , z > v - t + D and r=R, z=0, z=L,

dn

^ = 0 , for r=R, z=0, z=L;

dn

In these equations, φ(χ, t) is an continuous function which labels the two phases at the position x and the time t, with φ = 0 in the solid phase, φ = 1 in the liquid phase and 0 < φ < 1 in the interfacial region; M is the mobility coefficient, K_E and K_s are interfacial energy coefficients, W_E and W_S are barrier heights for double well potential function g(p) = φ²(\ -φ)² ; ρ(φ) = φ (3 - 2φ) ; L is the latent heat per unit volume and T_M is the melting point; T(x, t) is the temperature distribution at the position x and the time t; c_p is the heat capacity per unit volume, τ{φ) = ρ(φ) ; k^) permits different thermal conductivities k_L = k(0) and k_s = k(\) in the bulk liquid and solid phases, respectively; these governing equations being solved on a cylindrical domain x = (z,r) e [0, L] x [0, R], L is the length of phosphorus cylinder 2011/097961 The specification PCT/CN2011/000218 and R is the radius, z is the axial direction coordinate, r is the radial direction coordinate, D is the length of the heater, n is the outer normal vector of the boundary, h is the heat transfer coefficient, T₀ is the room temperature and T_h is the temperature of the heater.

3. The method of production of ultrapure phosphorus as mentioned in any one of the preceding claims, wherein the said phosphorus bar can be molded by a plastic mold, preferably a high density polyethylene or Teflon tube with a round hole in the middle.

4. The method of production of ultrapure phosphorus as mentioned in any one of the preceding claims, where the length of the heater is at least 2 times, preferably at least 3 times, more preferably at least 4 times, longer than the diameter of the phosphorus bar.

5. The method of production of ultrapure phosphorus as mentioned in any one of the preceding claims, wherein the temperature of heater is between 45°C and 80°C.

6. The method of production of ultrapure phosphorus as mentioned in any one of the preceding claims, wherein the temperature of heater is between 45°C and 60°C.

7. The method of production of ultrapure phosphorus as mentioned in any one of the preceding claims, wherein the temperature of heater is between 45°C and 50°C.

8. The method of production of ultrapure phosphorus as mentioned in any one of the preceding claims, wherein the velocity of the heater is between 0.02mm/s and O. lmm/s.

9. The method of production of ultrapure phosphorus as mentioned in any one of the preceding claims, wherein the velocity of the heater is between 0.05mm/s and

O. lmm/s.

10. A apparatus for the production of ultrapure phosphorus as mentioned in the

preceding claims, which comprises in a non-flammable environment:

g. heater;

h. a phosphorus bar fixed through the heater and between fixing units; i. a motor set up to control the heater moving along the phosphorus bar.