# Thermal Behavior of Ti-64 Primary Material in Electron Beam Melting Process

Jean-Pierre Bellot, Julien Jourdan, Jean-Sébastien Kroll-Rabotin, Thibault

Quatravaux, Alain Jardy

► To cite this version:

Jean-Pierre Bellot, Julien Jourdan, Jean-Sébastien Kroll-Rabotin, Thibault Quatravaux, Alain Jardy. Thermal Behavior of Ti-64 Primary Material in Electron Beam Melting Process. Materials, 2021, 14 (11), pp.2853. 10.3390/ma14112853. [hal-03845448](https://hal.science/hal-03845448)

HAL Id: hal-03845448

https://hal.science/hal-03845448

Submitted on 9 Nov 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Thermal Behavior of Ti-64 Primary Material in Electron Beam Melting Process

Jean-Pierre Bellot *ID*, Julien Jourdan, Jean-Sébastien Kroll-Rabotin ID, Thibault Quatravaux and Alain Jardy

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Bellot, J.-P.; Jourdan, J.; Kroll-Rabotin, J.-S.; Quatravaux, T.; Jardy, A. Thermal Behavior of Ti-64 Primary Material in Electron Beam Melting Process. *Materials* **2021**, *14*, 2853. https://doi.org/10.3390/ma14112853

Academic Editor: Katia Vutova

Received: 28 April 2021

Accepted: 23 May 2021

Published: 26 May 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Institut Jean Lamour—UMR CNRS 7198, LabEx DAMAS, Campus Artem, Université de Lorraine, 2 allée André Guinier, 54000 Nancy, France; julien.jourdan@univ-lorraine.fr (J.J.); jean-sebastien.kroll-rabotin@univ-lorraine.fr (J.-S.K.-R.); thibault.quatravaux@univ-lorraine.fr (T.Q.); alain.jardy@univ-lorraine.fr (A.J.)

* Correspondence: jean-pierre.bellot@univ-lorraine.fr; Tel.: +33-372-744-917

**Abstract:** The Electron Beam Melting (EBM) process has emerged as either an alternative or a complement to vacuum arc remelting of titanium alloys, since it is capable of enhancing the removal of exogenous inclusions by dissolution or sedimentation. The melting of the primary material is a first step of this continuous process, which has not been studied so far and is investigated experimentally and numerically in the present study. Experiments have been set up in a 100 kW laboratory furnace with the aim of analyzing the effect of melting rate on surface temperature of Ti-64 bars. It was found that melting rate is nearly proportional to the EB power while the overheating temperature remains roughly independent of the melting rate and equal to about 100 °C. The emissivity of molten Ti-64 was found to be 0.22 at an average temperature of about 1760 °C at the tip of the bar. In parallel, a mathematical model of the thermal behavior of the material during melting has been developed. The simulations revealed valuable results about the melting rate, global heat balance and thermal gradient throughout the bar, which agreed with the experimental values to a good extent. The modeling confirms that the overheating temperature of the tip of the material is nearly independent of the melting rate.

**Keywords:** melting; electron beam; melting temperature; numerical simulation

## 1. Introduction

Among the secondary remelting techniques, the EBM (Electron Beam Melting) process applies a high-power electron beam on the metallic material for melting, refining and controlling the casting and solidification stages. These operations ensure both the purification of the metal as it is gradually melted and the controlled solidification of the ingot in terms of structure and chemical homogeneity.

The EBM process can be used in two ways: Electron Beam Cold-Hearth Melting (EBCHM) and Electron Beam Drip Melting (EBDM). Schematic diagrams of these techniques are shown in Figure 1.

The main difference between these two methods is the inclusion of a cold hearth (circled in red in Figure 1a) in the EBCHM process as an intermediate refining stage between the melting and solidification steps. In the melting step, the electron beam heats and melts the bar tip. In the final solidification step, the molten metal is cast into a withdrawing water-cooled copper mold and solidified to form an ingot. All these operations are performed inside a vacuum chamber ($10^{-4}$ to $10^{-3}$ mbar) in order to guarantee proper operation of the electron guns and to avoid alloy contamination. The EBM methods have a number of advantages when compared to the classic vacuum arc remelting (VAR) process, that lead to a better ingot quality. In particular, melting is conducted in a higher vacuum and for a longer time, thus enabling more complete degassing and dissolution of exogenous inclusions such as low-density inclusions known as hard-alpha [1,2].
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Schematic diagrams of: (a) EBCHM; (b) EBDM.</figcaption></figure>

The EBM processes are commonly used for manufacturing refractory or reactive metals, such as niobium and titanium [2,3]. Since these techniques essentially concern very high-performance materials used in leading-edge applications demanding exceptional reliability, it is extremely important to choose the most appropriate operating parameters in order to achieve optimum product quality. Despite experimental and numerical works having been performed on the refining and solidification steps [2–5], the thermal behavior of the raw material during the melting stage has only been investigated in the case of alternative processes such as vacuum arc remelting in the literature [6,7] although it plays an important role on the process operation. The surface temperature of the tip of the load activates the volatilization mechanism, which controls the material refining and the losses in alloying elements such as aluminum in Ti-64. Furthermore, the melting rate influences the liquid pool shape and depth, which have a direct effect on the ingot solidification structure. From a modeling point of view, it is obvious that the numerical simulations must include the melting stage of the process so that numerical predictions can be strictly correlated to operating conditions such as electrical parameters without the need to specify the melt rate value as an input data of the simulation.

Thanks to the use of a well-dedicated electron beam laboratory furnace, a campaign of drip melting trials was performed with the aim of investigating the thermal behavior of a Ti-64 bar. Surface temperature of the bar tip was measured by pyrometer and infrared thermocamera, and the melting rate was recorded. As part of our modeling effort of the EBM processes, a mathematical model of the thermal behavior of the primary material (load) during melting has been developed. Details of this model are given in the present paper. The results of the model are compared to melt rate experimental data obtained with the laboratory EBM furnace, so that clear conclusions can be drawn.

## 2. Materials and Methods

### 2.1. *Experimental Setup and Procedure*

A set of experiments was conducted in a 100 kW laboratory Electron Beam Melting furnace (ALD Vacuum Technologies—Lab100, Hanau, Germany) shown in Figure 2. The EB gun was used at an electric power from 11 kW to 25 kW with a constant beam voltage of 40 kV. The beam path is controlled by the Escosys pilot program (ALD Vacuum Technologies, Hanau, Germany), which allows automatic control of beam displacement, as well as a choice of beam pattern shape and frequency. A FLIR X6540SC infrared (IR) camera (Teledyne FLIR Systems, Thousand Oaks, CA, USA) with a wavelength of 2–5 µm was set up at one side of the furnace to evaluate the temperature of the bar tip as well as qualitative IR pictures of the bar side. This camera provides temperature measurements from low ambient temperature to beyond 2000 °C with an uncertainty of ±2% of measured values. The view port was equipped
with a BaF₂ glass transparent to infrared radiation. Since the glass view was the object of
a possible quick coating with vaporized titanium and aluminum, a flow of argon in front
of the viewpoint shutter was applied to prevent deposition. A two-wavelength pyrometer
(IRCONE Modline Type-R, FLUKE Corp, Everett, WA, USA) was used to measure the surface
temperature and to calibrate the infrared camera.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** (a) Picture of the EB-Lab100 furnace; (b) schematic diagram of the Ti-64 bar melting.</figcaption></figure>

The Ti-64 bar (45 cm in length and with a section of 25.0 cm²) was placed in the bar feeder and the EB chamber vacuumed to below 5 × 10⁻⁴ mbar. The EB gun was then turned on, a first pattern being used to heat the bar tip while the ingot surface was heated by a second pattern. A few moments were spent in each trial on finding the exact required power of the beam for a steady melt according to the bar feeder speed. Initially, both the bar and ingot feeder speeds were changed manually; then, the speed of the ingot feeder was adjusted according to melting rate. Later, the bar and ingot feeder rates were set on automatic so that, when the quasi-stationary state is reached, the bar tip is maintained at the same location above the ingot crucible. During this quasi-stationary melting rate period temperature measurements using both pyrometer and IR camera were performed. After the consumption of 90% of the bar length, the beam was switched off and materials were cooled down. Finally, the post-experimental microstructure of the bar tips was observed and compared for different trials.

## 2.2. Numerical Modeling

This section describes the numerical modeling of the transient thermal behavior of
the consumable Ti-64 bar in the EBDM process. An original feature of the model is the
representation of liquid metal removal at the bar tip, which allows the model to predict the
evolution of the melt rate with time, for given operating parameters.

### 2.2.1. Heat Transfer Model
Model Equations

The modeling of the thermal behavior of the consumable bar is developed in a 2D
axisymmetric geometry, enabling a representation of the curvature of the bar tip. Within the
solid metal, heat is transferred by conduction. Under these conditions, the heat transport equation is written in the following form [8]:

$$ \rho C_p \frac{\partial T}{\partial t} = \nabla \times \left( \lambda \overleftrightarrow{\nabla} T \right) \qquad (1) $$

where T is the temperature, ρ the density, C_p the specific heat and λ the thermal conductivity. The temperature dependence of the alloy thermophysical properties is taken into account.

Dissipation of the latent heat of melting is represented using the equivalent specific heat method [9], which consists of replacing the specific heat in Equation (1) by

$$ C_p^* = C_p + \frac{\partial g_l}{\partial T} L \qquad (2) $$

where L is the latent heat of melting, and g_l is the liquid fraction.

The calculation of the equivalent specific heat requires the knowledge of the solidification path, i.e., the evolution of the liquid fraction with temperature. In the case of a multicomponent alloy (such as Ti-64), the solidification path is in general unknown. Under these conditions and because the interval between liquidus and solidus is very small, we consider a linear variation of the liquid fraction between the solidus (T_sol) and liquidus (T LIq) temperatures. Such a variation implies the assumption of a uniform dissipation of the latent heat L in the temperature interval [T_sol, T LIq], which is particularly narrow in the case of Ti-64. Equation (2) can thus be written as

$$ C_p^* = C_p + \frac{L}{T_{\text{liq}} - T_{\text{sol}}} \qquad (3) $$

Initially, the bar temperature is considered to be homogeneous, equal to the stub or pusher temperature (T = T_0) i.e., the room temperature. The symmetry condition on the bar axis is given by $(\partial T / \partial r)_{r=0} = 0$. The boundary conditions on the edge and the tip of the bar depend on the process itself.

## Boundary Conditions

At the bar tip, the kinetic energy of the electron beam is converted into thermal energy; however, the backscattering of electrons results in significant losses, which are in the range of 30% of the incident energy [10] for titanium. Thus, the heat power density applied to the tip surface of the bar is given by:

$$ \varphi_{\text{EB}} = \frac{\varepsilon_{\text{EB}} P_{\text{EB-tip}}}{\pi R^2} \qquad (4) $$

where P_EB-tip is the electrical power of the beam applied to the tip, and ε_EB is equal to 0.7. In addition, the net heat exchanged by radiation between the liquid film and the furnace wall must be taken into account. Since the surface ratio between the bar and the furnace wall is very small and in the approximation of the grey body, the thermal flux density transferred can be expressed as:

$$ \varphi_{\text{rad}}^{\text{tip}} = \sigma \varepsilon_1 (T^4 - T_w^4) \qquad (5) $$

where T_w is the temperature of the furnace wall.

At the bar lateral surface, heat transfer is controlled by thermal radiation between the material and the furnace walls with a similar expression as Equation (5):

$$ \varphi_{\text{rad}}^{\text{side}} = \sigma \varepsilon_s (T^4 - T_w^4) \qquad (6) $$
At the end of the bar, a contact resistance between the stub or pusher and the bar is taken into account using a heat transfer coefficient h<sub>sp</sub>. This coefficient is assigned a constant value of 500 W·m<sup>-2</sup>·K<sup>-1</sup>, describing a moderately effective thermal contact between the pusher and the bar [11,12].

$$ \varphi_{\text{push}} = h_{\text{sp}} (T - T_{\text{push}}) \qquad (7) $$

The integration of these flux densities (Equations (4)–(7)) over the respective surface leads to heat power defined as ε<sub>EB</sub>P<sub>EB-tip</sub>, Q<sub>tip</sub>, Q<sub>side</sub>, Q<sub>rad</sub>, Q<sub>push</sub> for respectively the effective EB heat power, heat lost by radiation at the tip, heat lost by radiation at the bar side and heat lost at the pusher contact.

The thermal properties of Ti-64 used in the model are reported in Table 1. Note that the thermal emissivity of the liquid titanium has been obtained by IR and pyrometer measurements (see Section 3.2).

<table><caption>**Table 1.** Thermal properties of Ti-64 used in the numerical simulations [13,14].</caption><thead><tr><th>Properties of Ti-64</th><th>Values or Expression</th></tr></thead><tbody><tr><td>ρ<sub>l</sub> liquid density (kg·m<sup>-3</sup>)</td><td>4100</td></tr><tr><td>ρ<sub>s</sub> solid density (kg·m<sup>-3</sup>)</td><td>4158</td></tr><tr><td>T<sub>liq</sub> liquidus (°C)</td><td>1670</td></tr><tr><td>T<sub>sol</sub> solidus (°C)</td><td>1650</td></tr><tr><td>L latent heat of melting (J·kg<sup>-1</sup>)</td><td>3.89 × 10<sup>5</sup></td></tr><tr><td>C<sub>p</sub> specific heat of solid (J·kg<sup>-1</sup>·K<sup>-1</sup>)</td><td>710</td></tr><tr><td>C<sub>p</sub> specific heat of liquid (J·kg<sup>-1</sup>·K<sup>-1</sup>)</td><td>794</td></tr><tr><td>λ<sub>s</sub> conductivity of solid (W·m<sup>-1</sup>·K<sup>-1</sup>)</td><td>18.4</td></tr><tr><td>λ<sub>l</sub> conductivity of liquid (W·m<sup>-1</sup>·K<sup>-1</sup>)</td><td>22.5</td></tr><tr><td>ε<sub>l</sub> emissivity of liquid</td><td>0.22</td></tr><tr><td>ε<sub>s</sub> emissivity of solid</td><td>0.43 [14]</td></tr></tbody></table>

### 2.2.2. Calculation Procedure

The heat transfer equation is solved using a finite volume method [15]. The numerical program called 'Ebmelting' is written in FORTRAN. A typical 160 × 1500 (r,z) orthogonal grid is applied for the spatial discretization. Furthermore, a fully implicit scheme is used for time discretization. During each time step (typically 0.5 s), we first calculate the temperature field in the material. Then, in order to simulate consumption of the bar associated to the fall of liquid metal droplets formed at the bar tip, the mesh cells whose temperature is greater than an "overheating temperature" T<sub>oh</sub> are removed from the computational domain. After mesh cell removal, the boundary conditions are set at the new bar tip for the next time step. Calculations are performed until full consumption of the bar or until a given processing time. Note that the overheating temperature T<sub>oh</sub> is an input variable for Ebmelting and is defined as the sum of the alloy liquidus temperature and a superheat.

Finally, Ebmelting calculates at each time increment the temperature field of the bar, the shape of the bar tip as well as the melting rate.

## 3. Results and Discussion

### 3.1. *Experimental Melting Rate and Overheating Temperature*

A set of eight trials has been set up in our laboratory electron beam furnace with associated melting rate and EB power as reported in Table 2. Notice that the total EB power P<sub>EB-tot</sub> has two contributions, one is the electrical power used to scan the pattern at the bar tip P<sub>EB-tip</sub> and the second is the power P<sub>EB-ing</sub> applied to the ingot top for solidification control as schematically described in Figure 2b.

$$ P_{\text{EB-tot}} = P_{\text{EB-tip}} + P_{\text{EB-ing}} \qquad (8) $$
<table><caption>**Table 2.** Experimental values for each individual run.</caption><thead><tr><th>Run #</th><th>Bar Feeder Speed (mm·min<sup>-1</sup>)</th><th>Melting Rate ṁ<sub>m</sub> (kg/h)</th><th>P<sub>EB-tip</sub> (kVA)</th><th>P<sub>EB-tot</sub> (kVA)</th><th>Measured Overheating Temperature (°C)</th></tr></thead><tbody><tr><td>1</td><td>8</td><td>5.15</td><td>5.81</td><td>11.6</td><td>1760</td></tr><tr><td>2</td><td>10</td><td>6.44</td><td>6.38</td><td>12.4</td><td>1770</td></tr><tr><td>3</td><td>12</td><td>7.73</td><td>7.65</td><td>13.6</td><td>1762</td></tr><tr><td>4</td><td>15</td><td>9.66</td><td>8.33</td><td>14.8</td><td>1760</td></tr><tr><td>5</td><td>20</td><td>12.88</td><td>9.68</td><td>17.2</td><td>1772</td></tr><tr><td>6</td><td>25</td><td>15.90</td><td>12.16</td><td>21.6</td><td>1772</td></tr><tr><td>7</td><td>27</td><td>18.58</td><td>13.06</td><td>23.2</td><td>1784</td></tr><tr><td>8</td><td>30</td><td>20.65</td><td>14.00</td><td>24.8</td><td>1779</td></tr></tbody></table>

Table 2 clearly reveals that the temperature of the bar tip is not correlated with the EB power applied, and its value remains in a range between 100 and 125 °C above the liquidus temperature. On the contrary, the melting rate increases almost linearly with the heating power as it is shown in Figure 3.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** EB power P<sub>EB-tip</sub> vs. melting rate. Slope = 1.8 kg/hr/kW.</figcaption></figure>

During the quasi-stationary state, a straightforward heat balance based on the schematic Figure 4 leads to

$$ (h_{\text{liq}}(T_{\text{oh}}) - h_{\text{sol}}(T_{0}))\dot{m}_{\text{m}} = \varepsilon_{\text{EB}}P_{\text{EB-tip}} - \left[ \dot{Q}_{\text{rad}}^{\text{tip}} + \dot{Q}_{\text{rad}}^{\text{side}} + \dot{Q}_{\text{push}} \right] \quad (9) $$

where h<sub>liq</sub> and h<sub>sol</sub> are respectively the specific enthalpy of liquid and solid, and ṁ<sub>m</sub> is the melting rate. The thermal power lost by radiation on the bar tip remains at a constant value (since the overheating temperature T<sub>oh</sub> does not change with the beam power), and the heat lost on the pusher side can be considered as negligible during the quasi-stationary regime.
<figure><img src="image_7.png" /><figcaption>Figure 4. Heat fluxes at the bar boundaries.</figcaption></figure>

If the heat lost by radiation on the bar side $\dot{Q}_{\text{rad}}^{\text{side}}$ was independent of the EB heat power applied on the bar tip, then the melting rate would be proportional to that EB power applied and the slope of the adjusted line would be equal to:

$$S = \varepsilon_{\text{EB}} (h_{\text{liq}}(T_{\text{oh}}) - h_{\text{sol}}(T_{0}))^{-1} \quad (10)$$

The analytical value of S (1.5 kg/hr/kW) does not match with the experimental slope (see Figure 3), which is equal to 1.8 kg/hr/kW. We will see later (Section 3.5) that the temperature profile on the bar side cannot be considered as independent of the beam power. Therefore, the assumption of an independent radiation heat flux on the bar side is wrong.

## 3.2. Emissivity Measurements and Temperature Profiles

In a first step, the temperature of the bar tip is continuously measured with the IR camera during a cooling test (EB is sharply switched off). The emissivity is then selected such as the measured temperature of the phase change ranges between 1650 °C ($T_{\text{sol}}$) and 1670 °C ($T_{\text{liq}}$). An emissivity of 0.22 was obtained. In a second step, the average surface temperature of the bar tip was monitored by the IR camera during a drip melting run and compared to the temperature measured with the pyrometer. Again, the value 0.22 was required to match IR and pyrometer temperatures. This value is in agreement with the findings of Choi [16] and Rai [17] (0.23 and 0.20, respectively).

An example of IR picture is presented in Figure 5 for a low (Run#1) and a high (Run#5) melting rate. It clearly emphasizes the steeper temperature gradient at high melting rate. The heat supplied by the EB diffuses more deeply into the bar at lower melting rate.

<figure><img src="image_8.png" /><figcaption>Figure 5. Temperature visualization for two conditions: (a) low melting rate—8 mm/min and (b) high melting rate—20 mm/min.</figcaption></figure>

## 3.3. Microstructure of the Bar Tip

In order to get a better understanding of the thermal treatment experienced by the bar during electron beam drip melting (EBDM), the microstructure of the remaining bar tips, after two heats performed at distinguishable melting rates, was studied. In both cases, the melt was sharply interrupted thanks to the EB switch-off. Two samples were
obtained from the experiments with a bar feeding speed of 8 and 20 mm/s (melting rate of 5.2 and 12.9 kg/h, respectively). The samples were prepared by cutting in the middle longitudinal section the first tens of millimeters from the front surface of the bar, then polished and etched by Kroll agent. Finally, they were examined under Zeiss-Axioplan 2 optical microscope (Zeiss, Iena, Germany).

Figure 6 illustrates the different grain structure in the first millimeters from the extremity of the bar. The tip of the bar is clearly distinguished on the right-hand side of each image.

<figure><img src="image_9.png" /><figcaption>**Figure 6.** Observed grain structure near the bar tip for (a) low melting rate—8 mm/min and (b) high melting rate—20 mm/min.</figcaption></figure>

The grain morphologies are equiaxed with a size noticeably larger for the lower melting rate. This is caused by the longer time this region was subjected to high temperature (above the beta transus—see below).

From these local images obtained by optical microscopy, mapping pictures of the whole bar tip could be built thanks to the automatic driven system option integrated in Axiovision software (Zeiss, Iena, Germany). Indeed, this system allows to stitch, with a controlled overlap percentage, a large number of images (each square, as the ones dashed in yellow, is a single optical microscope image). These microstructure maps are shown in the Figure 7 for both trials.

On Figure 7, the Heat Affected Zone (HAZ) can be easily unveiled because of microstructure change during heating. According to the Ti-64 phase diagram, beta grain growth occurs during heating at a temperature higher than the beta transus temperature $T_{\alpha\beta}$ (around 880 °C). Because of the higher mobility of beta Ti for grain growth due to its single-phase solution nature whereas the alpha Ti has a lower mobility due to its multiphase nature, this heating above the transus leads to a large grain structure [18] before quenching as electron beam power is turned off.

Consequently, HAZ boundaries have been sketched in dotted line on the images. On the left-hand side of this limit, primary Ti-64 that has not experienced the beta transus can be seen. On the other side, beta grain size enlargement can be noticed from the HAZ beginning to the bar tip. An average value for HAZ was approximately measured for two samples and determined to be around 21 and 11 mm at low and high melting rates, respectively.
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Assessment of the HAZ for (a) low melting rate—8 mm/min and (b) high melting rate—20 mm/min. Stitched optical microscope images for mapping (unit in dashed yellow).</figcaption></figure>

### 3.4. Aluminum Depletion by Volatilization

Because of high vacuum, the main drawback of the EBM processes is a significant loss by volatilization of alloying elements exhibiting a high vapor pressure (such as aluminum in Ti-64). This provokes a pollution of the chamber walls and makes it difficult to control the chemical composition of the as-cast ingot.

The significant temperature gradient in the liquid film makes the chemical analysis of the bar tip difficult and inaccurate. This is the reason why the volatilization losses have been obtained with a chemical analysis of the cast ingots using the Glow Discharge Optical Emission Spectrometry (GDOES) technique. Table 3 reports the aluminum content of the initial bar and the mean value of the radial ingot section for the low and high melting rate runs.

<table><caption>**Table 3.** Aluminum content in the raw material and in the final ingot.</caption><thead><tr><th></th><th>Al wt%—Initial Bar</th><th>Al wt%—Final Ingot</th></tr></thead><tbody><tr><td>Run#1—8 mm/min</td><td>6.0</td><td>4.67</td></tr><tr><td>Run#5—20 mm/min</td><td>6.0</td><td>5.39</td></tr></tbody></table>

The analyses emphasize that the EBDM at higher melting rate reduces the losses by volatilization. Since the overheating temperature does not change with the EB power, the volatilization flux on the bar tip remains roughly constant. It is therefore obvious that a higher melting rate reduces the residence time of the alloy in a liquid state and then the aluminum losses. This result agrees well with the literature [19,20] where the detrimental effect of low melting rate on aluminum depletion was established.

### 3.5. Results of the Numerical Simulation and Discussion

#### 3.5.1. Melting Rate and Overheating Temperature

The eight experimental runs have been simulated, and the excellent convergence of the heat transfer equation lead to a heat balance lower than 0.1%. One of the most interesting results provided by Ebmelting code is the melting rate profile as shown in Figure 8 for the Run#5 at 20 mm/min. Following the EB heating of the bar tip, the first droplets of liquid
metal are formed after 30 s. The melting rate increases very steeply and then levels off at a roughly constant value corresponding to the quasi-stationary state. Finally, the melting rate rises again at the end of the consumption of the bar. The oscillation of the melting rate can be readily explained by the ablation of mesh cells, whose temperatures are greater than $T_{oh}$. The amplitude of the oscillations is then strictly correlated to the mesh refinement.

<figure><img src="image_11.png" /><figcaption>**Figure 8.** Typical shape of the melting rate (example of Run#5; 20 mm/s)—Mesh refinement (r,z): 160 × 1500.</figcaption></figure>

For each experiment, a wide range of overheating temperature was set into the numerical simulation model as an input to calculate the corresponding melting rates as illustrated in Figure 8. The melting rate calculated by the model obviously decreases with the overheating temperature since a higher thermal energy is required to remove the liquid Ti-64 at the bar tip. Moreover, the heat lost by radiation $Q_{\text{rad}}^{\text{tip}}$ increases non-linearly with $T_{oh}$, which results in a lower melting rate.

Then, the measured (continuous line) and simulated (markers) melting rates were compared as shown as an example in Figure 9. The intersection of calculated and experimental melting rate gives an assessment of the overheating temperature required in the model.

<figure><img src="image_12.png" /><figcaption>**Figure 9.** Simulated melting rate versus overheating temperature $T_{oh}$ and measured temperature and melting rate (Run #2: 10 mm/min).</figcaption></figure>

The latter matches well with the measured tip temperature, as shown in Table 4, that gathers the simulated and experimental values for each individual run. A good agree-
ment between measured and indirectly calculated overheating temperatures is obtained in all cases except in Run#3 where the substantial difference could be attributed to an accumulation of experimental errors.

<table><caption>**Table 4.** Simulated and experimental values for each individual run.</caption><thead><tr><th>Run #</th><th>Bar Feeder Speed (mm/min)</th><th>Melting Rate m&#773;<sub>m</sub> (kg/hr)</th><th>P<sub>EB-tip</sub> (kVA)</th><th>Experimental Overheating Temperature (°C)</th><th>Simulated Liquid Temperature T<sub>oh</sub> (°C)</th></tr></thead><tbody><tr><td>1</td><td>8</td><td>5.15</td><td>5.81</td><td>1760</td><td>1790</td></tr><tr><td>2</td><td>10</td><td>6.44</td><td>6.38</td><td>1770</td><td>1775</td></tr><tr><td>3</td><td>12</td><td>7.73</td><td>7.65</td><td>1762</td><td>1822</td></tr><tr><td>4</td><td>15</td><td>9.66</td><td>8.33</td><td>1760</td><td>1798</td></tr><tr><td>5</td><td>20</td><td>12.88</td><td>9.68</td><td>1772</td><td>1750</td></tr><tr><td>6</td><td>25</td><td>15.90</td><td>12.16</td><td>1772</td><td>1770</td></tr><tr><td>7</td><td>27</td><td>18.58</td><td>13.06</td><td>1784</td><td>1762</td></tr><tr><td>8</td><td>30</td><td>20.65</td><td>14.00</td><td>1779</td><td>1781</td></tr></tbody></table>

Table 4 also confirms that the overheating temperature is not correlated with the melting rate, and its value remains in a range between 90 and 150 °C above the liquidus temperature ($T_{liq}$ = 1670 °C).

### 3.5.2. Thermal Profiles and HAZ

The numerical model provides at each time step of the simulation the temperature distribution in the bar. Figure 10 presents such computed thermal map under two different melting rates at $t = 1000$ s. The model predicts that the highest axial temperature gradient is confined close to the bar tip whereas the main part of the bar remains at room temperature. This thermal behavior is easily attributable to the competition between the consumption speed and the thermal diffusion velocity. Calculation of a dimensionless Péclet number allows assessment of the relative role of each of these two phenomena (ratio of advective to diffusive). It gives:

$$Pe_{\text{Run1}} = \frac{u_b L}{\alpha} = 7.3 \text{ and } Pe_{\text{Run5}} = 27 \qquad (11)$$

where L is the bar length, $u_b$ the bar speed and $\alpha$ the thermal diffusivity. Accordingly, these values testify that on the one hand the bar moves forward much faster than the heat diffuses and, on the other hand, the higher the melting rate, the steeper the temperature gradient.

<figure><img src="image_13.png" /><figcaption>**Figure 10.** Contour of the temperature after 1000 s of melting Run#1 at 8 mm/min and Run#5 at 20 mm/min.</figcaption></figure>
The concave shape of the bar tip results directly from the thermal losses by radiation along the edge of the bar. This shape remains essentially constant throughout the melt, but we can notice that the curvature of the tip is more pronounced at the beginning of the melt, when EB power is low. In that case, the effect of radiation losses is proportionally more important and Ti-64 melting is, from a thermal point of view, enhanced at the center compared to the periphery. However, the discrepancy with the real shape observed in Figure 7 can be readily explained by the fact that the model does not take into account the dynamic motion of the liquid film under forces such as wetting, buoyancy and thermocapillary [21].

As discussed above, the Heat-Affected Zone (HAZ) is the part of the bar that experiences a temperature over the beta transus $T_{\alpha\beta}$. Profiles of normalized temperature of the bar side are drawn in Figure 11 for the two melting rates (Run#1 and Run#5) and allow determination of the predicted HAZ. The normalized temperature is given by:

$$T^*(z) = \frac{T(z) - T_0}{T(z=0) - T_0} \qquad (12)$$

where z is here the distance from the bar tip.

<figure><img src="image_14.png" /><figcaption>Figure 11. Comparison of the normalized wall temperature profiles in two cases of low and high melting rates (Run#1—8 mm/min and Run#5—20 mm/min).</figcaption></figure>

The computed HAZ values for Run#1 and Run#5 (respectively 24 and 14 mm) match quite well with the measured values (respectively 21 and 11 mm) deduced from the micrographs on Figure 7.

### 3.5.3. Calculation of Thermal Balance

The global thermal balance applied to the bar was calculated during the quasi-stationary stage by the numerical model Ebmelting, and an example of the results is reported in Table 5. Definition of the variables are given in Section 2.2 except for the heat power requested to heat and melt the material:

$$\dot{Q}_{\text{melting}} = \dot{m}_{\text{m}}(h_{\text{liq}}(T_{\text{oh}}) - h_{\text{sol}}(T_0)) \qquad (13)$$
<table><caption>**Table 5.** Global thermal balance (W) of the bar calculated during quasi-stationary regime.</caption><thead><tr><th>Watt:</th><th>ε<sub>EB</sub>P<sub>EB-tip</sub></th><th>Q̇<sub>rad</sub><sup>tip</sup></th><th>Q̇<sub>rad</sub><sup>side</sup></th><th>Q̇<sub>push</sub></th><th>Q̇<sub>melting</sub></th><th>Balance (%)</th></tr></thead><tbody><tr><td>Run#1</td><td>4068</td><td>561</td><td>964</td><td>0</td><td>2542</td><td>0.02</td></tr><tr><td>Run#5</td><td>6780</td><td>529</td><td>450</td><td>0</td><td>5808</td><td>0.10</td></tr></tbody></table>

The excellent convergence of the simulation is proved by the small balance (lower than 0.1%). The heat lost at the pusher side remains negligible as long as the bar length is much larger than the HAZ. Only during the last 100 s of the melting, Q̇<sub>push</sub> is no longer negligible. Finally, as a consequence of the bar consumption, it is found that the heat lost by radiation clearly decreases when the melting rate increases, as discussed in Section 3.1.

## 4. Conclusions

Experiments of melting Ti-64 bars in a laboratory electron beam furnace were compared to a numerical model accounting for the various heat transfers in the titanium bar and at its surfaces. The experiments were designed to mimic the melting step in EBDM and EBCHM processes in order to investigate the effect of operating conditions in such processes. Melting rate and beam power heating the bar tip are adjusted altogether to achieve a quasi-steady melt. After experimentally determining the emissivity of liquid titanium, the overheating temperature at the bar tip was measured for 8 sets of operating parameters. These temperature measurements supported by the numerical results provided the undisputed evidence that the overheating temperature remains roughly independent of the melting rate (and therefore the EB power applied to the bar tip) and equal to about 100 °C above the liquidus temperature. This important finding indicates that high melting rate operations favor the control of the chemical composition of the re-melted material, notably in terms of aluminum depletion by volatilization. The heat affected zone in the bar was optically determined from the grain structure after quenching for two of these sets that differ significantly (melting rate more than doubled). While the overheating temperature is constant, the HAZ was found to strongly depend on operating parameters. The predictions of the numerical model for overheating temperature and HAZ are in good agreement with the experimental measurements, which validates the hypotheses on which the model is built.

Numerical simulations and experiments together demonstrate that the consumption speed of the raw material is much faster than the thermal diffusion velocity, and as a consequence, the higher the melting rate, the steeper the temperature gradient and the smaller the HAZ.

Overall, the experimental results on the melting of a titanium bar by EBM provide knowledge on this step that conditions industrial processes such as EBDM and EBCHM, but that was only scarcely studied in the scientific literature. The associated numerical model not only provides a way to quantify the respective contributions of several heat fluxes in the process, but since it has been validated against experiments, it also provides a relatively simple yet predictive model of EBM that can be used as input in further studies of the downstream steps that distinguish EBDM and EBCHM processes.

**Author Contributions:** Conceptualization, J.-P.B. and A.J.; methodology, J.-P.B.; software, A.J.; validation, J.-P.B., J.J. and A.J.; formal analysis, J.-P.B., J.-S.K.-R., T.Q. and A.J.; experimental investigation, J.J.; writing—original draft preparation, J.-P.B.; writing—review and editing, J.-P.B., A.J., J.-S.K.-R., and T.Q. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data presented in this study are available on request from the corresponding author.
**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses or interpretation of data; in the writing of the manuscript or in the decision to publish the results.

## References

1. Sandell, V.; Hansson, T.; Roychowdhury, S.; Mansson, T.; Delin, M.; Akerfeldt, P.; Antti, M.L. Defects in Electron Beam Melted Ti-6Al-4V: Fatigue Life Prediction Using Experimental Data and Extreme Value Statistics. *Materials* **2021**, *14*, 640. [CrossRef] [PubMed]

2. Bellot, J.P.; Foster, B.; Hans, S.; Hess, E.; Ablitzer, D.; Mitchell, A. Dissolution of hard-alpha inclusions in liquid titanium alloys. *Met. Trans. B* **1997**, *28*, 1001–1010. [CrossRef]

3. Cen, M.J.; Liu, Y.; Chen, X.; Zhang, H.W.; Li, Y.X. Calculation of flow, heat transfer and evaporation during the Electron Beam Cold Hearth Melting of Ti-6al-4V Alloy. *Rare Met. Mat. Eng.* **2020**, *49*, 833–841.

4. Bellot, J.P.; Dussoubs, B.; Reiter, G.; Flinspach, J. A comprehensive numerical modelling of electron beam cold hearth refining and ingot consolidation of Ti alloys. *Rare Met. Mat. Eng.* **2006**, *35* (Suppl. 1), 93–96.

5. 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. [CrossRef]

6. Bertram, L.A.; Zanner, F.J. Electrode tip melting simulation during vacuum arc remelting of Inconel 718. In Proceedings of the Modeling and Control of Casting and Welding Processes, Santa Barbara, CA, USA, 12 January 1986; p. 95.

7. Jardy, A.; Falk, L.; Ablitzer, D. Energy exchanges during vacuum arc remelting. *Ironmak. Steelmak.* **1992**, *19*, 226–232.

8. Carslaw, H.S.; Jaeger, J.C. *Conduction of Heat in Solids*; Oxford University Press: Oxford, UK, 1959.

9. Zanner, F.J.; Bertram, L.A. Vacuum arc remelting: An overview. In Proceedings of the 8th International Conference on Vacuum Metallurgy, Linz, Austria, 30 October–4 November 1985; Volume 1, pp. 512–552.

10. Shiller, S.; Heisig, U.; Panzer, S. *Electron Beam Technology*; Publ. John Wiley&Sons: New York, NY, USA, 1982.

11. Clites, P.G.; Beall, R.A. *A Study of Heat Transfer to Water-Cooled Copper Crucibles during Vacuum Arc Melting*; Report No. 7035; United States Bureau of Mines: Albany, OR, USA, 1967.

12. Khasin, G.A.; Bigashev, V.Z.; Ermanovich, N.A. Contact phenomena at bottom of an ingot in ESR and VAR processes. *Steel USSR* **1973**, *3*, 383–385.

13. Boivineau, M.; Cagran, C.; Doytier, D.; Eyraud, V.; Nadal, M.; Wilthan, B.; Pottlacher, G. Thermo- physical properties of solid and liquid TA6V titanium alloy. *Int. J. Thermophys.* **2006**, *27*, 507–529. [CrossRef]

14. Shur, B.A.; Peletskii, V.E. The Effect of Alloying Addition on the emissivity of Titanium in the Neighborhood of Polymorphous Transformation. *J. High Temp. 2004*, *42*, 414–420. [CrossRef]

15. Patankar, S.V. *Numerical Heat Transfer and Fluid Flow*; Hemisphere Publishing Corporation: Washington, DC, USA, 1980.

16. Choi, W.; Jourdan, J.; Matveichev, A.; Jardy, A.; Bellot, J.P. Kinetics of Evaporation of Alloying Elements under Vacuum: Application to Ti alloys in Electron Beam Melting. *High Temp. Mater. Process.* **2017**, *36*, 815–823. [CrossRef]

17. Rai, R.; Burgardt, P.; Milewski, J.O.; Lienert, T.J.; DebRoy, T. Heat transfer and fluid during electron beam welding of 21Cr-6Ni-9Mn steel and Ti-6Al-4V alloy. *J. Phys. D Appl. Phys.* **2009**, *42*, 1–12. [CrossRef]

18. Kobryn, P.A.; Semiatin, S.L. Microstructure and texture evolution during solidification processing of Ti-6Al-4V. *J. Mater. Process. Technol.* **2003**, *135*, 330–339. [CrossRef]

19. Bellot, J.P.; Hess, E.; Ablitzer, D. Aluminium volatilization and inclusion removal in the electron beam cold hearth melting of Ti alloys. *Met. Mater. Trans. B* **2000**, *31B*, 845–854. [CrossRef]

20. Westterberg, K.W.; McClelland, M.A. *Modeling of Material and Energy Flow in an EBCHR Casting System*; No. UCRL-JC-118172, CONF-9411237-1); Bakish, R., Ed.; Lawrence Livermore National Lab.: Livermore, CA, USA, 1994; pp. 185–201.

21. Bhar, R.; Jourdan, J.; Descotes, V.; Jardy, A. An experimental study of the inclusion behavior during maraging steel processing. *Met. Res. Technol.* **2019**, *116*, 517.