Next Article in Journal
CO2 Separation and Capture Properties of Porous Carbonaceous Materials from Leather Residues
Next Article in Special Issue
Simulations and Measurements of Human Middle Ear Vibrations Using Multi-Body Systems and Laser-Doppler Vibrometry with the Floating Mass Transducer
Previous Article in Journal
The Synthesis of α-MoO3 by Ethylene Glycol
Previous Article in Special Issue
Homogenized Elastic Properties of Graphene for Small Deformations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Electron Beam Melting and Refining of Metals: Computational Modeling and Optimization

Institute of Electronics, Bulgarian Academy of Sciences, 72 Tzarigradsko shosse, Sofia 1784, Bulgaria
*
Author to whom correspondence should be addressed.
Materials 2013, 6(10), 4626-4640; https://doi.org/10.3390/ma6104626
Submission received: 15 August 2013 / Revised: 7 October 2013 / Accepted: 9 October 2013 / Published: 18 October 2013
(This article belongs to the Special Issue Computational Modeling and Simulation in Materials Study)

Abstract

:
Computational modeling offers an opportunity for a better understanding and investigation of thermal transfer mechanisms. It can be used for the optimization of the electron beam melting process and for obtaining new materials with improved characteristics that have many applications in the power industry, medicine, instrument engineering, electronics, etc. A time-dependent 3D axis-symmetrical heat model for simulation of thermal transfer in metal ingots solidified in a water-cooled crucible at electron beam melting and refining (EBMR) is developed. The model predicts the change in the temperature field in the casting ingot during the interaction of the beam with the material. A modified Pismen-Rekford numerical scheme to discretize the analytical model is developed. These equation systems, describing the thermal processes and main characteristics of the developed numerical method, are presented. In order to optimize the technological regimes, different criteria for better refinement and obtaining dendrite crystal structures are proposed. Analytical problems of mathematical optimization are formulated, discretized and heuristically solved by cluster methods. Using important for the practice simulation results, suggestions can be made for EBMR technology optimization. The proposed tool is important and useful for studying, control, optimization of EBMR process parameters and improving of the quality of the newly produced materials.

1. Introduction

Electron beam melting and refining (EBMR) in a vacuum using an intense electron beam is a widely used, ecologically-friendly method in special electro metallurgy for new materials fabrication: the production of pure metals and special alloys [1,2,3,4,5,6]. The obtained materials have improved chemical compositions, structures, properties, and have many applications in medicine, aerospace engineering, nuclear industry, instrument engineering, etc. The electron beam melting and refining of metals is accomplished in a vacuum chamber [the working vacuum pressure is about (5–8) × 10 - 3 Pa] using electron beams as a heating source. The raw material is melted, refined and re-solidified in a water-cooled copper crucible, which is used due to the increasing demand for purity of the transformed materials (Figure 1). The electrons fall on the front side of the feeding rod (raw material) and heat it. Then drops of molten metal fall in the crucible. The top surface of the molten metal in the water-cooled crucible (G1) is also heated by the electrons (Figure 1). By using a water-cooled pulling mechanism (Figure 1), the operator can withdraw the bottom of the growing formed pure metal ingot. This allows the liquid pool surface to be maintained at a constant level suitable to be viewed by the operator. This method provides a high refining level (high level of impurities removal—gases, metals and non-metals), chemical composition homogeneity and optimal structure of the cast pure ingots [1,2,3,4,5,6].
Figure 1. Scheme of Electron beam melting and refining (EBMR) process. G1—the top surface of the formed pure ingot; G2—the interface molten ingot/ water-cooled crucible side wall; G3—the interface ingot/ vacuum; G4—the interface ingot/water-cooled puller.
Figure 1. Scheme of Electron beam melting and refining (EBMR) process. G1—the top surface of the formed pure ingot; G2—the interface molten ingot/ water-cooled crucible side wall; G3—the interface ingot/ vacuum; G4—the interface ingot/water-cooled puller.
Materials 06 04626 g001
Despite development of the technology, there are many unsolved problems concerning mechanisms and relationships of macroscopic heat and mass transfer during heating free liquid surface with intensive energy flow, that still exist [3,4,5,6,7]. For the EBMR process and equipment optimization and the price of the produced pure materials, the questions concerning the improved energy efficiency are crucial. The answer of these questions depends on the detailed study of the heat transfer processes taking place in the zone of interaction beam-material, on the processes and factors limiting the geometry of the molten pool and the precise evaluation of the temperature field dynamics and behavior of the metals and their compounds during the EBMR. Due to difficulties in acquiring real time information about the processes in the molten pool [1,2], the successful application and optimization of EBMR depends also on the adequate mathematical modeling of the processes that allow for studying the influence of many technological parameters and the various limiting factors. Stationary and quasi-stationary heat models in this field are already published [8,9,10,11,12,13,14,15] and reveal the importance of some process characteristics. In [16] a finite element based mathematical model to predict the evolution of the temperature and stress during EBMR of solar-grade silicon ingots is also presented. Applying a heuristic approach the flatness of the liquid/solid boundary is examined and studied [17,18].
A time-dependent thermal mathematical model of the electron beam drip melting and refining process is developed and described in this paper. The model is a continuation and an extension of our quasi-steady-state heat model [10,12,15]. The boundary conditions on the ingot interfaces are formulated assuming different mechanisms of the heat transfer. The temperature variations of the thermo-physical parameters (heat capacity, thermal conductivity, etc.) are taken into account in the heat model. An optimization scheme, based on the heat model, is proposed and applied for quality improvement of new materials obtained by EBMR.

2. Modeling and Mathematical Setting of the Problem

2.1. Main Equation

Let H be the ingot’s height, R—the ingot’s radius (Figure 1), and F—the maximal heating time that is chosen. The investigated temperature distributions are along cylindrical ingots (samples) and an assumption for angle symmetry is applied. Then in:
Ω = { x = ( r , z , t ) | 0 r R , 0 z H , 0 < t F }
T ( r , z , t ) denotes the temperature at time t at the points with height z and polar distance r (Figure 2). In the inner points of Ω, the temperature distribution T ( x ) is described by the heat Equation (2) in cylindrical coordinates:
1 r r r T r + 2 T z 2 + V a T z = ρ C p λ T t
where ρ [kg/m 3 ] is the density of the metal; C p [W · s/kg · K] is the heat capacity; λ [W/m · K] is the thermal conductivity; a [m 2 /s] is the thermal diffusivity and a = λ / ( ρ C p ) = 1 / k .
The last term in Equation (2) V a T z shows the casting, i.e., the heat added by the poured molten metal (from the melting raw material—Figure 1) into the crucible. It is given by the thermal energy transfer from the material moving with velocity V, coincident to the z-axis (Figure 1 and Figure 2). When V 0 , the model describes the solidification of the cast ingot in the water-cooled crucible—“drip melting process" (Figure 1), while the case V = 0 describes the “disks melting method" (no material is added to the ingot by pouring) [2,10]. The thermo-physical parameters λ and C p of the investigated metals are modeled as functions of the temperature (by linear regression method) using experimental data [19,20,21,22].
Figure 2. Cylindrical metal ingot, T ( r , z , t ) —the temperature at the moment t in the ingot’s point with coordinates ( r , z ) .
Figure 2. Cylindrical metal ingot, T ( r , z , t ) —the temperature at the moment t in the ingot’s point with coordinates ( r , z ) .
Materials 06 04626 g002

2.2. Equations for the Boundary Conditions

The boundaries G1, G2, G3, G4 (Figure 1) are described as follows:
Ω 1 = { ( r , H , t ) | 0 r R , t 0 }
Ω 2 = { ( R , z , t ) | H - Q z H , t 0 }
Ω 3 = { ( R , z , t ) | 0 z H - Q , t 0 }
Ω 4 = { ( r , 0 , t ) | 0 r R , t 0 }
The boundary Equations that correspond to these interfaces are:
λ ( T ( x ) ) T z ( x ) = + P s u r f ( r , t ) - α σ T 4 ( x ) - T r o o m 4 - C p ( T ( x ) ) · W v · T ( x ) , x Ω 1
λ ( T ( x ) ) T r ( x ) = T 1 r ( x ) , x Ω 2
λ ( T ( x ) ) T r ( x ) = - α σ T 4 ( x ) - T r o o m 4 , x Ω 3
λ ( T ( x ) ) T z ( x ) = T 2 z ( x ) , x Ω 4
The heat transfer through G1 boundary (Equation (7)), where the radiation losses predominate, is described by the Stefan-Boltzmann law [3,12]. In Equation (7) α is the metal’s emissivity, σ = 5 . 6704 × 10 - 8 [J/s · m 2 · K 4 ] is the Stefan-Boltzmann constant. P s u r f ( r , t ) is the e-beam power density function and usually Gaussian-like distribution is used [2], defined by the beam diameter’s value. The vapor losses are calculated by C p · W v · T (Equation (7)). The weight loss velocity W v [kg/s · m 2 ] is evaluated using experimental data for the material losses obtained during EBMR process; T r o o m is the room temperature (300 K).
Equation (8) describes Newton’s type of heat transfer through G2—the interface molten metal/water-cooled crucible side wall. The gradient T 1 / r is equal to λ 2 ( T w a t e r - T ( x ) ) / Δ w a l l ; λ 2 is the conductivity of the crucible’s material (copper). When the liquid metal pool does not reach the interface G2 and/or its volume is small this boundary condition can be neglected (i.e., Q = 0 , Figure 1). Δ wall is the width of the crucible side wall and T w a t e r is the mean water temperature in the water-cooling system, ( T w a t e r = 300 K).
The thermal transfer through the interface ingot/ vacuum (G3) is described by the Stefan-Boltzmann law and is represented by Equations (9) and (10) describes the heat transfer through the ingot bottom (G4)—area with an ideal heat contact. The gradient T 2 / z is equal to λ 2 ( T w a t e r - T ( x ) ) / Δ u n d e r , where Δ u n d e r is the width of the water-cooled puller (Figure 1).
The initial condition gives information about the temperature distribution at the moment t = 0 s. Usually, the temperature along the metal ingot at the first moment is a constant and is equal to the mean room temperature T r o o m . This condition is described by Equation (11):
T ( x ) = T 0 ( r , z ) , x Ω 0
Ω 0 = { ( r , z , 0 ) | 0 r R , 0 z H }

2.3. Heat Streams

The heat streams through the boundaries at the moment t are calculated as follows:

2.3.1. Energy Losses (Radiation and Vapor Losses) through G1:

S G 1 , r ( t ) = [ 0 , 2 π ] × [ 0 , R ] α σ ( T 4 ( r , H , t ) - T r o o m 4 ) r d r d ϕ
S G 1 , v ( t ) = [ 0 , 2 π ] × [ 0 , R ] C p ( T ( r , H , t ) ) · W v · T ( r , H , t ) r d r d ϕ

2.3.2. Streams through the Boundaries G2, G3 and G4:

S G 2 ( t ) = [ 0 , 2 π ] × [ H - Q , H ] λ ( T ( R , z , t ) ) T r ( R , z , t ) R d z d ϕ
S G 3 ( t ) = [ 0 , 2 π ] × [ O , H - Q ] λ ( T ( R , z , t ) ) T r ( R , z , t ) R d z d ϕ
S G 4 ( t ) = [ 0 , 2 π ] × [ 0 , R ] λ ( T ( r , 0 , t ) ) T z ( r , 0 , t ) r d r d ϕ

2.3.3. The Incoming Heat is Determined by the Heating Beam Energy Distribution P s u r f and by the Heat Added by the Poured Molten Drops P a d d :

S p ( t ) = [ 0 , 2 π ] × [ 0 , R ] P s u r f ( r , t ) r d r d ϕ
The additional heat due to the added liquid drops is calculated by:
P a d d = π R r o d 2 V ρ ( C p * T m e l t + q )
where R r o d is the radius of the feeding rod (Figure 1); q [W · s/kg] is the specific melting heat of the metal; T m e l t is the melting temperature of the investigated metal; C p * is an average of the heat capacity. The heat transfer processes during EBMR become stationary when the difference between the input heat and total energy losses through the boundaries becomes close to 0.

3. Numerical Method

A modified Pismen-Rekford method [23] is developed and applied for solving the problem Equations (2) and (7–11) and for calculating the temperature fields in the cast ingot. A regular net of points W h 1 , h 2 , τ in the domain Ω concerning Equations (2) and (7–11) is made:
W h 1 , h 2 , τ = { ( r i , z j , t n ) | r i = i h 1 , z j = j h 2 , t n = n τ ; i = 0 , N ¯ , j = 0 , M ¯ , n = 0 , P ¯ }
T i , j n denotes the approximated temperature in the point ( r i , z j , t n ) of the grid. Additionally, T n is the matrix of M × N points, corresponding to the temperature distribution at the moment t n = n τ . For the purposes of the algorithm intermediate planes T n + 1 / 2 (which are not included in the net W h 1 , h 2 , τ ) are built. Approximations T n of T ( r , z , τ n ) , continuously starting from n = 0 to n = P , are made. Geometrically, “jumps" are continuously made on rectangles Ω n = { ( r , z , τ n ) | 0 r R , 0 z H } , parallel to Ω 0 . For each n, two jumps are associated and also two mathematical problems ( A ) and ( B ) have to be solved:
( A ) k T n + 1 2 - T n 0 . 5 τ = Λ 1 T n + 1 2 + Λ 2 T n
( B ) k T n + 1 - T n + 1 2 0 . 5 τ = Λ 1 T n + 1 + Λ 2 T n + 1 2
Here Λ 1 and Λ 2 are difference operators, corresponding to 1 r r r T r and V a T z , respectively:
Λ 2 T i , j n = T i , j - 1 n - 2 T i , j n + T i , j + 1 n h 2 2
Λ 1 T i , j n = T i - 1 , j n - 2 T i , j n + T i + 1 , j n 2 h 1 2 + 1 i h 1 T i + 1 , j n - T i - 1 , j n 2 h 1
Several linear systems are solved in order to handle ( A ) and ( B ) . For the problem ( A ) each system consists of equations for the points { T i , j n + 1 2 } i = 0 , N ¯ along Ω n (j and n are fixed, Figure 3). For the problem ( B ) such a system consists of equations concerning the points { T i , j n + 1 } i = 0 , M ¯ (i and n are fixed, Figure 3). The obtained linear systems are three-diagonal and are solved via the Thomas method. The numerical scheme developed is absolutely stable and implicit in terms of τ due to the definition of problems ( A ) and ( B ) proposed in the non-stationary mathematical model.
Figure 3. Approximation patterns in the inner points for problems (A) and (B). In the left figure for (A)—the planes are Ω n , Ω n + 1 / 2 ; for (B)—the planes are Ω n + 1 / 2 , Ω n + 1 .
Figure 3. Approximation patterns in the inner points for problems (A) and (B). In the left figure for (A)—the planes are Ω n , Ω n + 1 / 2 ; for (B)—the planes are Ω n + 1 / 2 , Ω n + 1 .
Materials 06 04626 g003
The thermal conductivity λ and the specific heat capacity C p for each investigated metal are implemented in the presented model as functions of the temperature (see 2.1, Table 1). When solving (A) and (B), temperature distribution from the previous iteration is used for evaluation of λ and C p .

4. Optimization Model and Different Criteria

The knowledge of the geometry of the crystallization front (liquid/solid boundary) is very important for studying and optimizing the quality of the obtained pure metal after the EBMR process. The flatness of the liquid/solid contour is directly connected to the quality of the structure of the obtained metal. The flat crystallization front permits the formation of vertical dendrite structure and uniform impurities’ displacement toward the ingot top surface (Figure 4, Figure 5 and Figure 6). When the liquid/solid boundary is deep in the center part of the metal block and shallow at the periphery of the cast ingot, a non-uniform structure along the ingot radius will be formed. The flatness of the temperature lines in a vertical cross-section of the cylindrical metal ingot for a fixed moment of time depends on the values of T r . For fixed moment of the heating t = t f and height z = z f , T ( r ; z f ; t f ) is a decreasing function of the variable r and T r 0 . Our aim is to maximize the flatness of the contour liquid metal/ solid metal.

4.1. Different Criteria

In this paragraph we discuss different criteria aiming to obtain maximum flatness of the crystallization front shape. One approach is to minimize the average of - T r over an area that includes the liquid/ solid boundary. Then, the criterion is:
A ( P b , r b , V ) = - D μ ( r , z , t ) T r ( r , z , t ) d r d z d t m i n
where T is the temperature field determined by the solution of Equations (2) and (7–11) and μ ( r , z , t ) is a weight function. The control variables are the casting velocity V, the beam power P b and the beam radius r b · P b is connected with the heating beam energy distribution P s u r f . For a fixed heating time t f , the beam energy distribution P s u r f ( r , t f ) is a Gaussian-like function and:
0 R P s u r f ( r , t f ) d r = P b
Minimization of Equation (25) can be realized over an area D that includes the crystallization front (the molten pool contour, Figure 4a):
D = [ R 1 , R 2 ] × [ H 1 , H 2 ] × [ F 1 , F 2 ] Ω
In a more complicated case D (Figure 4b) can be defined as:
D = [ R 1 , R 2 ] × [ H 1 , H 2 ] × [ F 1 , F 2 ] { ( r , z , t ) | T ( r , z , t ) T m e l t }
Figure 4. Geometrical characteristics on a half of the vertical ingot cross-section, R 1 = 0 ; 1—solid metal; 2—liquid metal; 3—projection of D onto ( r , z ) plane; 4—liquid/solid metal contour. (a) D = [ R 1 , R 2 ] × [ H 1 , H 2 ] × [ F 1 , F 2 ] Ω ; (b) D = [ R 1 , R 2 ] × [ H 1 , H 2 ] × [ F 1 , F 2 ] { ( r , z , t ) | T ( r , z , t ) T m e l t } .
Figure 4. Geometrical characteristics on a half of the vertical ingot cross-section, R 1 = 0 ; 1—solid metal; 2—liquid metal; 3—projection of D onto ( r , z ) plane; 4—liquid/solid metal contour. (a) D = [ R 1 , R 2 ] × [ H 1 , H 2 ] × [ F 1 , F 2 ] Ω ; (b) D = [ R 1 , R 2 ] × [ H 1 , H 2 ] × [ F 1 , F 2 ] { ( r , z , t ) | T ( r , z , t ) T m e l t } .
Materials 06 04626 g004
Another approach is to minimize the mean curvature of the curve Γ ( t f ) , corresponding to the liquid/ solid boundary (Figure 4) at a fixed moment of heating t f :
0 S μ 2 ( s ) | | d Γ ( t f ) d s r ( s ) , z ( s ) | | d s m i n
where:
Γ ( t f ) = { ( r , z ) | T ( r , z , t f ) = T m e l t }
or to minimize the mean curvature over a heating time interval:
t 1 t 2 μ 1 ( t ) 0 S μ 2 ( s ) | | d Γ ( t ) d s r ( s ) , z ( s ) | | d s d t m i n
For Γ ( t ) , s is a natural parameter for the curve and μ 1 , μ 2 are weight functions. The weight functions in Equations (25) and (31) can be defined in different ways taking into account the peculiarities of the electron beam melting process of the investigated metal. The used weight function in the paper is constant.

4.2. Discretization of the Criteria

The proposed criteria for achieving flatness of the crystallization front shape can be discretized synchronically to the proposed numerical scheme for the heat model. For the criterion (25) and D, described by Equation (27) or Equation (28), the multidimensional trapezoidal rule can be used. In the case of Equation (27) lets denote:
S i , j n = μ ( t n ) T i , j n - T i - 1 , j n h 1 + O ( h 1 )
where { T i , j n } i = 1 , N ¯ , j = 1 , M ¯ n = 1 , P ¯ is the discrete temperature field calculated by the modified Pismen-Rekfort scheme. The subnet of W h 1 , h 2 , τ which is used to discretize the criterion is:
W h 1 , h 2 , τ 1 = { ( r i , z j , t n ) | r i = i h i , z j = j h 2 , t n = n τ ; i = N 1 , N 2 ¯ , j = M 1 , M 2 ¯ , n = P 1 , P 2 ¯ }
In Equation (33), N i corresponds to R i , M i corresponds to H i , and P i corresponds to F i , i = 1 , 2 . Formula ( 34 ) presents the approximation of Equation (25):
- 8 h 1 h 2 τ A = S N 1 , M 1 P 1 + S N 1 , M 1 P 2 + S N 1 , M 2 P 1 + S N 1 , M 2 P 2 + S N 2 , M 1 P 1 + S N 2 , M 1 P 2 + S N 2 , M 2 P 1 + S N 2 , M 2 P 2 + + 2 i = N 1 + 1 N 2 - 1 [ S i , M 1 P 1 + S i , M 2 P 1 + S i , M 1 P 2 + S i , M 2 P 2 ] + j = M 1 + 1 M 2 - 1 [ S N 1 , j P 1 + S N 2 , j P 1 + S N 1 , j P 2 + S N 2 , j P 2 ] + n = P 1 + 1 P 2 - 1 [ S N 1 , M 1 n + S N 1 , M 2 n + S N 2 , M 1 n + S N 2 , M 2 n ] + 4 i = N 1 + 1 N 2 - 1 j = M 1 + 1 M 2 - 1 [ S i , j P 1 + S i , j P 2 ] + j = M 1 + 1 M 2 - 1 n = P 1 + 1 P 2 - 1 [ S N 1 , j n + S N 2 , j n ] + n = P 1 + 1 P 2 - 1 i = N 1 + 1 N 2 - 1 [ S i , M 1 n + S i , M 2 n ] + 8 i = N 1 + 1 N 2 - 1 j = M 1 + 1 M 2 - 1 n = P 1 + 1 P 2 - 1 S i , j n
For discrete calculation of the criteria (29) and (31), the curve Γ is approximated by a set of points. Then the curve can be interpolated by linear regression method or some interpolation method.
After the criterion discretizing and choosing the optimization variables (the input power P b , the e-beam radius r b , the casting velocity V, etc.) for solving the defined optimization problem, an algorithm is needed. Heuristic methods (clustering optimization technique, genetic algorithms, etc.) are appropriate because the dependence of the criteria on the control variables is implicit.

5. Results and Discussion

The whole developed process tool of modeling, simulation and optimization of the process of EBMR of metals can be summarized as follows:
  • Theoretical mathematical modeling of the heating processes during EBMR of metals and alloys, describing the equations of the model;
  • Construction of a modified Pismen-Rekford numerical scheme of the model;
  • Development of a software for model simulation via the numerical scheme;
  • Verification of the model;
  • Development of analytical criteria for optimization of the quality of the obtained new materials by achieving flatness of the crystallization front shape;
  • Discretization of the criteria synchronically to the numerical discretization of the model and computer implementation;
  • Choice of heuristic optimization techniques and synchronization of all developed computer programs;
  • Tests using different combinations of criteria, control variables and optimization techniques.
A corresponding computer program, based on the proposed time–dependent heat model, is developed for the study of the thermal processes at EBMR of metals. In this paper, results for EBMR of copper and hafnium are presented and discussed. The temperature variations of the thermal conductivity λ and the heat capacity C p for Cu and Hf are estimated using experimental data [19,21]. The obtained dependencies are presented in Table 1.
Table 1. Material characteristics used in calculations.
Table 1. Material characteristics used in calculations.
ParameterHfCu
T m e l t , K25061356
λ , [ W m · K ] 25 . 17 , T < 293 K 26 . 05 - 1 . 54 10 2 T + 8 . 95 10 6 T 2 + 8 . 46 * 10 2 T , 293 K T 2000 K 31 . 3409 , T > 2000 K 439 . 1 - 0 . 0937 T , T < 1500 K 320 , T > 1500 K
C p , [ W · s g · K ] 0 . 14225 , T < 293 K 23 . 347 + 8 . 004 * 10 - 3 * T - 1 . 058 * 10 4 * T - 2 178 . 49 , 293 K T 2000 K 0 . 220474 , T > 2000 K 0 . 388 , T < 293 K 0 . 36 + 9 . 96 10 4 T , 293 K T 1356 K 0 . 494 , T > 1356 K
Simulation in copper (Cu) cylindrical ingots ( 2 R = 60 mm, H = 100 mm) for e-beam powers P b = 20 kW, beam radius r b = 12 mm and for 10 min heating time at V = 0 mm/min is made. The temperature fields are computed and the profiles of the molten metal pool are investigated. In Figure 5a calculated temperature field in a vertical cross-section of the ingot at the 10th minute of the heating is shown. The solid yellow curve represents the temperature level contour corresponding to the copper melting temperature T m e l t = 1356 K.
Figure 5. Temperature distribution in vertical cross-section of Cu cylindrical sample. The solid yellow curve ( T m e l t = 1356 K) is the liquid/solid boundary in the metal.
Figure 5. Temperature distribution in vertical cross-section of Cu cylindrical sample. The solid yellow curve ( T m e l t = 1356 K) is the liquid/solid boundary in the metal.
Materials 06 04626 g005
A series of experiments realizing different technological regimes are performed in our laboratory “Physical problems of electron beam technologies", Institute of electronics, Bulgarian Academy of Sciences, using ELIT 60 equipment for EBMR. Liquid pool shape variations (experimental and simulation data) in copper samples ( 2 R = 60 mm, H = 50 mm) vs. the refining time are shown in Table 2 at heating with different e-beam powers of 10 kW and 15 kW, r b = 12 mm and V = 0 mm/min. The obtained experimental results, investigated (after EMBR) by a metallographic etching method, are compared with simulation data, obtained by the presented heat model. The results about the geometry of the crystallization front shape—diameters d m and depths h m of the liquid metal pool are shown in Table 2. A good correspondence between calculated and experimentally obtained shapes of the crystallization front is observed.
Table 2. Experimental and simulation data about the geometry of the molten pool during EBMR of copper. E—experimental data; S—simulation data.
Table 2. Experimental and simulation data about the geometry of the molten pool during EBMR of copper. E—experimental data; S—simulation data.
P b , kWCharacteristicsat 10th minat 15th minat 20th minat 30th min
ESESESES
10 d m , mm5045606060606060
h m , mm1312141417151915
15 d m , mm6060606060606060
h m , mm2020222119202020
Using the developed simulation tool important data about the geometry of the liquid metal pool, heat streams through the boundaries, temperature fields in vertical and/or horizontal cross-sections of the cylindrical ingots during the EBMR process are obtained. Calculated and experimentally obtained crystallization front shapes are compared and a good correspondence is observed. Conclusions concerning the influence of a variety of regime parameters (such as e-beam radius, beam power, casting velocity, etc.) are made based on the simulation results [13,15,24]. They can be used for optimization of the heat streams according to the process requirements and for choosing proper process conditions. They are used in the developed optimization approach concerning the flatness of the crystallization front shape.
Using the developed tools and corresponding computer programs numerical experiments for EBMR of hafnium are made. The ingot’s dimensions are H = 50 mm and R = 30 mm, the total heating time is 6 min, V = 0 mm/min. The radius of the e-beam is 18 mm and the beam power is 18, 20, 22, 24 kW. Simulation results for the heat transfer processes are investigated. It is observed that for each beam power the molten pool is mostly flat at the moment before the liquid pool reaches G2 boundary for a first time (Table 3). For all the beam powers the molten pool periodically reaches and withdraws G2.
Table 3. First moment of contact between the molten pool and G2 boundary.
Table 3. First moment of contact between the molten pool and G2 boundary.
No. P b , kWTime of first contact, s
118300
220230
322180
424150
Figure 6. Temperature distribution in a vertical cross-section of Hf ingot ( H = 50 mm, R = 30 mm) at the 295th second of EBM process for P b = 18 kW, r b = 18 mm, T m e l t = 2506 K.
Figure 6. Temperature distribution in a vertical cross-section of Hf ingot ( H = 50 mm, R = 30 mm) at the 295th second of EBM process for P b = 18 kW, r b = 18 mm, T m e l t = 2506 K.
Materials 06 04626 g006
Using the obtained simulation results, one dimensional cluster optimization technique for EBMR of Hf is performed. The control variable is the beam power P b . The chosen criterion is Equation (25) with a constant weight function. The investigated range for the beam power is [18 kW, 24 kW] and the cluster analysis is made for three different areas D: [ 0 , R ] × [ 0 . 7 H , 0 . 9 H ] × [ 0 . 2 F , F ] ; [ 0 , R ] × [ 0 . 7 H , 0 . 9 H ] × [ 0 . 2 F , 0 . 6 F ] and [ 0 . 2 R , 0 . 5 R ] × [ 0 . 7 H , 0 . 9 H ] × [ 0 . 2 F , 0 . 6 F ] . The choice of the areas D is based on the results obtained by the heat model (Figure 4 and Figure 5). The results obtained by the one dimensional cluster optimization show that for all three different areas D, the flattest molten pool (according to the criterion (25)) is observed for P b = 18 kW (Figure 6). So a possible technological suggestion is that for the investigated technological regimes a lower e-beam power (18 kW) has to be kept.
The cluster optimization technique is applied to optimize the technological regime for electron beam melting of titanium. The criterion (25) is minimized for P b = 11 kW [25]. This result is compared and coincides to the value obtained by experimental data and statistical method [26]. An important advantage of the optimization scheme proposed is that experimental data for chemical analysis of the impurities’ concentrations for different technological regimes are not needed. The developed optimization approach can be applied to control and suggest proper technological parameters for electron beam melting of different metals for improving the quality of the obtained new materials.

6. Conclusions

A time-dependent thermal mathematical model, corresponding numerical method and computer program for simulation of heat transfer in metal ingots during EBMR are developed and presented. The heat model is 3D axis-symmetric and different heat transfer mechanisms through the boundaries are assumed. The developed numerical scheme is absolutely stable and implicit in terms of the time. The non-stationary heat model is programmed and corresponding computer software based on the model is also developed. The mathematical model allows acquiring important for the practice data and dependencies, which are otherwise difficult to obtain through EBMR experimental study (such as liquid pool geometry, energy losses, temperature distributions in metal ingots, etc.). The developed tool is employed to predict the evolution of the temperature, streams and crystallization front shape during electron beam melting and refining of different metals (hafnium, copper). Calculated results are compared to experimental data and good correspondence is observed. In order to optimize the EBMR process, different criteria, concerning the flatness of the crystallization front shape, are proposed and discussed. Analytical optimization problems are formulated, discretized using the heat model and solved by a cluster optimization method. The developed computational modeling and optimization tools and results obtained are important and give opportunity for better understanding, studying and optimizing the EBMR process and quality improvement of the purified materials produced by this expensive modern technology.

Acknowledgements

The authors express gratitude to V. Vassileva and R. Nikolov for the performed experiments. This research was funded by the National Fund for Scientific Research at the Ministry of Education, Youth and Science of Republic of Bulgaria under contract DO 02-127 (BIn-5/2009) and contract DNTS/Slovakia/ 01-10/2011.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Mitchell, A.; Wang, T. Electron Beam Melting Technology Review. In Proceedings of the Electron Beam Melting and Refining—State of the Art 2000, Reno, NV, USA, 29–31 October 2000; Bakish, R., Ed.; Bakish Materials Corporation: Englewood, NJ, USA, 2000; pp. 2–13. [Google Scholar]
  2. Mladenov, G. Electron and Ion Technologies (in Bulgarian); Academic Public House: Sofia, Bulgaria, 2009. [Google Scholar]
  3. Mladenov, G.; Koleva, E.; Vutova, K.; Vassileva, V. Experimental and theoretical studies of electron beam melting and refining. In Practical Aspects and Applications of Electron Beam Irradiation; Nemtanu, M., Brasoveanu, M., Eds.; Transword Research Network: Trivandrum, India, 2011; pp. 43–93. [Google Scholar]
  4. Murr, L.; Gaytan, S.; Ramirez, D.; Martinez, E.; Hernandez, J.; Amato, K.; Shindo, P.; Medina, F.; Wicker, R. Metal fabrication by additive manufacturing using laser and electron beam melting technologies. J. Mater. Sci. Technol. 2012, 28, 1–14. [Google Scholar] [CrossRef]
  5. Bakish, R. Electron Beam Melting 1995 to 2005. In Proceedings of the 7th International EBT Conference, Varna, Bulgaria, 5–10 June 2006; pp. 233–240.
  6. Choudhury, A.; Hengsberg, E. Electron beam melting and refining of metals and alloys. ISIJ Int. 1992, 32, 673–681. [Google Scholar] [CrossRef]
  7. Suhas, P. Computational modeling of flow and heat transfer in industrial applications. Int. J. Heat Fluid Flow 2002, 23, 222–231. [Google Scholar]
  8. Bellot, J.P.; Floris, E.; Jardy, A.; Ablitzer, D. Numerical Simulation of the E.B.C.H.R. Process. In Proceedings of the International Conference Electron Beam Melting and Refining—State of the Art 1993, Reno, NA, USA, 3–5 November 1993; Bakish, R., Ed.; Bakish Materials Corporation: Englewood, NJ, USA, 1993; pp. 139–152. [Google Scholar]
  9. Westerberg, K.; Clelland, M.; Finlayson, B. Finite element analysis of flow, heat transfer, and free interfaces in an electron-beam vaporization system for metals. Int. J. Numer. Methods Fluid. 1998, 26, 637–655. [Google Scholar] [CrossRef]
  10. Vutova, K.Zh.; Vassileva, V.I.; Mladenov, G.M. Simulation of the heat transfer process through treated metal, melted in a water-cooled crucible by an electron beam. Vacuum 1997, 48, 143–148. [Google Scholar] [CrossRef]
  11. Gupta, G.; Marwah, R. Estimation of Effective Thermal Conductivity of Hot Liquid Refractory Metal Generated by Electron Beam Heating. In Proceedings of the International DAE-BRNS Symposium on Electron Beam Technology and Applications, Mumbai, India, 28–30 September 2005; Mascarenhas, M., Kandaswamy, E., Sharma, A., Chakravarthy, D.P., Eds.; Allied Publishers PVT Limited: Mumbai, India, 2005; pp. 371–375. [Google Scholar]
  12. Vutova, K.Zh.; Koleva, E.G.; Mladenov, G.M. Simulation of thermal transfer process in cast ingots at electron beam melting and refining. Int. Rev. Mech. Eng. 2011, 5, 257–265. [Google Scholar]
  13. Vutova, K.Zh.; Donchev, V.D.; Vassileva, V.I.; Mladenov, G.M. Thermal transfer process through treated metal regenerated by electron beam melting and refining. J. Metal Sci. Heat Treat 2013, in press. [Google Scholar]
  14. Ou, J.; Chatterjee, A.; Reilly, C.; Maijer, D.M.; Cockcroft, S.L. Supplemental Proceedings: Materials Processing and Interfaces; John Wiley Sons, Inc.: Hoboken, NJ, USA, 2012; pp. 871–878. [Google Scholar]
  15. Vutova, K.; Donchev, V.; Vassileva, V.; Mladenov, G. Influence of Process and Thermo-physical Parameters on the Heat Transfer at Electron Beam Melting of Cu and Ta. In Supplemental Proceedings: Materials Processing and Interfaces; John Wiley Sons, Inc.: Hoboken, NJ, USA, 2012; pp. 125–132. [Google Scholar]
  16. Maijer, D.; Ikeda, T. Mathematical modeling of residual stress formation in electron beam remelting and refining of scrap silicon for the production of solar-grade silicon. Mater. Sci. Eng. A 2011, 390, 188–201. [Google Scholar] [CrossRef]
  17. Vutova, K.Zh.; Mladenov, G.M. Computer simulation of the heat transfer during electron beam melting and refning. Vacuum 1999, 53, 87–91. [Google Scholar] [CrossRef]
  18. Koleva, E.; Vutova, K.; Mladenov, G. The role of ingot crucible thermal contact in mathematical modelling of the heat transfer during electron beam melting. Vacuum 2001, 62, 189–196. [Google Scholar] [CrossRef]
  19. Larikov, L.; Jurchenko, J. Structure and Characteristics of Metals and Alloys. Heat Characteristics of Metals and Alloys (in Russian); Naukova Dumka: Kiev, Russia, 1987. [Google Scholar]
  20. Samsonov, G. Chemo-Physical Properties of Elements (in Russian); Naukova Dumka: Kiev, Russia, 1965. [Google Scholar]
  21. Onufriev, S.; Petukhov, V. The thermophysical properties of hafnium in the temperature range from 293 to 2000 K. High Temp. 2008, 46, 203–211. [Google Scholar] [CrossRef]
  22. Desal, P. Thermodynamic properties of manganese and molybdenum. J. Phys. Chem. Ref. Data 1987, 16, 91–108. [Google Scholar]
  23. Samarskii, A. Theory of Difference Schemes; Dekker, M., Ed.; CRC Press: New York, NY, USA, 2001. [Google Scholar]
  24. Vutova, K.; Donchev, V.; Vassileva, V.; Koleva, E.; Mladenov, G. Investigation of Electron Beam Drip Melting by a Time-dependent Heat Model. In Proceedings of the International Conference on High-Power Electron Beam Technology “EBEAM 2012", Reno, NA, USA, 14–16 October 2012; pp. 35–41.
  25. Donchev, V.; Vutova, K. Optimization method for electron beam melting and refining of metals. J. Phys. Conf. Ser. 2013, in press. [Google Scholar] [CrossRef]
  26. Vutova, K.; Vassileva, V.; Koleva, E.; Georgieva, E.; Mladenov, G.; Mollov, D.; Kardjiev, M. Investigation of electron beam melting and refining of titanium and tantalum scrap. J. Mater. Process. Technol. 2010, 210, 1089–1094. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Vutova, K.; Donchev, V. Electron Beam Melting and Refining of Metals: Computational Modeling and Optimization. Materials 2013, 6, 4626-4640. https://doi.org/10.3390/ma6104626

AMA Style

Vutova K, Donchev V. Electron Beam Melting and Refining of Metals: Computational Modeling and Optimization. Materials. 2013; 6(10):4626-4640. https://doi.org/10.3390/ma6104626

Chicago/Turabian Style

Vutova, Katia, and Veliko Donchev. 2013. "Electron Beam Melting and Refining of Metals: Computational Modeling and Optimization" Materials 6, no. 10: 4626-4640. https://doi.org/10.3390/ma6104626

APA Style

Vutova, K., & Donchev, V. (2013). Electron Beam Melting and Refining of Metals: Computational Modeling and Optimization. Materials, 6(10), 4626-4640. https://doi.org/10.3390/ma6104626

Article Metrics

Back to TopTop